Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,400,543 T→A 54.6% P86P (CCT→CCA yeiU → undecaprenyl pyrophosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,400,5430TA54.6% 12.8 / 71.8 66P86P (CCT→CCAyeiUundecaprenyl pyrophosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base T (26/4);  new base A (21/15);  total (47/19)
Fisher's exact test for biased strand distribution p-value = 1.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

TGTTCACTGCTGGCGATGGGTATGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTATTAA  >  minE/1400484‑1400601
                                                           |                                                          
tgttCACTGCTGGCGATGGGTATGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAgacga                                                   >  1:813228/1‑69 (MQ=255)
tgttCACTGCTGGCGATGGGTATGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAgacga                                                   >  1:446329/1‑69 (MQ=255)
tgttCACTGCTGGCGATGGGTATGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAgacga                                                   >  1:34814/1‑69 (MQ=255)
tgttCACTGCTGGCGATGGGTATGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAgacga                                                   >  1:3232775/1‑69 (MQ=255)
tgttCACTGCTGGCGATGGGTATGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAgacga                                                   >  1:1730065/1‑69 (MQ=255)
tgttCACTGCTGGCGATGGGTATGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAgacga                                                   >  1:1648086/1‑69 (MQ=255)
tgttCACTGCTGGCGATGGGTATGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAgacg                                                    >  1:2924579/1‑68 (MQ=255)
       tgctgGGGATGGGTATGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATc                                             >  1:238126/1‑68 (MQ=255)
       tgctgGCGATGGGTATGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATc                                             >  1:1353115/1‑68 (MQ=255)
       tgctgGCGATGGGTATGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATc                                             >  1:1762548/1‑68 (MQ=255)
           ggCGATGGGTATGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTg                                          <  1:3059179/67‑1 (MQ=255)
                     aTGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTGATTATtggtctg                              >  1:991979/1‑69 (MQ=255)
                     aTGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTGATTATtggtctg                              >  1:2830177/1‑69 (MQ=255)
                     aTGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTGATTATtggtctg                              >  1:3018066/1‑69 (MQ=255)
                     aTGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTGATTATtggtctg                              >  1:376333/1‑69 (MQ=255)
                     aTGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTGATTATtggtctg                              >  1:1445183/1‑69 (MQ=255)
                     aTGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTGATTATtggtctg                              >  1:463663/1‑69 (MQ=255)
                     aTGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTGATTATtggtctg                              >  1:556782/1‑69 (MQ=255)
                     aTGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTGATTATtggtctg                              >  1:901404/1‑69 (MQ=255)
                     aTGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTGATTATtggtct                               >  1:223389/1‑68 (MQ=255)
                       gTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTGATTATtggtctggt                            >  1:3186283/1‑69 (MQ=255)
                        ttGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTGat                                        >  1:2651483/1‑56 (MQ=255)
                                gAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTAtt                   <  1:1221278/69‑1 (MQ=255)
                                gAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTAtt                   <  1:255399/69‑1 (MQ=255)
                                   tttCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAAc                >  1:3108491/1‑69 (MQ=255)
                                       tGGCTGAAAGAAAACGCCCCAGGCAGACGACGTATCGTtattat                                     <  1:1438834/44‑1 (MQ=38)
                                       tGGCTGAAAGAAAACGCCCCAGGCAGACGACGTATCGTGattat                                     <  1:1596360/44‑1 (MQ=255)
                                       tGGCTGAAAGAAAACGCCCCAGGCAGACGACGTATCGTGattat                                     <  1:1087712/44‑1 (MQ=255)
                                       tGGCTGAAAGAAAACGCCCCAGGCAGACGACGTATCGTGattat                                     <  1:1494453/44‑1 (MQ=255)
                                       tGGCTGAAAGAAAACGCCCCAGGCAGACGACGTATCGTGattat                                     <  1:3422718/44‑1 (MQ=255)
                                       tGGCTGAAAGAAAACGCCCCAGGCAGACGACGTATCGTGattat                                     <  1:152932/44‑1 (MQ=255)
                                       tGGCTGAAAGAAAACGCCCCAGGCAGACGACGTATCGTGattat                                     <  1:3003860/44‑1 (MQ=255)
                                       tGGCTGAAAGAAAACGCCCCAGGCAGACGACGTATCGTGattat                                     <  1:2928773/44‑1 (MQ=255)
                                       tGGCTGAAAGAAAACGCCCCAGGCAGACGACGTATCGTGattat                                     <  1:1949802/44‑1 (MQ=255)
                                       tGGCTGAAAGAAAACGCCCCAGGCAGACGACGTATCGTGattat                                     <  1:2738290/44‑1 (MQ=255)
                                       tGGCTGAAAGAAAACGCCCCAGGCAGACGACGTATCGTGattat                                     <  1:1047542/44‑1 (MQ=255)
                                            gaaagaaaACGCCCCAGGCAGACGACGTATCGTGattat                                     <  1:2723192/39‑1 (MQ=255)
                                             aaagaaaACGCCCCAGGCAGACGACGTATCGTGattat                                     <  1:1740368/38‑1 (MQ=255)
                                                gaaaACGCCCCTGGCAGACGACGTATCGTGATTATtggtctggt                            <  1:2263585/44‑1 (MQ=255)
                                                gaaaACGCCCCAGGCAGCCGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:318812/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGattat                                     <  1:979338/35‑1 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGattat                                     <  1:886631/35‑1 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGTTCATGCTATTAACTGCAGTGGTAt     >  1:1788535/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:2171357/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:2429666/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:456437/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:1286101/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:463730/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:588307/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:2304721/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:839511/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:1136256/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:1468376/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:2410227/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:3362627/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:2197302/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:2488713/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:2563268/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:2599959/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:283202/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTAt     >  1:2106642/1‑67 (MQ=255)
                                                gaaaACGCCCCAGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGAAt     >  1:125169/1‑67 (MQ=255)
                                                   aaCGCCCCTGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTATTaa  >  1:3373381/1‑67 (MQ=255)
                                                   aaCGCCCCTGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTATTaa  >  1:3260770/1‑67 (MQ=255)
                                                   aaCGCCCCTGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTATTaa  >  1:60571/1‑67 (MQ=255)
                                                   aaCGACCCTGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTATTaa  >  1:2459089/1‑67 (MQ=255)
                                                           |                                                          
TGTTCACTGCTGGCGATGGGTATGTTGATGCTGAGTTTCTGGCTGAAAGAAAACGCCCCTGGCAGACGACGTATCGTGATTATTGGTCTGGTCATGCTATTAACTGCAGTGGTATTAA  >  minE/1400484‑1400601

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: