Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,396,005 C→A 5.6% I237I (ATC→ATA setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,396,0050CA 5.6% 287.4 / 2.3 88I237I (ATC→ATAsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base C (38/45);  new base A (2/3);  total (40/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

ACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGAT  >  minE/1395935‑1396067
                                                                      |                                                              
aCCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACatcatc                                                                <  1:726205/71‑1 (MQ=255)
  cGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCaa                                                              <  1:1385944/71‑1 (MQ=255)
  cGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCaa                                                              <  1:2807599/71‑1 (MQ=255)
  cGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCaa                                                              >  1:2711968/1‑71 (MQ=255)
           aCGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGc                                                      >  1:4964894/1‑70 (MQ=255)
            cGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCt                                                     <  1:5105210/70‑1 (MQ=255)
             gctgctgctGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCtat                                                   >  1:1000954/1‑71 (MQ=255)
               tgctgctgTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCtatt                                                  >  1:569779/1‑70 (MQ=255)
                 ctgctgTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACAt                                                           <  1:3880622/59‑1 (MQ=255)
                 ctgctgTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCtattt                                                 <  1:319693/69‑1 (MQ=255)
                    ctgTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATg                                                          >  1:1060381/1‑57 (MQ=255)
                       tttGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACAtcatca                                                               <  1:3329894/49‑1 (MQ=255)
                       tttGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCaa                                                              <  1:5269931/50‑1 (MQ=255)
                       tttGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATc                                           <  1:446060/69‑1 (MQ=255)
                       tttGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATc                                           <  1:552655/69‑1 (MQ=255)
                       tttGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATc                                           <  1:653063/69‑1 (MQ=255)
                       tttGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATc                                           <  1:704981/69‑1 (MQ=255)
                         tGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCaa                                         >  1:5250506/1‑69 (MQ=255)
                         tGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAAc                                        >  1:3959260/1‑70 (MQ=255)
                         tGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAAc                                        >  1:2972197/1‑70 (MQ=255)
                          gTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACg                                       <  1:1843657/70‑1 (MQ=255)
                           tCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCaa                                                              <  1:4411523/46‑1 (MQ=255)
                            cATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACAtcatca                                                               <  1:1508079/44‑1 (MQ=255)
                             aTTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACt                                   >  1:654201/1‑71 (MQ=255)
                             aTTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACt                                   >  1:1882355/1‑71 (MQ=255)
                              tttGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACg                                       <  1:5656636/66‑1 (MQ=255)
                              tttGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACg                                       <  1:557893/66‑1 (MQ=255)
                              tttGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTg                                  <  1:851785/71‑1 (MQ=255)
                              tttGTACATTGATGTGGGGCTCGAACAGCC‑‑TACATCATCAACATGCCGCTATTTATTATCAACg                                       <  1:5669471/64‑1 (MQ=255)
                               ttGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATttatta                                             <  1:4495391/59‑1 (MQ=255)
                               ttGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATttatta                                             <  1:4622908/59‑1 (MQ=255)
                               ttGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATttatta                                             <  1:1154904/59‑1 (MQ=255)
                               ttGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATttatta                                             <  1:3509749/59‑1 (MQ=255)
                               ttGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTg                                  >  1:5586422/1‑70 (MQ=255)
                               ctGTTCAATGATGTGGGCGGCGAATAATCTCTACATGATAAATATGCCGCTATTTATTATTGATGAACTg                                  >  1:2536813/2‑70 (MQ=255)
                                tGTTCAATGATGTGGGCGGCGAATAATCTCTACATGATAAATATGCCGCTATTtattatt                                           <  1:955745/60‑2 (MQ=255)
                                tGTTCAATGATGTGGGCGGCGAATAATCTCTACATGATAAATATGCCGCTATTTATTATTGATGAACTGCa                                >  1:4973423/1‑71 (MQ=255)
                                tGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCtatt                                                  >  1:5599928/1‑53 (MQ=255)
                                tGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCa                                >  1:197016/1‑71 (MQ=255)
                                 gTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCaa                                         >  1:3018465/1‑61 (MQ=255)
                                   aCATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCtattta                                                >  1:2934567/1‑52 (MQ=255)
                                   aCATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCtattta                                                >  1:1121831/1‑52 (MQ=255)
                                    cATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCtt                            <  1:3073421/71‑1 (MQ=255)
                                     aTTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAAc                                    >  1:4877722/1‑62 (MQ=255)
                                     aTTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTc                           <  1:1655155/71‑1 (MQ=255)
                                     aTTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTc                           <  1:60312/71‑1 (MQ=255)
                                       tGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCtatttat                                               >  1:4176034/1‑49 (MQ=255)
                                       tGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCtatttat                                               >  1:4668722/1‑49 (MQ=255)
                                       tGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCtatttat                                               >  1:5612835/1‑49 (MQ=255)
                                       tGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCtt                            >  1:2042681/1‑68 (MQ=255)
                                       tGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTccc                         >  1:3514209/1‑71 (MQ=255)
                                       tGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTccc                         >  1:3670382/1‑71 (MQ=255)
                                        gATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGaa                                     >  1:1339412/1‑58 (MQ=255)
                                         aTGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCga                       <  1:2398905/71‑1 (MQ=255)
                                          tgtgGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCgag                      <  1:3667849/71‑1 (MQ=255)
                                            tgGGGCTCGAACAGCCTGTACATCATCAACATGCCGc                                                      <  1:3925178/37‑1 (MQ=255)
                                            tgGGGCTCGAACAGCCTGTACATCATCAACATGCCGc                                                      <  1:2412291/37‑1 (MQ=255)
                                            tgGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCgaga                     <  1:3347261/70‑1 (MQ=255)
                                            tgGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCgaga                     <  1:2067995/70‑1 (MQ=255)
                                            tgGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGaa                    <  1:3575806/71‑1 (MQ=255)
                                            tgGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGaa                    <  1:4257263/71‑1 (MQ=255)
                                                gCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTc                           <  1:1034926/60‑1 (MQ=255)
                                                gCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTc                           <  1:2114125/60‑1 (MQ=255)
                                                  tCGAACAGCCTGTACATCATCAACATGCCGCtatt                                                  <  1:1338762/35‑1 (MQ=255)
                                                  tCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGc              >  1:4624053/1‑71 (MQ=255)
                                                    gAACAGCCTGTACATCATCAACATGCCGCTATTTATTATc                                           <  1:810116/40‑1 (MQ=255)
                                                    gAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGa                                      >  1:3867504/1‑45 (MQ=255)
                                                    gAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGa                                      >  1:5034800/1‑45 (MQ=255)
                                                    gAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGa                                      >  1:3130385/1‑45 (MQ=255)
                                                    gAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGa                                      >  1:1340580/1‑45 (MQ=255)
                                                    gAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCg            >  1:4434653/1‑71 (MQ=255)
                                                    gAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCg            <  1:4939742/71‑1 (MQ=255)
                                                     aaCAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGa                                      >  1:4155517/1‑44 (MQ=255)
                                                      aCAGCCTGTACATCATCAACATGCCGCTATttatta                                             >  1:5663360/1‑36 (MQ=255)
                                                      aCAGCCTGTACATCATCAACATGCCGCTATttatta                                             >  1:400335/1‑36 (MQ=255)
                                                      aCAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGgt          >  1:3574674/1‑71 (MQ=255)
                                                      aCAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGgt          >  1:2819425/1‑71 (MQ=255)
                                                        aGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCgg           <  1:3203946/68‑1 (MQ=255)
                                                         gCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCgg           <  1:1967018/67‑1 (MQ=255)
                                                         gCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGgtgt        <  1:5342752/70‑1 (MQ=255)
                                                         gCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGtgtg       <  1:4254711/71‑1 (MQ=255)
                                                           cTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGtgat     <  1:5405958/71‑1 (MQ=255)
                                                           cTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGtgat     <  1:2847149/71‑1 (MQ=255)
                                                              tACATGATAAATATGCCGCTATTTATTATTGATGAACTGc                                 <  1:4544882/40‑1 (MQ=12)
                                                              tACATCATCAACATGCCGCTATTTATTATCAACGAACTGc                                 <  1:3200847/40‑1 (MQ=255)
                                                              tACATCATCAACATGCCGCTATTTATTATCAACGAACTGc                                 <  1:2061786/40‑1 (MQ=255)
                                                              tACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGtgatga   <  1:172156/70‑1 (MQ=255)
                                                              tACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTgatgat  >  1:2900126/1‑71 (MQ=255)
                                                              tACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTgatgat  >  1:5058574/1‑71 (MQ=255)
                                                              tACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTgatgat  >  1:2076898/1‑71 (MQ=255)
                                                                catgatAAATATGCCGCTATTTATTATTGATGAACTGCATCTAACCGATAAACTg                <  1:507621/55‑1 (MQ=255)
                                                                  tcatcaACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGcc             <  1:1072108/56‑1 (MQ=255)
                                                                  tcatcaACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGcc             <  1:4310370/56‑1 (MQ=255)
                                                                   catcaACATGCCGCTATTTATTATCAACTAACTGCATCTTCCCGAGAAACTGGcc             <  1:3175233/55‑1 (MQ=255)
                                                                     tcaACATGCCGCTATTTATTATCAACGAACTGCATCTTccc                         <  1:1708762/41‑1 (MQ=255)
                                                                     tcaACATGCCGCTATTTATTATCAACGAACTGCATCTTccc                         <  1:3708/41‑1 (MQ=255)
                                                                     tcaACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCgaga                     <  1:5406828/45‑1 (MQ=255)
                                                                     tcaACATGCCGCTATTTATTATCAAAGAACTGCATCTTCCCGAGAAACTGGCCGGTGTgatgat  <  1:4699823/64‑1 (MQ=255)
                                                                      caACATGCCGCTATTTATTATCAACGAACTGCATCTTccc                         <  1:4090578/40‑1 (MQ=255)
                                                                      |                                                              
ACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGAT  >  minE/1395935‑1396067

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: