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_00091375,4470AG100.0% 113.6 / NA 33C12R (TGC→CGC) thiBthiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/15);  total (18/15)

CCACCAGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCACAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTATTTTCT  >  NC_000913/75315‑75573
                                                                                                                                    |                                                                                                                              
ccaccaGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagc                                                                                                                            >  2:268654/1‑137 (MQ=255)
ccaccaGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagc                                                                                                                            <  1:268654/137‑1 (MQ=255)
  accaGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagcaaca                                                                                                                        >  2:675319/1‑139 (MQ=255)
                                 cAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCa                                                                                         <  2:1447102/139‑1 (MQ=255)
                                        ttttttAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCaaa                                                                                  >  2:974930/1‑139 (MQ=255)
                                                   ccGGACCAGGCCCCCAGTCGGCTGCGAAGGAATCGTAGATATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagc                                                                                                                            <  1:975148/86‑1 (MQ=255)
                                                   ccGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagc                                                                                                                            >  2:975148/1‑86 (MQ=255)
                                                    cGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACAt                                                                                                          >  1:1139807/1‑103 (MQ=255)
                                                    cGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACAt                                                                                                          <  2:1139807/103‑1 (MQ=255)
                                                      gACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAg                                                                    >  1:1532672/1‑139 (MQ=255)
                                                            gCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagcaacag                                                                                                                       <  1:1048188/82‑1 (MQ=255)
                                                            gCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagcaacag                                                                                                                       >  2:1048188/1‑82 (MQ=255)
                                                              ccccAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTgag                                                            >  2:293200/1‑139 (MQ=255)
                                                                 cAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAg                                                         <  1:293200/139‑1 (MQ=255)
                                                                           cgAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAac                                                                                                 >  1:206370/1‑89 (MQ=255)
                                                                           cgAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAac                                                                                                 <  2:206370/89‑1 (MQ=255)
                                                                                     gTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGcc                                     >  2:717933/1‑139 (MQ=255)
                                                                                          tataAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAagag                                                                            >  2:713809/1‑95 (MQ=255)
                                                                                          tataAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAagag                                                                            <  1:713809/95‑1 (MQ=255)
                                                                                           ataAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTGAAAAAAGAGTGGCAAAGGACTTGAGAAGGa                                                       <  2:711348/115‑1 (MQ=255)
                                                                                           ataAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGt                                                                           <  1:959327/95‑1 (MQ=255)
                                                                                           ataAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGt                                                                           >  2:959327/1‑95 (MQ=255)
                                                                                           ataAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGa                                                       >  1:711348/1‑115 (MQ=255)
                                                                                            taAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTa                              <  2:1532672/139‑1 (MQ=255)
                                                                                                   tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagc                                                                                                                            <  2:1261083/38‑1 (MQ=255)
                                                                                                   tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagc                                                                                                                            >  1:1261083/1‑38 (MQ=37)
                                                                                                   tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTg                                                                                           <  2:1287741/71‑1 (MQ=255)
                                                                                                   tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTg                                                                                           >  1:1287741/1‑71 (MQ=255)
                                                                                                   tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCaa                                                                      <  1:397571/92‑1 (MQ=255)
                                                                                                   tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCaa                                                                      >  2:397571/1‑92 (MQ=255)
                                                                                                   tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGAt                       >  2:522704/1‑139 (MQ=255)
                                                                                                   tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGAt                       >  2:163917/1‑139 (MQ=255)
                                                                                                    cAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAgg                                                        >  2:523857/1‑105 (MQ=255)
                                                                                                    cAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAgg                                                        <  1:523857/105‑1 (MQ=255)
                                                                                                             tttAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTaa                                                                                                  >  2:342267/1‑54 (MQ=255)
                                                                                                             tttAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTaa                                                                                                  <  1:342267/54‑1 (MQ=255)
                                                                                                                 gCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGAc                                                                 <  1:1374847/83‑1 (MQ=255)
                                                                                                                 gCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGAc                                                                 >  2:1374847/1‑83 (MQ=255)
                                                                                                                 gCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGt         >  2:65579/1‑139 (MQ=255)
                                                                                                                      aaCGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCa                                                                       >  1:747391/1‑72 (MQ=255)
                                                                                                                      aaCGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCa                                                                       <  2:747391/72‑1 (MQ=255)
                                                                                                                        cGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTATTTtct  <  1:717933/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                              
CCACCAGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCACAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTATTTTCT  >  NC_000913/75315‑75573

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