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,4470AG98.2% 193.9 / ‑7.0 56C12R (TGC→CGC) thiBthiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (30/25);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ACCAGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCACAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTATTTTCTCAGCG  >  NC_000913/75317‑75578
                                                                                                                                  |                                                                                                                                   
accaGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagcaaca                                                                                                                             >  2:1636170/1‑139 (MQ=255)
 ccaGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagcaacag                                                                                                                            >  2:1482196/1‑139 (MQ=255)
         aGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCa                                                                                                                    >  2:1727042/1‑139 (MQ=255)
              gCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACAt                                                                                                               <  2:59022/139‑1 (MQ=255)
              gCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACAt                                                                                                               <  1:294074/139‑1 (MQ=255)
                     cAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCACAGCAACAGAAGGGGCAGACATTTTTTTa                                                                                                        <  2:1199336/139‑1 (MQ=255)
                         cGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTaa                                                                                                       >  2:1503263/1‑136 (MQ=255)
                         cGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTaa                                                                                                       <  1:1503263/136‑1 (MQ=255)
                          ggCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACt                                                                                                   >  2:1758693/1‑139 (MQ=255)
                            cTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACttt                                                                                                 >  2:375845/1‑139 (MQ=255)
                            cTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACttt                                                                                                 >  2:1494655/1‑139 (MQ=255)
                               cAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCa                                                                                              <  1:375845/139‑1 (MQ=255)
                                aaaGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCAc                                                                                             <  1:1636170/139‑1 (MQ=255)
                                       tttttAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACAt                                                                                                               >  2:1389570/1‑114 (MQ=255)
                                       tttttAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACAt                                                                                                               <  1:1389570/114‑1 (MQ=255)
                                       tttttAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCaaaa                                                                                      <  1:1482196/139‑1 (MQ=255)
                                                   ggACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACAt                                                                                                               >  1:1559391/1‑102 (MQ=255)
                                                   ggACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACAt                                                                                                               <  2:1559391/102‑1 (MQ=255)
                                                       cAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTg                                                                                                >  1:1171198/1‑113 (MQ=255)
                                                       cAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTg                                                                                                <  2:1171198/113‑1 (MQ=255)
                                                        aGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACt                                                                     >  2:1531072/1‑139 (MQ=255)
                                                          gCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAGACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTg                                                                   >  2:492332/1‑139 (MQ=255)
                                                                    ggcggcGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACAt                                                                                                               >  2:579554/1‑85 (MQ=255)
                                                                    ggcggcGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACAt                                                                                                               <  1:579554/85‑1 (MQ=255)
                                                                            aGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATcc                                                 >  2:289711/1‑139 (MQ=255)
                                                                                       gTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCg                                      <  1:1531072/139‑1 (MQ=255)
                                                                                          taAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGAc                                                                                                                 <  2:1018826/61‑1 (MQ=255)
                                                                                          taAACAGTCAGAAAGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGAc                                                                                                                 >  1:1018826/1‑61 (MQ=255)
                                                                                            aaCAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATc                                 <  1:289711/139‑1 (MQ=255)
                                                                                                gTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGa                             >  2:288502/1‑139 (MQ=255)
                                                                                                 tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTaa                                                                                                       >  1:373515/1‑64 (MQ=255)
                                                                                                 tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTaa                                                                                                       <  2:373515/64‑1 (MQ=255)
                                                                                                 tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCAc                                                                                             <  1:474906/74‑1 (MQ=255)
                                                                                                 tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCAc                                                                                             >  2:474906/1‑74 (MQ=255)
                                                                                                 tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCg                                      >  1:158490/1‑129 (MQ=255)
                                                                                                 tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCg                                      <  2:158490/129‑1 (MQ=255)
                                                                                                 tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGAt                            >  2:253933/1‑139 (MQ=255)
                                                                                                 tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGAt                            >  2:143082/1‑139 (MQ=255)
                                                                                                 tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGAt                            >  1:743518/1‑139 (MQ=255)
                                                                                                  cAGAACGGGTTTAGCGAAAACGGGCGCTGTGTGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCATCTCAAAAAtga                                                                                  <  1:1040212/84‑1 (MQ=255)
                                                                                                  cAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCa                                                                                              >  1:805812/1‑72 (MQ=255)
                                                                                                  cAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCa                                                                                              <  2:805812/72‑1 (MQ=255)
                                                                                                  cAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAaga                                                                                  >  2:1040212/1‑84 (MQ=255)
                                                                                                  cAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCaa                                                                           <  1:866274/91‑1 (MQ=255)
                                                                                                  cAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCaa                                                                           >  2:866274/1‑91 (MQ=255)
                                                                                                  cAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAgg                                                                        >  2:176204/1‑94 (MQ=255)
                                                                                                  cAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAgg                                                                        <  1:176204/94‑1 (MQ=255)
                                                                                                  cAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGAGTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATc                           >  1:556259/1‑139 (MQ=255)
                                                                                                         ggTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCg                    >  1:1758077/1‑139 (MQ=255)
                                                                                                               gCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACA‑TTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTa             >  1:115626/1‑139 (MQ=255)
                                                                                                                   aaaCGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCa                                                                                              <  1:69632/55‑1 (MQ=255)
                                                                                                                   aaaCGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCa                                                                                              >  2:69632/1‑55 (MQ=255)
                                                                                                                   aaaCGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTAttt          >  1:123958/1‑139 (MQ=255)
                                                                                                                    aaCGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGc                                           <  2:1460681/105‑1 (MQ=255)
                                                                                                                    aaCGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGc                                           >  1:1460681/1‑105 (MQ=255)
                                                                                                                      cGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTATCACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTATTTtct       <  1:1727042/139‑1 (MQ=255)
                                                                                                                           gcTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGTGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTATTTTCTCAgcg  <  2:123958/139‑1 (MQ=255)
                                                                                                                           gcTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTATTTTCTCAgcg  <  1:288502/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                   
ACCAGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCACAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTATTTTCTCAGCG  >  NC_000913/75317‑75578

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: ‑