Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,727,747 A→T 46.3% V161V (GTT→GTA purL ← phosphoribosylformyl‑glycineamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,727,7470AT46.3% 3.6 / 43.4 41V161V (GTT→GTApurLphosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base A (13/9);  new base T (3/16);  total (16/25)
Fisher's exact test for biased strand distribution p-value = 9.25e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

AGCCAAGACGCAGGTTAGCGTCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCAAAAAA  >  minE/1727692‑1727812
                                                       |                                                                 
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                        <  1:1638554/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                        <  1:2687414/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                        <  1:2881304/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                        <  1:3073610/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                        <  1:3368557/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                        <  1:200277/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                        <  1:395604/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                        <  1:14080/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                        <  1:1398471/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                        <  1:788299/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                        <  1:1340355/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                        <  1:1312948/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                        <  1:101536/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCCGCCCCAGCAAATCTACGCTGGTGAc                                                        <  1:3005027/66‑1 (MQ=255)
 gccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                        <  1:3208896/66‑1 (MQ=255)
         gCAGGTTAGCGTCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGtt                                             <  1:3388933/69‑1 (MQ=255)
          cAGGTTAGCGTCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGtt                                             <  1:2194700/68‑1 (MQ=255)
           aGGTTAGCGTCAATTAGAGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                        <  1:1227277/56‑1 (MQ=255)
                    tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa                                  >  1:772448/1‑69 (MQ=255)
                    tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa                                  >  1:83341/1‑69 (MQ=255)
                    tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa                                  >  1:3355250/1‑69 (MQ=255)
                    tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa                                  >  1:2104785/1‑69 (MQ=255)
                    tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa                                  >  1:1233591/1‑69 (MQ=255)
                    tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa                                  >  1:1580735/1‑69 (MQ=255)
                    tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa                                  >  1:2621650/1‑69 (MQ=255)
                    tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa                                  >  1:2236198/1‑69 (MQ=255)
                                          gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg                          <  1:231309/55‑1 (MQ=255)
                                          gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg                          <  1:2036553/55‑1 (MQ=255)
                                          gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg                          <  1:1719437/55‑1 (MQ=255)
                                          gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg                          <  1:694959/55‑1 (MQ=255)
                                          gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg                          <  1:2778854/55‑1 (MQ=255)
                                          gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg                          <  1:1356538/55‑1 (MQ=255)
                                                 cAAATCTACGCTGGTGACCGGAGTCGGTTGATGGTGTGCAAACAACTGCTCTGCATCATCTAAAGCaa      >  1:3026406/1‑68 (MQ=255)
                                                 cAAATCTACGCTGGTGACCGGAGTCGGTTGATGGTGTGCAAACAACTGCTCTGCATCATCTAAAGCaa      >  1:3030001/1‑68 (MQ=255)
                                                 cAAATCTACGCTGGTGACCGGAGTCGGTTGATGGTGTGCAAACAACTGCTCTGCATCATCTAAAGCaa      >  1:1594409/1‑68 (MQ=255)
                                                  aaaTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTg                     >  1:2578648/1‑52 (MQ=255)
                                                  aaaTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCaa      >  1:2254551/1‑67 (MQ=255)
                                                  aaaTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCaa      >  1:2845279/1‑67 (MQ=255)
                                                  aaaTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCaa      >  1:465013/1‑67 (MQ=255)
                                                  aaaTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCaa      >  1:2740670/1‑67 (MQ=255)
                                                    aTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCaaaaaa  <  1:558475/69‑1 (MQ=255)
                                                       |                                                                 
AGCCAAGACGCAGGTTAGCGTCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCAAAAAA  >  minE/1727692‑1727812

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