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%	140.6 / NA	39	T29A (ACG→GCG)	rspA	bifunctional D‑altronate/D‑mannonate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/18);  total (21/18)

ACCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTG  >  NC_000913/1654923‑1655190                                                                                                                                       |                                                                                                                                      aCCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATa                                                                                                                                   <  1:56139/139‑1 (MQ=255)   cTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATAcc                                                                                                                                 >  2:704038/1‑139 (MQ=255)               aCTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAgtggt                                                                                                                     <  1:70244/139‑1 (MQ=255)                  gCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGAt                                                                                                                  >  2:38883/1‑139 (MQ=255)                                    ggTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAAtt                                                                                                >  2:1089416/1‑139 (MQ=255)                                        cgcATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGc                                                                                            >  1:173803/1‑139 (MQ=255)                                               cgACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGAc                                                                                     <  2:173803/139‑1 (MQ=255)                                               cgACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACG‑CCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACa                                                                                    >  2:75212/1‑139 (MQ=255)                                                  ccAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAgg                                                                                  >  2:19755/1‑139 (MQ=255)                                                        aGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACaa                                                                            <  1:1089416/139‑1 (MQ=255)                                                            gCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAAc                                                                        >  1:947510/1‑139 (MQ=255)                                                              ggACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCt                                                                                                                            <  2:68603/85‑1 (MQ=255)                                                              ggACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCt                                                                                                                            >  1:68603/1‑85 (MQ=255)                                                               gACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTc                                                                     >  1:129970/1‑139 (MQ=255)                                                                            cTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAAtt                                                                                                >  1:546163/1‑99 (MQ=255)                                                                            cTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAAtt                                                                                                <  2:546163/99‑1 (MQ=255)                                                                              gCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATa                                                                                                                                   >  1:808064/1‑62 (MQ=255)                                                                              gCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATa                                                                                                                                   <  2:808064/62‑1 (MQ=255)                                                                              gCAAATAAGAGGCCACGGAAAGCTCACG‑CCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTc                                                     <  1:75212/139‑1 (MQ=255)                                                                               cAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTc                                                                                                                              >  2:776980/1‑66 (MQ=255)                                                                               cAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTc                                                                                                                              <  1:776980/66‑1 (MQ=255)                                                                                    aaGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGtt                                                <  1:794467/139‑1 (MQ=255)                                                                                     agagGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTc                                               >  1:441016/1‑139 (MQ=255)                                                                                                    cTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATa                                                                                                                                   <  1:938379/40‑1 (MQ=255)                                                                                                    cTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATa                                                                                                                                   >  2:938379/1‑40 (MQ=255)                                                                                                    cTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGc                                     >  2:611225/1‑134 (MQ=255)                                                                                                    cTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGc                                <  1:19755/139‑1 (MQ=255)                                                                                                    cTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACCAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGc                                     <  1:611225/134‑1 (MQ=255)                                                                                                         gTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTc                                                                                                                              <  2:33325/40‑1 (MQ=255)                                                                                                         gTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTc                                                                                                                              >  1:33325/1‑40 (MQ=255)                                                                                                                  gggTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCAt                  <  2:947510/139‑1 (MQ=255)                                                                                                                   ggTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATg                 >  2:921810/1‑139 (MQ=255)                                                                                                                     tGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGaa               >  1:465600/1‑139 (MQ=255)                                                                                                                       gCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGcc                                                                                           >  2:1131003/1‑61 (MQ=255)                                                                                                                       gCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGcc                                                                                           <  1:1131003/61‑1 (MQ=255)                                                                                                                       gCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAAt             >  2:76646/1‑139 (MQ=255)                                                                                                                                 ggCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGa                       >  2:362117/1‑119 (MQ=255)                                                                                                                                 ggCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGa                       <  1:362117/119‑1 (MQ=255)                                                                                                                                  gCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTg  <  1:76646/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      ACCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTG  >  NC_000913/1654923‑1655190

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