https://bioinformaticsonline.com/snippets/view/44437/perl-script-to-calculate-gc-content? https://bioinformaticsonline.com/snippets/view/44437/perl-script-to-calculate-gc-content Tue, 30 Jan 2024 05:20:10 -0600 https://bioinformaticsonline.com/snippets/view/44437/perl-script-to-calculate-gc-content #!/usr/bin/perl sub calculate_gc_content { my ($sequence) = @_; $sequence = uc($sequence); # Convert the sequence to uppercase my $gc_count = () = $sequence =~ /[GC]/g; my $total_bases = length($sequence); my $gc_content = ($gc_count / $total_bases) * 100; return $gc_content; } # Example usage: my $dna_sequence = "ATGCGCTAAAGCGAGCGAAGCGCTAGATCGATCGATCGATCGATCGATCGATCGATCGATCG"; my $gc_content = calculate_gc_content($dna_sequence); printf "GC content: %.2f%%\n", $gc_content;]]> BioStar