Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,078 G→C 51.5% V265V (GTC→GTG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,0780GC51.5% ‑5.1 / 80.5 66V265V (GTC→GTGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (17/15);  new base C (19/15);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

CCGTTAAACTCTGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTC  >  minE/312025‑312143
                                                     |                                                                 
ccGTTAAACTCTGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACgcagcag                                                     <  1:1576463/68‑1 (MQ=255)
ccGTTAAACTCTGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACgcagcag                                                     <  1:2100040/68‑1 (MQ=255)
     aaaCTCTGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCAcgc                                                >  1:3373956/1‑68 (MQ=255)
     aaaCTCTGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCAcg                                                 >  1:2173572/1‑67 (MQ=255)
        ctctGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGa                                            <  1:2103654/69‑1 (MQ=255)
        ctctGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGa                                            <  1:2318535/69‑1 (MQ=255)
        ctctGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGa                                            <  1:1677583/69‑1 (MQ=255)
        ctctGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGa                                            <  1:139572/69‑1 (MQ=255)
        ctctGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGa                                            <  1:3010153/69‑1 (MQ=255)
        ctctGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGa                                            <  1:3228847/69‑1 (MQ=255)
        ctctGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGa                                            <  1:1285871/69‑1 (MQ=255)
        ctctGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGGCACGCAGCAGCACGCGGGa                                            <  1:471545/69‑1 (MQ=255)
                          gTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAAcc                                         >  1:312660/1‑54 (MQ=255)
                          gTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACt                            >  1:3331587/1‑67 (MQ=255)
                          gTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACt                            >  1:1839509/1‑67 (MQ=255)
                          gTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACt                            >  1:2393231/1‑67 (MQ=255)
                          gTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACt                            >  1:1653997/1‑67 (MQ=255)
                          gTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCATGATTCACt                            >  1:1573893/1‑67 (MQ=255)
                             ctcACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTc                         >  1:1800483/1‑67 (MQ=255)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATGCTGATTCACTTTcagcag                    <  1:3233701/67‑1 (MQ=255)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagtag                    <  1:2544803/67‑1 (MQ=37)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                    <  1:206128/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                    <  1:2526230/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                    <  1:1110998/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                    <  1:2762564/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                    <  1:624887/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                    <  1:1416420/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                    <  1:1373969/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                    <  1:1378122/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                    <  1:1352721/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                    <  1:1294240/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                    <  1:1637512/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                    <  1:1274927/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                    <  1:1922523/67‑1 (MQ=38)
                                          tCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCg          <  1:3238717/69‑1 (MQ=255)
                                          tCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCg          <  1:2736906/69‑1 (MQ=255)
                                            aaTTTTCGCCACATCGCGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGa         >  1:2782144/1‑68 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTg                                   >  1:2254960/1‑42 (MQ=21)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTg             >  1:2588721/1‑64 (MQ=255)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa        >  1:443847/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa        >  1:33225/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa        >  1:2943572/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa        >  1:3250340/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa        >  1:446922/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa        >  1:1135944/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa        >  1:553757/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa        >  1:2873955/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa        >  1:2012804/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa        >  1:2710901/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa        >  1:2589232/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa        >  1:1641630/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa        >  1:1889840/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa        >  1:1719141/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGa         >  1:2712103/1‑68 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGa         >  1:618401/1‑68 (MQ=25)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa        <  1:724151/67‑1 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa        <  1:1099736/67‑1 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAAct      >  1:751086/1‑69 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAAct      >  1:57438/1‑69 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAAct      >  1:453572/1‑69 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAAct      >  1:1691461/1‑69 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAAct      >  1:2262731/1‑69 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAAct      >  1:3293436/1‑69 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAAct      >  1:2656429/1‑69 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAAct      >  1:1353835/1‑69 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTAAGCAGGATTTTGCCGaa        <  1:151253/67‑1 (MQ=255)
                                                  cgcgACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTc  >  1:1195304/1‑69 (MQ=255)
                                                     |                                                                 
CCGTTAAACTCTGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTC  >  minE/312025‑312143

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: