Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,090 C→T 47.9% L261L (CTG→CTA acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,0900CT47.9% 6.1 / 60.2 71L261L (CTG→CTAacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (18/19);  new base T (19/15);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

CCGTTAAACTCTGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAG  >  minE/312025‑312148
                                                                 |                                                          
ccGTTAAACTCTGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACgcagcag                                                          <  1:2100040/68‑1 (MQ=255)
ccGTTAAACTCTGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACgcagcag                                                          <  1:1576463/68‑1 (MQ=255)
     aaaCTCTGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCAcgc                                                     >  1:3373956/1‑68 (MQ=255)
     aaaCTCTGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCAcg                                                      >  1:2173572/1‑67 (MQ=255)
        ctctGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGa                                                 <  1:139572/69‑1 (MQ=255)
        ctctGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGa                                                 <  1:2103654/69‑1 (MQ=255)
        ctctGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGa                                                 <  1:1677583/69‑1 (MQ=255)
        ctctGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGa                                                 <  1:2318535/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:1839509/1‑67 (MQ=255)
                          gTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACt                                 >  1:2393231/1‑67 (MQ=255)
                          gTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACt                                 >  1:1653997/1‑67 (MQ=255)
                          gTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACt                                 >  1:3331587/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:2526230/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                         <  1:1922523/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                         <  1:206128/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:1378122/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                         <  1:1274927/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                         <  1:1416420/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                         <  1:1373969/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                         <  1:1294240/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                         <  1:1352721/67‑1 (MQ=38)
                                  caccCAGCTCAATTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTcagcag                         <  1:1637512/67‑1 (MQ=38)
                                          tCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCg               <  1:2736906/69‑1 (MQ=255)
                                          tCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCg               <  1:3238717/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:446922/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa             >  1:3250340/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa             >  1:2943572/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa             >  1:2873955/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa             >  1:33225/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa             >  1:2710901/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa             >  1:1135944/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa             >  1:2589232/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa             >  1:553757/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa             >  1:1641630/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa             >  1:1719141/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa             >  1:2012804/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa             >  1:1889840/1‑69 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGa              >  1:2712103/1‑68 (MQ=25)
                                            aaTTTTCGCCACATCACGCAGTAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGa              >  1:618401/1‑68 (MQ=25)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa             <  1:1099736/67‑1 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGaa             <  1:724151/67‑1 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAAct           >  1:751086/1‑69 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAAct           >  1:453572/1‑69 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAAct           >  1:57438/1‑69 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAAct           >  1:3293436/1‑69 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAAct           >  1:1353835/1‑69 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAAct           >  1:2656429/1‑69 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAAct           >  1:2262731/1‑69 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAAct           >  1:1691461/1‑69 (MQ=255)
                                              tCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTAAGCAGGATTTTGCCGaa             <  1:151253/67‑1 (MQ=255)
                                                  cgcgACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTc       >  1:1195304/1‑69 (MQ=255)
                                                       cGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTag  <  1:2952417/69‑1 (MQ=255)
                                                       cGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTag  <  1:2785909/69‑1 (MQ=255)
                                                       cGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTag  <  1:466893/69‑1 (MQ=255)
                                                       cGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTag  <  1:1525744/69‑1 (MQ=255)
                                                       cGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTag  <  1:1680373/69‑1 (MQ=255)
                                                       cGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTag  <  1:1881060/69‑1 (MQ=255)
                                                                 |                                                          
CCGTTAAACTCTGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAG  >  minE/312025‑312148

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: