[EMBOSS] notseq and fasta definition headers
Andres Pinzon
andrespinzon at gmail.com
Tue Jun 17 19:47:17 UTC 2008
Hi,
Im using notseq to obtain a subset of fasta seqs from a multiple fasta file:
notseq -junkoutseq 1000-1.fasta -sequence 7135seqs.fasta -exclude
@xaa.list.fasta -outseq leftSeqs.fast
The output is correct, but notseq changes the definition in the fasta
headers, so if the fasta header in "xaa.list.fasta" was:
lcl|29855|ORF26673_6
the corresponding fasta header in sequence in 1000-1.fasta is:
29855
Is there a way to tell "notseq" to keep the original fasta headers intact?
Thanks in advance,
--
Andrés Pinzón cPhD
http://bioinf.ibun.unal.edu.co/~apinzon/
Bioinformatics Center, Colombia EMBnet node
http://bioinf.ibun.unal.edu.co
Tel +57 3165000 ext 16961 Fax +571 3165415
Micology and Phytopathology Laboratory - Los Andes University.
http://bioinf.uniandes.edu.co
Tel +571 3394949 ext. 2768
More information about the EMBOSS
mailing list