Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,587,708 T→C 46.2% S251G (AGC→GGC)  mdtE ← multidrug resistance efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,587,7080TC46.2% 3.7 / 54.2 65S251G (AGC→GGC) mdtEmultidrug resistance efflux transporter
Reads supporting (aligned to +/- strand):  ref base T (8/27);  new base C (20/10);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 4.81e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

GTCACGGAGCCCGTGGTTTCATCCACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGATTTGCC  >  minE/2587647‑2587773
                                                             |                                                                 
gTCACGGAGCCCGTGGTTTCATCCACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCg                                                             <  1:1204937/68‑1 (MQ=255)
gTCACGGAGCCCGTGGTTTCATCCACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCg                                                             <  1:21138/68‑1 (MQ=255)
              ggTTTCATCCACTGTCGGGTTGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAg                                              <  1:2344097/69‑1 (MQ=255)
              ggTTTCATCCACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAg                                              <  1:2542163/69‑1 (MQ=255)
              ggTTTCATCCACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAg                                              <  1:2596977/69‑1 (MQ=255)
              ggTTTCATCCACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAg                                              <  1:1418123/69‑1 (MQ=255)
              ggTTTCATCCACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAg                                              <  1:1430114/69‑1 (MQ=255)
              ggTTTCATCCACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAg                                              <  1:2138441/69‑1 (MQ=255)
              ggTTTCATCCACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAg                                              <  1:3178914/69‑1 (MQ=255)
              ggTTTCATCCACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTACCATTTTCCAg                                              <  1:1443334/68‑1 (MQ=255)
               gTTTCATCCACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGa                                             >  1:2354922/1‑69 (MQ=255)
               gTTTCATCCACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGa                                             >  1:2702675/1‑69 (MQ=255)
               gTTTCATCCACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGa                                             >  1:3129476/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGGCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:1546154/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTcc                                                >  1:1346458/1‑61 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGa                                         >  1:3279920/1‑68 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGa                                         >  1:3390262/1‑68 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:3192152/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:2995316/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:2969939/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:1112156/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:2836769/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:53074/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:667615/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:669185/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:2476623/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:2407440/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:790031/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:978791/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:2292017/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:2202871/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:2248026/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                        >  1:2264476/1‑69 (MQ=255)
                       cACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAg                                              >  1:2013957/1‑60 (MQ=255)
                       cACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTg                                     <  1:3046563/69‑1 (MQ=255)
                       cACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTg                                     <  1:2057577/69‑1 (MQ=255)
                         cTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTa                                   <  1:757554/69‑1 (MQ=255)
                                gtcggAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCg                             <  1:3132426/68‑1 (MQ=255)
                                gtcggAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCg                             <  1:1377712/68‑1 (MQ=255)
                                gtcggAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCg                             <  1:1848933/68‑1 (MQ=255)
                                gtcggAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCg                             <  1:1912672/68‑1 (MQ=255)
                                gtcggAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCg                             <  1:3009437/68‑1 (MQ=255)
                                gtcggAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCg                             <  1:2920119/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                        <  1:88257/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                        <  1:1326692/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                        <  1:1823752/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                        <  1:339121/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                        <  1:2472215/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                        <  1:701585/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                        <  1:2291813/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                        <  1:1340359/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                        <  1:825021/68‑1 (MQ=255)
                                         tttCAGCGTGCCGGTCTGGCTGTAGCTTTTACCATTTTCCAGATTGAGCTg                                     <  1:1871559/51‑1 (MQ=255)
                                         tttCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCg                             <  1:1361297/59‑1 (MQ=255)
                                                   ccGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAg                                        >  1:538107/1‑38 (MQ=255)
                                                   ccGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTg         >  1:247945/1‑69 (MQ=255)
                                                   ccGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTg         >  1:2171403/1‑69 (MQ=255)
                                                   ccGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTg         >  1:2015536/1‑69 (MQ=255)
                                                     ggTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGa        <  1:3246172/68‑1 (MQ=255)
                                                     ggTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGa        <  1:66385/68‑1 (MQ=255)
                                                     ggTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGa        <  1:3372011/68‑1 (MQ=255)
                                                     ggTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGa        <  1:3244596/68‑1 (MQ=255)
                                                     ggTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGa        <  1:2295251/68‑1 (MQ=255)
                                                     ggTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGa        <  1:224992/68‑1 (MQ=255)
                                                         tGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                        <  1:20886/48‑1 (MQ=255)
                                                             tGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGATTTGcc  <  1:294737/66‑1 (MQ=255)
                                                             |                                                                 
GTCACGGAGCCCGTGGTTTCATCCACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGATTTGCC  >  minE/2587647‑2587773

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: