Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,395,408 C→A 12.6% L38L (CTC→CTA setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,395,4080CA12.6% 195.1 / 9.8 63L38L (CTC→CTAsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base C (24/31);  new base A (4/4);  total (28/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 = 1.00e+00

ACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAGCGCTG  >  minE/1395340‑1395478
                                                                    |                                                                      
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:2077634/71‑1 (MQ=255)
           gTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTa                                                           >  1:704505/1‑71 (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: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:3050815/70‑1 (MQ=255)
                                          gCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtgg                             <  1:118666/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:2050547/71‑1 (MQ=255)
                                                 ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc                     <  1:2458615/71‑1 (MQ=255)
                                                 ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc                     <  1:2890552/71‑1 (MQ=255)
                                                 ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc                     <  1:23755/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: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: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)
                                                                  cTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTcc                                     <  1:3320866/38‑1 (MQ=255)
                                                                    cAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAGCGCTg  >  1:3149683/1‑71 (MQ=255)
                                                                    |                                                                      
ACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAGCGCTG  >  minE/1395340‑1395478

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: