Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	75,447	A→G	C12R (TGC→CGC)	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	75,447	0	A	G	100.0%	119.0 / NA	33	C12R (TGC→CGC)	thiB	thiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/12);  total (21/12)

CGCCACCAGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCACAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTATTTTCTCAG  >  NC_000913/75313‑75576                                                                                                                                       |                                                                                                                                  cgcCACCAGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagc                                                                                                                               >  2:732811/1‑139 (MQ=255)                   gCAATTACAGTCGGCTTCAAAGGCTTTTTCAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACAt                                                                                                             <  1:793500/139‑1 (MQ=255)                             tCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAaca                                                                                                   >  2:185395/1‑139 (MQ=255)                                gCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACtt                                                                                                >  2:874024/1‑139 (MQ=255)                                 cTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACttt                                                                                               >  2:185148/1‑139 (MQ=255)                                    cAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCa                                                                                            >  1:1012384/1‑139 (MQ=255)                                      aaGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCAcc                                                                                          <  1:949953/139‑1 (MQ=255)                                                       cGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCaa                                                                         <  2:1012384/139‑1 (MQ=255)                                                         gACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAg                                                                       >  2:166618/1‑139 (MQ=255)                                                                           cggcgAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTc                                                     >  1:926907/1‑139 (MQ=255)                                                                                            gTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCg                                    <  1:732811/139‑1 (MQ=255)                                                                                             tataAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGt                                   >  2:1049865/1‑139 (MQ=255)                                                                                               taAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACAt                                                                                                             >  1:335834/1‑63 (MQ=255)                                                                                               taAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACAt                                                                                                             <  2:335834/63‑1 (MQ=255)                                                                                               taAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTa                                 <  1:1049865/139‑1 (MQ=255)                                                                                                aaaCAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTAt                                >  2:342371/1‑139 (MQ=255)                                                                                                  aCAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATcc                              >  2:653128/1‑139 (MQ=255)                                                                                                    aGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCgg                            >  2:647205/1‑139 (MQ=255)                                                                                                      tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAg                                                         >  1:382010/1‑108 (MQ=255)                                                                                                      tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAg                                                         <  2:382010/108‑1 (MQ=255)                                                                                                      tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGAt                          >  1:94604/1‑139 (MQ=255)                                                                                                      tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGAt                          >  1:38037/1‑139 (MQ=255)                                                                                                      tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGAt                          >  1:740947/1‑139 (MQ=255)                                                                                                      tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGAt                          <  1:655711/139‑1 (MQ=255)                                                                                                      tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGAt                          >  2:317063/1‑139 (MQ=255)                                                                                                        aGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCa                        >  2:514214/1‑139 (MQ=255)                                                                                                           aCGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGt                     <  2:926907/139‑1 (MQ=255)                                                                                                           aCGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGt                     >  1:779339/1‑139 (MQ=255)                                                                                                               gTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGa                 <  1:514214/139‑1 (MQ=255)                                                                                                                tttAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGAc                >  1:363581/1‑139 (MQ=255)                                                                                                                 ttAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACg               <  1:862426/139‑1 (MQ=255)                                                                                                                         aaCGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACg               <  2:459602/131‑1 (MQ=255)                                                                                                                         aaCGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACg               >  1:459602/1‑131 (MQ=255)                                                                                                                              gcgcTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTATTTTCTCAg  >  2:1022827/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                  CGCCACCAGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCACAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTATTTTCTCAG  >  NC_000913/75313‑75576

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