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 _______________________________________________ EMBOSS mailing list [email protected] http://lists.open-bio.org/mailman/listinfo/emboss
