Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,258,394 +G 6.5% coding (983/1350 nt) yjcD → predicted permease

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,258,3941.→G 6.5% 4.0 ‑0.7 ‑0.231coding (983/1350 nt)yjcDpredicted permease
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (16/13):  total (18/13)
Fisher's exact test strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.41e-01

GGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTATCACC  >  REL606/4258359‑4258428
                                    |                                  
ggCAGCGGTATACATCGAGTCTGCGGCGGGTAcggc                                     <  1:7211439/36‑1 (MQ=255)
ggCAGCGGTATACATCGAGTCTGCGGCGGGCAcggc                                     >  1:6262362/1‑36 (MQ=255)
 gCAGCGGTATACATCGAGTCTGCGGCGGGTAcggc.g                                   <  1:5881924/36‑1 (MQ=255)
 gCAGCGGTATACATCGAGTCTGCGGCGGGTAcggc.g                                   <  1:2519688/36‑1 (MQ=255)
   aGCGGTATACATCGAGTCTGCGGCGGGTAcggc.ggg                                 >  1:7207378/1‑35 (MQ=255)
    gCGGTATACATCGAGTCTGCGGCGGGTAcggc.gggg                                >  1:1534921/1‑34 (MQ=255)
    gCGGTATACATCGAGTCTGCGGCGGGTAcggc.gggg                                >  1:250693/1‑34 (MQ=255)
      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:7570779/36‑1 (MQ=255)
          taCATCGAGTCTGCGGCGGGTACGGC.GGCGGGCGGt                          <  1:3347204/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:5573814/36‑1 (MQ=255)
                    cTGCGGCGGGTACGGC.GGCGGGCGGTAAAACCGGtt                <  1:5008001/36‑1 (MQ=255)
                    cTGCGGCGGGTACGGC.GGCGGGCGGTAAAACCGGtt                <  1:1225732/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:2096701/36‑1 (MQ=255)
                           gggTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGc         >  1:768597/1‑36 (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)
                                    |                                  
GGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC.GGCGGGCGGTAAAACCGGTTTGACGGCTATCACC  >  REL606/4258359‑4258428

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