Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,395,399 C→T 11.7% T35T (ACC→ACT setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,395,3990CT11.7% 184.8 / 12.8 68T35T (ACC→ACTsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base T (4/4);  total (32/36)
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.28e-01

GTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTAC  >  minE/1395330‑1395467
                                                                     |                                                                    
gTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAccc                                                                     >  1:2887438/1‑71 (MQ=255)
gTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAccc                                                                     >  1:2012424/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: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)
                         ttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGac                                                                 >  1:1698816/1‑50 (MQ=255)
                           aaTCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGa                                                                  >  1:3798978/1‑47 (MQ=255)
                           aaTCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGa                                                                  >  1:2647527/1‑47 (MQ=255)
                           aaTCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGa                                                                  >  1:350033/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:2890552/71‑1 (MQ=255)
                                                           ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc          <  1:2458615/71‑1 (MQ=255)
                                                           ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc          <  1:23755/71‑1 (MQ=255)
                                                           ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc          >  1:3719713/1‑71 (MQ=255)
                                                           ctctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc          <  1:2050547/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)
                                                                     |                                                                    
GTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTAC  >  minE/1395330‑1395467

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: