Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	436,026	G→A	T131T (ACG→ACA)	thiL →	thiamine monophosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	436,026	0	G	A	100.0%	207.9 / NA	60	T131T (ACG→ACA)	thiL	thiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/34);  total (26/34)

GCTTGAGTCCTTCAGCGACAGTTTGTTTGATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACGTTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGACCGGTACACCGGGC  >  NC_000913/435930‑436119                                                                                                 |                                                                                              gCTTGAGTCCTTCAGCGACAGTTTGTTTGATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTg                                                                                            <  2:1020363/100‑1 (MQ=255)   ttGAGTCCTTCAGCGACAGTTTGTTTGATCTTCTTAATTATTACGATATGCAACCCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTggg                                                                                          <  2:1542229/100‑1 (MQ=255)   ttGAGTCCTTCAGCGACAGTTTGTTTGATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTggg                                                                                          <  1:893421/100‑1 (MQ=255)        tCCTTCAGCGACAGTTTGTTTGATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATcc                                                                                     <  1:1337181/100‑1 (MQ=255)         ccTTCAGCGACAGTTTGTTTGATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCa                                                                                    >  2:533132/1‑100 (MQ=255)          cTTCAGCGACAGTTTGTTTGATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCAc                                                                                   >  2:174869/1‑100 (MQ=255)             cAGCGACAGTTTGTTTGATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGc                                                                                <  2:225548/100‑1 (MQ=255)              aGCGACAGTTTGTTTGATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCt                                                                               >  1:1655953/1‑100 (MQ=255)                   cAGTTTGTTTGATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGtt                                                                          >  2:322796/1‑100 (MQ=255)                     gtttgtttGATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCTCGCGTGAGCCATTATCAATGACATTGGGTATACACGGATTTGtt                                                                          >  2:560876/1‑98 (MQ=255)                     gtttgtttGATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTcc                                                                        <  2:613520/100‑1 (MQ=255)                     gtttgtttGATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTcc                                                                        <  1:655860/100‑1 (MQ=255)                         gtttgATCTTCTCAATTATTACGATATGCAACTCAGTGGCGGCGATACCACGCGTGGGCCATTATCAATGACAGTGGGTATCCACGGCTTTGTTCCGATg                                                                    <  2:984274/100‑1 (MQ=255)                          tttgATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATg                                                                    <  2:277408/99‑1 (MQ=255)                          tgtgATCTATTCAATTATTACGATATGTATTTTATTAGCGGCGATACCACGCGTGGACGATTATCATTGACATTGTGTATCCACGTCATTTTTCCAATgg                                                                   <  1:121709/98‑1 (MQ=255)                            tgATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGa                                                                 <  1:836200/100‑1 (MQ=255)                                  tctcAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGcc                                                           >  1:1223734/1‑100 (MQ=255)                                  tctcAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGcc                                                           <  1:1414130/100‑1 (MQ=255)                                   ctcAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCt                                                          >  2:1506051/1‑100 (MQ=255)                                     cAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTa                                                        >  2:448761/1‑100 (MQ=255)                                     cAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTa                                                        >  1:746434/1‑100 (MQ=255)                                        ttatttCGTTATGCAATTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAAcg                                                     <  2:1121598/100‑1 (MQ=255)                                        ttattaCGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAAcg                                                     <  1:533128/100‑1 (MQ=255)                                        ttattaCGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAAcg                                                     <  2:1241239/100‑1 (MQ=255)                                         tattaCGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAAcgc                                                    >  1:1654667/1‑100 (MQ=255)                                          attaCGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAAcgcg                                                   >  2:226355/1‑100 (MQ=255)                                            taCGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGAGct                                                 >  1:1748966/1‑100 (MQ=255)                                                 tatGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTggg                                            <  1:1016064/100‑1 (MQ=255)                                                    gCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCg                                         <  1:322794/100‑1 (MQ=255)                                                     cAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGa                                        >  2:1398121/1‑100 (MQ=255)                                                       aCTCATTGGTGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGaaa                                      >  1:931683/1‑100 (MQ=255)                                                       aCTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGaaa                                      >  2:797327/1‑100 (MQ=255)                                                        cTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAAc                                     <  1:1389629/100‑1 (MQ=255)                                                        cTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAAc                                     >  1:1585192/1‑100 (MQ=255)                                                          cATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCg                                   <  2:1342747/100‑1 (MQ=255)                                                           aTTGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCgg                                  <  2:1654678/100‑1 (MQ=255)                                                             tGGCGGCGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGt                                <  2:820456/100‑1 (MQ=255)                                                              ggcggcGATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTg                               >  1:454009/1‑100 (MQ=255)                                                                   cgATACCACGCGTGGGCTATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTgg                          <  1:140563/100‑1 (MQ=255)                                                                    gATACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGa                         <  1:746856/100‑1 (MQ=255)                                                                     aTACCACGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGAt                        >  2:650598/1‑100 (MQ=255)                                                                          aCGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTAtg                   <  2:1081108/100‑1 (MQ=255)                                                                          aCGCGTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTAtg                   <  2:746438/100‑1 (MQ=255)                                                                           cgtgTGGGACATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTTACTGGATCTAtg                   <  2:1194809/96‑1 (MQ=255)                                                                            gcgTGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGTGCTCTGGGGCGAAACCGGGTGACTGGATCTAtgtg                 >  2:1053975/1‑100 (MQ=255)                                                                             cgTGGGCCATTATCAATGACATTGGGTATCAACGGCTTTGTTCCGATAGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGa                >  1:581617/1‑100 (MQ=255)                                                                               tGGGCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGAcc              <  1:307840/100‑1 (MQ=255)                                                                                gggCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGACCg             <  2:1424590/100‑1 (MQ=255)                                                                                 ggCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGACCgg            <  2:1382598/100‑1 (MQ=255)                                                                                 ggCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGACCgg            >  1:823847/1‑100 (MQ=255)                                                                                  gCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGACCGGt           >  2:316136/1‑100 (MQ=255)                                                                                  gCCATTATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGACCGGt           >  2:1131247/1‑100 (MQ=255)                                                                                       tATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGACCGGTACAcc      <  2:1318204/100‑1 (MQ=255)                                                                                       tATCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGACCGGTACAcc      <  2:1300596/100‑1 (MQ=255)                                                                                        aTCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGACCGGTACACCg     >  1:1218372/1‑100 (MQ=255)                                                                                        aTCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGACCGGTACACCg     <  1:1702599/100‑1 (MQ=255)                                                                                         tCAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGACCGGTACACCgg    <  1:1683695/100‑1 (MQ=255)                                                                                          cAATGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGACCGGTACACCggg   >  2:723223/1‑100 (MQ=255)                                                                                           aaTGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGACCGGTACACCGGGc  >  2:431937/1‑100 (MQ=255)                                                                                           aaTGACATTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGACCGGTACACCGGGc  <  2:282753/100‑1 (MQ=255)                                                                                                 |                                                                                              GCTTGAGTCCTTCAGCGACAGTTTGTTTGATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACGTTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGCGCTCTGGGGCGAAACCGGGTGACTGGATCTATGTGACCGGTACACCGGGC  >  NC_000913/435930‑436119

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: ‑