Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,059 C→G 7.8% V406V (GTG→GTC mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,0590CG 7.8% 280.1 / 10.1 90V406V (GTG→GTCmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base C (36/47);  new base G (3/4);  total (39/51)
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.29e-01

Only 100 of 101 total aligned reads displayed.

CCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCC  >  minE/2585992‑2586125
                                                                   |                                                                  
ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                                  <  1:1032638/70‑1 (MQ=255)
ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                                  <  1:3623113/70‑1 (MQ=255)
ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                                  >  1:2138510/1‑70 (MQ=255)
ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTATCCACCACCACGATGGCGTCATccac                                                                  >  1:984765/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)
        gggaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGccc                                                         <  1:626892/68‑1 (MQ=9)
         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)
          tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                                       <  1:3260103/71‑1 (MQ=255)
          gaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGcccga                                                       >  1:222344/2‑69 (MQ=255)
          gaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGcccga                                                       <  1:4053857/70‑3 (MQ=255)
           aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAg                                                               >  1:1719769/1‑62 (MQ=255)
           aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCt                                                        >  1:4113399/1‑69 (MQ=255)
           aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTAt                                                      >  1:1670558/1‑71 (MQ=255)
            gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCGCGATGGCGTCATcca                                                                   >  1:1881326/1‑57 (MQ=255)
            gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcca                                                                   >  1:1924630/1‑57 (MQ=255)
            gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcca                                                                   >  1:579385/1‑57 (MQ=255)
             cTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAc                                                           <  1:1695424/64‑1 (MQ=255)
               tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                                   >  1:920585/1‑70 (MQ=255)
                aTCTTCCGCAATGACACGCTCGACGTTCTCCACCCCCACGATGGCGTCATccacca                                                                >  1:837820/1‑56 (MQ=255)
                aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGt                                                              >  1:1107565/1‑58 (MQ=255)
                aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                                   >  1:1404749/1‑69 (MQ=255)
                aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                                   >  1:899383/1‑69 (MQ=255)
                 tCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCt                                                        <  1:3825497/63‑1 (MQ=255)
                  cTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAg                                                >  1:2156977/1‑70 (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:3473991/67‑1 (MQ=255)
                          aTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCa                                           <  1:1663049/67‑1 (MQ=255)
                           tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATccc                                       >  1:1348189/1‑70 (MQ=255)
                           tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCg                                      <  1:814906/71‑1 (MQ=255)
                             acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                                    >  1:2630334/4‑71 (MQ=2)
                             acacGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca                                                                >  1:2704051/1‑43 (MQ=255)
                             acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                             >  1:3560893/1‑62 (MQ=255)
                             acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGaa                                    >  1:3412148/1‑71 (MQ=255)
                             acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGaa                                    >  1:501919/1‑71 (MQ=255)
                               acGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca                                                                >  1:3612529/1‑41 (MQ=255)
                               acGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca                                                                >  1:4037653/1‑41 (MQ=255)
                               acGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                  <  1:3927935/71‑1 (MQ=255)
                               acGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                  <  1:2006510/71‑1 (MQ=255)
                                   tCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAAc                                   <  1:379082/66‑1 (MQ=9)
                                    cGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGt                              >  1:3853016/1‑70 (MQ=2)
                                     gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                             <  1:3155042/70‑1 (MQ=255)
                                     gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                             <  1:2081560/70‑1 (MQ=255)
                                     gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                             <  1:1334357/70‑1 (MQ=255)
                                     gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCa                            <  1:115142/71‑1 (MQ=255)
                                     gACGTTCTCCAACACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                             <  1:1914800/70‑1 (MQ=255)
                                         ttCTCCACCACCACGATGGCGTCATCCACCAGTAAccc                                                         <  1:1429061/38‑1 (MQ=255)
                                         ttCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATa                                <  1:920801/63‑1 (MQ=255)
                                          tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                        >  1:4138693/1‑54 (MQ=255)
                                          tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                        >  1:932832/1‑54 (MQ=255)
                                           ctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAGCATACCGAACATAGTCAACgtgt                      >  1:3743712/1‑71 (MQ=255)
                                           ctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                          >  1:2643675/1‑67 (MQ=255)
                                           ctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                          >  1:3828203/1‑67 (MQ=255)
                                           ctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgtg                       <  1:2930515/70‑1 (MQ=255)
                                              caccaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                                   <  1:4077617/55‑1 (MQ=255)
                                              caccaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGtt                     <  1:2970651/69‑1 (MQ=255)
                                              caccaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTg                    <  1:1637659/70‑1 (MQ=255)
                                              caccaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGa                   <  1:667836/71‑1 (MQ=255)
                                               accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                  <  1:1880151/71‑1 (MQ=255)
                                                ccaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                  <  1:3462353/54‑1 (MQ=255)
                                                ccaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCa                            <  1:591886/60‑1 (MQ=255)
                                                 caccacGATTGCGTCATCCACAAGTAACCCTATCGCCAGCACCATCCCGAACa                                  <  1:1477405/53‑1 (MQ=21)
                                                  accacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGa                   >  1:1766881/1‑67 (MQ=255)
                                                   ccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGa                                     <  1:3112642/48‑1 (MQ=255)
                                                   ccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGa                   >  1:4011644/1‑66 (MQ=255)
                                                   ccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTg              <  1:3833700/71‑1 (MQ=255)
                                                    cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTg              >  1:1881580/1‑70 (MQ=255)
                                                    cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa             <  1:2403952/71‑1 (MQ=255)
                                                    cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa             <  1:1635053/71‑1 (MQ=255)
                                                    cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa             <  1:3836052/71‑1 (MQ=255)
                                                     acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa             >  1:3312876/1‑70 (MQ=255)
                                                     acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa             >  1:2046148/1‑70 (MQ=255)
                                                     acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACAATCCCGAACATAGTCAACGTGTTGATGGTGa             >  1:3819739/1‑70 (MQ=255)
                                                     acGATGGCGTCATCCACCACTAACCCTATCGACAGCAAGATCCCGAACATAGTCAACGTGTAGAAGGt               >  1:924041/1‑68 (MQ=25)
                                                         tGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa            <  1:3964587/67‑1 (MQ=255)
                                                         tGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa            <  1:35588/67‑1 (MQ=255)
                                                         tGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAacc         >  1:2920476/1‑70 (MQ=255)
                                                           gCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAaccgac      >  1:3746943/1‑71 (MQ=255)
                                                            cGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                          >  1:4051349/1‑50 (MQ=255)
                                                            cGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAaccgacc     >  1:1517052/1‑71 (MQ=255)
                                                            cGTCATCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAaccgacc     <  1:3706088/70‑1 (MQ=255)
                                                              tCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGtt                     <  1:2191374/53‑1 (MQ=255)
                                                               cATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc  >  1:1595167/1‑71 (MQ=255)
                                                               cATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc  >  1:2909912/1‑71 (MQ=255)
                                                                aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAaccgac      >  1:338247/1‑66 (MQ=255)
                                                                aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc  >  1:3148884/1‑70 (MQ=255)
                                                                 tCCACCAGTAACCCTATCGCTAGCACCATCCCGAACATAGTc                             <  1:2491892/42‑1 (MQ=255)
                                                                 tCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgtg                       <  1:2125051/48‑1 (MQ=255)
                                                                   |                                                                  
CCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCC  >  minE/2585992‑2586125

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: