Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,727,756 A→G 46.4% V158V (GTT→GTC purL ← phosphoribosylformyl‑glycineamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,727,7560AG46.4% ‑0.8 / 32.7 41V158V (GTT→GTCpurLphosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base A (13/9);  new base G (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 = 9.59e-01

AGCCAAGACGCAGGTTAGCGTCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCAAAAAAGACCGTTTCCAT  >  minE/1727692‑1727824
                                                                |                                                                    
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                                    <  1:101536/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:200277/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                                    <  1:3368557/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                                    <  1:1638554/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:1340355/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                                    <  1:1312948/66‑1 (MQ=255)
ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc                                                                    <  1:788299/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:3355250/1‑69 (MQ=255)
                    tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa                                              >  1:83341/1‑69 (MQ=255)
                    tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa                                              >  1:2621650/1‑69 (MQ=255)
                    tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa                                              >  1:1580735/1‑69 (MQ=255)
                    tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa                                              >  1:2236198/1‑69 (MQ=255)
                    tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa                                              >  1:2104785/1‑69 (MQ=255)
                    tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa                                              >  1:1233591/1‑69 (MQ=255)
                                          gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg                                      <  1:231309/55‑1 (MQ=255)
                                          gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg                                      <  1:694959/55‑1 (MQ=255)
                                          gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg                                      <  1:1719437/55‑1 (MQ=255)
                                          gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg                                      <  1:1356538/55‑1 (MQ=255)
                                          gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg                                      <  1:2036553/55‑1 (MQ=255)
                                          gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg                                      <  1:2778854/55‑1 (MQ=255)
                                                 cAAATCTACGCTGGTGACCGGAGTCGGTTGATGGTGTGCAAACAACTGCTCTGCATCATCTAAAGCaa                  >  1:3030001/1‑68 (MQ=255)
                                                 cAAATCTACGCTGGTGACCGGAGTCGGTTGATGGTGTGCAAACAACTGCTCTGCATCATCTAAAGCaa                  >  1:3026406/1‑68 (MQ=255)
                                                 cAAATCTACGCTGGTGACCGGAGTCGGTTGATGGTGTGCAAACAACTGCTCTGCATCATCTAAAGCaa                  >  1:1594409/1‑68 (MQ=255)
                                                  aaaTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTg                                 >  1:2578648/1‑52 (MQ=255)
                                                  aaaTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCaa                  >  1:2845279/1‑67 (MQ=255)
                                                  aaaTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCaa                  >  1:2740670/1‑67 (MQ=255)
                                                  aaaTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCaa                  >  1:2254551/1‑67 (MQ=255)
                                                  aaaTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCaa                  >  1:465013/1‑67 (MQ=255)
                                                    aTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCaaaaaa              <  1:558475/69‑1 (MQ=255)
                                                                aaCCGGAGTCGGTTGATGGTGGGCGAACAACTGCTCTGCATCATCTAAAGCAAAAAAGACCGTTTCcat  <  1:1734617/69‑1 (MQ=255)
                                                                aaCCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCAAAAAAGACCGTTTCcat  <  1:1078344/69‑1 (MQ=255)
                                                                aaCCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCAAAAAAGACCGTTTCcat  <  1:501575/69‑1 (MQ=255)
                                                                aaCCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCAAAAAAGACCGTTTCcat  <  1:362829/69‑1 (MQ=255)
                                                                aaCCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCAAAAAAGACCGTTTCcat  <  1:3289500/69‑1 (MQ=255)
                                                                aaCCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCAAAAAAGACCGTTTCcat  <  1:3071167/69‑1 (MQ=255)
                                                                aaCCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCAAAAAAGACCGTTTCcat  <  1:2213500/69‑1 (MQ=255)
                                                                |                                                                    
AGCCAAGACGCAGGTTAGCGTCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCAAAAAAGACCGTTTCCAT  >  minE/1727692‑1727824

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: