Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,341,659 T→A 100% T115T (ACA→ACT wcaA ← predicted glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,341,6590TA77.2% 54.7 / 12.8 35T115T (ACA→ACTwcaApredicted glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (7/1);  new base A (19/8);  total (26/9)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTC  >  minE/1341616‑1341718
                                           |                                                           
cTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCg                                     >  1:111669/1‑68 (MQ=255)
cTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCg                                     >  1:2488445/1‑68 (MQ=255)
cTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCg                                     >  1:295009/1‑68 (MQ=255)
          gCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACg                            >  1:1522388/1‑67 (MQ=255)
          gCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACg                            >  1:2262721/1‑67 (MQ=255)
          gCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACg                            >  1:2364468/1‑67 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGtt                 >  1:1371967/1‑67 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg               >  1:935436/1‑69 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg               >  1:347257/1‑69 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg               >  1:3411165/1‑69 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg               >  1:3242030/1‑69 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg               >  1:314112/1‑69 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg               >  1:2988958/1‑69 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg               >  1:960683/1‑69 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg               >  1:2897721/1‑69 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg               >  1:2895094/1‑69 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg               >  1:2679188/1‑69 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg               >  1:989384/1‑69 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg               >  1:1447229/1‑69 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg               >  1:1699186/1‑69 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg               >  1:2035548/1‑69 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg               >  1:2034128/1‑69 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTg                >  1:347211/1‑68 (MQ=255)
                     cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTAgg               >  1:2452965/1‑69 (MQ=255)
                     cGTTAGCGTACAAAAAGCCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg               >  1:2587561/1‑69 (MQ=255)
                            gTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGGGTCCa        <  1:2523805/69‑1 (MQ=255)
                            gTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTGGGTGTcc         <  1:3197606/68‑1 (MQ=255)
                            gTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTGGGTGTcc         <  1:12258/68‑1 (MQ=255)
                            gTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTGGGTGTcc         <  1:145253/68‑1 (MQ=255)
                            gTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTGGGTGTcc         <  1:2636382/68‑1 (MQ=255)
                            gTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTGGGTGTcc         <  1:2372719/68‑1 (MQ=255)
                            gTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTGGGTGTcc         <  1:2646389/68‑1 (MQ=255)
                            gTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTACGGGCGAGGAAGACGCTCAGGCGGTTGGGTGTcc         <  1:36516/68‑1 (MQ=255)
                               cAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtc     >  1:2969983/1‑69 (MQ=255)
                                     ggCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTGGGTGTcc         <  1:2886680/59‑1 (MQ=255)
                                         tgtgtgACCAGTTGCTGTTTATTGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgt   >  1:630451/1‑61 (MQ=255)
                                         tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTTTCCATtcgtc  >  1:452599/1‑62 (MQ=255)
                                         tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtc  >  1:1011528/1‑62 (MQ=255)
                                         tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtc  >  1:972908/1‑62 (MQ=255)
                                         tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtc  >  1:2712766/1‑62 (MQ=255)
                                         tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtc  >  1:3386243/1‑62 (MQ=255)
                                         tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtc  >  1:3357022/1‑62 (MQ=255)
                                         tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtc  >  1:2300057/1‑62 (MQ=255)
                                         tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtc  >  1:1815452/1‑62 (MQ=255)
                                         tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtc  >  1:3087243/1‑62 (MQ=255)
                                         tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGATCAGACGGTTGGGAGTCCATtcgtc  >  1:165489/1‑62 (MQ=255)
                                           |                                                           
CTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTC  >  minE/1341616‑1341718

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: