BOL: Raku script to calculate GC content !

BioScripts
LEGE
Raku script to calculate GC content !

Raku script to calculate GC content !

By LEGE 325 days ago

sub calculate-gc-content(Str $sequence) {
    my $gc-count = $sequence.comb(/<[GCgc]>/).elems;
    my $total-bases = $sequence.chars;

    return $gc-count / $total-bases * 100;
}

my $dna_sequence = "ATGCGCTAAAGCGCGCGCCTTACGCGCGCGCGC";
my $gc_content = calculate-gc-content($dna_sequence);

say "DNA Sequence: $dna_sequence";
say "GC Content: $gc_content%";