Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 241,319 C→T 78.4% S968S (TCG→TCA sbcC ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE241,3190CT78.4% 58.7 / 10.6 46S968S (TCG→TCAsbcCexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (6/4);  new base T (7/29);  total (13/33)
Fisher's exact test for biased strand distribution p-value = 1.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CAGCGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCG  >  minE/241256‑241387
                                                               |                                                                    
cAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAag                                                                 <  1:3170734/69‑1 (MQ=255)
cAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAag                                                                 <  1:2674692/69‑1 (MQ=255)
cAGTGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAag                                                                 <  1:2184802/69‑1 (MQ=255)
             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:1130134/65‑1 (MQ=255)
             ccTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAagcgccagtgc                                                        <  1:128314/65‑1 (MQ=255)
             ccTTCATCAAGGAACAGCGAGTCAATACGCGTTTTATGGCTGACCAGATCTGAAagcgccagtgc                                                        <  1:1254014/65‑1 (MQ=255)
                   tcaAGGAACAGCGAGTCAATACGTGTTTTATGGCTGACCACATCCGAAAGcgccagcgccagcgccag                                               <  1:1010051/68‑1 (MQ=255)
                     aaGGAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAagcgcc                                                             >  1:2178865/1‑52 (MQ=255)
                     aaGGAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAagcgcc                                                             >  1:2329744/1‑52 (MQ=255)
                     aaGGAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAagcgc                                                              >  1:260596/1‑51 (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:700222/68‑1 (MQ=39)
                                                aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                  <  1:2057757/68‑1 (MQ=39)
                                                aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                  <  1:3252170/68‑1 (MQ=39)
                                                aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                  <  1:1422102/68‑1 (MQ=39)
                                                aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                  <  1:3077584/68‑1 (MQ=39)
                                                aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                  <  1:2202276/68‑1 (MQ=39)
                                                aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                  <  1:1499451/68‑1 (MQ=39)
                                                aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                  <  1:2254444/68‑1 (MQ=39)
                                                aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                  <  1:2701259/68‑1 (MQ=39)
                                                aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                  <  1:243388/68‑1 (MQ=39)
                                                aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                  <  1:1483184/68‑1 (MQ=39)
                                                aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                  <  1:2164540/68‑1 (MQ=39)
                                                    cTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                  <  1:2455489/64‑1 (MQ=38)
                                                      gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                         >  1:343557/1‑55 (MQ=38)
                                                      gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                         >  1:1805095/1‑55 (MQ=38)
                                                      gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                         >  1:603413/1‑55 (MQ=38)
                                                      gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                         >  1:865774/1‑55 (MQ=38)
                                                      gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                         >  1:1459480/1‑55 (MQ=38)
                                                      gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                         >  1:2891849/1‑55 (MQ=38)
                                                      gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                         >  1:1814076/1‑55 (MQ=38)
                                                               cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGTGTACGGGTATcgcg  <  1:3141985/69‑1 (MQ=255)
                                                               cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATcgcg  <  1:2288762/69‑1 (MQ=255)
                                                               cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATcgcg  <  1:631475/69‑1 (MQ=255)
                                                               cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATcgcg  <  1:1242502/69‑1 (MQ=255)
                                                               cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATcgcg  <  1:2875911/69‑1 (MQ=255)
                                                               cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTACGCCGCCGGAAAGGGTACGGGTATcgcg  <  1:1182503/69‑1 (MQ=255)
                                                               |                                                                    
CAGCGTGCCAAAACCTTCATCAAGGAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCG  >  minE/241256‑241387

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: