Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,316,032 +G 9.4% coding (54/1968 nt) yhcP ← p‑hydroxybenzoic acid efflux system component

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,316,0321.→G 9.4% 3.2 ‑2.1 0.032coding (54/1968 nt)yhcPp‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (9/20):  total (12/20)
Fisher's exact test strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTC.GCCAGTTTTACCGCAAAGCGAATATGTTGATTAGC  >  REL606/3315997‑3316067
                                    |                                   
tGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTc                                      <  1:7282890/36‑1 (MQ=255)
tGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTc                                      <  1:6078117/36‑1 (MQ=255)
 gAAAGCCAACAAACAGCGCCAGTACAATGGCGGGC.g                                    >  1:4859775/1‑36 (MQ=255)
  aaaGCCAACAAACAGCGCCAGTACAATGGCGGGC.gg                                   >  1:6549148/1‑35 (MQ=255)
   aatccAACAAACAGCGCCAGTACAATGGCGGTC.Gcc                                  <  1:2771699/33‑1 (MQ=255)
   aaGCCAACAAACAGCGCCAGTACAATGGCGGGC.ggc                                  >  1:3113358/1‑34 (MQ=37)
   aaGCCAACAAACAGCGCCAGTACAATGGCGGGC.Gcc                                  >  1:3116835/1‑36 (MQ=255)
    aGCCAACAAACAGCGCCAGTACAATGGCGGTC.GCCa                                 <  1:1813903/36‑1 (MQ=255)
    aGCCAACAAACAGCGCCAGTACAATCGCGGTC.ggcc                                 >  1:3297385/1‑33 (MQ=37)
     gCCAACAAACAGCGCCAGTACAATGGCGGTC.GCCAg                                <  1:5761263/36‑1 (MQ=255)
      ccAACAAACATCGCCAGTACAAAGGCGGTC.GCCAGt                               <  1:89120/36‑1 (MQ=37)
        aacaaacaGCGCCAGTACAATGGCGGTG.GGCCGttt                             >  1:2352210/1‑36 (MQ=25)
        aacaaacaGCGCCAGTACAATGGCGGTC.GCCAGttt                             <  1:7335281/36‑1 (MQ=255)
         acaaacaGCGCCAGTACAATGGCGGTC.GCCAGtttt                            <  1:3467642/36‑1 (MQ=255)
          caaacaGCGCCAGTACAATGGCGGGCGGCCGGtttt                            >  1:7964772/1‑36 (MQ=19)
            aacaGCGCCAGTACAATGGCGGTG.GCCAGgtttttc                         >  1:6042671/1‑33 (MQ=25)
             acaGCGCCAGTACAATGGCGGTCTGCCGGTTTTAcc                         >  1:1328495/1‑36 (MQ=25)
             acaGCGCCAGTACAATGGCGGTCGGCCGGtttttcg                         >  1:6899042/1‑33 (MQ=19)
              caGCGCCAGTACAATGGCGGTC.GCCAGTTTTACCGc                       <  1:4524348/36‑1 (MQ=255)
              caGCGCCAGTACAATGGCGGTC.GCCAGTTTTACCGc                       <  1:2303635/36‑1 (MQ=255)
               aGCGCCAGTACAATGGCGGTC.GCCAGTTTTACCGCa                      <  1:1299452/36‑1 (MQ=255)
               aGCGCCAGTACAATGGCGGTC.GCCAGTTTTACCGCa                      <  1:5506276/36‑1 (MQ=255)
               aGCGCCAGTACAATGGCGGTC.GCCAGTTTTACCGCa                      <  1:2984980/36‑1 (MQ=255)
               aGCGCCAGTACAATGGCGGGC.GCCCGTTTTTTCGCa                      >  1:85819/1‑36 (MQ=16)
                 cgcCAGTACAATGGCGGTC.GCCAGTTTTACCGCaaa                    <  1:6698598/36‑1 (MQ=255)
                  gcCAGTACAATGGCGGTC.GCCAGTTTTACCGCAAAg                   <  1:7718821/36‑1 (MQ=255)
                   ccAGTACAATGGCGGTC.GCCAGTTTTACCGCAAAGc                  <  1:7261866/36‑1 (MQ=255)
                        aCAATGGCGGTC.GCCAGTTTTACCGCAAAGCGAata             <  1:5847966/36‑1 (MQ=255)
                        aCAATGGCGGTC.GCCAGTTTTACCGCAAAGCGAata             <  1:130402/36‑1 (MQ=255)
                            tGGCGGTC.GCCAGTTTTACCGCAAAGCGAATATGtt         <  1:7831700/36‑1 (MQ=255)
                             ggCGGTC.GCCAGTTTTACCGCAAAGCGAATATGTTg        >  1:7941517/1‑36 (MQ=255)
                               cGGTC.GCCAGTTTTACCGCAAAGCGAATATGTTGAt      >  1:1036877/1‑36 (MQ=255)
                               cGGTC.GCCAGTTTTACCGCAAAGCGAATATGTTGAt      <  1:7667266/36‑1 (MQ=255)
                               cGGTC.GCCAGTTTTACCGCAAAGCGAATATGTTGAt      >  1:3868335/1‑36 (MQ=255)
                                ggTC.GCCAGTTTTACCGCAAAGCGAATATGTTGAtt     >  1:7363031/1‑36 (MQ=255)
                                ggTC.GCCAGTTTTACCGCAAAGCGAATATGTTGAtt     >  1:8260704/1‑36 (MQ=255)
                                 gTC.GCCAGTTTTACCGCAAATAGAATATGTTGATTa    <  1:808361/36‑1 (MQ=37)
                                  tC.GCCAGTTTTACCGCAAAGCGAATATGTTGATTAg   <  1:4630944/36‑1 (MQ=255)
                                   c.GCCAGTTTTACCGCAAAGCGTATATGTTGATTAGc  <  1:2905065/36‑1 (MQ=255)
                                    |                                   
TGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTC.GCCAGTTTTACCGCAAAGCGAATATGTTGATTAGC  >  REL606/3315997‑3316067

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG