[EMBOSS] notseq and fasta definition headers

Andres Pinzon andrespinzon at gmail.com
Tue Jun 17 19:49:59 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/<http://bioinf.ibun.unal.edu.co/%7Eapinzon/>
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