Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,620,003 C→G 17.5% A353G (GCG→GGG)  slt → lytic murein transglycosylase, soluble

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,620,0030C→G17.5% 3.7 ‑4.3 ‑0.0281058sltlytic murein transglycosylase, soluble
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (8/15):  total (13/15)
Fisher's exact test strand distribution p-value = 1.31e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.40e-01

GGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTACTGCTGGAACGCGGACGTGAAGCTGAA  >  REL606/4619968‑4620037
                                   |                                  
ggCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGc                                    <  1:889621/36‑1 (MQ=255)
   gaaagagaaagaTGAATGGCGTTACTGGCAGGGggg                                 >  1:5642007/1‑35 (MQ=255)
   gaaagagaaagaTGAATGGCGTTACTGGCAGGGggg                                 >  1:3047800/1‑35 (MQ=255)
    aaagagaaagaTGAATGGCGTTACTGGCAGGGgggg                                >  1:1391529/1‑34 (MQ=37)
     aagagaaagaTGAATGGCGTTACTGGTAGGCGGAtt                               <  1:4434923/36‑1 (MQ=255)
     aagagaaagaTGAATGGCGTTACTGGCAGGGggggt                               >  1:4232906/1‑33 (MQ=37)
       tagaaataTGAATGGCGTTACTGGCAGGCGGATTTa                             <  1:4081602/35‑1 (MQ=255)
        agaaagaTGAATGGCGTTACTGGCAGGCGGATTTAc                            <  1:1895144/36‑1 (MQ=255)
        agaaagaTGAATGGCGTTACTGGCAGGCGGATTTAc                            <  1:4158163/36‑1 (MQ=255)
         gaaagaTGAATGGCGTTACTGGCAGGCGGATTTAct                           <  1:622074/36‑1 (MQ=255)
         gaaagaTGAATGGCGTTACTGGCAGGCGGATTTAct                           <  1:1159693/36‑1 (MQ=255)
            agaTGAATGGCGTTACTGGCAGGCGGATTTActgct                        <  1:3691193/36‑1 (MQ=255)
               tGAATGGCGTTACTGGCAGGCGGATTTACTGCTGGa                     <  1:3954523/36‑1 (MQ=255)
               tGAATGGCGTTACTGGCAGGCGGATTTACTGCTGGa                     <  1:4560970/36‑1 (MQ=255)
                   tGGCGTTACTGGCAGGCGGATTTACTGCTGGAAcgc                 <  1:4637804/36‑1 (MQ=255)
                     gCGTTACTGGCAGGGGGGGTTATTGGTGGGAAGCgg               >  1:3432925/1‑36 (MQ=0)
                     gCGTTACTGGCAGGCGGGTTTTTTGCTTGGAAGCgg               >  1:637551/1‑36 (MQ=9)
                       gTTACTGGCAGGCGGGTTTTTTGCTTGGACGCgggc             >  1:1197602/1‑34 (MQ=12)
                            tGGCAGGCGGATTTACTGCTGGAACGCGGACGtgaa        <  1:4614899/36‑1 (MQ=255)
                             ggCAGGCGGATTTACTGCTGGAACGCGGACGtgaag       >  1:1702836/1‑36 (MQ=255)
                              gCAGGCGGATTTACTGCTGGAACGCGGACGtgaagc      <  1:5454994/36‑1 (MQ=255)
                              gCAGGCGGATTTACTGCTGGAAAGCGGACGtgaagc      >  1:5706767/1‑36 (MQ=255)
                               cAGGCGGATTTACTGCTGGAACGCGGACGtgaagct     <  1:983523/36‑1 (MQ=255)
                               cAGGCGGATTTACTGCTGGAACGCGGACGtgaagct     >  1:540083/1‑36 (MQ=255)
                               cAGGCGGATTTACTGCTGGAACGCGGACGtgaagct     >  1:3150462/1‑36 (MQ=255)
                                aGGCGGATTTACTGCTGGAACGCGGACGtgaagctg    >  1:5458380/1‑36 (MQ=255)
                                 ggCGGATTTCCTGCTGGAACGCGGACGtgaagctga   <  1:1836790/36‑1 (MQ=255)
                                  gCGGATTTACTTCGGGAACGCGGACGtgaagctgaa  <  1:5083883/36‑1 (MQ=38)
                                  gCGGATTTACTGCTGGAACGCGGACGtgaagctgaa  >  1:728387/1‑36 (MQ=255)
                                   |                                  
GGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTACTGCTGGAACGCGGACGTGAAGCTGAA  >  REL606/4619968‑4620037

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