Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,011 T→C 9.3% E422E (GAA→GAG mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,0110TC 9.3% 299.0 / 19.0 107E422E (GAA→GAGmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base T (57/40);  new base C (2/8);  total (59/48)
Fisher's exact test for biased strand distribution p-value = 4.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 113 total aligned reads displayed.

AATACCGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCA  >  minE/2585942‑2586081
                                                                     |                                                                      
aaTACCGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCtt                                                                       <  1:2468813/71‑1 (MQ=255)
   aCCGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCg                                                                    >  1:2518016/1‑71 (MQ=255)
    ccGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCg                                                                    >  1:4009889/1‑70 (MQ=255)
    ccGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCg                                                                    >  1:2012078/1‑70 (MQ=255)
    ccGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGc                                                                   <  1:3648048/71‑1 (MQ=255)
       aCCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAAt                                                                <  1:3550423/71‑1 (MQ=255)
       aCCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAAt                                                                <  1:412882/71‑1 (MQ=255)
            cgcACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGacac                                                           <  1:2927509/71‑1 (MQ=255)
            cgcACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGacac                                                           <  1:3478591/71‑1 (MQ=255)
             gcACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGa                                                              >  1:3380668/1‑67 (MQ=255)
             gcACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCAGCAATGACACg                                                          >  1:446923/1‑71 (MQ=255)
              cACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGa                                                              >  1:818546/1‑66 (MQ=255)
                 gTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTc                                                       <  1:3039609/70‑1 (MQ=255)
                   tGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGAc                                                    >  1:1116533/1‑71 (MQ=255)
                    ggATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACg                                                   >  1:3187793/1‑71 (MQ=255)
                    ggATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACg                                                   >  1:284226/1‑71 (MQ=255)
                    ggATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTAACTTCCGCAATGACACGCTCGACg                                                   >  1:86847/1‑71 (MQ=255)
                       tCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCaa                                                                 >  1:3157120/1‑54 (MQ=255)
                         tGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTtctc                                              >  1:1540611/1‑71 (MQ=255)
                         tGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTtctc                                              >  1:2667679/1‑71 (MQ=255)
                         tGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTtctc                                              >  1:2822200/1‑71 (MQ=255)
                         tGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTtctc                                              >  1:3710093/1‑71 (MQ=255)
                          gCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTtctc                                              >  1:2567313/1‑70 (MQ=255)
                           cgacATCGATTTTTCCGTCGCTTCTTTTGGCGGGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTcca                                            <  1:3377/68‑1 (MQ=255)
                           ccccATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTtctc                                              <  1:2347520/69‑1 (MQ=255)
                              cATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCtt                                                                       <  1:828214/41‑1 (MQ=255)
                               aTCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccac                                        >  1:104508/1‑71 (MQ=255)
                               aTCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccac                                        <  1:2292232/71‑1 (MQ=255)
                                tCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCtt                                                                       <  1:406211/39‑1 (MQ=255)
                                tCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGa                                                              <  1:2468640/48‑1 (MQ=255)
                                tCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGAGACGCTCGACGTTCTccaccac                                        >  1:646283/1‑70 (MQ=255)
                                tCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACg                                                   <  1:2909084/59‑1 (MQ=255)
                                tCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccacc                                       <  1:3466649/71‑1 (MQ=255)
                                tCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccac                                        >  1:781318/1‑70 (MQ=255)
                                tCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccac                                           <  1:1987355/67‑1 (MQ=255)
                                 cGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGAAATGACACGCTCGACGTTCTccaccacca                                      >  1:3596741/1‑71 (MQ=255)
                                  gATTTTTCCGTCGCTTCTTTTGGCGGGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTccaccacca                                      >  1:1034524/1‑70 (MQ=255)
                                  gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGacac                                                           >  1:1658697/1‑49 (MQ=255)
                                  gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGa                                                              >  1:3665277/1‑46 (MQ=255)
                                  gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGa                                                              >  1:1817004/1‑46 (MQ=255)
                                  gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccacca                                      >  1:2664385/1‑70 (MQ=255)
                                  gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTcc                                             >  1:3998611/1‑63 (MQ=255)
                                  gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                     >  1:3514972/1‑71 (MQ=255)
                                  gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                     <  1:780664/71‑1 (MQ=255)
                                  gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                     >  1:2474361/1‑71 (MQ=255)
                                  gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                     >  1:2179504/1‑71 (MQ=255)
                                  gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                     >  1:3083339/1‑71 (MQ=255)
                                  gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                     <  1:3123705/71‑1 (MQ=255)
                                  gATTTATGCGTAGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                     <  1:4086823/71‑1 (MQ=255)
                                   aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGa                                                     <  1:3385717/54‑1 (MQ=255)
                                   aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGa                                                     <  1:2230086/54‑1 (MQ=255)
                                   aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                     <  1:2486509/70‑1 (MQ=255)
                                   aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                     <  1:2686758/70‑1 (MQ=255)
                                   aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                     <  1:2636548/70‑1 (MQ=255)
                                    tttATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGtt                                                 >  1:1567068/1‑57 (MQ=255)
                                    tttATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGa                                   >  1:2428704/1‑71 (MQ=255)
                                       aTGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATg                                 <  1:3102437/70‑1 (MQ=255)
                                        tGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTcca                                            >  1:932963/1‑58 (MQ=255)
                                          cGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCg                              <  1:1850986/70‑1 (MQ=255)
                                          cGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCg                              <  1:2174545/70‑1 (MQ=255)
                                            tCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATTACACGCTCGACGTTCTCCACCACCACGATGGCGTCa                           >  1:3425377/1‑71 (MQ=255)
                                            tCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCa                           >  1:548634/1‑71 (MQ=255)
                                              gCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATc                         >  1:2163682/1‑71 (MQ=255)
                                               cttcttttGGCGGGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCatcatc                         <  1:3797960/70‑1 (MQ=12)
                                                  ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                      <  1:1032638/70‑1 (MQ=255)
                                                  ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                      <  1:3623113/70‑1 (MQ=255)
                                                  ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                      >  1:2138510/1‑70 (MQ=255)
                                                      cggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGt                  <  1:2144167/70‑1 (MQ=255)
                                                        gcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAc               <  1:3349906/71‑1 (MQ=255)
                                                          gggaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGccc             <  1:3071193/68‑1 (MQ=9)
                                                          gggaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGccc             <  1:626892/68‑1 (MQ=9)
                                                           gTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGa                                   <  1:661694/48‑1 (MQ=255)
                                                           gTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAccc             >  1:790588/1‑70 (MQ=255)
                                                            tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCCCCAGTAACCCTa           >  1:1636305/1‑71 (MQ=255)
                                                            tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCt            <  1:52052/70‑1 (MQ=255)
                                                            tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa           >  1:767643/1‑71 (MQ=255)
                                                            gaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGcccga           <  1:4053857/70‑3 (MQ=255)
                                                            gaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGcccga           >  1:222344/2‑69 (MQ=255)
                                                             aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAg                   >  1:1719769/1‑62 (MQ=255)
                                                             aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCt            >  1:4113399/1‑69 (MQ=255)
                                                             aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTAt          >  1:1670558/1‑71 (MQ=255)
                                                              gCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccac                                        >  1:1035014/1‑40 (MQ=255)
                                                              gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCGCGATGGCGTCATcca                       >  1:1881326/1‑57 (MQ=255)
                                                              gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcca                       >  1:1924630/1‑57 (MQ=255)
                                                              gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcca                       >  1:579385/1‑57 (MQ=255)
                                                              gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATc                         >  1:1506712/1‑55 (MQ=255)
                                                               cTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAc               <  1:1695424/64‑1 (MQ=255)
                                                               cTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCatcatcg                        <  1:3377221/55‑2 (MQ=255)
                                                               cTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCatcatcg                        <  1:3896454/55‑2 (MQ=255)
                                                                 tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc       >  1:920585/1‑70 (MQ=255)
                                                                  aTCTTCCGCAATGACACGCTCGACGTTCTCCACCCCCACGATGGCGTCATccacca                    >  1:837820/1‑56 (MQ=255)
                                                                  aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGc                               <  1:3594923/45‑1 (MQ=255)
                                                                  aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGc                               <  1:2556031/45‑1 (MQ=255)
                                                                  aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGt                  >  1:1107565/1‑58 (MQ=255)
                                                                  aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc       >  1:1404749/1‑69 (MQ=255)
                                                                  aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc       >  1:899383/1‑69 (MQ=255)
                                                                   tCTTCCGCAATGACACGCTCGACGTTCTccaccac                                        <  1:3735670/35‑1 (MQ=255)
                                                                   tCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCt            <  1:3825497/63‑1 (MQ=255)
                                                                     ttCCGCAATGACACGCTCGACGTTCTccaccacct                                      >  1:1109879/1‑34 (MQ=255)
                                                                     ttCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa  >  1:3505798/1‑71 (MQ=255)
                                                                     |                                                                      
AATACCGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCA  >  minE/2585942‑2586081

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: