Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,655,057	T→C	T29A (ACG→GCG)	rspA ←	bifunctional D‑altronate/D‑mannonate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,655,057	0	T	C	100.0%	79.7 / NA	23	T29A (ACG→GCG)	rspA	bifunctional D‑altronate/D‑mannonate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/13);  total (10/13)

AAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAAC  >  NC_000913/1654935‑1655188                                                                                                                             |                                                                                                                                    aaaCTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAgtg                                                                                                                     >  2:821485/1‑139 (MQ=255)   aCTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAgtggt                                                                                                                   <  1:754737/139‑1 (MQ=255)      gCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATa                                                                                                                                 <  2:655749/122‑1 (MQ=255)      gCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATa                                                                                                                                 >  1:655749/1‑122 (MQ=255)         aGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGAtttt                                                                                                             <  2:18106/139‑1 (MQ=255)              tCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAAtg                                                                                                        <  1:821485/139‑1 (MQ=255)                     tACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCt                                                                                                                          <  2:476231/114‑1 (MQ=255)                     tACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCt                                                                                                                          >  1:476231/1‑114 (MQ=255)                        ggTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGa                                                                                                     <  1:1278507/132‑1 (MQ=255)                        ggTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGa                                                                                                     >  2:1278507/1‑132 (MQ=255)                        ggTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAAtt                                                                                              <  1:116218/139‑1 (MQ=255)                             gcATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGcc                                                                                         <  1:437950/139‑1 (MQ=255)                                                 cGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACt                                                                     >  2:477252/1‑139 (MQ=255)                                                          gTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCt                                                              <  2:500501/139‑1 (MQ=255)                                                                    aaaTAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATa                                                                                                                                 >  2:1101984/1‑60 (MQ=255)                                                                    aaaTAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATa                                                                                                                                 <  1:1101984/60‑1 (MQ=255)                                                                          gagGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTa                                                           <  2:1191746/124‑1 (MQ=255)                                                                          gagGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTa                                                           >  1:1191746/1‑124 (MQ=255)                                                                                           aCGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTg                           >  1:1330768/1‑139 (MQ=255)                                                                                             gTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTc                         >  2:462285/1‑139 (MQ=255)                                                                                                             gCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAAt           >  2:174640/1‑139 (MQ=255)                                                                                                                      aGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAAc  <  2:105993/139‑1 (MQ=255)                                                                                                                      aGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAAc  <  1:477252/139‑1 (MQ=255)                                                                                                                             |                                                                                                                                    AAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATT‑‑GAGGGTGGCATCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAAC  >  NC_000913/1654935‑1655188

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑