[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