Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,901,169 G→A 56.4% N624N (AAC→AAT recB ← exonuclease V (RecBCD complex), beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,901,1690GA56.4% 10.1 / 40.9 55N624N (AAC→AATrecBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base G (14/10);  new base A (14/17);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACC  >  minE/1901111‑1901235
                                                          |                                                                  
ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                                          <  1:1152314/69‑1 (MQ=255)
ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                                          <  1:848687/69‑1 (MQ=255)
ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                                          <  1:1957973/69‑1 (MQ=255)
ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                                          <  1:702369/69‑1 (MQ=255)
ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                                          <  1:1280235/69‑1 (MQ=255)
ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                                          <  1:582955/69‑1 (MQ=255)
ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                                          <  1:2185511/69‑1 (MQ=255)
ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                                          <  1:2321438/69‑1 (MQ=255)
ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                                          <  1:2745888/69‑1 (MQ=255)
                 cACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGa                                                   >  1:2622591/1‑59 (MQ=255)
                 cACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGAc                                                  >  1:2476637/1‑60 (MQ=255)
                 cACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGAc                                                  >  1:1642098/1‑60 (MQ=255)
                 cACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGAc                                                  >  1:2115174/1‑60 (MQ=255)
                     cATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCTCGTTCAGCCCCATCATTGACGTTGCCAGCGCAc                                     <  1:121856/69‑1 (MQ=255)
                     cATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCAc                                     <  1:2547072/69‑1 (MQ=255)
                     cATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCAc                                     <  1:2801361/69‑1 (MQ=255)
                     cATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCAc                                     <  1:496769/69‑1 (MQ=255)
                     cATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCAc                                     <  1:1853653/69‑1 (MQ=255)
                      aTGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCAc                                     <  1:1621848/68‑1 (MQ=255)
                      aTGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCAc                                     <  1:2023993/68‑1 (MQ=255)
                         ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGc                                                              <  1:986201/40‑1 (MQ=255)
                         ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGc                                                              <  1:154365/40‑1 (MQ=255)
                         ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGc                                                              <  1:1627853/40‑1 (MQ=255)
                         ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                   >  1:2318095/1‑67 (MQ=255)
                         ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                   >  1:3020365/1‑67 (MQ=255)
                         ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                   >  1:305474/1‑67 (MQ=255)
                         ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                   >  1:92702/1‑67 (MQ=255)
                         ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                   >  1:3166948/1‑67 (MQ=255)
                         ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                   >  1:551093/1‑67 (MQ=255)
                         ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                   >  1:2050756/1‑67 (MQ=255)
                         ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                   >  1:2046787/1‑67 (MQ=255)
                         ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                   >  1:1814983/1‑67 (MQ=255)
                         ttCGTCATTGTTCAGCGTTTCGATATCCAGAGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                   >  1:2400026/1‑67 (MQ=255)
                                 tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                         >  1:1404426/3‑69 (MQ=255)
                                 tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                         >  1:1324516/3‑69 (MQ=255)
                                 tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                         >  1:1207288/3‑69 (MQ=255)
                                 tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                         >  1:1275721/3‑69 (MQ=255)
                                 tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                         >  1:2967814/3‑69 (MQ=255)
                                 tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                         >  1:2907740/3‑69 (MQ=255)
                                 tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                         >  1:1513974/3‑69 (MQ=255)
                                 tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                         >  1:816600/3‑69 (MQ=255)
                                 tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                         >  1:1848835/3‑69 (MQ=255)
                                 tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                         >  1:886310/3‑69 (MQ=255)
                                 tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                         >  1:2565046/3‑69 (MQ=255)
                                 tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                         >  1:1555806/3‑69 (MQ=255)
                                 tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                         >  1:889485/3‑69 (MQ=255)
                                 tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                         >  1:1659770/3‑69 (MQ=255)
                                                      gcgcATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGtcatca    <  1:2739752/69‑1 (MQ=255)
                                                      gcgcATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGtcatca    <  1:285945/69‑1 (MQ=255)
                                                      gcgcATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGtcatca    <  1:3072191/69‑1 (MQ=255)
                                                      gcgcATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGtcatca    <  1:1562450/69‑1 (MQ=255)
                                                      gcgcATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGtcatca    <  1:487007/69‑1 (MQ=255)
                                                      gcgcATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGtcatca    <  1:1965675/69‑1 (MQ=255)
                                                      gcgcATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGtcatca    <  1:2116892/69‑1 (MQ=255)
                                                      gcgcATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACATTCGGGTGtcatca    <  1:651576/69‑1 (MQ=255)
                                                       cgcaTTCAGCCCCATCATTGACGTTGCCAGCGCACTGCg                                 >  1:186074/1‑39 (MQ=37)
                                                       cgcaTTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGtcatca    >  1:2847672/1‑68 (MQ=255)
                                                       cgcaTTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGTCATCAc   >  1:932365/1‑69 (MQ=255)
                                                       cgcaTTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGTCATCAc   >  1:2815581/1‑69 (MQ=255)
                                                       cgcaTTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGTCATCAc   >  1:1552229/1‑69 (MQ=255)
                                                       cgcaTTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGTCATCAc   >  1:2513352/1‑69 (MQ=255)
                                                       cgcaTTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGTCATCAc   >  1:1557313/1‑69 (MQ=255)
                                                       cgcaTTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGTCATCAc   >  1:2146854/1‑69 (MQ=255)
                                                        gcgTTCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCAcc  <  1:1155636/69‑1 (MQ=255)
                                                          |                                                                  
CTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACC  >  minE/1901111‑1901235

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