Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 241,346 A→G 47.6% V959V (GTT→GTC sbcC ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE241,3460AG47.6% ‑2.7 / 34.6 42V959V (GTT→GTCsbcCexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (5/17);  new base G (7/13);  total (12/30)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATC  >  minE/241280‑241411
                                                                  |                                                                 
gAACAGCGAGTCAATACGTGTTTTATGGCTGCCCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTa                                                                   >  1:2581714/1‑67 (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:2202276/68‑1 (MQ=39)
                        aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                                          <  1:2164540/68‑1 (MQ=39)
                        aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                                          <  1:2057757/68‑1 (MQ=39)
                        aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                                          <  1:243388/68‑1 (MQ=39)
                        aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                                          <  1:2701259/68‑1 (MQ=39)
                        aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                                          <  1:1499451/68‑1 (MQ=39)
                        aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                                          <  1:1483184/68‑1 (MQ=39)
                        aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                                          <  1:1422102/68‑1 (MQ=39)
                        aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                                          <  1:3077584/68‑1 (MQ=39)
                        aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                                          <  1:3252170/68‑1 (MQ=39)
                        aTGGCTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                                          <  1:700222/68‑1 (MQ=39)
                            cTGACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTCGCCGCCGGaaa                                          <  1:2455489/64‑1 (MQ=38)
                              gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                                                 >  1:2891849/1‑55 (MQ=38)
                              gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                                                 >  1:343557/1‑55 (MQ=38)
                              gACCAGATCTGAAAGCGCCAGTGCCAGCGCCAGACTGACGAGGAAACTTTcgccg                                                 >  1:603413/1‑55 (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)
                                       cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGTGTACGGGTATcgcg                          <  1:3141985/69‑1 (MQ=255)
                                       cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATcgcg                          <  1:2288762/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)
                                       cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTACGCCGCCGGAAAGGGTACGGGTATcgcg                          <  1:1182503/69‑1 (MQ=255)
                                             cgccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGc                    <  1:1167534/69‑1 (MQ=255)
                                             cgccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGc                    <  1:639583/69‑1 (MQ=255)
                                              gccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCa                   <  1:2320207/69‑1 (MQ=255)
                                              gccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCa                   <  1:421646/69‑1 (MQ=255)
                                              gccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCa                   <  1:417642/69‑1 (MQ=255)
                                              gccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCa                   <  1:339306/69‑1 (MQ=255)
                                              gccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCa                   <  1:3092188/69‑1 (MQ=255)
                                              gccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCa                   <  1:1523155/69‑1 (MQ=255)
                                              gccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCa                   <  1:1831790/69‑1 (MQ=255)
                                                  gcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATctg               <  1:2654947/69‑1 (MQ=255)
                                                  gcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATctg               <  1:1033587/69‑1 (MQ=255)
                                                                cTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATc  >  1:1358326/1‑68 (MQ=255)
                                                                cTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATc  >  1:1410214/1‑68 (MQ=255)
                                                                cTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATc  >  1:1735281/1‑68 (MQ=255)
                                                                  |                                                                 
GAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATC  >  minE/241280‑241411

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: