Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	435,879	T→A	D82E (GAT→GAA)	thiL →	thiamine monophosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	435,879	0	T	A	100.0%	130.9 / NA	36	D82E (GAT→GAA)	thiL	thiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (16/20);  total (16/20)

CCTGATATCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGATCCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAGTCCTTCAGCGACAGTTTGT  >  NC_000913/435802‑435955                                                                              |                                                                             ccTGATATCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAAcc                                                                            <  1:92653/80‑1 (MQ=255) ccTGATATCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCg                                                                           <  1:92655/81‑1 (MQ=255)    gATATCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGcc                                                                        >  1:92691/1‑81 (MQ=255)     atatCGATCCTGCTGTTCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGcc                                                                        <  2:92658/80‑1 (MQ=255)     atatCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGcc                                                                        <  2:92657/80‑1 (MQ=255)     atatCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCt                                                                       >  1:92704/1‑81 (MQ=255)     atatCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCt                                                                       <  1:92656/81‑1 (MQ=255)       atCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTgg                                                                     >  2:92698/1‑81 (MQ=255)        tCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTgg                                                                     <  2:92659/80‑1 (MQ=255)          gATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTGTCt                                                                   >  2:92703/1‑80 (MQ=255)          gATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTGGCt                                                                   >  2:92684/1‑80 (MQ=255)          gATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTGGCt                                                                   <  2:92660/80‑1 (MQ=255)           aTCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTGGCTGa                                                                 >  2:92692/1‑81 (MQ=255)            tCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGTACCGGCCTGGCTGAc                                                                <  1:92661/81‑1 (MQ=255)                 ctgATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTGGCTGACGCt                                                             <  1:92662/79‑1 (MQ=255)                               aaaGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCg                                             >  1:92715/1‑81 (MQ=255)                                aaGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCg                                             <  2:92663/80‑1 (MQ=255)                                  gCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGAc                                          >  2:92696/1‑81 (MQ=255)                                    aCTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGt                                        <  1:92664/81‑1 (MQ=255)                                        gCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAg                                      <  1:92665/79‑1 (MQ=255)                                          ggTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACg                                   >  1:92701/1‑80 (MQ=255)                                             gAACCTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGc                               >  2:92718/1‑81 (MQ=255)                                                 cTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTgg                           >  1:92711/1‑81 (MQ=255)                                                 cTAAGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTgg                           >  1:92702/1‑81 (MQ=255)                                                   aaGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCt                         >  1:92726/1‑81 (MQ=255)                                                    aGCGATCTGGCAGCGATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCtt                        <  2:92666/81‑1 (MQ=255)                                                       gATCTGGCAGCGATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAg                     >  2:92713/1‑81 (MQ=255)                                                             gCAGCGATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAg                     <  2:92667/75‑1 (MQ=255)                                                             gCAGCGATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAGTcc                  <  2:92669/78‑1 (MQ=255)                                                             gCAGCGATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAGTCCt                 <  1:92668/79‑1 (MQ=255)                                                                 cGATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAGTCCTTCAGCg           <  1:92670/81‑1 (MQ=255)                                                                  gATGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAGTCCTTCAGCGa          >  2:92700/1‑81 (MQ=255)                                                                   aTGGGGGCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAGTCCTTCAGCGa          >  2:92716/1‑80 (MQ=255)                                                                       gggCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAGTCCTTCAGCGACAgtt     <  2:92671/81‑1 (MQ=255)                                                                       gggCCGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAGTCCTTCAGCGACAgtt     <  1:92672/81‑1 (MQ=255)                                                                          ccGAACCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAGTCCTTCAGCGACAgtttgt  <  1:92673/81‑1 (MQ=255)                                                                              |                                                                             CCTGATATCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGATCCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAGTCCTTCAGCGACAGTTTGT  >  NC_000913/435802‑435955

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