Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,258,397 C→G 18.3% A329G (GCG→GGG)  yjcD → predicted permease

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,258,3970C→G18.3% 2.1 ‑3.7 0.030986yjcDpredicted permease
Reads supporting (aligned to +/- strand):  new base (6/0):  ref base (11/13):  total (17/13)
Fisher's exact test strand distribution p-value = 2.37e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GGTATACATCGAGTCTGCGGCGGGTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTG  >  REL606/4258365‑4258432
                                 |                                   
ggTATACATCGAGTCTGCGGCGGGGACGGC.gggggg                                  >  1:847094/1‑36 (MQ=38)
 gTATACATCGAGTCTGCGGCGGGTACGGC.ggggggg                                 >  1:6733159/1‑35 (MQ=255)
   ataCATCGAGTCTGCGGCGGGCACGGC.Ggggggggg                               >  1:1910861/1‑36 (MQ=25)
    taCATCGAGTCTGCGGCGGGTACGGC.GGGGGGGGGt                              >  1:5866390/1‑36 (MQ=38)
    taCATCGAGTCTGCGGCGGGTACGGC.GGCGGGCGGt                              <  1:3347204/36‑1 (MQ=255)
    taCATCGAGTCTGCGGCGGGTACGGC.GGCGGGCGGt                              <  1:7570779/36‑1 (MQ=255)
       aTCGAGTCTGCGGCGGGTACGGC.GGGGGGGGGTgaa                           >  1:3012536/1‑36 (MQ=37)
         cGAGTCTGCGGCGGGTACGGC.GGCGGGCGGTAAAAc                         <  1:2567719/36‑1 (MQ=255)
          gAGTCTGCGGCGGGTACGGCGGGGGGGGGGTAAAAc                         >  1:4441891/1‑36 (MQ=19)
          gAGTCTGCGGCGGGTACGGC.GGGGGGGGGTtaaaac                        >  1:1165789/1‑34 (MQ=25)
          gAGTCTGCGGCGGGTACGGC.GGCGGGGGGGAAAAcc                        >  1:31872/1‑36 (MQ=37)
           aGTCTGCGGCGGGTACGGCGGGGGGGGGGTAAAAcc                        >  1:454124/1‑36 (MQ=19)
            gTCTGCGGCGGGTACGGC.GGCGGGGGGGAAAAACgg                      >  1:3375419/1‑36 (MQ=25)
             tCTGCGGCTGGTACGGC.TGCGGGCGGTAAAACCGGt                     <  1:234476/36‑1 (MQ=37)
             tCTGCGGCGGGTACGGC.GGCGGGCGGTAAAACCGGt                     <  1:2550432/36‑1 (MQ=255)
              cTGCGGCGGGTACGGC.GGCGGGGGGTGAAAACgggt                    >  1:6046438/1‑34 (MQ=21)
              cTGCGGCGGGTACGGC.GGCGGGCGGTAAAACCGGtt                    <  1:1225732/36‑1 (MQ=255)
              cTGCGGCGGGTACGGC.GGCGGGCGGTAAAACCGGtt                    <  1:5008001/36‑1 (MQ=255)
              cTGCGGCGGGTACGGC.GGCGGGCGGTAAAACCGGtt                    <  1:5573814/36‑1 (MQ=255)
               tGCGGCGGGTACGGC.GGCGGGCGGTAAAACCGGttt                   <  1:2710373/36‑1 (MQ=255)
                  ggcggGTACGGC.GGGGGGGGGTAAAACCGGGTTggg                >  1:4348592/1‑34 (MQ=21)
                  ggcggGTACGGC.GGCGGGCGGTAAAACCGGTTTGAc                <  1:6421999/36‑1 (MQ=255)
                     gggTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGc             <  1:4623459/36‑1 (MQ=255)
                     gggTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGc             >  1:1232430/1‑36 (MQ=255)
                     gggTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGc             >  1:768597/1‑36 (MQ=255)
                     gggTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGc             <  1:2096701/36‑1 (MQ=255)
                     gggTACGGC.GGCGGGCGGTAAAACCGGTTTGAAGGc             >  1:1154312/1‑36 (MQ=255)
                      ggTACGGC.GGCGGGGGGTAAAACCGGTTTGACgggt            >  1:3689972/1‑34 (MQ=37)
                      ggTACGGC.GGCGGGCGGTAAAACCGGTTTGACgggt            >  1:5354824/1‑34 (MQ=255)
                      ggTACGGC.GGCGGGCGGTAAAACCGGGTTGACgggt            >  1:4093776/1‑34 (MQ=25)
                       gTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTa           <  1:4308670/36‑1 (MQ=255)
                         aCGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTATc         >  1:4167857/1‑36 (MQ=255)
                         aCGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTATc         >  1:1006850/1‑36 (MQ=255)
                          cggc.ggcggGCGGTAAAACCGGTTTGACGGCTAGCa        >  1:4808740/1‑36 (MQ=255)
                           ggc.ggcggGCGGTAAAACCGGTTTTACGGCTAttgc       >  1:2340452/1‑33 (MQ=37)
                            gc.ggcggGCGGTAAAACCGGTTTGACGGCTATCAcc      >  1:6168679/1‑36 (MQ=255)
                               ggcgggcgGTAAAACCGGTTTGACGGCTATCACCGt    >  1:3937159/1‑36 (MQ=255)
                               ggcgggcgGTAAAACCGGTTTGACGGCTATCACCGt    >  1:329349/1‑36 (MQ=255)
                                gcgggcggTAAAACCGGTTTGACGGCTATCACCGtt   >  1:6564767/1‑36 (MQ=255)
                                gcgggcggTAAAACCGGTTTGACGGCTATCACCGtt   <  1:3511313/36‑1 (MQ=255)
                                gcgggcggTAAAACCGGTTTGACGGCTATCACCGtt   <  1:6803787/36‑1 (MQ=255)
                                 cgggcggTAAAACCGGTTTGACGGCTATCACCGTTg  <  1:4810757/36‑1 (MQ=255)
                                 |                                   
GGTATACATCGAGTCTGCGGCGGGTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTG  >  REL606/4258365‑4258432

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 16 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG