Skip to main content
3 of 10
deleted 20 characters in body

How to properly annotate errors and variations present in sequences in a GenBank file format

I am maintaining a database of sequence files describing a collection of bacterial strains and plasmids that we have in our lab.

I am confronted to two sorts of issues:

  • some sequence files contain errors
  • other sequence files describe molecules that themselves contain variations (mutations, insertion or deletions) that occurred during the molecular cloning process.

I want to annotate the sequences to mark those errors and variations using the GenBank format. However, from the feature keys of the GenBank format described on the INSDC website, it is not clear to me which feature key should be used in each case.

To better explain my case, here below are what seems to me the most relevant feature keys, and my comment of the problem I have with them:

Feature key           variation
Definition            a related strain contains stable mutations from the same
                      gene (e.g., RFLPs, polymorphisms, etc.) which differ
                      from the presented sequence at this location (and
                      possibly others);
...
Comment               used to describe alleles, RFLP's,and other naturally 
                      occurring mutations and  polymorphisms;
                  >>> variability arising as a result of genetic 
                      manipulation (e.g. site directed mutagenesis) should 
                      be described with the misc_difference feature;

According to the line marked >>> above, the feature key misc_difference should be used...

Feature Key           misc_difference
Definition            feature sequence is different from that presented in the
                      entry and cannot be described by any other difference
                      key (old_sequence, variation, or modified_base);
...
Comment               the misc_difference feature key should be used to
                      describe variability that arises as a result of 
                      genetic manipulation (e.g. site directed mutagenesis);
                  >>> use /replace="" to annotate deletion, e.g.
                      misc_difference 412..433                                                     
                                      /replace=""  

But now, according to the example marked >>> above, it means the variation (here a deletion) is NOT in the considered entry... So how to indicate variations (mutation, insertion, deletion, errors) that are actually present in an entry ?