Shulou(Shulou.com)06/01 Report--

Xiaobian to share with you how bioperl directly read the fasta sequence of the output gz compression format, I believe most people still do not know how to share this article for your reference, I hope you have a lot of harvest after reading this article, let's go to understand it together!

bioperl directly reads the fasta sequence of the output gz compressed format, code example:

die "perl $0 " unless(@ARGV==3);use Math::BigFloat;use Bio::SeqIO;use Bio::Seq;use Data::Dumper;use PerlIO::gzip;open $FI, ":gzip", "$ARGV[2]" or die "$! ";$in = Bio::SeqIO->new(-fh => $FI , -format => 'Fasta');$out = Bio::SeqIO->new(-fh => $FO , -format => 'Fasta');

Another way of writing is also possible:

use Bio::SeqIO;use Bio::Seq;use Data::Dumper;use PerlIO::gzip;open FI, ":gzip", "$ARGV[2]" or die "$! ";$in = Bio::SeqIO->new(-fh => \*FI , -format => 'Fasta');$out = Bio::SeqIO->new(-fh => \*FO , -format => 'Fasta'); above is "bioperl how to directly read the fasta sequence of output gz compressed format" All the contents of this article, thank you for reading! I believe that everyone has a certain understanding, hope to share the content to help everyone, if you still want to learn more knowledge, welcome to pay attention to the industry information channel!

Welcome to subscribe "Shulou Technology Information " to get latest news, interesting things and hot topics in the IT industry, and controls the hottest and latest Internet news, technology news and IT industry trends.