Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,258,401 C→G 21.7% G330G (GGC→GGG yjcD → predicted permease

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,258,4010C→G21.7% 2.8 ‑3.4 ‑0.031990yjcDpredicted permease
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (13/11):  total (20/11)
Fisher's exact test strand distribution p-value = 3.31e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.55e-01

GTATACATCGAGTCTGCGGCGGGTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGT  >  REL606/4258366‑4258436
                                    |                                   
gTATACATCGAGTCTGCGGCGGGTACGGC.ggcgggc                                     <  1:459457/36‑1 (MQ=255)
gTATACATCGAGTCTGCGGCGGGTACGGC.ggcgggc                                     <  1:1237444/36‑1 (MQ=255)
 tataCATCGAGTCTGCGGCGGGTACGGC.ggcgggcg                                    <  1:3396368/36‑1 (MQ=255)
  ataCATCGAGTCTGCGGCGGGTACGGC.Ggcgggcgg                                   <  1:2751417/36‑1 (MQ=255)
   taCATCGAGTCTGCGGCGGGTACGGC.GGGGGGGGGt                                  >  1:4992205/1‑36 (MQ=38)
   taCATCGAGTCTGCGGCGGGTACGGC.GGCGGGCGGt                                  <  1:6770554/36‑1 (MQ=255)
    aCATCGAGTCTGCGGCGGGTACGGC.Ggggggagggg                                 >  1:2127201/1‑34 (MQ=37)
    aCATCGAGTCTGCGGCGGGTACGGC.GGGGGGGGGtt                                 >  1:6474325/1‑35 (MQ=37)
     cATCGAGTCTGCGGCGGGTACCGA.GGCGGGGGGtta                                >  1:5689296/1‑34 (MQ=25)
      aTCGAGTCTGCGGCGGGTACGGC.GGGGGGGGGTtaa                               >  1:6760379/1‑36 (MQ=37)
      aTCGAGTCTGCGGCGGGTACGGC.GGGGGGCGGTgaa                               >  1:4121442/1‑36 (MQ=37)
      aTCGAGTCTGCGGCGGGTACGGC.GGCGGGCGGTaaa                               >  1:3428270/1‑36 (MQ=255)
         gAGTCTGCGGCGGGTACGGC.GGGGGGCGGTaaaaac                            >  1:518850/1‑34 (MQ=38)
            tCTGCGGCGGGTACGGC.GGCGGGCGGTAAAACCGGt                         <  1:3345093/36‑1 (MQ=255)
             cTGCGGCGGGTACGGC.GGCGGGCGGTCAAACCgggt                        >  1:1239698/1‑34 (MQ=37)
              tGCGGCGGGTACGGCGGGGGGGGGGGAAAACCGGtt                        >  1:6904346/1‑36 (MQ=16)
              tGCGGCGGGTACGGC.GGGGGGGGGGGAAAACGGgtt                       >  1:2619389/1‑36 (MQ=12)
                cggcggGTACGGC.GGCGGGCGGTAAAACCGGTTTGa                     <  1:266521/36‑1 (MQ=255)
                 ggcggGTACGGC.GGCGGGGGGTAAAACCGGTTTggc                    >  1:4053422/1‑34 (MQ=37)
                    gggTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGc                 <  1:7260964/36‑1 (MQ=255)
                    gggTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGc                 >  1:6878771/1‑36 (MQ=255)
                    gggTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGc                 >  1:4345657/1‑36 (MQ=255)
                     ggTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGCt                >  1:6948390/1‑36 (MQ=255)
                      gTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTa               >  1:727101/1‑36 (MQ=255)
                      gTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTa               >  1:6639261/1‑36 (MQ=255)
                       tACGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTAt              <  1:819816/36‑1 (MQ=255)
                       tACGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTAt              <  1:6928010/36‑1 (MQ=255)
                       tACGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTAt              <  1:4854260/36‑1 (MQ=255)
                       tACGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTAt              <  1:5662650/36‑1 (MQ=255)
                       tACGGC.GGAGGGCGGTAAAACCGGTTTGACGGCTAt              <  1:1984006/36‑1 (MQ=255)
                        aCGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTAtt             >  1:6304843/1‑35 (MQ=255)
                        aCGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTATc             >  1:353450/1‑36 (MQ=255)
                        aCGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTATc             <  1:1966418/36‑1 (MQ=255)
                        aCGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTATc             >  1:3731765/1‑36 (MQ=255)
                        aCGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTATc             >  1:2429241/1‑36 (MQ=255)
                            c.ggcggGCGGTAAAACCGGTTTGACGGCTATCACCg         <  1:3704423/36‑1 (MQ=255)
                              ggcgggcgGTAAAACCGGTTTGACGGCTATCACCGt        <  1:975560/36‑1 (MQ=255)
                              ggcgggcgGTAAAACCGGTTTGACGGCTATCACCGt        >  1:7015850/1‑36 (MQ=255)
                               gcgggcggTAAAACCGGTTTGACGGCTATCACCGtt       >  1:2265231/1‑36 (MQ=255)
                                 gggcggTAAAACCGGTTTGACGGCTATCACCGTTgg     >  1:3787865/1‑36 (MQ=255)
                                 gggcggTAAAACCGGTTTGACGGCTATCACCGTTgg     >  1:2145293/1‑36 (MQ=255)
                                 gggcggTAAAACCGGTTTGACGGCTATCACCGTTgg     >  1:5677993/1‑36 (MQ=255)
                                 gggcggTAAAACCGGTTTGACGGCTATCACCGTTgg     <  1:944876/36‑1 (MQ=255)
                                 gggcggTAAAACCGGTTTGACGGCTATCACCGTTgg     >  1:5432914/1‑36 (MQ=255)
                                  ggcggTAAAACCGGTTTGACGGCTATCACCGTTGGc    >  1:2080388/1‑36 (MQ=255)
                                  ggcggTAAAACCGGTTTGACGGATATTAGCGTTggg    >  1:5839869/1‑35 (MQ=21)
                                   gcggTAAAACCGGTTTGACGGCTATCACCGTTGGCg   >  1:4150958/1‑36 (MQ=255)
                                   gcggTAAAACCGGTTTGACGGCTATCACCGTTGGCg   >  1:1776533/1‑36 (MQ=255)
                                   gcggTAAAACCGGTTTGACGGCTATCACCGTTGGCg   >  1:924846/1‑36 (MQ=255)
                                    cggTAAAGCCGGTTTGACGGCTATCACCATTGGCGt  <  1:3494792/36‑1 (MQ=37)
                                    |                                   
GTATACATCGAGTCTGCGGCGGGTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGT  >  REL606/4258366‑4258436

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 15 ≤  ATCG  < 26 ≤  ATCG  < 37 ≤  ATCG