Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,395,405 A→C 13.1% T37T (ACA→ACC setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,395,4050AC13.1% 126.3 / 18.2 61T37T (ACA→ACCsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base A (23/30);  new base C (4/4);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

CCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCA  >  minE/1395336‑1395472
                                                                     |                                                                   
cctgacctCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGacac                                                                    >  1:3348220/1‑71 (MQ=255)
    acctCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAg                                                                <  1:2421926/71‑1 (MQ=255)
     cctCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAg                                                                >  1:1699721/1‑70 (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:704505/1‑71 (MQ=255)
               gTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTa                                                     <  1:2077634/71‑1 (MQ=255)
               gTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTa                                                     >  1:182198/1‑71 (MQ=255)
                tttttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCACGACACTCAGTATTTTTCTTa                                                     <  1:2188535/70‑1 (MQ=255)
                 ttttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCtcg                                                     >  1:3049004/1‑67 (MQ=255)
                          ttGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTAc                                                    <  1:1059526/61‑1 (MQ=255)
                          ttGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGaa                                             >  1:1456976/1‑68 (MQ=255)
                           tGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTa                                          >  1:2385743/1‑70 (MQ=255)
                            gCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTAc                                         >  1:1693384/1‑70 (MQ=255)
                             ccTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACAt                                       <  1:2230401/71‑1 (MQ=255)
                              cTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATg                                      <  1:3904721/71‑1 (MQ=255)
                                      ccGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTa                                          >  1:1866052/1‑59 (MQ=255)
                                      ccGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTa                                          >  1:2163567/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:3719713/1‑71 (MQ=255)
                                                     ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc               <  1:2890552/71‑1 (MQ=255)
                                                     ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc               <  1:23755/71‑1 (MQ=255)
                                                     ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc               <  1:2050547/71‑1 (MQ=255)
                                                     ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc               <  1:2458615/71‑1 (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:2067071/69‑1 (MQ=255)
                                                        tGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttct              <  1:453422/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:224769/66‑1 (MQ=255)
                                                           aaaCCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttct              <  1:1306072/66‑1 (MQ=255)
                                                             aCCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTAc       >  1:3405782/1‑71 (MQ=255)
                                                             aCCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTAc       >  1:2556724/1‑71 (MQ=255)
                                                                 cGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCgg    >  1:604785/1‑70 (MQ=255)
                                                                  gACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCa  <  1:1698872/71‑1 (MQ=255)
                                                                   acacTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTa        >  1:1237545/1‑64 (MQ=255)
                                                                     |                                                                   
CCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCA  >  minE/1395336‑1395472

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: