Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,395,414 T→A 10.6% I40I (ATT→ATA setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,395,4140TA10.6% 176.7 / 12.1 66I40I (ATT→ATAsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base A (2/5);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

GACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAGCGCTGTCATT  >  minE/1395345‑1395483
                                                                     |                                                                     
gACGGCGTTTATAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTAttt                                                                      >  1:345318/1‑71 (MQ=255)
  cGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTAtattc                                                                    <  1:1329297/71‑2 (MQ=255)
    gcatttttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCt                                                                  <  1:5433845/68‑1 (MQ=255)
    gcatttttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCt                                                                  <  1:433957/68‑1 (MQ=255)
      gTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTa                                                                >  1:3885576/1‑71 (MQ=255)
      atttttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCtcg                                                                >  1:3109744/2‑69 (MQ=255)
        ttttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTAcc                                                              >  1:2495962/1‑71 (MQ=255)
        ttttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTAcc                                                              >  1:5049856/1‑71 (MQ=255)
         tttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTAcc                                                              >  1:4926033/1‑70 (MQ=255)
          ttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGa                                                            >  1:4581485/1‑71 (MQ=255)
             aTCGTTGCCTTTCTCACCGGTATTGCGGGTGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATg                                                          <  1:2188262/70‑1 (MQ=255)
             aTCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGa                                                         <  1:4649244/71‑1 (MQ=255)
              tCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCtt                                                                 <  1:5064535/62‑1 (MQ=255)
                 ttGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTa                                                     <  1:3752649/71‑1 (MQ=255)
                 ttGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTa                                                     <  1:5577743/71‑1 (MQ=255)
                   gCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGt                                                      >  1:5173362/1‑68 (MQ=255)
                    ccTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACAt                                                  <  1:3398826/71‑1 (MQ=255)
                       ttCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGaa                                                        >  1:1215099/1‑62 (MQ=255)
                       ttCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGcc                                               >  1:796732/1‑71 (MQ=255)
                       ttCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGcc                                               >  1:4008181/1‑71 (MQ=255)
                        tctcACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTGTTACCGATGAAGTACATGccc                                              >  1:2692158/1‑71 (MQ=255)
                        tctcACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGccc                                              >  1:1987298/1‑71 (MQ=255)
                         ctcACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGccc                                              >  1:557475/1‑70 (MQ=255)
                            aCCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTc                                           <  1:4786186/70‑1 (MQ=255)
                                 tATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGa                                                         <  1:5633364/51‑1 (MQ=255)
                                 tATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGa                                                         <  1:2040793/51‑1 (MQ=255)
                                  aTTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTcctat                                       <  1:4691839/68‑4 (MQ=255)
                                  aTTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTcctat                                       <  1:1824549/68‑4 (MQ=255)
                                  aTTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTcctat                                       <  1:5156775/68‑4 (MQ=255)
                                   ttGCGGGCGCTCTGCAAACCCCGACACTCAGTAttt                                                                      <  1:1616264/36‑1 (MQ=255)
                                   ttGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAg                                                       <  1:4564215/51‑1 (MQ=255)
                                   ttGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtg                                   <  1:3216434/71‑1 (MQ=255)
                                   ttGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtg                                   <  1:4153850/71‑1 (MQ=255)
                                   ttGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACGATGAAGTAc                                                    >  1:261351/1‑53 (MQ=255)
                                     gCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTa                                                                >  1:3901986/1‑40 (MQ=255)
                                     gCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTa                                                                >  1:4262473/1‑40 (MQ=255)
                                     gCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtggt                                 >  1:4023571/1‑71 (MQ=255)
                                      cGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtggtg                                >  1:4463386/1‑71 (MQ=255)
                                       gggCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtggtgg                               <  1:2389281/71‑1 (MQ=255)
                                          cgcTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAt                            >  1:5515954/1‑71 (MQ=255)
                                              ctGCAAACCCCGACACTCAGTATTTTTCTTACCGATg                                                          <  1:4276862/37‑1 (MQ=255)
                                              ctGCAAACCCCGACACTCAGTATTTTTCTTACCGATg                                                          <  1:3362466/37‑1 (MQ=255)
                                              ctGCAAACCCCGACACTCAGTATTTTTCTTACCGATg                                                          <  1:2538815/37‑1 (MQ=255)
                                              ctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGcc                                               >  1:4824744/1‑48 (MQ=255)
                                              ctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtgg                                  <  1:5537199/61‑1 (MQ=255)
                                              ctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtgg                                  <  1:2533110/61‑1 (MQ=255)
                                               tGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttcttc                       >  1:2153550/1‑71 (MQ=255)
                                               tGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttcttc                       >  1:459481/1‑71 (MQ=255)
                                               tGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttcttc                       >  1:5036779/1‑71 (MQ=255)
                                               tGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttct                         <  1:5132943/69‑1 (MQ=255)
                                                gCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAtt                           >  1:2390881/1‑66 (MQ=255)
                                                 cAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACa                                                   >  1:2636017/1‑41 (MQ=255)
                                                 cAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACa                                                   >  1:2365508/1‑41 (MQ=255)
                                                 cAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttctt                        <  1:3241071/68‑1 (MQ=255)
                                                 cAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTcttctt                     >  1:696417/1‑71 (MQ=255)
                                                 cAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTcttctt                     >  1:2719088/1‑71 (MQ=255)
                                                  aaaCCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATtcttct                      <  1:1268916/69‑1 (MQ=255)
                                                    aCCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc                          >  1:4026626/1‑63 (MQ=255)
                                                      cccGACACTCAGTATTTTTCTTACCGATGAAGTACa                                                   >  1:1671558/1‑36 (MQ=255)
                                                       ccGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCgg               <  1:1553053/71‑1 (MQ=255)
                                                            acTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGa                                     >  1:2832268/1‑44 (MQ=255)
                                                            acTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAgcg          <  1:4657633/71‑1 (MQ=255)
                                                               cAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCg                                      <  1:1141473/40‑1 (MQ=255)
                                                                aGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAGCGCTg       >  1:3821947/1‑70 (MQ=255)
                                                                  tATTTTTCTTACCGATGAAGTACATGCCCGTCCGGc                                       <  1:5234548/36‑1 (MQ=255)
                                                                  tATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAt                                    <  1:4165423/39‑1 (MQ=255)
                                                                  tATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAGCGCTGTCa    <  1:2977169/71‑1 (MQ=255)
                                                                  tATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAGCGCTGTCa    >  1:1913715/1‑71 (MQ=255)
                                                                   aTTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAGCGCTGTCa    <  1:1774201/70‑1 (MQ=255)
                                                                    tttttCTTACCGATGAAGTACATGCCCGTCCGGCGa                                     >  1:2240864/1‑36 (MQ=255)
                                                                    tttttCTTACCGATGAAGTACATGCCCGTCCGGCGa                                     >  1:4204129/1‑36 (MQ=255)
                                                                    tttttCTTACCGATGAAGTACATGCCCGTCCGGCGa                                     >  1:2495543/1‑36 (MQ=255)
                                                                    tttttCTTACCGATGAAGTACATGCCCGTCCGGCGa                                     >  1:5504989/1‑36 (MQ=255)
                                                                    tttttCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAGCGCTGTCAtt  <  1:2018490/71‑1 (MQ=255)
                                                                    tttttCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAGCGCTGTCAtt  >  1:2979538/1‑71 (MQ=255)
                                                                    tttttCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAGCGCTGTCAtt  >  1:3624727/1‑71 (MQ=255)
                                                                    tttttCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAGCGCTGTCAt   >  1:3838273/1‑70 (MQ=255)
                                                                    tttttCTTACCGATGAAGTACATGCACGTCCGGCGa                                     >  1:5568328/1‑36 (MQ=255)
                                                                    tttttCTTACCGATGAAGTACATGCACGTCCGGCGa                                     >  1:5574928/1‑36 (MQ=255)
                                                                     ttttCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAt                            <  1:1920099/44‑1 (MQ=255)
                                                                     |                                                                     
GACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAGCGCTGTCATT  >  minE/1395345‑1395483

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: