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%	97.7 / NA	29	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 (15/14);  total (15/14)

CCACCAGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCACAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTATTTT  >  NC_000913/75315‑75571                                                                                                                                     |                                                                                                                             ccaccaGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagcaa                                                                                                                        >  2:200901/1‑139 (MQ=255)   accaGTTTCAGTTCGCAATTACGGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagcaaca                                                                                                                      >  2:768372/1‑139 (MQ=255)              ttCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGCTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGa                                                                                                           >  2:769425/1‑139 (MQ=255)                                  cAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCa                                                                                       <  2:673206/139‑1 (MQ=255)                                  cAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCa                                                                                       <  1:515796/139‑1 (MQ=255)                                  cAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCa                                                                                       <  2:745129/139‑1 (MQ=255)                                    aaGGCTTTTTTAACCACCGGACCAGGCACCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCAcc                                                                                     <  2:696080/139‑1 (MQ=255)                                          tttttAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCaaaa                                                                               >  2:705818/1‑139 (MQ=255)                                           ttttAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCaaaaa                                                                              >  1:1062461/1‑139 (MQ=255)                                                     cGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagc                                                                                                                          <  2:964838/85‑1 (MQ=255)                                                     cGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagc                                                                                                                          >  1:964838/1‑85 (MQ=255)                                                     cGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCaa                                                                    >  2:654147/1‑139 (MQ=255)                                                     cGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCaa                                                                    <  2:1062461/139‑1 (MQ=255)                                                       gACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAg                                                                  <  1:200901/139‑1 (MQ=255)                                                        aCCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAgg                                                                 <  1:768372/139‑1 (MQ=255)                                                                                      gTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGcc                                   >  1:916728/1‑139 (MQ=255)                                                                                          gTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTATGCGCAGCAACAGCAGGGGCAGACATCTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAGTCCCTTCGCCGGCg                               <  1:654147/139‑1 (MQ=255)                                                                                            ataAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTa                                                                                                 <  1:265808/71‑1 (MQ=255)                                                                                            ataAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTa                                                                                                 >  2:265808/1‑71 (MQ=255)                                                                                                  aGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCaa                                                                    <  1:819804/94‑1 (MQ=255)                                                                                                  aGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCaa                                                                    >  2:819804/1‑94 (MQ=255)                                                                                                  aGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCgg                       >  1:154575/1‑139 (MQ=255)                                                                                                   gTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCa                                                                     <  1:1060800/92‑1 (MQ=255)                                                                                                   gTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCa                                                                     >  2:1060800/1‑92 (MQ=255)                                                                                                    tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagc                                                                                                                          >  2:560168/1‑38 (MQ=37)                                                                                                    tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagc                                                                                                                          <  1:560168/38‑1 (MQ=37)                                                                                                    tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACAt                                                                                                        >  1:992990/1‑56 (MQ=255)                                                                                                    tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACAt                                                                                                        <  2:992990/56‑1 (MQ=255)                                                                                                    tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCa                                                                                       >  2:650335/1‑73 (MQ=255)                                                                                                    tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCa                                                                                       <  1:650335/73‑1 (MQ=255)                                                                                                     cAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATc                    >  1:1015005/1‑139 (MQ=255)                                                                                                         aCGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCg                               <  2:264664/124‑1 (MQ=255)                                                                                                         aCGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCg                               >  1:264664/1‑124 (MQ=255)                                                                                                                       aaCGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTAtttt  >  2:578302/1‑139 (MQ=255)                                                                                                                       aaCGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTAtttt  >  2:721012/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                             CCACCAGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCACAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTATTTT  >  NC_000913/75315‑75571

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