Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,395,402 G→T 11.5% P36P (CCG→CCT setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,395,4020GT11.5% 218.3 / 9.3 69P36P (CCG→CCTsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base G (29/32);  new base T (4/4);  total (33/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.35e-01

TGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCA  >  minE/1395333‑1395472
                                                                     |                                                                      
tgacctgaccTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGa                                                                       >  1:112253/1‑71 (MQ=255)
   cctgacctCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGacac                                                                    >  1:3348220/1‑71 (MQ=255)
       acctCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAg                                                                <  1:2421926/71‑1 (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:182198/1‑71 (MQ=255)
                  gTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTa                                                     >  1:704505/1‑71 (MQ=255)
                  gTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTa                                                     <  1:2077634/71‑1 (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:350033/1‑47 (MQ=255)
                        aaTCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGa                                                                       >  1:2647527/1‑47 (MQ=255)
                        aaTCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGa                                                                       >  1:3798978/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)
                               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: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: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:2050547/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: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)
                                                                     |                                                                      
TGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCA  >  minE/1395333‑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: