Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,987,528 T→C 26.0% K65K (AAA→AAG ansB ← periplasmic L‑asparaginase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,987,5280TC26.0% 103.6 / 35.6 92K65K (AAA→AAGansBperiplasmic L‑asparaginase II
Reads supporting (aligned to +/- strand):  ref base T (48/20);  new base C (22/2);  total (70/22)
Fisher's exact test for biased strand distribution p-value = 5.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

CAGACATTATCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACTGTGTAGTTAGATTTG  >  minE/1987460‑1987594
                                                                    |                                                                  
cAGACATTATCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCt                                                                    <  1:1245805/69‑1 (MQ=255)
     aTTATCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAg                                                               >  1:2434283/1‑69 (MQ=255)
     aTTATCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAg                                                               >  1:843496/1‑69 (MQ=255)
     aTTATCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAg                                                               >  1:3204337/1‑69 (MQ=255)
     aTTATCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAg                                                               >  1:1532010/1‑69 (MQ=255)
      ttATCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGt                                                              <  1:2665741/69‑1 (MQ=255)
      ttATCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGt                                                              <  1:2853530/69‑1 (MQ=255)
      ttATCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGt                                                              <  1:1486874/69‑1 (MQ=255)
        aTCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCtttt                                                                 >  1:408731/1‑64 (MQ=255)
        aTCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGtt                                                             >  1:2620362/1‑68 (MQ=255)
        aTCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTg                                                            >  1:996662/1‑69 (MQ=255)
        aTCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTg                                                            >  1:2215833/1‑69 (MQ=255)
        aTCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTg                                                            >  1:2139675/1‑69 (MQ=255)
        aTCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTg                                                            >  1:1886598/1‑69 (MQ=255)
        aTCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTg                                                            >  1:3208390/1‑69 (MQ=255)
        aTCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTg                                                            >  1:1123337/1‑69 (MQ=255)
        aTCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTg                                                            >  1:507910/1‑69 (MQ=255)
        aTCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGAAATGTCTTTTAGTTg                                                            >  1:1208726/1‑69 (MQ=255)
              cATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCAc                                                      <  1:3050328/69‑1 (MQ=255)
              cATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCAc                                                      <  1:2453616/69‑1 (MQ=255)
                           ccGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGa                                         <  1:3395006/69‑1 (MQ=255)
                                atTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTc                                    <  1:156288/69‑1 (MQ=255)
                                atTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTc                                    <  1:222840/69‑1 (MQ=255)
                                atTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTc                                    <  1:3049464/69‑1 (MQ=255)
                                atTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTc                                    <  1:3337327/69‑1 (MQ=255)
                                 ttCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTc                                    <  1:1440053/68‑1 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCtttt                                                                 <  1:244699/37‑1 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCtttt                                                                 <  1:1988819/37‑1 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTa                                                                >  1:2689205/1‑38 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:2644251/1‑69 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:680615/1‑69 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:612432/1‑69 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:2818299/1‑69 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:366330/1‑69 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:2872762/1‑69 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:3206587/1‑69 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:1105818/1‑69 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:2399513/1‑69 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:1448924/1‑69 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:1274082/1‑69 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:1872285/1‑69 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:1719500/1‑69 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:1944416/1‑69 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:1269966/1‑69 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:234892/1‑69 (MQ=255)
                                   cACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:1597039/1‑69 (MQ=255)
                                   cACTACCAGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTAc                                 >  1:3133770/1‑69 (MQ=255)
                                     cTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTc                                    <  1:1272214/64‑1 (MQ=255)
                                                    tAACGTTCGCAATGTCTTTTAGTTGCGGCAGCGCATTAACCAGATTTTCTACGCCAACTTTACCCACtg                >  1:345563/1‑69 (MQ=255)
                                                    tAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACtg                >  1:2683782/1‑69 (MQ=255)
                                                    tAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACtg                >  1:1740392/1‑69 (MQ=255)
                                                    tAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACtg                >  1:545449/1‑69 (MQ=255)
                                                    tAACGTTCGCAATGTCCTTCAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCaa                            <  1:2067361/57‑1 (MQ=255)
                                                    tAACGTTCGCAATGTCCTTCAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCaa                            <  1:2628781/57‑1 (MQ=255)
                                                     aaCGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACtg                >  1:687544/1‑68 (MQ=255)
                                                     aaCGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACtg                >  1:2806792/1‑68 (MQ=255)
                                                        gTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACTGTGTa            <  1:2226694/69‑1 (MQ=255)
                                                           cGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACt                          >  1:38681/1‑52 (MQ=255)
                                                           cGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACt                          >  1:2365500/1‑52 (MQ=255)
                                                           cGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACt                          >  1:2276528/1‑52 (MQ=255)
                                                           cGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACTGTGTa            >  1:560734/1‑66 (MQ=255)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:41301/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:545743/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:709168/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:836192/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:1396192/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:867009/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:3379114/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:936675/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:1697325/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:3261777/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:2387882/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:3204049/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:2424349/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:2192905/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:3202565/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:2011172/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:1982762/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:2847171/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:2936026/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:17880/1‑39 (MQ=37)
                                                             cAATGTCCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:1827816/1‑39 (MQ=37)
                                                             cAATGACCTTCAGTTGCGGCACCGCATTAACCAGAtttt                                     >  1:1509455/1‑39 (MQ=25)
                                                              gaTGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACTGTGTAGTTAg       <  1:132769/67‑1 (MQ=255)
                                                              aaTGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACTGTGTAGTTAg       <  1:2702240/68‑1 (MQ=255)
                                                              aaTGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACTGTGTAGTTAg       <  1:3424011/68‑1 (MQ=255)
                                                              aaTGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACTGTGTAGTTAg       <  1:2231626/68‑1 (MQ=255)
                                                              aaTGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACTGTGTAGTTAg       <  1:88299/68‑1 (MQ=255)
                                                              aaTGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACTGTGTAGTTAg       <  1:10188/68‑1 (MQ=255)
                                                                 gTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACTGTGTAGTTAg       <  1:1177671/65‑1 (MQ=255)
                                                                  tCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCGCTGTGTAGTTAGATTTg  >  1:2566263/1‑69 (MQ=255)
                                                                  tCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACTGTGTAGTTAGATTTg  >  1:317675/1‑69 (MQ=255)
                                                                  tCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACTGTGTAGTTAGATTTg  >  1:3194385/1‑69 (MQ=255)
                                                                  tCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACTGTGTAGTTAGATTTg  >  1:1436163/1‑69 (MQ=255)
                                                                  tCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACTGTGTAGTTAGATTTg  >  1:3252358/1‑69 (MQ=255)
                                                                  tCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACTGTGTAGTTAGATTTg  >  1:977032/1‑69 (MQ=255)
                                                                    |                                                                  
CAGACATTATCGTTCATGTCCTGGGAGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGATTTTCTACGCCAACTTTACCCACTGTGTAGTTAGATTTG  >  minE/1987460‑1987594

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: