Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,163,714	T→C	L152P (CTG→CCG)	thiK →	thiamine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,163,714	0	T	C	100.0%	244.8 / NA	69	L152P (CTG→CCG)	thiK	thiamine kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (42/27);  total (42/27)

CGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCTGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCACGCCGGAAATTTAGTGCATAGCGCGTCAGGGTTAAAAC  >  NC_000913/1163623‑1163809                                                                                            |                                                                                                cGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAg                                                                                         >  2:1288967/1‑100 (MQ=255)     tAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGcgcg                                                                                     >  1:1592296/1‑100 (MQ=255)      aaCGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCgg                                                                                    <  1:1678916/100‑1 (MQ=255)       aCGGTGTTGCCGGTACTGGAACTGTACTTGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAGTTTTAAAACGTCCGCGCAAAGAGCggg                                                                                   <  1:738067/100‑1 (MQ=255)       aCGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCggg                                                                                   >  1:1224879/1‑100 (MQ=255)       aCGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCggg                                                                                   >  2:1654593/1‑100 (MQ=255)          cTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGTGCAAAGCGCGGGAAc                                                                                >  2:1794535/1‑100 (MQ=255)           tGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAAcc                                                                               <  2:1138974/100‑1 (MQ=255)              tGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACg                                                                            >  1:878146/1‑100 (MQ=255)                ccGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGc                                                                          <  1:1178374/100‑1 (MQ=255)                ccGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGc                                                                          >  2:985146/1‑100 (MQ=255)                  gTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCt                                                                        <  1:1817343/100‑1 (MQ=255)                   ttACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCtt                                                                       >  1:1317601/1‑100 (MQ=255)                        ggAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGc                                                                  >  2:512234/1‑100 (MQ=255)                         gAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCt                                                                 <  2:668511/100‑1 (MQ=255)                         gAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCt                                                                 <  1:62618/100‑1 (MQ=255)                         gAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCt                                                                 <  1:62608/100‑1 (MQ=255)                          aaCTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCtt                                                                <  1:1615092/100‑1 (MQ=255)                            cTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTaa                                                              >  2:994244/1‑100 (MQ=255)                             tGTACTGGCAGCAAAGCTATCCGGCGCGGCGGACAGTGGGTTAGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTaa                                                              <  1:1441182/99‑1 (MQ=255)                             tGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAg                                                             >  1:697237/1‑100 (MQ=255)                             tGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAg                                                             <  2:125175/100‑1 (MQ=255)                              gTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAg                                                             >  2:744232/1‑99 (MQ=255)                              gTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGt                                                            <  2:251221/100‑1 (MQ=255)                                aCTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTcc                                                          >  2:728160/1‑100 (MQ=255)                                 cTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCa                                                         >  1:1371526/1‑100 (MQ=255)                                  tGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCAt                                                        <  1:308636/100‑1 (MQ=255)                                   ggCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCAtt                                                       >  2:50450/1‑100 (MQ=255)                                   ggCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCAtt                                                       >  2:259643/1‑100 (MQ=255)                                   ggCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCAtt                                                       <  2:682523/100‑1 (MQ=255)                                   ggCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCAtt                                                       >  1:549457/1‑100 (MQ=255)                                   ggCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCAAtt                                                       >  2:148563/1‑100 (MQ=255)                                     cagcaAAGCGATCAGGCGCGGCGGACAGTGGATTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGc                                                     >  2:1225768/1‑100 (MQ=255)                                       gcaAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCat                                                   >  2:1787538/1‑100 (MQ=255)                                        caAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCata                                                  >  1:328713/1‑100 (MQ=255)                                        caAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCata                                                  >  1:748384/1‑100 (MQ=255)                                        caAAGCGATCCGGCGCGGCGGACAGTGGGTAGGGTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGGCCATTGCata                                                  >  1:1195487/1‑100 (MQ=255)                                         aaaGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCatat                                                 >  2:537053/1‑100 (MQ=255)                                           aGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATgg                                               <  2:284159/100‑1 (MQ=255)                                           aGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATgg                                               <  1:985138/100‑1 (MQ=255)                                           aGCGATCCGGCGCGGCCGACAGTTGGTTGGCTGCGAATGTTAAAACGTCCTCGCAAAGCGCGGTAACCACGGCCTTTACGCTTAAGTCCATTGCATATg                                                <  1:1684467/99‑1 (MQ=255)                                            gCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGa                                              <  1:126350/100‑1 (MQ=255)                                              gATCCGGCTCTGCGGATAGTGGGTTTGTTTCGAATGTTAAAACGTCCGTGCAAAGCGTGGGAACTTCGGTTTTTACGCTTAAGTCCATTACATATgg                                               <  1:1395381/97‑1 (MQ=255)                                              gATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATg                                            >  2:797716/1‑100 (MQ=255)                                              gATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATg                                            >  1:510757/1‑100 (MQ=255)                                              gATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATg                                            >  1:1687345/1‑100 (MQ=255)                                              gATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATg                                            <  2:903238/100‑1 (MQ=255)                                               aTCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATg                                            <  2:1518933/99‑1 (MQ=255)                                               aTCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGt                                           >  2:1238043/1‑100 (MQ=255)                                                tCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTc                                          >  2:1273519/1‑100 (MQ=255)                                                 ccGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTcc                                         >  1:889252/1‑100 (MQ=255)                                                   ggcgcggcgGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCAc                                       >  2:1369987/1‑100 (MQ=255)                                                   ggcgcggcgGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCAc                                       >  1:168995/1‑100 (MQ=255)                                                     cgcggcgGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCACGc                                     <  1:486758/100‑1 (MQ=255)                                                      gcggcgGACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCACGcc                                    <  2:978949/100‑1 (MQ=255)                                                           ggACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCACGCCGGaaa                               >  2:1774415/1‑100 (MQ=255)                                                           ggACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCACGCCGGaaa                               >  2:820359/1‑100 (MQ=255)                                                            gACAGTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGAGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCACGCCGGAAAt                              <  2:878152/100‑1 (MQ=255)                                                               tgTGGGTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCACGCCGGAAATTTa                           >  1:228058/2‑100 (MQ=255)                                                                  gggTTGGCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCACGCCGGAAATTTAGTg                        >  1:885460/1‑100 (MQ=255)                                                                   ggTTGGCTGCAAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGACTTTACGCTTAAGTCCATTGCATATGGATGTCCACGCCGGAAATTTAGTGc                       <  1:95209/100‑1 (MQ=255)                                                                       ggCTGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCACGCCGGAAATTTAGTGCATa                    >  2:845592/1‑99 (MQ=255)                                                                          tGCGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCACGCCGGAAATTTAGTGCATAgcg                 <  1:1581754/99‑1 (MQ=255)                                                                            cGAATGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCACGCCGGAAATTTAGTGCATAGCGCGt              >  2:269311/1‑100 (MQ=255)                                                                              aaTGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCACGCCGGAAATTTAGTGCATAGCGCGTCa            >  2:232608/1‑100 (MQ=255)                                                                               aTGTTAAAACGTCCGCGCAAAGCGCGGGAACCACGGCTTTTACGCTTAAGTCCATTGCATATGGATGTCCACGCCGGAAATTTAGTGCATAGCGCGTCAg           <  2:508288/100‑1 (MQ=255)                                                                                  ntAAAACGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCACGCCGGAAATTTAGTGCATAGCGCGTCAGGGt        >  2:1332380/2‑100 (MQ=255)                                                                                    aaaaCGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCACGCCGGAAATTTAGTGCATAGCGCGTCAGGGTTa      >  2:1263402/1‑100 (MQ=255)                                                                                        cGTCCGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGACGTCCACGCCGGAAATTTAGTGCATAGCGCGTCAGGGTTAAAAc  >  1:1158515/1‑100 (MQ=255)                                                                                            |                                                                                                CGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCTGCGCAAAGCGCGGGAACCACGGCCTTTACGCTTAAGTCCATTGCATATGGATGTCCACGCCGGAAATTTAGTGCATAGCGCGTCAGGGTTAAAAC  >  NC_000913/1163623‑1163809

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑