Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	74,829	Δ1 bp	coding (652/984 nt)	thiB ←	thiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	74,829	0	A	.	100.0%	227.2 / NA	48	coding (652/984 nt)	thiB	thiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (19/29);  total (19/29)

ATAAACTGGAGGAATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTT  >  NC_000913/74684‑74966                                                                                                                                                  |                                                                                                                                          aTAAACTGGAGGAATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑ggtgg                                                                                                                                      >  1:585769/1‑150 (MQ=255)    aaCTGGAGGAATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTa                                                                                                                                  <  1:257286/151‑1 (MQ=255)         gaggaATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCa                                                                                                                             <  1:103107/151‑1 (MQ=255)          aggaATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTc                                                                                                                              >  1:386218/1‑149 (MQ=255)           ggaATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGt                                                                                                                           >  1:395614/1‑151 (MQ=255)            gaATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTa                                                                                                                          <  1:836286/151‑1 (MQ=255)               tttttGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTAcc                                                                                                                        >  2:332339/1‑150 (MQ=255)                ttttGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCa                                                                                                                       <  1:500349/150‑1 (MQ=255)                      gcCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCt                                                                                                                >  2:183532/1‑151 (MQ=255)                          gCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTc                                                                                                             >  2:777280/1‑150 (MQ=255)                              cGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCtt                                                                                                        >  1:779035/1‑151 (MQ=255)                               ggCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCttt                                                                                                       <  2:829407/151‑1 (MQ=255)                                   gcttgctGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTaaa                                                                                                   <  1:127134/151‑1 (MQ=255)                                      tgctGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTaaaa                                                                                                  <  2:567202/149‑1 (MQ=255)                                      tgctGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAAc                                                                                                >  2:667713/1‑151 (MQ=255)                                                 gTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGc                                                                                     <  1:468877/151‑1 (MQ=255)                                                   gcgggcggCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCtt                                                                                   >  2:24836/1‑151 (MQ=255)                                                   gcgggcggCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCt                                                                                    <  2:555283/150‑1 (MQ=255)                                                    cgggcggCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTc                                                                                  <  1:523140/151‑1 (MQ=255)                                                                  cACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTgg                                                                    >  2:753773/1‑151 (MQ=255)                                                                      tGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGgtg                                                                  <  1:669591/149‑1 (MQ=255)                                                                       gCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGgtgacc                                                               >  2:779032/1‑151 (MQ=255)                                                                         aGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGgtgaccgt                                                             >  1:658016/1‑151 (MQ=255)                                                                            tAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGgtgaccgtga                                                           <  1:667671/150‑1 (MQ=255)                                                                             aGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGgtgaccgtgacc                                                         >  2:724833/1‑151 (MQ=255)                                                                              gTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGCGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGgtgaccgtgtcc                                                         >  2:817287/1‑150 (MQ=255)                                                                                 aCCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGttt                                                     <  1:694745/151‑1 (MQ=255)                                                                                 aCCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGttt                                                     <  1:825053/151‑1 (MQ=255)                                                                                 aCCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGttt                                                     >  2:757065/1‑151 (MQ=255)                                                                                     tCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTttct                                                  <  2:467903/150‑1 (MQ=255)                                                                                     tCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTttct                                                  <  1:524918/150‑1 (MQ=255)                                                                                      cGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTttctt                                                 <  2:31890/150‑1 (MQ=255)                                                                                       gCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCg                                               <  1:466845/151‑1 (MQ=255)                                                                                        cTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTttcttc                                                >  1:611531/1‑149 (MQ=255)                                                                                                cgcgGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCt                                      <  1:278096/151‑1 (MQ=255)                                                                                                         gtcgTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTACCTCAGTACCCGATCGCTTTCTCCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGAGACCGAGACCATTTTCTTCGCCAGGTTCTGCCAGGCtt                             >  1:209026/4‑151 (MQ=255)                                                                                                            gTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTtgtg                          <  2:179695/151‑1 (MQ=255)                                                                                                               aTCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTggg                        <  2:457690/150‑1 (MQ=255)                                                                                                               aTCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTggg                        <  1:292097/150‑1 (MQ=255)                                                                                                                   ttcttcTCTTCGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCATCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGc                      <  2:175412/148‑1 (MQ=255)                                                                                                                     cttctcttcGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTCAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGTGGCGTCAt                 >  1:61661/1‑151 (MQ=255)                                                                                                                     cttctcttcGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCTCCAGTTTCTGCCAGGCTTGTGGGGCGTc                   >  1:336097/1‑149 (MQ=255)                                                                                                                       tctcttcGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATc                <  1:618753/150‑1 (MQ=255)                                                                                                                             cGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATa          <  2:448563/150‑1 (MQ=255)                                                                                                                             cGAGAATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATa          <  1:531296/150‑1 (MQ=255)                                                                                                                              gagaATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCAt           >  2:221117/1‑148 (MQ=255)                                                                                                                               agaATGTGATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGAc       <  1:49492/151‑1 (MQ=255)                                                                                                                                    gtgATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACttttt  <  1:531351/151‑1 (MQ=255)                                                                                                                                    gtgATAAGCCGGAG‑GGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACtttt   <  2:262398/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                          ATAAACTGGAGGAATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTT  >  NC_000913/74684‑74966

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