Results for "my"
Bio-Scripts
Perl script to find overlaps between two bed files !
#!/usr/bin/perl use strict; use warnings; # Check if the correct number of arguments are provided if (@ARGV != 2) { die "Usage: $0 file1.bed file2.bed\n"; } # Read the contents of the two BED files my $file1 = shift @ARGV; my $file2 = shift @ARGV; open my $fh1, '151 days ago
Perl script to parse VCF file !
...t.vcf die "Usage: ./parse_vcf.pl input.vcf\n" unless @ARGV; my $vcf_file = shift @ARGV; my @vcf_entries = parse_vcf($vcf...t "---------------------------\n"; my %chromosome_counts; for my $entry (@vcf_entries) {...151 days ago
Perl script to parse blast results and plot basic stats !
...st.pl blast_result.txt\n" unless @ARGV; my $blast_file = shift @ARGV; my @blast_entries = parse_blast(...t detailed information for each entry for my $entry (@blast_entries) {...; } # Calculate additional statistics my $avg_percent_identity = calcu...151 days ago
Perl script to calculate the basic stats of the assembled genome !
...bject to read the FASTA file my $seqio = Bio::SeqIO->new(-fil...stics my $total_length = 0; my $num_contigs = 0; my @contig...culate additional statistics my $min_contig_length = $contig_...length distribution foreach my $length (@contig_lengths) {...151 days ago
