bioperl可以读取fastq文件,对其进行处理,而且一般是双端的序列,可以同时读取两个文件,方法如下:
while ( my $obj1=$fq1->next_seq() and my $obj2=$fq2->next_seq() ) {
my ($id1,$id2)=($obj1->id,$obj2->id);
....
}
这样就可以同时读取read1、read2中的对应序列。
更多perl语言知识可学习 Perl语言高级编程 !
如果觉得我的文章对您有用,请随意打赏。你的支持将鼓励我继续创作!