[EMBOSS] Bug in extractfeat ?
Charles Plessy
charles-listes-emboss at plessy.org
Tue Nov 19 12:26:45 UTC 2013
Dear EMBOSS developers,
I have used extractfeats in EMBOSS 6.4 to extract some features of GenBank
files with commands like the following.
$ extractfeat TRG at .gb -filter -join TRUE -tag 'standard_name' -value "TRGC1" -describe standard_name -osformat2 fasta
>NG_007033_36562_39077 [C_region] (standard_name="TRGC1") Mus musculus T cell receptor gamma chain (Tcrg) on chromosome 13.
acaaaaggcttgatgcagacatttcccccaagcccactattttccttccttctgttgctg
aaacaaatctccataagactgggacatacctttgtctccttgaaaagttctttccagatg
ttataagggtgtattggaaagaaaaggatggcaatactatcctggactcccaggaagggg
atacgctgaagactaacgacacatacatgaagtttagctggcttacagtgcctgaaaggg
caatggggaaagagcacagatgtattgtcaaacatgagaacaacaaaggaggagcagatc
aagagattttctttccttcaataaagaaagttgctgtgagtaccaagcctacaacttgct
ggcaagataaaaatgatgtgctgcagctccagttcacgatcacctctgcctactacacct
accttctcctgctcctcaagagtgtgatctacttggccatcatcagcttctctctgctta
gaagaacatctgtctgtggcaatgagaagaaatcctaa
However, with EMBOSS 6.6, there is no output anymore.
$ extractfeat TRG at .gb -filter -join TRUE -tag 'standard_name' -value "TRGC1" -describe standard_name -osformat2 fasta
Warning: No sequences written to output file 'stdout'
The TRG at .gb file was downloaded from NCBI's website and a copy is available at
the following URL.
http://source.clonotyper.branchable.com/?p=source.git;a=blob;f=references/TRG@.gb;h=9d48fafbb7e4209a462dbb1ab552195c00e920e3;hb=HEAD
The commands ran on a Debian system with the Debian packages 6.4.0-2 or
6.6.0-1.
Is that a bug in EMBOSS > 6.4 ?
Have a nice day,
Charles Plessy
--
Charles Plessy, Ph.D. – RIKEN Center for Life Science Technologies
Division of Genomic Technologies – Genomics Miniaturization Technology Unit
1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa 230-0045 Japan
http://population-transcriptomics.org/
More information about the EMBOSS
mailing list