Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,377,653 G→A 20.7% R245R (CGC→CGT ubiB ← 2‑octaprenylphenol hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,377,6530GA20.7% 134.4 / 23.1 91R245R (CGC→CGTubiB2‑octaprenylphenol hydroxylase
Reads supporting (aligned to +/- strand):  ref base G (41/31);  new base A (8/11);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 3.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

CAGCAATTTCATGTTAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTCCGGGATGTAGAGC  >  minE/2377587‑2377713
                                                                  |                                                            
cAGCAATTTCATGTTAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATgcg                                                            <  1:1145573/69‑1 (MQ=255)
cAGCAATTTCATGTTAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATgcg                                                            <  1:281103/69‑1 (MQ=255)
cAGCAATTTCATGTTAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATgcg                                                            <  1:19449/69‑1 (MQ=255)
cAGCAATTTCATGTTAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATgcg                                                            <  1:1458299/69‑1 (MQ=255)
cAGCAATTTCATGTTAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATgcg                                                            <  1:1442682/69‑1 (MQ=255)
             tttgTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTccatcaccatc                                               >  1:286635/4‑69 (MQ=255)
             ttAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTccatcaccatc                                               >  1:58394/1‑69 (MQ=255)
             ttAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTccatcaccatc                                               >  1:1601649/1‑69 (MQ=255)
             ttAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTccatcaccatc                                               >  1:326249/1‑69 (MQ=255)
             ttAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTccatcaccatc                                               >  1:2557997/1‑69 (MQ=255)
             ttAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTccatcaccatc                                               >  1:1504535/1‑69 (MQ=255)
             ttAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTccatcaccatc                                               >  1:2469356/1‑69 (MQ=255)
             ttAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTccatcaccatc                                               >  1:1902921/1‑69 (MQ=255)
             ttAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTccatcaccatc                                               >  1:547035/1‑69 (MQ=255)
             ttAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTccatcaccatc                                               >  1:1236162/1‑69 (MQ=255)
             ttAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTccatcaccatc                                               >  1:1218131/1‑69 (MQ=255)
             ttAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTccatcaccatc                                               >  1:534357/1‑69 (MQ=255)
             ttAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTccatcaccatc                                               >  1:456924/1‑69 (MQ=255)
                  gctgTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATccc                                          <  1:1466400/66‑1 (MQ=255)
                   ccGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATccc                                          >  1:1023231/1‑68 (MQ=255)
                   ccGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATccc                                          >  1:254767/1‑68 (MQ=255)
                   ccGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATccc                                          >  1:306253/1‑68 (MQ=255)
                   ccGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATccc                                          >  1:3142068/1‑68 (MQ=255)
                       ttttCTCCAGCGTTGCAACATCAGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCAc                                     >  1:2058511/1‑69 (MQ=255)
                       ttttCTCCAGCGTCGCAACATCAGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCa                                      >  1:681017/1‑68 (MQ=255)
                       ttttCTCCAGCGTCGCAACATCAGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCa                                      >  1:2899432/1‑68 (MQ=255)
                       ttttCTCCAGCGTCGCAACATCAGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCAc                                     >  1:821217/1‑69 (MQ=255)
                       ttttCTCCAGCGTCGCAACATCAGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCAc                                     >  1:3327206/1‑69 (MQ=255)
                       ttttCTCCAGCGTCGCAACATCAGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCAc                                     >  1:1849843/1‑69 (MQ=255)
                       ttttCTCCAGCGTCGCAACATCAGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCAc                                     >  1:1242161/1‑69 (MQ=255)
                       ttttCTCCAGCGTCGCAACATCAGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCAc                                     >  1:2267879/1‑69 (MQ=255)
                          tctcCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATgcg                                                            <  1:1563111/43‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:2481329/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:1894779/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:3370967/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:2993215/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:2875238/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:2829641/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:2598033/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:2775032/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:3252956/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:2025060/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:115354/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:1280910/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:823424/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:1505744/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:784290/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:1663063/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:1667021/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:1704787/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:1881861/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:659698/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:1873934/69‑1 (MQ=255)
                                       aaCATCAGACACCGGAATGCCGTAAATGCGCTCCATAACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:618577/69‑1 (MQ=255)
                                           tCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:3104470/65‑1 (MQ=255)
                                           tCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:1263341/65‑1 (MQ=255)
                                           tCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:86090/65‑1 (MQ=255)
                                           tCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:1603097/65‑1 (MQ=255)
                                           tCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:3319444/65‑1 (MQ=255)
                                           tCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     <  1:2684561/65‑1 (MQ=255)
                                           tCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAAc                  >  1:2970424/1‑68 (MQ=255)
                                           tCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACt                 >  1:1728677/1‑69 (MQ=255)
                                           tCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACt                 >  1:2628627/1‑69 (MQ=255)
                                             aGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTc               <  1:89011/69‑1 (MQ=255)
                                             aGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTc               <  1:921376/69‑1 (MQ=255)
                                             aGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTc               <  1:3129882/69‑1 (MQ=255)
                                             aGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTc               <  1:38834/69‑1 (MQ=255)
                                             aGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTc               <  1:1595965/69‑1 (MQ=255)
                                             aGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTc               <  1:2797829/69‑1 (MQ=255)
                                             aGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTc               <  1:2757761/69‑1 (MQ=255)
                                             aGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTc               <  1:2547058/69‑1 (MQ=255)
                                             aGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTc               <  1:2224963/69‑1 (MQ=255)
                                             aGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTc               <  1:1332600/69‑1 (MQ=255)
                                             aGACACCGGAATGCCGTAAATACGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTc               <  1:2091225/69‑1 (MQ=255)
                                                cacCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTCCgg            >  1:1517767/1‑69 (MQ=255)
                                                cacCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTCCgg            >  1:344376/1‑69 (MQ=255)
                                                    ggAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTc                           >  1:367856/1‑50 (MQ=255)
                                                    ggAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTCCGGGa          >  1:2923424/1‑67 (MQ=255)
                                                    ggAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTCCGGGa          >  1:1776078/1‑67 (MQ=255)
                                                    ggAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTCCGGGa          >  1:2066765/1‑67 (MQ=255)
                                                    ggAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTCCGGGa          >  1:488150/1‑67 (MQ=255)
                                                    ggAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTCCGGGa          >  1:280756/1‑67 (MQ=255)
                                                          ccGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTCCGGGATGTAGAGc  <  1:2153293/69‑1 (MQ=255)
                                                          ccGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTCCGGGATGTAGAGc  <  1:1960327/69‑1 (MQ=255)
                                                               aaTGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGt                      >  1:134298/1‑44 (MQ=255)
                                                               aaTGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     >  1:869270/1‑45 (MQ=255)
                                                               aaTGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     >  1:849886/1‑45 (MQ=255)
                                                               aaTGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     >  1:2969024/1‑45 (MQ=255)
                                                               aaTGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     >  1:1366813/1‑45 (MQ=255)
                                                               aaTGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     >  1:1096081/1‑45 (MQ=255)
                                                               aaTGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     >  1:1462694/1‑45 (MQ=255)
                                                               aaTGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     >  1:1771984/1‑45 (MQ=255)
                                                               aaTGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     >  1:1851436/1‑45 (MQ=255)
                                                               aaTGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     >  1:492691/1‑45 (MQ=255)
                                                               aaTGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     >  1:415497/1‑45 (MQ=255)
                                                               aaTGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     >  1:3392464/1‑45 (MQ=255)
                                                               aaTGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTa                     >  1:328727/1‑45 (MQ=255)
                                                                  |                                                            
CAGCAATTTCATGTTAGTGCCGTTTTTCTCCAGCGCCGCAACATCAGACACCGGAATGCCGTAAATGCGCTCCATCACCATCATCCCTTCACTACAATAGTCAGGGTAAACTTCCGGGATGTAGAGC  >  minE/2377587‑2377713

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: