5
seqretsplit - split a sequence (or stream) into a single file per sequence
9
seqretsplit file1 file2 ..
15
The script will split all sequences from fasta file(s) (or stdin) to
16
individual files. The filename is the sequence ID (everything before
17
the 1st whitespace in a FASTA header). Currently it doesn't check to
18
see that it isn't overwriting an existing file so IDs should be unique
20
This is inspired by EMBOSS seqretsplit tool.
26
User feedback is an integral part of the evolution of this and other
27
Bioperl modules. Send your comments and suggestions preferably to
28
the Bioperl mailing list. Your participation is much appreciated.
30
bioperl-l@bioperl.org - General discussion
31
http://bioperl.org/wiki/Mailing_lists - About the mailing lists
35
Report bugs to the Bioperl bug tracking system to help us keep track
36
of the bugs and their resolution. Bug reports can be submitted via
39
https://redmine.open-bio.org/projects/bioperl/
43
Jason Stajich E<lt>jason_AT_bioperl_DOT_orgE<gt>
49
my $in = Bio::SeqIO->new(-format => 'fasta',
51
while( my $s = $in->next_seq ) {
52
Bio::SeqIO->new(-format => 'fasta',
53
-file => ">".$s->id.".fa")->write_seq($s);