Off topic:Re-Perl How To Isolate Fasta Sequences With A Specific Keyword
0
0
Entering edit mode
13.1 years ago
Raghul ▴ 200

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 • 1.8k views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 3381 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6