Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,193,996 T→G D70A (GAT→GCT)  thiE ← thiamine phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,193,9960TG93.5% 93.9 / ‑3.4 33D70A (GAT→GCT) thiEthiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base G (18/11);  minor base A (2/0);  total (21/12)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCT  >  NC_000913/4193760‑4194234
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                
gCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAAt                                                                                                                                                                                                                                    >  1:177113/1‑251 (MQ=255)
          ccccTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTtat                                                                                                                                                                                                                          <  1:62833/251‑1 (MQ=255)
                    gCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGcc                                                                                                                                                                                                                >  1:123785/1‑251 (MQ=255)
                      agcaGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACATCGGTTTCCGTATCGTCATGTGTCGAAACGCCCAGCCGCGGGCCGGCCGCGCGGATGGCATTGCGATCGGCGGCTTGCACATCTTCCTGCCCCACATGGACGCCCGACGCCTGATGCTTGATCGCCAGCCACCTGCAAGCGGTGAGAAACAATCGCGCGGTAGAgcgccgtcca                                                                                                                                                                                                               >  2:5699/4‑243 (MQ=255)
                       gcaGGCATCTGTTTGGATTGCGTCGGGAACACATGTCCCAGCGAGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCCTGTGACGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCCTTGAGATCGGTGGCTTGCAACTCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCACGCCGCCCGTAAACGTTGATAAACAATCCCACGTTATAacggagcacc                                                                                                                                                                                                               >  2:247937/3‑243 (MQ=255)
                                           cGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACAt                                                                                                                                                                                         <  1:313845/251‑1 (MQ=255)
                                                 gAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCtt                                                                                                                                                                                   >  1:173777/1‑251 (MQ=255)
                                                 gAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCtt                                                                                                                                                                                   >  1:174282/1‑251 (MQ=255)
                                                                                     tgcCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCtt                                                                                                                                                                                   >  1:164351/1‑215 (MQ=255)
                                                                                      gcCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCTGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCCGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTg                                                                                                                                              >  1:126657/1‑251 (MQ=255)
                                                                                           cgcgACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGACGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCCTACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCAACTCTTc                                                                                                                                                                         >  1:115390/1‑219 (MQ=255)
                                                                                                          cATATCGTCATGTGTCGAAACGCCCAGCCGCATGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCCAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGCTACACAATcgcgcg                                                                                                                                                                                                                              >  2:125397/1‑151 (MQ=255)
                                                                                                          cATATCGTAATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATcgcgcg                                                                                                                                                                                                                              <  1:125397/151‑1 (MQ=255)
                                                                                                           atatCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATc                                                                                                                         >  1:263292/1‑251 (MQ=255)
                                                                                                                         cGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCg                                                                                                            >  1:186959/1‑250 (MQ=255)
                                                                                                                          gAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAgcggcg                                                                                                                                                                                                                   >  1:218573/1‑146 (MQ=255)
                                                                                                                                  cAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGt                                                                                                   >  1:157598/1‑250 (MQ=255)
                                                                                                                                   aGCCGCAGGCCTGCCGCGCGGCTGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCAGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCCTCATCGCGCCGATCTTTTATTGGTAGATGGAGAGTACGGACGCCTGCATCCAACAggcggtccataccatgtcc                                                                                                 >  2:288956/1‑233 (MQ=255)
                                                                                                                                           gCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTcca                                                                                         <  1:91089/251‑1 (MQ=255)
                                                                                                                                                ccGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTACACCTCTTCATCGCGCCGATCTTTGATGCGTCACTGGTGAGGACGGACGCCTGCATCCAACAGACGTTCGATACCCAGTAAGCTGTcccccccc                                                                                    >  2:263134/1‑245 (MQ=255)
                                                                                                                                                                          gCTTGCAAATCTTCCCGCCCCAAATGGACGCCATACGCCTGATGATATATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTAc                                                         <  1:1910/251‑1 (MQ=255)
                                                                                                                                                                                 aaccTTCCTGCCCCAAATGGACGCCATACGCCGTCTCCTGGAACGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAggagga                                                  <  1:167895/248‑1 (MQ=255)
                                                                                                                                                                                           ccccAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACAACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCt                                       >  1:123713/1‑251 (MQ=255)
                                                                                                                                                                                            cccaaagggaccacactcGCCTGATGCTTGATCGCCAGCCGCCAGTAACCGTTGATAAACAACCGCGCGTTATAGCGGCGGCCCGGCGCAATTGCCGCCACGACATCGGCTTCCCCCGCTTCAGCGCTCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTg                                      <  2:44861/234‑1 (MQ=255)
                                                                                                                                                                                             ccAAATGGCCGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGACGCCAAGACATCGGCTTCCACCTCTTCATcgc                                                                                                                                                                    >  2:205210/1‑126 (MQ=255)
                                                                                                                                                                                             ccAAATGGACGCCAGACGCCTTATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATATCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATcgc                                                                                                                                                                    <  1:205210/126‑1 (MQ=255)
                                                                                                                                                                                                          aTACGCTTGATGCTTGATGCCCAGCTGCCAGTAAGCGTTGATAAACAATGGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCGCTTCATCGCGCCGATCTTTGATGCGGAGCTGGAGAGTACGTACGCCTGCAGCCAACAGGCGTTCGATCCACGGTACGCTGTCCACCTCCGGTTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTtc                        <  1:18576/251‑1 (MQ=255)
                                                                                                                                                                                                            ccgctagatgttggagagCCACCCGCCAGTAAGCGGGGAGAAACAATCGCGCGTTATAGCGGGGGCCAGGCGCAATTGCCGCCACGACTGCGGCTTCCACCTCTTCATCGCGCCGATCGTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGTTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTtcct                      <  2:318872/234‑1 (MQ=255)
                                                                                                                                                                                                            aCGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCAc                             >  1:9072/1‑244 (MQ=255)
                                                                                                                                                                                                                   aTGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTGCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACg               <  1:155220/251‑1 (MQ=255)
                                                                                                                                                                                                                         gATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTa         <  1:290543/251‑1 (MQ=255)
                                                                                                                                                                                                                          aTCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCCATTGCCGCCACGACATCGGCTTACACCTCTTCATCGCGCCGACCTTTGATGCGTAGCTGGAGCGTACGGAAGCCTGCATCCCACAGACGTTAGATCCACTTTACGCTGTCCCACACCGGGTACAGTCCTGAACGAAACAGGACAGGAGGAAAATACAGCTGCAACATCACGCTTCATCCTTACgacggg          >  2:230155/1‑244 (MQ=255)
                                                                                                                                                                                                                                aGCCGCCAGAAAACGTTGAAAAACAATCGCGCGGTATAGCGCCGGCCCAGCGCACTTGCCGCCACGACCTCGGCTTACACCTCTGCATCGCCCCGATCTTTGATGCACAGCTGGAGGGTACGTACGCCTGCATCCAACAGACGTTCGATCCCATGTACGATGTCCACCACCGGGTACAGGCCTAAACGACAAAGTCCAGGAGGAAAATCAGGCTGGTCCTGCACGCTTCCCCCTTCAGCAGTTAGATTtct  >  2:249897/1‑251 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                
GCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCT  >  NC_000913/4193760‑4194234

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: