Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,395,396 A→G 11.6% Q34Q (CAA→CAG setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,395,3960AG11.6% 181.8 / 13.5 69Q34Q (CAA→CAGsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base G (4/4);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

AATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTC  >  minE/1395326‑1395462
                                                                      |                                                                  
aaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCaa                                                                    <  1:648864/71‑1 (MQ=255)
  tCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAc                                                                  >  1:3296002/1‑71 (MQ=255)
    gTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAccc                                                                >  1:2012424/1‑71 (MQ=255)
    gTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAccc                                                                >  1:2887438/1‑71 (MQ=255)
       tgacctgaccTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGa                                                             >  1:112253/1‑71 (MQ=255)
          cctgacctCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGacac                                                          >  1:3348220/1‑71 (MQ=255)
              acctCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAg                                                      <  1:2421926/71‑1 (MQ=255)
               cctCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAccc                                                                >  1:3802522/1‑60 (MQ=255)
               cctCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAg                                                      >  1:1699721/1‑70 (MQ=255)
                 tCGACGGCGTTTTCAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGa                                                             <  1:178586/61‑1 (MQ=255)
                   gACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTAtt                                                  <  1:213680/70‑1 (MQ=255)
                     cGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTAtattc                                               >  1:491273/1‑70 (MQ=255)
                      ggCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTAtattcc                                              <  1:666666/71‑3 (MQ=255)
                       gcatttttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGt                                                     >  1:3578279/4‑63 (MQ=255)
                         gTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTa                                           <  1:2077634/71‑1 (MQ=255)
                         gTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTa                                           >  1:182198/1‑71 (MQ=255)
                         gTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTa                                           >  1:704505/1‑71 (MQ=255)
                          tttttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCACGACACTCAGTATTTTTCTTa                                           <  1:2188535/70‑1 (MQ=255)
                           ttttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCtcg                                           >  1:3049004/1‑67 (MQ=255)
                            tttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCaa                                                                    <  1:984101/43‑1 (MQ=255)
                            tttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCaa                                                                    <  1:1047866/43‑1 (MQ=255)
                             ttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGac                                                            >  1:1698816/1‑50 (MQ=255)
                               aaTCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGa                                                             >  1:350033/1‑47 (MQ=255)
                               aaTCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGa                                                             >  1:3798978/1‑47 (MQ=255)
                               aaTCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGa                                                             >  1:2647527/1‑47 (MQ=255)
                                    ttGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTAc                                          <  1:1059526/61‑1 (MQ=255)
                                    ttGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGaa                                   >  1:1456976/1‑68 (MQ=255)
                                     tGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTa                                >  1:2385743/1‑70 (MQ=255)
                                      gCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAcccc                                                               <  1:3530642/38‑1 (MQ=255)
                                      gCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTAc                               >  1:1693384/1‑70 (MQ=255)
                                       ccTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACAt                             <  1:2230401/71‑1 (MQ=255)
                                        cTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATg                            <  1:3904721/71‑1 (MQ=255)
                                         tttCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGa                                                             <  1:2548102/37‑1 (MQ=255)
                                         tttCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGa                                                             <  1:2382104/37‑1 (MQ=255)
                                                ccGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTa                                >  1:2163567/1‑59 (MQ=255)
                                                ccGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTa                                >  1:1866052/1‑59 (MQ=255)
                                                 gggTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTccta                   <  1:1933850/70‑3 (MQ=255)
                                                 gggTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTc                      <  1:1432052/67‑1 (MQ=255)
                                                 cGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTAc                               >  1:1954311/1‑59 (MQ=255)
                                                  ggTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTcctat                  <  1:3361635/71‑4 (MQ=255)
                                                    tATTGCGGGCGCTCAGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTa                                <  1:3766127/55‑1 (MQ=255)
                                                      ttGCGGGCGCTCTTCAGACTCCTACCCTAAGTAta                                                  >  1:1284793/1‑34 (MQ=255)
                                                      ttGCGGGCGCTCTGCAAACCCCGACACTCAGTAttttt                                               <  1:2068429/38‑1 (MQ=255)
                                                      ttGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTc                                              <  1:2622723/39‑1 (MQ=255)
                                                      ttGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtg              <  1:214826/71‑1 (MQ=255)
                                                      ttGCGGGCGCTCTGCAAACCACGACACTCAGTATTTTTCTTACCGATg                                     <  1:3534199/48‑1 (MQ=255)
                                                        gCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTc                      >  1:1289580/1‑61 (MQ=255)
                                                        gCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtgg             <  1:118666/70‑1 (MQ=255)
                                                        gCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtgg             <  1:3050815/70‑1 (MQ=255)
                                                           ggCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGaa                                   <  1:3137898/45‑1 (MQ=255)
                                                           ggCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtggtgg          <  1:2748935/70‑1 (MQ=255)
                                                            gcgcTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTcc                     <  1:1331092/58‑1 (MQ=255)
                                                             cgcTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAt       >  1:3148345/1‑71 (MQ=255)
                                                               ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGa                                       <  1:1116466/37‑1 (MQ=255)
                                                               ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc     <  1:2050547/71‑1 (MQ=255)
                                                               ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc     <  1:23755/71‑1 (MQ=255)
                                                               ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc     <  1:2890552/71‑1 (MQ=255)
                                                               ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc     <  1:2458615/71‑1 (MQ=255)
                                                               ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc     >  1:3719713/1‑71 (MQ=255)
                                                                tctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGc                           >  1:2068475/1‑48 (MQ=255)
                                                                tctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGa        <  1:2726225/67‑1 (MQ=255)
                                                                tctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGa        <  1:1563046/67‑1 (MQ=255)
                                                                tctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttct    >  1:1808239/1‑71 (MQ=255)
                                                                 ctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACAt                             >  1:239995/1‑45 (MQ=255)
                                                                  tGCAAACCCCGACACTCAGTATTTTTCTTACCGATg                                     <  1:2726315/36‑1 (MQ=255)
                                                                  tGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGt                                 >  1:1391204/1‑40 (MQ=255)
                                                                  tGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttcttc  >  1:3081368/1‑71 (MQ=255)
                                                                  tGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttct    <  1:453422/69‑1 (MQ=255)
                                                                  tGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttct    <  1:2067071/69‑1 (MQ=255)
                                                                    cAAACCGCGACACTCAGTATTTTTCTTACCGATGAAGTACATGc                           >  1:3908/1‑44 (MQ=255)
                                                                    cAAACCGCGACACTCAGTATTTTTCTTACCGATGAAGTACATGc                           >  1:8561/1‑44 (MQ=255)
                                                                    cAAACCCCGACACTCAGTATTTTTCTTACCGATGaa                                   >  1:578147/1‑36 (MQ=255)
                                                                     aaaCCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttct    <  1:1306072/66‑1 (MQ=255)
                                                                     aaaCCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttct    <  1:224769/66‑1 (MQ=255)
                                                                      |                                                                  
AATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTC  >  minE/1395326‑1395462

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: