Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,395,393 G→T 11.2% L33L (CTG→CTT setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,395,3930GT11.2% 221.1 / 9.2 71L33L (CTG→CTTsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base G (29/34);  new base T (4/4);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

AAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTC  >  minE/1395325‑1395462
                                                                    |                                                                     
aaaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGc                                                                      >  1:980850/1‑70 (MQ=255)
 aaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCaa                                                                    <  1:648864/71‑1 (MQ=255)
 aaTCGTTTGAACTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCa                                                                     >  1:1825302/1‑70 (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:2382104/37‑1 (MQ=255)
                                          tttCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGa                                                             <  1:2548102/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: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:2890552/71‑1 (MQ=255)
                                                                ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc     <  1:2050547/71‑1 (MQ=255)
                                                                ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc     <  1:2458615/71‑1 (MQ=255)
                                                                ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc     <  1:23755/71‑1 (MQ=255)
                                                                 tctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGc                           >  1:2068475/1‑48 (MQ=255)
                                                                 tctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGa        <  1:1563046/67‑1 (MQ=255)
                                                                 tctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGa        <  1:2726225/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)
                                                                    |                                                                     
AAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTC  >  minE/1395325‑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: