Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,396,008 C→T 5.3% N238N (AAC→AAT setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,396,0080CT 5.3% 287.9 / 3.4 94N238N (AAC→AATsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base C (39/50);  new base T (2/3);  total (41/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

ACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCC  >  minE/1395946‑1396077
                                                              |                                                                     
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)
            tttGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATc                                                     <  1:704981/69‑1 (MQ=255)
            tttGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATc                                                     <  1:653063/69‑1 (MQ=255)
            tttGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATc                                                     <  1:446060/69‑1 (MQ=255)
            tttGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATc                                                     <  1:552655/69‑1 (MQ=255)
              tGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCaa                                                   >  1:5250506/1‑69 (MQ=255)
              tGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAAc                                                  >  1:2972197/1‑70 (MQ=255)
              tGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAAc                                                  >  1:3959260/1‑70 (MQ=255)
               gTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACg                                                 <  1:1843657/70‑1 (MQ=255)
                  aTTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACt                                             >  1:1882355/1‑71 (MQ=255)
                  aTTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACt                                             >  1:654201/1‑71 (MQ=255)
                   tttGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACg                                                 <  1:557893/66‑1 (MQ=255)
                   tttGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACg                                                 <  1:5656636/66‑1 (MQ=255)
                   tttGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTg                                            <  1:851785/71‑1 (MQ=255)
                   tttGTACATTGATGTGGGGCTCGAACAGCC‑‑TACATCATCAACATGCCGCTATTTATTATCAACg                                                 <  1:5669471/64‑1 (MQ=255)
                    ttGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATttatta                                                       <  1:4622908/59‑1 (MQ=255)
                    ttGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATttatta                                                       <  1:4495391/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:60312/71‑1 (MQ=255)
                          aTTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTc                                     <  1:1655155/71‑1 (MQ=255)
                            tGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCtatttat                                                         >  1:4176034/1‑49 (MQ=255)
                            tGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCtatttat                                                         >  1:5612835/1‑49 (MQ=255)
                            tGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCtatttat                                                         >  1:4668722/1‑49 (MQ=255)
                            tGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCtt                                      >  1:2042681/1‑68 (MQ=255)
                            tGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTccc                                   >  1:3670382/1‑71 (MQ=255)
                            tGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTccc                                   >  1:3514209/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:2412291/37‑1 (MQ=255)
                                 tgGGGCTCGAACAGCCTGTACATCATCAACATGCCGc                                                                <  1:3925178/37‑1 (MQ=255)
                                 tgGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCgaga                               <  1:2067995/70‑1 (MQ=255)
                                 tgGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCgaga                               <  1:3347261/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:3130385/1‑45 (MQ=255)
                                         gAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGa                                                >  1:1340580/1‑45 (MQ=255)
                                         gAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGa                                                >  1:3867504/1‑45 (MQ=255)
                                         gAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGa                                                >  1:5034800/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:400335/1‑36 (MQ=255)
                                           aCAGCCTGTACATCATCAACATGCCGCTATttatta                                                       >  1:5663360/1‑36 (MQ=255)
                                           aCAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGgt                    >  1:2819425/1‑71 (MQ=255)
                                           aCAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGgt                    >  1:3574674/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:2847149/71‑1 (MQ=255)
                                                cTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGtgat               <  1:5405958/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)
                                                            aaCATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGAccgc   >  1:2834346/1‑71 (MQ=255)
                                                             aCATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGAccgcc  <  1:4852869/71‑1 (MQ=255)
                                                             aCATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGAccgcc  >  1:2268127/1‑71 (MQ=255)
                                                              |                                                                     
ACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCC  >  minE/1395946‑1396077

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: