Re-Perl How To Isolate Fasta Sequences With A Specific Keyword
Hi to neilfws & all others Thanks for the response. I made-up the code, please correct it. Can you people suggest me tutorial links or textbooks for Bioperl.
#!/usr/bin/perl -w
use strict;
use Bio::SeqIO;
my $seqin = Bio::SeqIO->new(-file => "euplotes.txt", -format => "fasta");
my $seqout = Bio::SeqIO->new(-file => ">outfile.txt", -format => "fasta");
while(my $seq = $seqin->next_seq)
{
if($seq->desc) =~/^FULL-LENGTH\s+/ {
$seqout->write_seq($seq);
}
}
• 1,965 views
•
link
0 answers
No answers yet.
Log in to answer this question.
please post follows up on your original questions, if you have a new and specific question make sure to detail the problem
Ok thanks, sorry for the confusion created