Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,901,157 G→A 35.8% I628I (ATC→ATT 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,1570GA35.8% 42.2 / 35.7 64I628I (ATC→ATTrecBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base G (31/10);  new base A (14/9);  total (45/19)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

AACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGT  >  minE/1901099‑1901212
                                                          |                                                       
aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAgcg                                               >  1:715777/1‑69 (MQ=255)
aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAgcg                                               >  1:1941883/1‑69 (MQ=255)
aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAgcg                                               >  1:1972807/1‑69 (MQ=255)
aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAgcg                                               >  1:2136316/1‑69 (MQ=255)
aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAgcg                                               >  1:2937659/1‑69 (MQ=255)
aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAgcg                                               >  1:1710810/1‑69 (MQ=255)
aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAgcg                                               >  1:1658908/1‑69 (MQ=255)
aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAgcg                                               >  1:2750194/1‑69 (MQ=255)
aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAgcg                                               >  1:2938034/1‑69 (MQ=255)
aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAgcg                                               >  1:3410116/1‑69 (MQ=255)
aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAgcg                                               >  1:735064/1‑69 (MQ=255)
aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAgcg                                               >  1:3167896/1‑69 (MQ=255)
aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAgcg                                               >  1:971237/1‑69 (MQ=255)
aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAgcg                                               >  1:1242357/1‑69 (MQ=255)
aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAgcg                                               >  1:2603492/1‑69 (MQ=255)
aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAgcg                                               >  1:3273643/1‑69 (MQ=255)
aaCCATCGAACTCTTCGACTACCACATCCAACGCATGTTCGTCATTGTTCAGCGTATCGATATCCAgcg                                               >  1:1415872/1‑69 (MQ=255)
            ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                   <  1:582955/69‑1 (MQ=255)
            ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                   <  1:2321438/69‑1 (MQ=255)
            ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                   <  1:1152314/69‑1 (MQ=255)
            ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                   <  1:2185511/69‑1 (MQ=255)
            ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                   <  1:2745888/69‑1 (MQ=255)
            ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                   <  1:1957973/69‑1 (MQ=255)
            ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                   <  1:702369/69‑1 (MQ=255)
            ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                   <  1:848687/69‑1 (MQ=255)
            ctTCGACTACCGCATCCCACGCATGTTCGTCATTATTCAGCGTTTCAATATCCAGCGCATTCAGCCCca                                   <  1:1280235/69‑1 (MQ=255)
                             cACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGa                            >  1:2622591/1‑59 (MQ=255)
                             cACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGAc                           >  1:2115174/1‑60 (MQ=255)
                             cACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGAc                           >  1:1642098/1‑60 (MQ=255)
                             cACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGAc                           >  1:2476637/1‑60 (MQ=255)
                                 cATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCTCGTTCAGCCCCATCATTGACGTTGCCAGCGCAc              <  1:121856/69‑1 (MQ=255)
                                 cATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCAc              <  1:1853653/69‑1 (MQ=255)
                                 cATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCAc              <  1:496769/69‑1 (MQ=255)
                                 cATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCAc              <  1:2801361/69‑1 (MQ=255)
                                 cATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCAc              <  1:2547072/69‑1 (MQ=255)
                                  aTGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCAc              <  1:2023993/68‑1 (MQ=255)
                                  aTGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCAc              <  1:1621848/68‑1 (MQ=255)
                                     ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGc                                       <  1:986201/40‑1 (MQ=255)
                                     ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGc                                       <  1:1627853/40‑1 (MQ=255)
                                     ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGc                                       <  1:154365/40‑1 (MQ=255)
                                     ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa            >  1:551093/1‑67 (MQ=255)
                                     ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa            >  1:2318095/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)
                                     ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa            >  1:3020365/1‑67 (MQ=255)
                                     ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa            >  1:305474/1‑67 (MQ=255)
                                     ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa            >  1:3166948/1‑67 (MQ=255)
                                     ttCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa            >  1:92702/1‑67 (MQ=255)
                                     ttCGTCATTGTTCAGCGTTTCGATATCCAGAGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTa            >  1:2400026/1‑67 (MQ=255)
                                             tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt  >  1:2967814/3‑69 (MQ=255)
                                             tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt  >  1:2907740/3‑69 (MQ=255)
                                             tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt  >  1:2565046/3‑69 (MQ=255)
                                             tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt  >  1:1848835/3‑69 (MQ=255)
                                             tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt  >  1:1659770/3‑69 (MQ=255)
                                             tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt  >  1:1555806/3‑69 (MQ=255)
                                             tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt  >  1:1513974/3‑69 (MQ=255)
                                             tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt  >  1:816600/3‑69 (MQ=255)
                                             tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt  >  1:1404426/3‑69 (MQ=255)
                                             tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt  >  1:886310/3‑69 (MQ=255)
                                             tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt  >  1:889485/3‑69 (MQ=255)
                                             tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt  >  1:1324516/3‑69 (MQ=255)
                                             tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt  >  1:1275721/3‑69 (MQ=255)
                                             tattCAGCGTTTCAATATCCAGCGCATTCAGCCCCATCATTGACGTTGCCAGCGCACTGCGCAGGgtgt  >  1:1207288/3‑69 (MQ=255)
                                                          |                                                       
AACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGT  >  minE/1901099‑1901212

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: