Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,396,047 G→T 11.5% E251D (GAG→GAT setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,396,0470GT11.5% 161.6 / 5.3 52E251D (GAG→GATsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base G (25/21);  new base T (2/4);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATT  >  minE/1395979‑1396117
                                                                    |                                                                      
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)
    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)
          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:3319209/1‑70 (MQ=255)
                tGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGtgat                                                       >  1:3267809/1‑70 (MQ=255)
                 gTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTgg                                                                 >  1:1235573/1‑59 (MQ=255)
                 gTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTgg                                                                 >  1:74692/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)
                           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: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:1376619/70‑1 (MQ=255)
                                         tattatCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGaa                              <  1:2062447/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)
                                                                    |                                                                      
TGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATT  >  minE/1395979‑1396117

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