Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,396,059 T→A 13.3% G255G (GGT→GGA setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,396,0590TA13.3% 111.1 / 10.5 45G255G (GGT→GGAsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base T (23/16);  new base A (2/4);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

ACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAAC  >  minE/1395989‑1396126
                                                                      |                                                                   
aCAGCCAGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGgt                                                                     >  1:505103/1‑71 (MQ=255)
  aGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGgtgt                                                                   <  1:1420683/71‑1 (MQ=255)
     cTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGgt                                                                     >  1:3257240/1‑66 (MQ=255)
     cTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGtgat                                                                <  1:3587145/71‑1 (MQ=255)
      tGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGtgat                                                                >  1:3319209/1‑70 (MQ=255)
      tGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGtgat                                                                >  1:3267809/1‑70 (MQ=255)
      tGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGtgat                                                                >  1:2801997/1‑70 (MQ=255)
        tACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTgatgat                                                             <  1:1041592/71‑1 (MQ=255)
        tACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTgatgat                                                             <  1:3755553/71‑1 (MQ=255)
          catcatCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGatgatg                                                            <  1:1046174/70‑1 (MQ=255)
          catcatCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATgg                                                           >  1:97344/1‑71 (MQ=255)
          catcatCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATgg                                                           >  1:2594023/1‑71 (MQ=255)
          catcatCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATgg                                                           <  1:1777031/71‑1 (MQ=255)
           atcatcAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTTATGATggg                                                          <  1:242835/71‑1 (MQ=255)
           atcatcAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATggg                                                          <  1:590192/71‑1 (MQ=255)
            tcatcaACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGgt                                                                     <  1:1497129/59‑1 (MQ=255)
            tcatcaACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATggg                                                          >  1:1562378/1‑70 (MQ=255)
             catcaACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGgtg                                                                    >  1:1169678/1‑59 (MQ=255)
                   cATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGTGAccgcc                                                   >  1:606880/1‑70 (MQ=255)
                     tGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGtga                                                                 <  1:3938387/54‑1 (MQ=255)
                      gCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCAGGTGTgatgat                                                             <  1:3237665/57‑1 (MQ=255)
                        cGCTATTTATTATTGATGAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCggt                                             <  1:3109331/71‑2 (MQ=255)
                        cGCTATTTATTATTGATGAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCggt                                             <  1:3437316/71‑2 (MQ=255)
                          cTATTTATTATTGATGAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCGGTCt                                           <  1:1336594/71‑1 (MQ=255)
                          cTATTTATTATTGATGAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCGGTCt                                           <  1:2095777/71‑1 (MQ=255)
                             tttattATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATgggg                                                         <  1:3233829/54‑1 (MQ=255)
                             tttattATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTgg                                         >  1:1163031/1‑70 (MQ=255)
                               tattatCCACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGAccgccgc                                                 >  1:1807884/1‑60 (MQ=255)
                               tattatCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGatgatg                                                            <  1:1271501/49‑1 (MQ=255)
                               tattatCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGaa                                       <  1:2062447/70‑1 (MQ=255)
                               tattatCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGaa                                       <  1:1376619/70‑1 (MQ=255)
                                    tCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAGATCCcg                                 <  1:1542191/71‑1 (MQ=255)
                                       atgAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTAtcgctgccg                                               >  1:2349334/3‑54 (MQ=255)
                                       aCGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGAccg                                                     >  1:116538/1‑48 (MQ=255)
                                        cGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCTcgacga                             <  1:506966/71‑1 (MQ=255)
                                         gAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCGGTCTGGAAATTCCgatgat                            >  1:1785542/1‑71 (MQ=255)
                                           aCTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGt                          >  1:2811283/1‑71 (MQ=255)
                                              gCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGa                       >  1:458160/1‑71 (MQ=255)
                                               cATCTTCCCGAGAAACTGGCCGGTGTGATGATggg                                                          <  1:602895/35‑1 (MQ=255)
                                               cATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTg                        >  1:2621546/1‑69 (MQ=255)
                                                 tCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGA‑‑‑CCGCCGGGCTGGAAATcc                                   >  1:3070076/1‑53 (MQ=255)
                                                  cTTCCCGAGAAACTGGCCGGTGTGATGATGGGGAc                                                       >  1:2547902/1‑35 (MQ=255)
                                                  cTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGc                   >  1:2723217/1‑71 (MQ=255)
                                                         agaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATc                                    >  1:929507/1‑47 (MQ=255)
                                                          gaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGatat            >  1:3039409/1‑70 (MQ=255)
                                                          gaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATAtt           >  1:2691213/1‑71 (MQ=255)
                                                          gaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATAtt           >  1:2215856/1‑71 (MQ=255)
                                                               tGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGc       >  1:2007992/1‑70 (MQ=255)
                                                               tGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGc       >  1:954128/1‑70 (MQ=255)
                                                                   cGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAAc  >  1:1135815/1‑71 (MQ=255)
                                                                      tgtgATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGa               >  1:3357020/1‑55 (MQ=255)
                                                                      tgtgATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGa               >  1:557733/1‑55 (MQ=255)
                                                                      |                                                                   
ACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAAC  >  minE/1395989‑1396126

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: