Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 241,331 C→T 62.5% A964A (GCG→GCA sbcC ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE241,3310CT62.5% 13.1 / 17.2 32A964A (GCG→GCAsbcCexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (3/9);  new base T (7/13);  total (10/22)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTCATCAAGGAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTG  >  minE/241269‑241398
                                                              |                                                                   
ccTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATTGCTGACCAGATCTGAAagcgccagtgc                                                                   <  1:344115/65‑1 (MQ=255)
ccTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGTTCTGAAagcgccagtgc                                                                   <  1:367465/65‑1 (MQ=255)
ccTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAagcgccagtgc                                                                   <  1:2670966/65‑1 (MQ=255)
ccTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAagcgccagtgc                                                                   <  1:959968/65‑1 (MQ=255)
ccTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAagcgccagtgc                                                                   <  1:914926/65‑1 (MQ=255)
ccTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAagcgccagtgc                                                                   <  1:575114/65‑1 (MQ=255)
ccTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAagcgccagtgc                                                                   <  1:541697/65‑1 (MQ=255)
ccTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAagcgccagtgc                                                                   <  1:3119/65‑1 (MQ=255)
ccTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAagcgccagtgc                                                                   <  1:2873036/65‑1 (MQ=255)
ccTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAagcgccagtgc                                                                   <  1:2867860/65‑1 (MQ=255)
ccTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAagcgccagtgc                                                                   <  1:1254014/65‑1 (MQ=255)
ccTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAagcgccagtgc                                                                   <  1:1130134/65‑1 (MQ=255)
ccTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAagcgccagtgc                                                                   <  1:128314/65‑1 (MQ=255)
      tcaAGGAACAGCGAGTCAATACGTGTTTTATGGCTGACCACATCCGAAAGcgccagcgccagcgccag                                                          <  1:1010051/68‑1 (MQ=255)
        aaGGAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACt                                                       <  1:579790/69‑1 (MQ=255)
        aaGGAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACt                                                       <  1:1523701/69‑1 (MQ=255)
           gAACAGCGAGTCAATACGTGTTTTATGGCTGCCCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTa                                                      >  1:2581714/1‑67 (MQ=255)
             aCAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACt                                                       <  1:606175/64‑1 (MQ=255)
                           cGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTcgccg                                    >  1:113528/1‑69 (MQ=255)
                           cGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTcgccg                                    >  1:2837329/1‑69 (MQ=255)
                                   aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                             <  1:2254444/68‑1 (MQ=39)
                                   aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                             <  1:1483184/68‑1 (MQ=39)
                                   aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                             <  1:3252170/68‑1 (MQ=39)
                                   aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                             <  1:1499451/68‑1 (MQ=39)
                                   aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                             <  1:3077584/68‑1 (MQ=39)
                                   aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                             <  1:1422102/68‑1 (MQ=39)
                                   aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                             <  1:2701259/68‑1 (MQ=39)
                                   aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                             <  1:2057757/68‑1 (MQ=39)
                                   aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                             <  1:2164540/68‑1 (MQ=39)
                                   aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                             <  1:243388/68‑1 (MQ=39)
                                   aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                             <  1:700222/68‑1 (MQ=39)
                                   aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                             <  1:2202276/68‑1 (MQ=39)
                                       cTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                             <  1:2455489/64‑1 (MQ=38)
                                         gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                                    >  1:865774/1‑55 (MQ=38)
                                         gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                                    >  1:1814076/1‑55 (MQ=38)
                                         gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                                    >  1:1805095/1‑55 (MQ=38)
                                         gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                                    >  1:1459480/1‑55 (MQ=38)
                                         gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                                    >  1:2891849/1‑55 (MQ=38)
                                         gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                                    >  1:603413/1‑55 (MQ=38)
                                         gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                                    >  1:343557/1‑55 (MQ=38)
                                                  cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGTGTACGGGTATcgcg             <  1:3141985/69‑1 (MQ=255)
                                                  cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATcgcg             <  1:631475/69‑1 (MQ=255)
                                                  cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATcgcg             <  1:2875911/69‑1 (MQ=255)
                                                  cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATcgcg             <  1:1242502/69‑1 (MQ=255)
                                                  cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATcgcg             <  1:2288762/69‑1 (MQ=255)
                                                  cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTACGCCGCCGGAAAGGGTACGGGTATcgcg             <  1:1182503/69‑1 (MQ=255)
                                                        cgccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGc       <  1:1167534/69‑1 (MQ=255)
                                                        cgccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGc       <  1:639583/69‑1 (MQ=255)
                                                         gccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCa      <  1:417642/69‑1 (MQ=255)
                                                         gccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCa      <  1:421646/69‑1 (MQ=255)
                                                         gccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCa      <  1:339306/69‑1 (MQ=255)
                                                         gccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCa      <  1:3092188/69‑1 (MQ=255)
                                                         gccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCa      <  1:2320207/69‑1 (MQ=255)
                                                         gccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCa      <  1:1831790/69‑1 (MQ=255)
                                                         gccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCa      <  1:1523155/69‑1 (MQ=255)
                                                             gcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATctg  <  1:2654947/69‑1 (MQ=255)
                                                             gcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATctg  <  1:1033587/69‑1 (MQ=255)
                                                              |                                                                   
CCTTCATCAAGGAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTG  >  minE/241269‑241398

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: