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,4470AG96.4% 190.5 / ‑2.6 56C12R (TGC→CGC) thiBthiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (28/26);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

CACCAGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCACAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTATTTT  >  NC_000913/75316‑75571
                                                                                                                                   |                                                                                                                            
caccaGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagcaac                                                                                                                       >  2:816427/1‑139 (MQ=255)
   caGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCgcagcaacagc                                                                                                                    <  2:280572/139‑1 (MQ=255)
             tCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGAc                                                                                                          <  1:763877/139‑1 (MQ=255)
               gCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACAt                                                                                                        <  2:1154506/139‑1 (MQ=255)
                      cAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTa                                                                                                 <  2:658116/139‑1 (MQ=255)
                         tCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAaca                                                                                              <  2:118421/139‑1 (MQ=255)
                                   aGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCt                                                                                    >  1:1278704/1‑139 (MQ=255)
                                               ccaccGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTg                                                                        >  1:1059393/1‑139 (MQ=255)
                                                  ccGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGNAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCa                                                                     >  1:555646/1‑139 (MQ=255)
                                                  ccGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCa                                                                     <  2:1059393/139‑1 (MQ=255)
                                                   cGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCaa                                                                    >  1:1013336/1‑139 (MQ=255)
                                                   cGGACCAGGCCCCCAGTCGGCGGCGAAGAAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCaa                                                                    <  2:804607/139‑1 (MQ=255)
                                                      aCCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAgg                                                                 <  1:946517/139‑1 (MQ=255)
                                                           gCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTg                                                            >  2:454380/1‑139 (MQ=255)
                                                                     ggcggcGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTa                                                                                                 >  1:177289/1‑92 (MQ=255)
                                                                     ggcggcGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTa                                                                                                 <  2:177289/92‑1 (MQ=255)
                                                                             aGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTaa                                                                                                <  2:153492/85‑1 (MQ=255)
                                                                             aGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTaa                                                                                                >  1:153492/1‑85 (MQ=255)
                                                                                      aGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTATGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCgg                                 <  1:454380/139‑1 (MQ=255)
                                                                                       ggTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGAc                                                                                                          >  1:515925/1‑65 (MQ=255)
                                                                                       ggTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGAc                                                                                                          <  2:515925/65‑1 (MQ=255)
                                                                                       ggTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAATGAGCCTCAAATCCCTTCGCCGGc                                >  1:1268209/1‑139 (MQ=255)
                                                                                        gTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCa                                                                                       <  1:1068614/83‑1 (MQ=255)
                                                                                        gTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCa                                                                                       >  2:1068614/1‑83 (MQ=255)
                                                                                        gTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCg                               <  2:1312865/139‑1 (MQ=255)
                                                                                          ataAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCa                                                                                       <  2:323155/81‑1 (MQ=255)
                                                                                          ataAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCa                                                                                       >  1:323155/1‑81 (MQ=255)
                                                                                          ataAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGtt                             >  2:246435/1‑139 (MQ=255)
                                                                                           taAACAGTCAGAACGGGTTTAGCGAAAACGGACGCTGTGCGCAGCAACAGCAGTGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAAGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTa                            >  1:547983/1‑139 (MQ=255)
                                                                                                  tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCaaa                                                                                >  1:1071948/1‑80 (MQ=255)
                                                                                                  tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCaaa                                                                                <  2:1071948/80‑1 (MQ=255)
                                                                                                  tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCg                               <  2:414269/129‑1 (MQ=255)
                                                                                                  tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCg                               >  1:414269/1‑129 (MQ=255)
                                                                                                  tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGAt                     >  2:1240401/1‑139 (MQ=255)
                                                                                                  tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGAt                     >  2:907253/1‑139 (MQ=255)
                                                                                                  tCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGAt                     >  1:902101/1‑139 (MQ=255)
                                                                                                    aGAACGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCa                   <  2:749812/139‑1 (MQ=255)
                                                                                                      aaCGGGTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAgg                 >  1:1272924/1‑139 (MQ=255)
                                                                                                          ggTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAg                                                                  <  2:33469/86‑1 (MQ=255)
                                                                                                          ggTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAg                                                                  >  1:33469/1‑86 (MQ=255)
                                                                                                          ggTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCg             <  1:1270976/139‑1 (MQ=255)
                                                                                                           gTTTAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAATAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGa            <  1:816427/139‑1 (MQ=255)
                                                                                                            tttAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTa                                                                                                 >  1:744817/1‑53 (MQ=255)
                                                                                                            tttAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTa                                                                                                 <  2:744817/53‑1 (MQ=255)
                                                                                                            tttAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATc                          <  1:64273/124‑1 (MQ=255)
                                                                                                            tttAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATc                          >  2:64273/1‑124 (MQ=255)
                                                                                                            tttAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGAc           >  1:941539/1‑139 (MQ=255)
                                                                                                              tAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGAc                                                               >  1:586834/1‑85 (MQ=255)
                                                                                                              tAGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGAc                                                               <  2:586834/85‑1 (MQ=255)
                                                                                                               aGCGAAAACGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACggg        <  1:1228630/139‑1 (MQ=255)
                                                                                                                     aaCGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGc                                    >  1:668749/1‑105 (MQ=255)
                                                                                                                     aaCGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGc                                    <  2:668749/105‑1 (MQ=255)
                                                                                                                     aaCGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTAtttt  >  1:115668/1‑139 (MQ=255)
                                                                                                                     aaCGGGCGCTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTAtttt  >  1:725927/1‑139 (MQ=255)
                                                                                                                          gcgcTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTg                                                                                         <  2:422909/47‑1 (MQ=255)
                                                                                                                          gcgcTGTGCGCAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTg                                                                                         >  1:422909/1‑47 (MQ=255)
                                                                                                                                 gCACAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATc                          <  1:32483/103‑1 (MQ=255)
                                                                                                                                 gCACAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATc                          >  2:32483/1‑103 (MQ=255)
                                                                                                                                   |                                                                                                                            
CACCAGTTTCAGTTCGCAATTACAGTCGGCTTCAAAGGCTTTTTTAACCACCGGACCAGGCCCCCAGTCGGCGGCGAAGGAATCGTAGGTATAAACAGTCAGAACGGGTTTAGCGAAAACGGGCGCTGTGCACAGCAACAGCAGGGGCAGACATTTTTTTAACACTTTGCACCTCAAAAAAGAGTGGCAAAGGACTTGAGAAGGAGCCTCAAATCCCTTCGCCGGCGTTATCCGGATCAGGTTCGACGGGTATTTT  >  NC_000913/75316‑75571

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