Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,668,051 G→A 7.4% Q378Q (CAG→CAA xseA → exonuclease VII, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,668,0510GA 7.4% 218.1 / 5.8 67Q378Q (CAG→CAAxseAexonuclease VII, large subunit
Reads supporting (aligned to +/- strand):  ref base G (39/23);  new base A (2/3);  total (41/26)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CGAAGATTCATCGCGCGCAAACGCGCATTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCGTAAGCCCACTGTCAACGC  >  minE/1667981‑1668121
                                                                      |                                                                      
cGAAGATTCATCGCGCGCAAACGCGCATTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAg                                                                        <  1:1834460/71‑1 (MQ=255)
cGAAGATTCATCGCGCGCAAACGCGCATTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAg                                                                        >  1:2776228/1‑71 (MQ=255)
    gATTCATCGCGCGCAAACGCGCATTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTa                                                                    <  1:1328964/71‑1 (MQ=255)
       tCATCGCGCGCAAACGCGCATTCAGCAATTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAg                                                                   >  1:2827011/1‑69 (MQ=255)
       tCATCGCGCGCAAACGCGCATTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAgcg                                                                 >  1:2726229/1‑71 (MQ=255)
       tCATCGCGCGCAAACGCGCATTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAgcg                                                                 >  1:1599255/1‑71 (MQ=255)
       tCATCGCGCGCAAACGCGCATTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAg                                                                   >  1:2880410/1‑69 (MQ=255)
        cATCGCGCGCAAACGCGCATTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAgcgc                                                                >  1:2538778/1‑71 (MQ=255)
        cATCGCGCGCAAACGCGCATTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAgcgc                                                                >  1:3342998/1‑71 (MQ=255)
         aTCGCGCGCAAACGCGCATTCAGCAACTGCAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAgcgc                                                                <  1:1545338/70‑1 (MQ=255)
           cgcgcgCAAACGCGCATTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCAc                                                             >  1:1217911/1‑71 (MQ=255)
            gcgcgcAAACGCGCATTCAGCAACTGGAATATGGTTTAGCAGAAACCCTGCGCGCACAACTTAGCGCCAc                                                             <  1:998507/70‑1 (MQ=255)
            gcgcgcAAACGCGCATTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCAcg                                                            >  1:2356337/1‑71 (MQ=255)
                     cgcgCATTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGt                                                   <  1:750265/71‑1 (MQ=255)
                       cgcATTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAg                                                                        <  1:1358129/48‑1 (MQ=255)
                        gcATTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTc                                                <  1:2076749/71‑1 (MQ=255)
                         cATTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCg                                               <  1:1883326/71‑1 (MQ=255)
                          aTTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAACTTAGCGCCACGCGTGAACGTTTCgg                                              >  1:1940892/1‑71 (MQ=255)
                            tCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGt                                             <  1:32139/70‑1 (MQ=255)
                            tCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGt                                             <  1:672631/70‑1 (MQ=255)
                             cAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGaa                                                      >  1:509871/1‑60 (MQ=255)
                             cAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTaa                                           >  1:374521/1‑71 (MQ=255)
                             cAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTaa                                           >  1:1207715/1‑71 (MQ=255)
                              aGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAAt                                          <  1:2662899/71‑1 (MQ=255)
                               gCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTaa                                           <  1:2715041/69‑1 (MQ=255)
                                 aaCTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATg                                         >  1:1138619/1‑69 (MQ=255)
                                    tGGAATATCGTTTAGCAGAAACCCTGCGCGCACGGCTTAGCGCCAcgc                                                           >  1:503265/1‑48 (MQ=255)
                                    tGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCAcgc                                                           >  1:994969/1‑48 (MQ=255)
                                    tGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTaa                                           <  1:2734045/64‑1 (MQ=255)
                                     ggAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTaa                                   >  1:528320/1‑71 (MQ=255)
                                      gAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAAc                                  >  1:1506945/1‑71 (MQ=255)
                                      gAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAAc                                  >  1:661386/1‑71 (MQ=255)
                                      gAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAAc                                  >  1:1838380/1‑71 (MQ=255)
                                           tCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGaa                                                      >  1:3165921/1‑46 (MQ=255)
                                           tCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAg                                      >  1:3012595/1‑62 (MQ=255)
                                           tCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCAcc                             >  1:117551/1‑71 (MQ=255)
                                           tCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCAc                              <  1:83656/70‑1 (MQ=255)
                                            cGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCt                            >  1:2094414/1‑71 (MQ=255)
                                            cGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCt                            >  1:278823/1‑71 (MQ=255)
                                             gTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTc                           >  1:3293850/1‑71 (MQ=255)
                                             gTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTc                           >  1:2211774/1‑71 (MQ=255)
                                             gTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTc                           >  1:3297865/1‑71 (MQ=255)
                                               ttAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTaa                                   >  1:1810441/1‑61 (MQ=255)
                                                tAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAAc                                                     >  1:2279537/1‑42 (MQ=255)
                                                tAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAAc                                                     >  1:1613533/1‑42 (MQ=255)
                                                tAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGa                         <  1:2453562/70‑1 (MQ=255)
                                                  gCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGc                      <  1:988467/71‑1 (MQ=255)
                                                    aGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCg                    <  1:880239/71‑1 (MQ=255)
                                                    aGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCg                    <  1:2766424/71‑1 (MQ=255)
                                                     gAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAAc                                                     >  1:1590019/1‑37 (MQ=255)
                                                     gAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCg                    >  1:1277627/1‑70 (MQ=255)
                                                     gAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCg                    >  1:354026/1‑70 (MQ=255)
                                                     gAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCg                    >  1:67204/1‑70 (MQ=255)
                                                      aaaCCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCGt                   >  1:1545598/1‑70 (MQ=255)
                                                       aaCCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATg                                         <  1:2807527/47‑1 (MQ=255)
                                                       aaCCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATg                                         <  1:705929/47‑1 (MQ=255)
                                                       aaCCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATg                                         <  1:2084581/47‑1 (MQ=255)
                                                       aaCCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATg                                         <  1:2623776/47‑1 (MQ=255)
                                                        aCCCTGCGCGCACAACTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCa                               >  1:2141022/1‑56 (MQ=255)
                                                         cccTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCg                          <  1:1184429/60‑1 (MQ=255)
                                                           cTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCGTAAGccc             <  1:319577/71‑1 (MQ=255)
                                                           cTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCGTAAGccc             <  1:1798415/71‑1 (MQ=255)
                                                              cgcgcACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCGTAAGCCCACt          >  1:1750000/1‑71 (MQ=255)
                                                                  cacaACTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCGTAAGCCCACTGTCa      <  1:2414169/71‑1 (MQ=255)
                                                                  cacaACTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCGTAAGCCCACTGTCa      <  1:2198468/71‑1 (MQ=255)
                                                                   acaGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCGTAAGCCCACTGTCaa     >  1:3124719/1‑71 (MQ=255)
                                                                   acaGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCGTAAGCCCACTGTCaa     >  1:280801/1‑71 (MQ=255)
                                                                    caGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGa                         <  1:16424/50‑1 (MQ=255)
                                                                    caGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCg                    >  1:2954890/1‑55 (MQ=255)
                                                                     aGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTaa                                   >  1:112272/1‑39 (MQ=255)
                                                                      gCTTATCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCGTAAGCCCACTGTCAACGc  <  1:1300291/71‑1 (MQ=255)
                                                                      gCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCGTAAGCCCACTGTCAACGc  >  1:544314/1‑71 (MQ=255)
                                                                      gCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCGTAAGCCCACTGTCAACGc  <  1:2630732/71‑1 (MQ=255)
                                                                      gCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCGTAAGCCCACTGTCAACGc  >  1:2535297/1‑71 (MQ=255)
                                                                      gCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCGTAAGCCCACTGTCAACGc  >  1:1902854/1‑71 (MQ=255)
                                                                      gCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCGTAAGCCCACTGTCAACGc  <  1:2248021/71‑1 (MQ=255)
                                                                      |                                                                      
CGAAGATTCATCGCGCGCAAACGCGCATTCAGCAACTGGAATATCGTTTAGCAGAAACCCTGCGCGCACAGCTTAGCGCCACGCGTGAACGTTTCGGTAATGCAGTAACGCACCTCGAAGCCGTAAGCCCACTGTCAACGC  >  minE/1667981‑1668121

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: