[Bioperl-l] Fasta headers

Sean Davis sdavis2 at mail.nih.gov
Mon Jan 3 10:27:05 EST 2005


Rasa,

You can parse the fasta file with seqio, but if you only want the headers 
"as-is", something like this from the command line might do:

cat fastafile.fa > perl -e 'while (<>) {print "$_\n" if ($_ =~ /^>/)}'

Sorry, I didn't test this....

Sean

----- Original Message ----- 
From: "Rasa Gulbinaite" <rasa at obj.hopto.org>
To: <bioperl-l at portal.open-bio.org>
Sent: Monday, January 03, 2005 9:41 AM
Subject: [Bioperl-l] Fasta headers


> Hello,
>
> i'm new to bioperl and a bit confused with fasta file headers. I'm working
> with SNPs and i would like to get only the fasta headers form the fasta
> file, not the sequences. What would be the best way to do this? Thank you.
>
> Rasa
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at portal.open-bio.org
> http://portal.open-bio.org/mailman/listinfo/bioperl-l
> 




More information about the Bioperl-l mailing list