Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,341,698 A→G 36.9% R102R (CGT→CGC wcaA ← predicted glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,341,6980AG36.9% 39.3 / 44.6 73R102R (CGT→CGCwcaApredicted glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (26/20);  new base G (19/8);  total (45/28)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

CGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGATATATTCCCCTTGC  >  minE/1341637‑1341752
                                                             |                                                      
cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGtt                                                   >  1:1371967/1‑67 (MQ=255)
cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg                                                 >  1:3242030/1‑69 (MQ=255)
cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg                                                 >  1:2679188/1‑69 (MQ=255)
cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg                                                 >  1:2895094/1‑69 (MQ=255)
cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg                                                 >  1:2897721/1‑69 (MQ=255)
cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg                                                 >  1:2035548/1‑69 (MQ=255)
cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg                                                 >  1:2034128/1‑69 (MQ=255)
cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg                                                 >  1:2988958/1‑69 (MQ=255)
cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg                                                 >  1:314112/1‑69 (MQ=255)
cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg                                                 >  1:3411165/1‑69 (MQ=255)
cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg                                                 >  1:1699186/1‑69 (MQ=255)
cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg                                                 >  1:347257/1‑69 (MQ=255)
cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg                                                 >  1:935436/1‑69 (MQ=255)
cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg                                                 >  1:989384/1‑69 (MQ=255)
cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg                                                 >  1:960683/1‑69 (MQ=255)
cGTTAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTgg                                                 >  1:1447229/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:12258/68‑1 (MQ=255)
       gTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTGGGTGTcc                                           <  1:145253/68‑1 (MQ=255)
       gTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTGGGTGTcc                                           <  1:2646389/68‑1 (MQ=255)
       gTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTGGGTGTcc                                           <  1:2636382/68‑1 (MQ=255)
       gTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTGGGTGTcc                                           <  1:2372719/68‑1 (MQ=255)
       gTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTGGGTGTcc                                           <  1:3197606/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:972908/1‑62 (MQ=255)
                    tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtc                                    >  1:3087243/1‑62 (MQ=255)
                    tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtc                                    >  1:2712766/1‑62 (MQ=255)
                    tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtc                                    >  1:3357022/1‑62 (MQ=255)
                    tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtc                                    >  1:3386243/1‑62 (MQ=255)
                    tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtc                                    >  1:1011528/1‑62 (MQ=255)
                    tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtc                                    >  1:2300057/1‑62 (MQ=255)
                    tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtc                                    >  1:1815452/1‑62 (MQ=255)
                    tgtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGATCAGACGGTTGGGAGTCCATtcgtc                                    >  1:165489/1‑62 (MQ=255)
                          aCCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgt                                     >  1:2045563/1‑55 (MQ=255)
                          aCCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgt                                     >  1:2003202/1‑55 (MQ=255)
                          aCCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgt                                     >  1:936255/1‑55 (MQ=255)
                          aCCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgt                                     >  1:646910/1‑55 (MQ=255)
                          aCCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgt                                     >  1:529970/1‑55 (MQ=255)
                          aCCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgt                                     >  1:1607085/1‑55 (MQ=255)
                          aCCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgt                                     >  1:1901758/1‑55 (MQ=255)
                          aCCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCAAtcgt                                     >  1:2650620/1‑55 (MQ=255)
                                tGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTg                 <  1:1977865/69‑1 (MQ=255)
                                tGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTg                 <  1:453922/69‑1 (MQ=255)
                                tGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTg                 <  1:49371/69‑1 (MQ=255)
                                 gCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGgtgta                                            <  1:1256483/41‑2 (MQ=255)
                                 gCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTc                                            <  1:1062394/41‑1 (MQ=255)
                                 gCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTc                                            <  1:907011/41‑1 (MQ=255)
                                 gCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTc                                            <  1:3020485/41‑1 (MQ=255)
                                 gCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTc                                            <  1:1971455/41‑1 (MQ=255)
                                 gCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGa                <  1:1848446/69‑1 (MQ=255)
                                 gCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGa                <  1:1046048/69‑1 (MQ=255)
                                 gCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGa                <  1:872098/69‑1 (MQ=255)
                                 gCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGa                <  1:430835/69‑1 (MQ=255)
                                 gCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGa                <  1:1831644/69‑1 (MQ=255)
                                 gCTGTTTATGGGCGAGGAAGACGCTCAGACGGGTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGa                <  1:514618/69‑1 (MQ=255)
                                       tATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGatatat           <  1:1364812/68‑1 (MQ=255)
                                       tATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGatatat           <  1:658621/68‑1 (MQ=255)
                                       tATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGatatat           <  1:1719137/68‑1 (MQ=255)
                                       tATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGatatat           <  1:333489/68‑1 (MQ=255)
                                       tATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGATATAtt          >  1:1349685/1‑69 (MQ=255)
                                       tATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGATATAtt          >  1:3204552/1‑69 (MQ=255)
                                       tATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGATATAtt          >  1:102471/1‑69 (MQ=255)
                                       tATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATCCGTCATCGTCATCAATCCCGGTGATATAtt          >  1:1895783/1‑69 (MQ=255)
                                           ggCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGATATATTccc       <  1:1932441/68‑1 (MQ=255)
                                                 gAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGATATATTCCCCTTgc  >  1:661292/1‑67 (MQ=255)
                                                 gAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGATATATTCCCCTTgc  >  1:977093/1‑67 (MQ=255)
                                                             |                                                      
CGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCAATCCCGGTGATATATTCCCCTTGC  >  minE/1341637‑1341752

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: