To get fasta header
0
0
Entering edit mode
7.8 years ago

I am trying to open the fasta file with sequences. My bio perl script opens the sequence but not with fasta header. How can I get fasta header with bio perl. My script:

#!/bin/perl -w
use Bio::SeqIO;

$seqio_obj = Bio::SeqIO->new(-file=> "no_plasmid.fasta", -format=>"fasta");
$seq_obj = $seqio_obj->next_seq;
$acc = $so->accession_number;

   while($seq_obj = $seqio_obj->next_seq){

             print $seq_obj-> seq, "\n";
}
RNA-Seq perl bio perl • 1.9k views
ADD COMMENT
1
Entering edit mode
  • What is $so?
  • Have you tried $seq_obj->id?
ADD REPLY
0
Entering edit mode

Hi bandanaschapagain,

Please edit your post to clarify the questions below or add a comment to reflect its status.

ADD REPLY

Login before adding your answer.

Traffic: 1590 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