Create genome scaffolding with Perl
...return(scalar(reverse($seq))); } sub getConsensus { my ($b1, $b2) = @_; i...)){ ## equal bases, or absent bases, so consensus is easy return($b1);...lacementSeqs{$sName}{score}) && ($consensusLength > $replacementSeqs{$sNa...2447 days ago
Perl script to create a consensus of nucleotide sequences !
...warnings; my @instances = qw ( AAAAA ATCGA ATAAA ); my @instances2 = qw ( AAAAA AACGA ATAAA AGAAA AGAAA); print consensus(@instances),"\n"; # ATAAA print consensus(@instances2),"\n"; # AT...2170 days ago