Ler uma sequências fasta e separar por tamanho [Bioinformática]

Descrição

Esconder código-fonte

#!/usr/bin/perl

open in,"<$ARGV[0]";
$i=0;
while($line = <in>){
   chomp $line;
   if($line =~ m/>/){
         $key = $line;      
   }else{
      $fasta{$key} .= $line;
   }
}

$i=0;
while( ($key,$seq) = each %fasta){
      if( $i < 9000){
         open out,">>1000/$key_9000.fasta";
         print out $key,"\n";
         print out $seq,"\n";
         close out;
      }elsif($i < 18000){
         open out,">>1000/$key_9000_18000.fasta";
         print out $key,"\n";
         print out $seq,"\n";
         close out;         
         
      }else{
         open out,">>1000/$key_18000_20504.fasta";
         print out $key,"\n";
         print out $seq,"\n";
         close out;   
         
      }
   
}
close in;