Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 241,259 C→T 55.3% T988T (ACG→ACA sbcC ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE241,2590CT55.3% 3.2 / 27.7 38T988T (ACG→ACAsbcCexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (8/9);  new base T (18/3);  total (26/12)
Fisher's exact test for biased strand distribution p-value = 1.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATCAAGGGCGGTATCCAGCGTTTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAG  >  minE/241224‑241324
                                   |                                                                 
gcATCAAGGGCGGTATCCAGCGTTTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGCGAGTc                                   >  1:1398402/1‑68 (MQ=255)
gcATCAAGGGCGGTATCCAGCGTTTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGCGAGTc                                   >  1:1645330/1‑68 (MQ=255)
gcATCAAGGGCGGTATCCAGCGTTTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGCGAGTc                                   >  1:1105083/1‑68 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCaa                                                            >  1:2153170/1‑38 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATAc                              >  1:999405/1‑68 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:266452/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:459326/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:440841/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:381284/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:354098/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:3425701/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:2858370/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:2813448/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:281087/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:266612/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:2225851/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:1117661/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:2168190/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:1396678/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:1525109/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:170689/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGc                                        <  1:2768952/58‑1 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGc                                        <  1:2152267/58‑1 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGc                                        <  1:1652382/58‑1 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGCGAGTc                                   <  1:2448024/63‑1 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGCGAGTc                                   <  1:607089/63‑1 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGCGAGTc                                   <  1:1313682/63‑1 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGCGAGTc                                   <  1:2013470/63‑1 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGCGAGTc                                   <  1:2527112/63‑1 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:2896419/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:2641117/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATCCAGCGTGCCAAAACCATCATCAAGGAACAGCGAGTCAATACg                             >  1:882720/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTTTCGCTATACAGCGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACg                             >  1:861211/1‑69 (MQ=255)
     aaGGGCGGTATCCAGCGTCTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGCGa                                      >  1:1663323/1‑60 (MQ=255)
              aTCCAGCGTTTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGAGAGTc                                   <  1:342532/54‑1 (MQ=255)
                                cAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAag  <  1:3170734/69‑1 (MQ=255)
                                cAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAag  <  1:2674692/69‑1 (MQ=255)
                                cAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAag  <  1:2184802/69‑1 (MQ=255)
                                   |                                                                 
GCATCAAGGGCGGTATCCAGCGTTTCGCTATCCAGCGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAG  >  minE/241224‑241324

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: