Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,396,026 C→T 8.4% I244I (ATC→ATT setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,396,0260CT 8.4% 210.1 / 7.6 71I244I (ATC→ATTsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base T (2/4);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCG  >  minE/1395959‑1396095
                                                                   |                                                                     
ttGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAAc                                                                    >  1:3388241/1‑71 (MQ=255)
 tGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAAc                                                                    <  1:809476/70‑1 (MQ=255)
   tCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGaa                                                                 >  1:2855214/1‑71 (MQ=255)
     aTTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACt                                                               >  1:646595/1‑71 (MQ=255)
     aTTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACt                                                               <  1:1159838/71‑1 (MQ=255)
     aTTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGATATTTATTATCAACGAACt                                                               >  1:1169728/1‑71 (MQ=255)
      tttGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTg                                                              >  1:1600476/1‑71 (MQ=255)
       ttGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTg                                                              >  1:803080/1‑70 (MQ=255)
        tGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACg                                                                   >  1:808224/1‑64 (MQ=255)
        tGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGc                                                             <  1:2249331/70‑1 (MQ=255)
         gTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCAt                                                           >  1:2009725/1‑71 (MQ=255)
          tACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCAt                                                           <  1:135227/70‑1 (MQ=255)
          tACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCAt                                                           <  1:3691508/70‑1 (MQ=255)
           aCATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCt                                                         >  1:1395095/1‑71 (MQ=255)
            caatGATGTGGGCGGCGAATAATCTCTACATGATAAATATGCCGCTATTTATTATTGATGAACTGCATCt                                                         >  1:2888310/4‑70 (MQ=255)
            caatGATGTGGGCGGCGAATAATCTCTACATGATAAATATGCCGCTATTTATTATTGATGAACTGCATCt                                                         >  1:1082490/4‑70 (MQ=255)
            cATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCt                                                         <  1:1044361/70‑1 (MQ=255)
             aTTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTc                                                       >  1:845882/1‑71 (MQ=255)
             aTTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTc                                                       <  1:1315720/71‑1 (MQ=255)
              ttGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAAc                                                                <  1:2934615/61‑1 (MQ=255)
               tGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATc                                                                       <  1:2244516/53‑1 (MQ=255)
               tGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAAc                                                                >  1:2326790/1‑60 (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:829317/71‑1 (MQ=255)
                        gCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAAAGAACTGCATCTTCCCGAGAAACTg                                            <  1:2930472/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)
                                     gTACATCATCAACATGCCGCTATTTATTATCAACg                                                                   <  1:3194278/35‑1 (MQ=255)
                                     gTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTgg                                           >  1:1235573/1‑59 (MQ=255)
                                     gTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTgg                                           >  1:74692/1‑59 (MQ=255)
                                      tACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTgatgat                              <  1:1041592/71‑1 (MQ=255)
                                      tACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTgatgat                              <  1:3755553/71‑1 (MQ=255)
                                        catcatCAACATGCCGCTATTTATTATCAACGAACTGc                                                             >  1:2327996/1‑38 (MQ=255)
                                        catcatCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGatgatg                             <  1:1046174/70‑1 (MQ=255)
                                        catcatCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATgg                            >  1:97344/1‑71 (MQ=255)
                                        catcatCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATgg                            <  1:1777031/71‑1 (MQ=255)
                                        catcatCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATgg                            >  1:2594023/1‑71 (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:3437316/71‑2 (MQ=255)
                                                      cGCTATTTATTATTGATGAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCggt              <  1:3109331/71‑2 (MQ=255)
                                                        cTATTTATTATTGATGAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCGGTCt            <  1:2095777/71‑1 (MQ=255)
                                                        cTATTTATTATTGATGAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCGGTCt            <  1:1336594/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:1376619/70‑1 (MQ=255)
                                                             tattatCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGaa        <  1:2062447/70‑1 (MQ=255)
                                                                  tCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAGATCCcg  <  1:1542191/71‑1 (MQ=255)
                                                                   |                                                                     
TTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCG  >  minE/1395959‑1396095

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: