Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,587,712 G→A 45.4% R249R (CGC→CGT mdtE ← multidrug resistance efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,587,7120GA45.4% 12.0 / 51.0 66R249R (CGC→CGTmdtEmultidrug resistance efflux transporter
Reads supporting (aligned to +/- strand):  ref base G (8/28);  new base A (20/10);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 4.12e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTCACGGAGCCCGTGGTTTCATCCACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGATTTGCCCACTGGC  >  minE/2587647‑2587780
                                                                 |                                                                    
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:2264476/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                               >  1:2248026/1‑69 (MQ=255)
                    aTCCACCGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAg                                               >  1:2202871/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:1912672/68‑1 (MQ=255)
                                gtcggAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCg                                    <  1:1848933/68‑1 (MQ=255)
                                gtcggAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCg                                    <  1:3132426/68‑1 (MQ=255)
                                gtcggAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCg                                    <  1:3009437/68‑1 (MQ=255)
                                gtcggAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCg                                    <  1:1377712/68‑1 (MQ=255)
                                gtcggAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCg                                    <  1:2920119/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                               <  1:1823752/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                               <  1:701585/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                               <  1:339121/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                               <  1:1340359/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                               <  1:825021/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                               <  1:88257/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                               <  1:1326692/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                               <  1:2291813/68‑1 (MQ=255)
                                     agaATTTCAGCGTGCCGGTCTGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                               <  1:2472215/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:2015536/1‑69 (MQ=255)
                                                   ccGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTg                >  1:247945/1‑69 (MQ=255)
                                                   ccGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTg                >  1:2171403/1‑69 (MQ=255)
                                                     ggTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGa               <  1:66385/68‑1 (MQ=255)
                                                     ggTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGa               <  1:2295251/68‑1 (MQ=255)
                                                     ggTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGa               <  1:224992/68‑1 (MQ=255)
                                                     ggTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGa               <  1:3244596/68‑1 (MQ=255)
                                                     ggTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGa               <  1:3372011/68‑1 (MQ=255)
                                                     ggTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGa               <  1:3246172/68‑1 (MQ=255)
                                                         tGGCCGTAACGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTg                               <  1:20886/48‑1 (MQ=255)
                                                             tGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGATTTGcc         <  1:294737/66‑1 (MQ=255)
                                                                 gCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGATTTGCCCACTGGc  >  1:3181361/1‑69 (MQ=255)
                                                                 |                                                                    
GTCACGGAGCCCGTGGTTTCATCCACTGTCGGGTCGGAGAATTTCAGCGTGCCGGTCTGGCTGTAGCGTTTACCATTTTCCAGATTGAGCTGTACTGGCGTACTGCCCTGAACCTGTTTGATTTGCCCACTGGC  >  minE/2587647‑2587780

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: