Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,860,429 C→T 44.7% D452D (GAC→GAT mutL → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,860,4290CT44.7% 13.2 / 52.0 67D452D (GAC→GATmutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (3/34);  new base T (22/8);  total (25/42)
Fisher's exact test for biased strand distribution p-value = 2.91e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

TTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCAC  >  minE/2860377‑2860484
                                                    |                                                       
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGc                                                      >  1:1947663/1‑56 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACa                                                  >  1:1620093/1‑60 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACa                                                  >  1:2128535/1‑60 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:2303976/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:178659/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:1836082/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:1916601/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:1920778/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:2067972/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:2099653/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:1077786/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:2597839/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:2920261/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:3127185/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:3138428/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:3316540/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:545512/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:617553/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:858870/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:1657237/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:1196012/1‑61 (MQ=255)
ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:1611910/1‑61 (MQ=255)
   gTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACtttt                                        <  1:1698779/67‑1 (MQ=255)
    tCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTAt                                      <  1:1249398/68‑1 (MQ=255)
    tCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTAt                                      <  1:2209277/68‑1 (MQ=255)
    tCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTAt                                      <  1:658208/68‑1 (MQ=255)
    tCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTAt                                      <  1:3224927/68‑1 (MQ=255)
         tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc                                <  1:1560742/69‑1 (MQ=255)
         tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc                                <  1:3129408/69‑1 (MQ=255)
         tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc                                <  1:1099131/69‑1 (MQ=255)
         tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc                                <  1:1593080/69‑1 (MQ=255)
         tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc                                <  1:2351107/69‑1 (MQ=255)
         tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc                                <  1:2294070/69‑1 (MQ=255)
         tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                <  1:484351/69‑1 (MQ=255)
         tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                <  1:538611/69‑1 (MQ=255)
         tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                <  1:640577/69‑1 (MQ=255)
         tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                <  1:3327439/69‑1 (MQ=255)
         tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                <  1:1849690/69‑1 (MQ=255)
         tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                <  1:206960/69‑1 (MQ=255)
         tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                <  1:682065/69‑1 (MQ=255)
          aCTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGcc                               >  1:2159661/1‑69 (MQ=255)
                      ccATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAAc                    <  1:1454079/68‑1 (MQ=255)
                         ttCCGACTGTCCGTTCTTGGAGCGCGACGGCAACAATTCACTTTTATCCTTGc                                <  1:1558319/53‑1 (MQ=255)
                          tCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc                                <  1:865534/52‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:598319/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:537392/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:494510/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:493841/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:70436/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:1023358/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:928653/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:3309405/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:1202238/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:2980272/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:120397/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:2479329/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:2472657/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:1937381/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:1571837/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:1900866/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:1895417/46‑1 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTgg               >  1:3048764/1‑68 (MQ=255)
                           ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGc              >  1:2009470/1‑69 (MQ=255)
                            cGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:1168025/45‑1 (MQ=255)
                                   gCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc                                <  1:1237948/43‑1 (MQ=38)
                                     gTTGCTGGAGCGCGACGTCAACATTTCACTTTTATCCTTGCCAGTGGCAGAAc                    <  1:1558147/53‑1 (MQ=255)
                                           ggAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGcac  <  1:1203891/65‑1 (MQ=255)
                                                    |                                                       
TTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCAC  >  minE/2860377‑2860484

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