Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,023 A→G 9.6% R418R (CGT→CGC mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,0230AG 9.6% 280.5 / 19.0 104R418R (CGT→CGCmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base A (50/44);  new base G (2/8);  total (52/52)
Fisher's exact test for biased strand distribution p-value = 9.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 115 total aligned reads displayed.

CGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACA  >  minE/2585954‑2586093
                                                                     |                                                                      
cgcACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGacac                                                                       <  1:2927509/71‑1 (MQ=255)
cgcACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGacac                                                                       <  1:3478591/71‑1 (MQ=255)
 gcACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCAGCAATGACACg                                                                      >  1:446923/1‑71 (MQ=255)
     gTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTc                                                                   <  1:3039609/70‑1 (MQ=255)
       tGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGAc                                                                >  1:1116533/1‑71 (MQ=255)
        ggATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGAc                                                                <  1:148213/70‑1 (MQ=255)
        ggATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACg                                                               >  1:284226/1‑71 (MQ=255)
        ggATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTAACTTCCGCAATGACACGCTCGACg                                                               >  1:86847/1‑71 (MQ=255)
             tGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTtctc                                                          >  1:2667679/1‑71 (MQ=255)
             tGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTtctc                                                          >  1:2822200/1‑71 (MQ=255)
             tGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTtctc                                                          >  1:1540611/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)
                   aTCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccac                                                    >  1:104508/1‑71 (MQ=255)
                   aTCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccac                                                    <  1:2292232/71‑1 (MQ=255)
                    tCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGAGACGCTCGACGTTCTccaccac                                                    >  1:646283/1‑70 (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)
                      gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccacca                                                  >  1:2664385/1‑70 (MQ=255)
                      gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTcc                                                         >  1:3998611/1‑63 (MQ=255)
                      gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                                 >  1:3514972/1‑71 (MQ=255)
                      gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                                 >  1:2474361/1‑71 (MQ=255)
                      gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                                 <  1:780664/71‑1 (MQ=255)
                      gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                                 >  1:2179504/1‑71 (MQ=255)
                      gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                                 >  1:2975242/1‑71 (MQ=255)
                      gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                                 >  1:3083339/1‑71 (MQ=255)
                      gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                                 <  1:3123705/71‑1 (MQ=255)
                       aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGa                                                                 <  1:2230086/54‑1 (MQ=255)
                       aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGa                                                                 <  1:3385717/54‑1 (MQ=255)
                       aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                                 <  1:2486509/70‑1 (MQ=255)
                       aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                                 <  1:2686758/70‑1 (MQ=255)
                       aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                                 <  1:347909/70‑1 (MQ=255)
                        tttATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGtt                                                             >  1:1567068/1‑57 (MQ=255)
                           aTGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATg                                             <  1:3102437/70‑1 (MQ=255)
                            tGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTcca                                                        >  1:932963/1‑58 (MQ=255)
                              cGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCg                                          <  1:2174545/70‑1 (MQ=255)
                                tCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATTACACGCTCGACGTTCTCCACCACCACGATGGCGTCa                                       >  1:3425377/1‑71 (MQ=255)
                                tCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCa                                       >  1:548634/1‑71 (MQ=255)
                                   cttcttttGGCGGGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCatcatc                                     <  1:3797960/70‑1 (MQ=12)
                                      ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                  <  1:3623113/70‑1 (MQ=255)
                                      ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                  <  1:1032638/70‑1 (MQ=255)
                                      ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                  >  1:2138510/1‑70 (MQ=255)
                                          cggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGt                              <  1:2144167/70‑1 (MQ=255)
                                           ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTaa                            <  1:2388022/71‑1 (MQ=255)
                                            gcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAc                           <  1:3349906/71‑1 (MQ=255)
                                              gggaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGccc                         <  1:3071193/68‑1 (MQ=9)
                                               gTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGa                                               <  1:661694/48‑1 (MQ=255)
                                               gTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAccc                         >  1:790588/1‑70 (MQ=255)
                                                tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCt                        <  1:52052/70‑1 (MQ=255)
                                                tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                       >  1:767643/1‑71 (MQ=255)
                                                tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                       <  1:3260103/71‑1 (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:579385/1‑57 (MQ=255)
                                                  gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcca                                   >  1:1924630/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:176843/45‑1 (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:899383/1‑69 (MQ=255)
                                                      aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                   >  1:1404749/1‑69 (MQ=255)
                                                       tCTTCCGCAATGACACGCTCGACGTTCTccaccac                                                    <  1:3735670/35‑1 (MQ=255)
                                                        cTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGt                                         >  1:4047199/1‑45 (MQ=255)
                                                        cTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAg                >  1:2156977/1‑70 (MQ=255)
                                                         ttCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATc                                     <  1:562592/48‑1 (MQ=255)
                                                         ttCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa              >  1:3505798/1‑71 (MQ=255)
                                                           ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAg                <  1:1017283/67‑1 (MQ=255)
                                                           ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAg                <  1:2466172/67‑1 (MQ=255)
                                                           ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa              <  1:2726256/69‑1 (MQ=255)
                                                              cAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCAt          <  1:3209622/70‑1 (MQ=255)
                                                              cAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATc         <  1:2465905/71‑1 (MQ=255)
                                                               aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATc         <  1:644711/70‑1 (MQ=255)
                                                                aTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCt                        <  1:1102717/54‑1 (MQ=255)
                                                                aTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCa           <  1:1663049/67‑1 (MQ=255)
                                                                aTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCa           <  1:3473991/67‑1 (MQ=255)
                                                                 tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATccc       >  1:1348189/1‑70 (MQ=255)
                                                                 tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCg      <  1:814906/71‑1 (MQ=255)
                                                                   acacGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca                                >  1:2704051/1‑43 (MQ=255)
                                                                   acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGaa    >  1:501919/1‑71 (MQ=255)
                                                                     acGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca                                >  1:4037653/1‑41 (MQ=255)
                                                                     acGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca                                >  1:3612529/1‑41 (MQ=255)
                                                                     acGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa  <  1:3927935/71‑1 (MQ=255)
                                                                     acGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa  <  1:2006510/71‑1 (MQ=255)
                                                                     |                                                                      
CGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACA  >  minE/2585954‑2586093

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: