Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,395,332 +G 14.3% coding (226/750 nt) yjfP → predicted hydrolase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,395,3321.→G14.3% 9.2 ‑2.8 ‑0.221coding (226/750 nt)yjfPpredicted hydrolase
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (6/12):  total (9/12)
Fisher's exact test strand distribution p-value = 6.32e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.72e-01

ATCACGGTAGTCGTTTTAGTGGCGACGCAGCGCGGC.GGTTAAATCAATTCTGGCAAATCTTGCTACAAAGT  >  REL606/4395297‑4395367
                                    |                                   
aTCACGGTAGTCGTTTTAGTGGCGACGCAGCgcggc                                      >  1:851736/1‑36 (MQ=255)
   aCGTTAGTCCTTTTAGTGGCGACGCAGCGCGGC.GGt                                  <  1:2194347/36‑1 (MQ=38)
   aCGGTAGTCGTTTTAGTGGCGACGCCCCGcggc.ggg                                  >  1:1802633/1‑35 (MQ=37)
     ggTAGTCGTTTTAGTGGCGACGCAGCGCGGC.GGTTa                                <  1:5083092/36‑1 (MQ=255)
         gTCGTTTTAGTGGCGAGGCAGCGAGGCGGGTTaaaa                             >  1:1496332/1‑35 (MQ=16)
         gTCGTTTTAGTGGCGACGCAGCGCGGC.GGTTAAATc                            <  1:1743666/36‑1 (MQ=255)
            gTTTTAGTGGCGACGCAGCGCGGCGGGTTAATTCat                          >  1:722978/1‑35 (MQ=21)
            gTTTTAGTGGCGACGCAGCGCGGCGGGTTAAAtaat                          >  1:4809268/1‑33 (MQ=25)
            gTTTTAGTGGCGACGCAGCGCGGC.GGGTTAATTAAt                         >  1:3210590/1‑36 (MQ=25)
             ttttAGTGGCGACGCAGCGCGGC.GGTTAAATCAAtt                        <  1:1098702/36‑1 (MQ=255)
             ttttAGTGGCGACGCAGCGCGGC.GGGTATATAaaat                        >  1:2937929/1‑34 (MQ=25)
              tttAGTGGCGACGCAGCGCGGC.GGTTAAATCAATTc                       <  1:5270529/36‑1 (MQ=255)
              tttAGTGGCGACGCAGCGCGGC.GGTTAAATCAATTc                       <  1:4987639/36‑1 (MQ=255)
               ttAGTGGCGACGCAGCGCGGC.GTTTAAATCAATTCt                      <  1:1453248/36‑1 (MQ=255)
                 aGTGGCGACGCAGCGCGGC.GGGTTAATTAAATCttg                    >  1:1643970/1‑34 (MQ=16)
                  gTGGCGACGCAGCGCGGC.GGTTAAATCAATTCTGGc                   <  1:56583/36‑1 (MQ=255)
                  gTGGCGACGCAGCGCGGC.GGTTAAATCAATTCTGGc                   <  1:4885904/36‑1 (MQ=255)
                  gTGGCGACGCAGCGCGGC.GGTTAAATCAATTCTGGc                   <  1:4047621/36‑1 (MQ=255)
                  gTGGCGACGCAGCGCGGC.GGTTAAATCAATTCTGGc                   <  1:3868397/36‑1 (MQ=255)
                       gACGCAGCGCGGC.GGGTAAAACAAATCTTGCgaaac              >  1:5780157/1‑34 (MQ=9)
                         cGCAGCGCGGC.GGTTTAATTAATTCTTGCAAATCtt            >  1:1978255/1‑36 (MQ=25)
                          gCAGCGCGGC.GGTTTAATCAAATCTGGCAAAACttt           >  1:1817090/1‑35 (MQ=21)
                            aGCGCGGC.GGTTAAATCAATTCTGGCAAATCTTGCt         <  1:1390620/36‑1 (MQ=255)
                               gctgc.ggTTCAATCAATTCTGGCAAATCTTGCTACa      <  1:4264401/33‑1 (MQ=37)
                               gcggc.gGTTAAATCAATTCTGGCAAATCTTGCTACa      <  1:2602171/36‑1 (MQ=255)
                                cggc.ggTTAAATCAATTCTGGCAAATCTTGCTACaa     <  1:3531589/36‑1 (MQ=255)
                                 ggc.ggTTAAATCAATTCTTGCCAGTCTTGCTACaaa    <  1:4743382/36‑1 (MQ=25)
                                 ggc.ggTTAAATCAATTCTGGCAAATCTTTCTACaaa    >  1:496186/1‑36 (MQ=255)
                                 ggc.ggTTAAATCAATTCTGGCAAATCTTGCTACaaa    <  1:5127009/36‑1 (MQ=255)
                                 ggc.ggTTAAATCAATTCTGGCAAATCTTGCTACaaa    <  1:1298204/36‑1 (MQ=255)
                                  gc.ggTTAAATCAATTCTGGCAAATCTTGCTACAAAg   >  1:529184/1‑36 (MQ=255)
                                   c.ggTTAAATCAATTCTGGCAAATCTTGCTACAAAGt  <  1:5447905/36‑1 (MQ=255)
                                    |                                   
ATCACGGTAGTCGTTTTAGTGGCGACGCAGCGCGGC.GGTTAAATCAATTCTGGCAAATCTTGCTACAAAGT  >  REL606/4395297‑4395367

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 13 ≤  ATCG  < 26 ≤  ATCG  < 37 ≤  ATCG