Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,023 A→G 7.7% 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 7.7% 215.6 / 7.7 78R418R (CGT→CGCmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base A (31/41);  new base G (4/2);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

GTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCC  >  minE/2585959‑2586088
                                                                |                                                                 
gTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTc                                                              >  1:2087638/1‑70 (MQ=255)
gTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTc                                                              >  1:3465481/1‑70 (MQ=255)
gTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCg                                                             <  1:2147322/71‑1 (MQ=255)
  tGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGAc                                                           >  1:3050303/1‑71 (MQ=255)
  tGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGATCGa                                                            >  1:737124/1‑70 (MQ=255)
          ccccATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTtctc                                                     <  1:1335659/69‑1 (MQ=255)
          ccccATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTcca                                                   <  1:1763130/71‑1 (MQ=255)
            ccATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccacc                                                 <  1:3363532/71‑1 (MQ=255)
             cATCGATTTTTCCGTCGCTTCTTTTGGCGGGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTccacca                                                >  1:1872414/1‑71 (MQ=255)
              aTCGATTTTTCCGTCGCTTCTTTTGGCGGGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTcc                                                    <  1:3054131/66‑1 (MQ=255)
              aTCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccac                                               >  1:20320/1‑71 (MQ=255)
              aTCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccac                                               >  1:3949271/1‑71 (MQ=255)
              aTCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccac                                               >  1:1533686/1‑71 (MQ=255)
               tCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccac                                               >  1:3540943/1‑70 (MQ=255)
               tCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccac                                               >  1:1977280/1‑70 (MQ=255)
                 gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                            >  1:2219701/1‑71 (MQ=255)
                  aTTTTTCCGTCGCTTCTTTTGGCGGGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTccaccacca                                             >  1:1671793/1‑69 (MQ=255)
                  aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                            >  1:3886668/1‑70 (MQ=255)
                  aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACg                                           <  1:2485590/71‑1 (MQ=255)
                  aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACg                                           <  1:3494392/71‑1 (MQ=255)
                  aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACg                                           <  1:559040/71‑1 (MQ=255)
                   tttATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCg                                                             <  1:437012/52‑1 (MQ=255)
                   tttATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGa                                          <  1:2510934/71‑1 (MQ=255)
                   tttATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGa                                          <  1:2038118/71‑1 (MQ=255)
                    ttATGCGTCGCTTCCTTCGGTGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGa                                          >  1:225237/1‑70 (MQ=255)
                    ttATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGAt                                         >  1:1175403/1‑71 (MQ=255)
                     tATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGAt                                         <  1:407216/70‑1 (MQ=255)
                      aTGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATg                                        <  1:1426324/70‑1 (MQ=255)
                      aTGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATg                                        <  1:3945938/70‑1 (MQ=255)
                       tGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATgg                                       <  1:274622/70‑1 (MQ=255)
                       tGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGc                                      <  1:2282291/71‑1 (MQ=255)
                           tCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACg                                                          <  1:689279/47‑1 (MQ=255)
                           tCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccac                                                  <  1:184344/55‑1 (MQ=255)
                              cttccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcc                               <  1:2360006/71‑1 (MQ=255)
                              cttccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcc                               <  1:576957/71‑1 (MQ=255)
                                tccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                             >  1:625026/1‑71 (MQ=255)
                                tccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                             >  1:1285240/1‑71 (MQ=255)
                                 ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTTATccac                             <  1:511122/70‑1 (MQ=255)
                                 ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                             <  1:2412808/70‑1 (MQ=255)
                                    ttggcggGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAg                          <  1:2026811/68‑1 (MQ=12)
                                    tcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                                            <  1:1894718/52‑1 (MQ=255)
                                     cggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCGCCAGt                         >  1:1686169/1‑70 (MQ=255)
                                      ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTaa                       >  1:3302608/1‑71 (MQ=255)
                                      ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTa                        <  1:1783259/70‑1 (MQ=255)
                                      ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTa                        <  1:1024032/70‑1 (MQ=255)
                                      ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTa                        <  1:3740850/70‑1 (MQ=255)
                                       gcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAc                      >  1:3736467/1‑71 (MQ=255)
                                         ggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca                           <  1:1971060/64‑1 (MQ=255)
                                          ggaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGccc                    >  1:1364286/3‑70 (MQ=255)
                                           tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                  <  1:290303/71‑1 (MQ=255)
                                           tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                  <  1:3595756/71‑1 (MQ=255)
                                            aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTAt                 >  1:3489502/1‑71 (MQ=255)
                                             gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                <  1:1975963/71‑1 (MQ=255)
                                             gCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGAtg                >  1:3295347/1‑70 (MQ=255)
                                               ttATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc              >  1:2878324/1‑71 (MQ=255)
                                                tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc              >  1:3234193/1‑70 (MQ=255)
                                                tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc              >  1:1274685/1‑70 (MQ=255)
                                                 aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAccc                    <  1:1540280/63‑1 (MQ=255)
                                                 aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCg               <  1:2119080/68‑1 (MQ=255)
                                                 aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGcc             <  1:3554745/70‑1 (MQ=255)
                                                  tCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACg                                           <  1:2427728/39‑1 (MQ=255)
                                                  tCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCa            <  1:2884450/70‑1 (MQ=255)
                                                    ttCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATACACCAGTAACCCTATCGCCAGCa         >  1:3488866/1‑71 (MQ=255)
                                                     tCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc        >  1:521560/1‑71 (MQ=255)
                                                      ccGCAATGACACTCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa         <  1:1871226/69‑1 (MQ=255)
                                                      ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa         <  1:3080823/69‑1 (MQ=255)
                                                       cGCAATTACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAcc       <  1:2019901/70‑1 (MQ=255)
                                                       cGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAcc       <  1:59734/70‑1 (MQ=255)
                                                       cGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCa      <  1:3776617/71‑1 (MQ=255)
                                                        gCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTAt                 >  1:179953/1‑59 (MQ=255)
                                                          aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTaa                       >  1:2727636/1‑51 (MQ=255)
                                                          aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                  <  1:1279990/56‑1 (MQ=255)
                                                          aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc   >  1:3364156/1‑71 (MQ=255)
                                                          aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc   >  1:253318/1‑71 (MQ=255)
                                                            tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                >  1:3748009/1‑56 (MQ=255)
                                                            tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                >  1:1467828/1‑56 (MQ=255)
                                                            tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATccc  >  1:416224/1‑70 (MQ=255)
                                                             gACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc        <  1:786780/63‑1 (MQ=255)
                                                                |                                                                 
GTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCC  >  minE/2585959‑2586088

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: