[EMBOSS] extractfeat problem with complement() in CDS
Anne Pajon
ap12 at sanger.ac.uk
Thu Nov 12 17:34:41 UTC 2009
Dear,
I have an embl file containing these lines
FT CDS complement(2232788..2233492)
FT /systematic_id="SARLGA251_20000"
when I use
extractfeat -sequence embl::LGA251_reduced.embl -type CDS
on emboss version 6.1, the output contains these lines for the
corresponding CDS:
>SARLGA251_2232788_2233492 [CDS]
t
and on emboss version 6.0, the output contains these lines:
>SARLGA251_2232788_2233492 [CDS]
atggctagtgttttatatcaacatggtacattaggcacattaatggcaggtttattaaaa
ggaactgcatcaataaatgaattattgcaacatggtgatttaggtatcgctacactaaca
ggttcaaatggtgaggtcatctttttagatggaaaagcctatcatgcaaatgaacataaa
gaatttgtggaattaaaaggtgatgagttaacaccatatgcaactgtaactaaatttgta
gcagatacaagctatgaaacgaaagataaatcttctgaagcagtttttgcagaaattaag
gaaaagatgttgagtgaaaatttattttctgcagtaaaaatttcaggtgtgtttaaaaaa
atgcatgtacgtatgatgccggcacaagagccaccttatacacgtttaattgattcagct
agaagacaacctgaacaaactgaaacatacgtcaaaggttcagtcgttgggttctttaca
ccagaactattccatggtatcgggtcagcaggattccatgtacactttgcaaatgatgat
cgtaactttggtggacatgtcttagattttgaagtagaagatgttaaggtagaaatccaa
aatatagaaacatttgaacagcattttccaattcaagatgaagatttcataaatgcaaat
attgactataaagatattgcagacgaaattagagaagctgaataa
Does anyone know why it is not extracting the correct sequence for the
CDS features with the latest version of emboss? I am running it on Mac
OSX and I've installed it with fink.
Thanks in advance for your help.
Kind regards,
Anne.
--
Dr Anne Pajon - Pathogen Genomics
Sanger Institute, Wellcome Trust Genome Campus, Hinxton
Cambridge CB10 1SA, United Kingdom
+44 (0)1223 494 798 (office) | +44 (0)7958 511 353 (mobile)
--
The Wellcome Trust Sanger Institute is operated by Genome Research
Limited, a charity registered in England with number 1021457 and a
company registered in England with number 2742969, whose registered
office is 215 Euston Road, London, NW1 2BE.
More information about the EMBOSS
mailing list