This is a test version of Biostars. For the public version, visit https://www.biostars.org.
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);
  }
}
perl fasta sequence retrieval

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

0 answers

No answers yet.

Log in to answer this question.