Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,901,202 A→G 48.4% R613R (CGT→CGC 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,2020AG48.4% ‑4.3 / 55.4 62R613R (CGT→CGCrecBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base A (20/12);  new base G (22/8);  total (42/20)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCT  >  minE/1901136‑1901265
                                                                  |                                                               
ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                                                 >  1:2046787/1‑67 (MQ=255)
ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                                                 >  1:305474/1‑67 (MQ=255)
ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                                                 >  1:92702/1‑67 (MQ=255)
ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                                                 >  1:2050756/1‑67 (MQ=255)
ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                                                 >  1:3020365/1‑67 (MQ=255)
ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                                                 >  1:3166948/1‑67 (MQ=255)
ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                                                 >  1:1814983/1‑67 (MQ=255)
ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                                                 >  1:2318095/1‑67 (MQ=255)
ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                                                 >  1:551093/1‑67 (MQ=255)
ttCGTCATTGTTCAGCGTTTCGATATCCAGAGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa                                                                 >  1:2400026/1‑67 (MQ=255)
        tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                                                       >  1:2565046/3‑69 (MQ=255)
        tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                                                       >  1:1848835/3‑69 (MQ=255)
        tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                                                       >  1:2967814/3‑69 (MQ=255)
        tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                                                       >  1:2907740/3‑69 (MQ=255)
        tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                                                       >  1:1659770/3‑69 (MQ=255)
        tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                                                       >  1:816600/3‑69 (MQ=255)
        tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                                                       >  1:1555806/3‑69 (MQ=255)
        tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                                                       >  1:886310/3‑69 (MQ=255)
        tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                                                       >  1:1513974/3‑69 (MQ=255)
        tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                                                       >  1:1404426/3‑69 (MQ=255)
        tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                                                       >  1:1324516/3‑69 (MQ=255)
        tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                                                       >  1:889485/3‑69 (MQ=255)
        tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                                                       >  1:1275721/3‑69 (MQ=255)
        tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt                                                       >  1:1207288/3‑69 (MQ=255)
                             gcgcATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGtcatca                                  <  1:285945/69‑1 (MQ=255)
                             gcgcATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGtcatca                                  <  1:2116892/69‑1 (MQ=255)
                             gcgcATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGtcatca                                  <  1:487007/69‑1 (MQ=255)
                             gcgcATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGtcatca                                  <  1:3072191/69‑1 (MQ=255)
                             gcgcATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGtcatca                                  <  1:1965675/69‑1 (MQ=255)
                             gcgcATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGtcatca                                  <  1:1562450/69‑1 (MQ=255)
                             gcgcATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGtcatca                                  <  1:2739752/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:2513352/1‑69 (MQ=255)
                              cgcaTTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGTCATCAc                                 >  1:2815581/1‑69 (MQ=255)
                              cgcaTTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGTCATCAc                                 >  1:2146854/1‑69 (MQ=255)
                              cgcaTTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGTCATCAc                                 >  1:1557313/1‑69 (MQ=255)
                              cgcaTTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGTCATCAc                                 >  1:1552229/1‑69 (MQ=255)
                              cgcaTTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGGTGTTCTCACGTTCGGGTGTCATCAc                                 >  1:932365/1‑69 (MQ=255)
                               gcgTTCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCAcc                                <  1:1155636/69‑1 (MQ=255)
                                  ttCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGcc                             <  1:2483153/69‑1 (MQ=255)
                                  ttCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGcc                             <  1:1779058/69‑1 (MQ=255)
                                  ttCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGcc                             <  1:1276975/69‑1 (MQ=255)
                                  ttCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGcc                             <  1:3360472/69‑1 (MQ=255)
                                  ttCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGcc                             <  1:1863878/69‑1 (MQ=255)
                                  ttCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGcc                             <  1:2027138/69‑1 (MQ=255)
                                  ttCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGcc                             <  1:2869556/69‑1 (MQ=255)
                                  ttCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGcc                             <  1:124085/69‑1 (MQ=255)
                                    cAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCt                            >  1:1590702/1‑68 (MQ=255)
                                    cAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCt                            >  1:2092940/1‑68 (MQ=255)
                                    cAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCt                            >  1:2046358/1‑68 (MQ=255)
                                    cAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCt                            >  1:3065833/1‑68 (MQ=255)
                                    cAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCt                            >  1:2292200/1‑68 (MQ=255)
                                    cAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCt                            >  1:324432/1‑68 (MQ=255)
                                    cAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCt                            >  1:248281/1‑68 (MQ=255)
                                    cAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCt                            >  1:4872/1‑68 (MQ=255)
                                              atTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCa                  >  1:797897/1‑68 (MQ=255)
                                                   cGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCa             <  1:760844/68‑1 (MQ=255)
                                                   cGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCa             <  1:2343696/68‑1 (MQ=255)
                                                   cGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCa             <  1:679651/68‑1 (MQ=255)
                                                       gCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTcc        >  1:1538342/1‑69 (MQ=255)
                                                       gCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTcc        >  1:1672356/1‑69 (MQ=255)
                                                       gCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTcc        >  1:365893/1‑69 (MQ=255)
                                                       gCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTcc        >  1:3221103/1‑69 (MQ=255)
                                                       gCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTcc        >  1:3077478/1‑69 (MQ=255)
                                                       gCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTcc        >  1:2169611/1‑69 (MQ=255)
                                                             gcACTACGCAGGGTGCTCTCACGTTCGGGCGTCATCACCGCCTGCAACACCCAAAGCATTTCCCGCGCt  >  1:2143112/1‑69 (MQ=255)
                                                             gcACTACGCAGGGTGCTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTgcgc   >  1:1115169/1‑68 (MQ=255)
                                                             gcACTACGCAGGGTGCTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCt  >  1:1785831/1‑69 (MQ=255)
                                                             gcACTACGCAGGGTGCTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCt  >  1:2705822/1‑69 (MQ=255)
                                                             gcACTACGCAGGGTGCTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACAAAAGCATTTCCTGCGCt  >  1:2333187/1‑69 (MQ=255)
                                                                  |                                                               
TTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCT  >  minE/1901136‑1901265

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: