Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,860,420 G→T 42.1% L449L (CTG→CTT mutL → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,860,4200GT42.1% 49.7 / 74.1 83L449L (CTG→CTTmutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (3/45);  new base T (23/12);  total (26/57)
Fisher's exact test for biased strand distribution p-value = 6.18e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCAC  >  minE/2860360‑2860484
                                                            |                                                                
ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAgcgc                                                            <  1:825617/67‑1 (MQ=255)
ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAgcgc                                                            <  1:1878135/67‑1 (MQ=255)
ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAgcgc                                                            <  1:1867123/67‑1 (MQ=255)
ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAgcgc                                                            <  1:502237/67‑1 (MQ=255)
ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc                                                            <  1:317424/67‑1 (MQ=255)
ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc                                                            <  1:1791249/67‑1 (MQ=255)
ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc                                                            <  1:1809246/67‑1 (MQ=255)
ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc                                                            <  1:1494761/67‑1 (MQ=255)
ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc                                                            <  1:1420487/67‑1 (MQ=255)
ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc                                                            <  1:1402120/67‑1 (MQ=255)
ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc                                                            <  1:1330371/67‑1 (MQ=255)
ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc                                                            <  1:3208810/67‑1 (MQ=255)
ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc                                                            <  1:3434730/67‑1 (MQ=255)
ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc                                                            <  1:17328/67‑1 (MQ=255)
ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGGTGCTGGAgcgc                                                            <  1:2065550/67‑1 (MQ=255)
                 ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAgc                                                              >  1:2998873/1‑48 (MQ=255)
                 ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGc                                                      >  1:1947663/1‑56 (MQ=255)
                 ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACa                                                  >  1:2128535/1‑60 (MQ=255)
                 ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACa                                                  >  1:1620093/1‑60 (MQ=255)
                 ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:1916601/1‑61 (MQ=255)
                 ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:2067972/1‑61 (MQ=255)
                 ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:1920778/1‑61 (MQ=255)
                 ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:1077786/1‑61 (MQ=255)
                 ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:2099653/1‑61 (MQ=255)
                 ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:2303976/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:1196012/1‑61 (MQ=255)
                 ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:1611910/1‑61 (MQ=255)
                 ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:1657237/1‑61 (MQ=255)
                 ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:178659/1‑61 (MQ=255)
                 ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt                                                 >  1:1836082/1‑61 (MQ=255)
                    gTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACtttt                                        <  1:1698779/67‑1 (MQ=255)
                     tCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTAt                                      <  1:1249398/68‑1 (MQ=255)
                     tCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTAt                                      <  1:3224927/68‑1 (MQ=255)
                     tCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTAt                                      <  1:2209277/68‑1 (MQ=255)
                     tCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTAt                                      <  1:658208/68‑1 (MQ=255)
                          tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc                                <  1:3129408/69‑1 (MQ=255)
                          tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc                                <  1:2294070/69‑1 (MQ=255)
                          tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc                                <  1:1593080/69‑1 (MQ=255)
                          tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc                                <  1:2351107/69‑1 (MQ=255)
                          tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc                                <  1:1099131/69‑1 (MQ=255)
                          tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc                                <  1:1560742/69‑1 (MQ=255)
                          tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                <  1:3327439/69‑1 (MQ=255)
                          tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                <  1:538611/69‑1 (MQ=255)
                          tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                <  1:640577/69‑1 (MQ=255)
                          tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                <  1:682065/69‑1 (MQ=255)
                          tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                <  1:206960/69‑1 (MQ=255)
                          tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                <  1:484351/69‑1 (MQ=255)
                          tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                <  1:1849690/69‑1 (MQ=255)
                           aCTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGcc                               >  1:2159661/1‑69 (MQ=255)
                             tGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc                                                            <  1:661426/38‑1 (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:537392/46‑1 (MQ=255)
                                            ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:1202238/46‑1 (MQ=255)
                                            ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:70436/46‑1 (MQ=255)
                                            ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:598319/46‑1 (MQ=255)
                                            ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:928653/46‑1 (MQ=255)
                                            ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:1023358/46‑1 (MQ=255)
                                            ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:1900866/46‑1 (MQ=255)
                                            ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:494510/46‑1 (MQ=255)
                                            ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:493841/46‑1 (MQ=255)
                                            ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:120397/46‑1 (MQ=255)
                                            ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:3309405/46‑1 (MQ=255)
                                            ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:2980272/46‑1 (MQ=255)
                                            ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                     <  1:1571837/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: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)
                                                            |                                                                
GGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCAC  >  minE/2860360‑2860484

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: