Raku script to calculate GC content !
sub calculate-gc-content(Str $sequence) { my $gc-count = $sequence.comb(//).elems;...my $dna_sequence = "ATGCGCTAAAGCGCGCGCCTTACGCGCGCGCGC"; my $gc_content = calculate-gc-content($dna_sequence);...105 days ago
Perl script to calculate GC content !
#!/usr/bin/perl sub calculate_gc_content { my ($sequence) = @_; $sequence = uc($sequence);...CGAGCGAAGCGCTAGATCGATCGATCGATCGATCGATCGATCGATCGATCG"; my $gc_content = calculate_gc_content($dna_sequence); p...89 days ago
Bash script to calculate difference between column !
Space Separated awk '{ $5 = $2 - $4 } 1' inputput.txt > outdiff.txt Or with tab separation: awk 'BEGIN { OFS = "\t" } { $5 = $2 - $4 } 1' inputput.txt > outdiff.txt87 days ago
Perl script to parse blast results and plot basic stats !
...print "Bit Score: ", $entry->{BITSCORE}, "\n"; print "---------------------------\n"; } # Calculate additional statistics my $avg_percent_identity = calculate_average(\@blast_entries, 'PER...87 days ago
Perl script to calculate the basic stats of the assembled genome !
...gths = sort { $b $a } @contig_lengths; # Calculate additional statistics my $mi.../ $num_contigs; my $median_contig_length = calculate_median(\@contig_lengths);...print "$length\n"; } # Subroutine to calculate N50 sub calculate_n50 {...87 days ago
Python script for basic stats of the assembled genome !
...rder contig_lengths.sort(reverse=True) # Calculate additional statistics min_co...gth = statistics.median(contig_lengths) # Calculate N50 def calculate_n50(length...gth: {median_contig_length}") print(f"N50: {calculate_n50(contig_lengths)}") print...87 days ago