Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,396,041 T→A 9.8% L249L (CTT→CTA setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,396,0410TA 9.8% 163.6 / 9.6 61L249L (CTT→CTAsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base T (27/28);  new base A (2/4);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

ATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGC  >  minE/1395972‑1396109
                                                                     |                                                                    
aTTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTc                                                                     >  1:845882/1‑71 (MQ=255)
aTTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTc                                                                     <  1:1315720/71‑1 (MQ=255)
       tgGGGCTCGAACAGCCTGTATATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCgaga                                                               <  1:248083/70‑1 (MQ=255)
       tgGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCgaga                                                               >  1:1443450/1‑70 (MQ=255)
       tgGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGaa                                                              <  1:685368/71‑1 (MQ=255)
         gggCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAAc                                                            >  1:144944/1‑71 (MQ=255)
          ggCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTcc                                                                    <  1:296029/62‑1 (MQ=255)
           gCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAAAGAACTGCATCTTCCCGAGAAACTg                                                          <  1:2930472/71‑1 (MQ=255)
           gCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAAAGAACTGCATCTTCCCGAGAAACTg                                                          <  1:829317/71‑1 (MQ=255)
             tCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGc                                                        >  1:1492950/1‑71 (MQ=255)
               gAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCg                                                      >  1:846309/1‑71 (MQ=255)
                aaCAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTccc                                                                   <  1:3022143/57‑1 (MQ=255)
                 aCAGCCAGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGgt                                                    >  1:505103/1‑71 (MQ=255)
                   aGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGgtgt                                                  <  1:1420683/71‑1 (MQ=255)
                     ccTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCgag                                                                <  1:905819/55‑1 (MQ=255)
                      cTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGgt                                                    >  1:3257240/1‑66 (MQ=255)
                      cTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGtgat                                               <  1:3587145/71‑1 (MQ=255)
                       tGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCtag                                                                <  1:1943649/53‑1 (MQ=255)
                       tGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGtgat                                               >  1:2801997/1‑70 (MQ=255)
                       tGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGtgat                                               >  1:3267809/1‑70 (MQ=255)
                       tGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGtgat                                               >  1:3319209/1‑70 (MQ=255)
                        gTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTgg                                                         >  1:74692/1‑59 (MQ=255)
                        gTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTgg                                                         >  1:1235573/1‑59 (MQ=255)
                         tACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTgatgat                                            <  1:3755553/71‑1 (MQ=255)
                         tACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTgatgat                                            <  1:1041592/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)
                               atcaACATGCCGCTATTTATTATCAACGAACTGCATCTTccc                                                                   <  1:2780872/42‑1 (MQ=255)
                                  aaCATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTg                                                          >  1:3747023/1‑48 (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)
                                                                     |                                                                    
ATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGC  >  minE/1395972‑1396109

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: