Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,453,694 +G 6.1% coding (2142/2553 nt) mrcA → fused penicillin‑binding protein 1a: murein transglycosylase/murein transpeptidase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,453,6941.→G 6.1% 3.7 ‑1.6 ‑0.733coding (2142/2553 nt)mrcAfused penicillin‑binding protein 1a: murein transglycosylase/murein transpeptidase
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (13/18):  total (15/18)
Fisher's exact test strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.09e-01

CGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC.GGGAAAACCGGTACCACTAACAGTTCGAAAGATG  >  REL606/3453659‑3453728
                                    |                                  
cGTGCAGGTCGTGATTTGCGGCGTCGCGATATcggc                                     >  1:4597351/1‑36 (MQ=255)
cGTGCAGGTCGTGATTTGCAGCGTCGCGATATcggc                                     >  1:4029498/1‑36 (MQ=255)
cGTGCAGGTCGTGATTTGCAGCGTCGCGATATcggc                                     >  1:3773878/1‑36 (MQ=255)
 gTGCAGGTCGTGATTTGCAGCGTCGCGATATcggc.g                                   >  1:6843295/1‑36 (MQ=255)
 gTGCAGGTCGTGATTTGCAGCGTCGCGATATcggc.g                                   <  1:2502913/36‑1 (MQ=255)
 gTGCAGGTCGTGATTTGCAGCGTCGCGATATcggc.g                                   >  1:2123351/1‑36 (MQ=255)
 gTGCAGGTCGTGATTTGCAGCGTCGCGATATcggc.g                                   <  1:1464704/36‑1 (MQ=255)
 gTGCAGGTCGTGATTTGCAGCGTCGCGATATcggc.g                                   >  1:1544550/1‑36 (MQ=255)
  tGCAGGTCGTGATTTGCAGCGTCGCGATATcggc.gg                                  >  1:2751056/1‑36 (MQ=255)
   gCAGGTCGTGATTTGCAGCGTCGCGATATCGGC.ggg                                 >  1:1611453/1‑36 (MQ=255)
      ggTCGTGATTTGCAGCGTCGCGATATCGGC.ggggga                              >  1:6504121/1‑33 (MQ=255)
       gTCTTGATTTGCACCGCCGCGATATCGGC.GGGaaaa                             <  1:5680035/36‑1 (MQ=25)
       gTCGTGATTTGCAGCGTCGCGATATCGGC.GGGgaaa                             >  1:2455643/1‑36 (MQ=255)
        tCGTGATTTGCAGCGTCGCGATATCGGC.GGGGAAAc                            >  1:8140820/1‑36 (MQ=255)
        tCGTGATTTGCAGCGTCGCGATATCGGC.GGGAAAAc                            <  1:8336320/36‑1 (MQ=255)
         cgggATTGGCAGTGTCGCGATATCTGC.GGGAAAAcc                           <  1:6329713/33‑1 (MQ=19)
         cGTGATTTGCATCGTCTCGATATCGGC.GGGAAAAcc                           <  1:3309908/36‑1 (MQ=37)
         cGTGATTTGCAGCGTCGCGATATCGGC.GGGGAAAcc                           >  1:4383416/1‑36 (MQ=255)
          gTGATTTGCAGCGTCGCGATATCGGCGGGGAAAAcc                           >  1:5549352/1‑36 (MQ=37)
          gTGATTTGCAGCGTCGCGATATCGGC.GGGAAAACCg                          <  1:3265508/36‑1 (MQ=255)
           tGATTTGCAGCGTCGCGATATCGGCGGGGAAAACCg                          >  1:809379/1‑36 (MQ=37)
           tGATTTGCAGCGTCGCGATATCGGC.GGGGAAACCgg                         >  1:3081895/1‑36 (MQ=255)
             aTTTGCAGCGTCGCGATATCGGC.GGGGAAACCCGTa                       >  1:6559921/1‑36 (MQ=37)
             aTTTGCAGCGTCGCGATATCGGC.GGGAAAACCGGTa                       <  1:4723115/36‑1 (MQ=255)
                 gCAGCGTCGCGATATCGGC.GGGAAAACCGGTACCAc                   <  1:2649419/36‑1 (MQ=255)
                  cAGCGTCGCGATATCGGC.GGGAAAACCGGTACCACt                  <  1:5227/36‑1 (MQ=255)
                  cAGCGTCGCGATATCGGC.GGGAAAACCGGTACCACt                  <  1:5018652/36‑1 (MQ=255)
                    gCGTCGCGATATCGGC.GGGGAAACCGGTACCACTaa                >  1:8443617/1‑36 (MQ=255)
                     cGTCGCGATATCGGC.GGGAAAACCGGTACCACTAAc               >  1:718048/1‑36 (MQ=255)
                     cGTCGCGATATCGGC.GGGAAAACCGGTACCACTAAc               <  1:7515203/36‑1 (MQ=255)
                     cGTCGCGATATCGGC.GGGAAAACCGGGACCCCTCAc               >  1:3600675/1‑36 (MQ=25)
                      gTCGCGATATCGGC.GGGGAAACCGGTACCACTCACa              >  1:6927227/1‑36 (MQ=37)
                      gTCGCGATATCGGC.GGGAAAACCGGTACCACTAACa              <  1:5770365/36‑1 (MQ=255)
                      gTCGCGATATCGGC.GGGAAAACCGGTACCACTAACa              <  1:8133561/36‑1 (MQ=255)
                      gTCGCGATATCGGC.GGGAAAAACGGTACCACTAAcc              >  1:8195788/1‑35 (MQ=255)
                       tCGCGATATCGGC.GGGAAAGCCGGTACCACTAACAg             <  1:2325150/36‑1 (MQ=255)
                       tCGCGATATCGGC.GGGAAAACCGGTACCACTAACAg             <  1:1470525/36‑1 (MQ=255)
                          cgATATCGGC.GGGAAAACCGGTACCACTAACAGTTc          <  1:4826498/36‑1 (MQ=255)
                          cgATATCGGC.GGGAAAACCGGTACCACTAACAGTTc          <  1:4634778/36‑1 (MQ=255)
                           gATATCGGC.GGGAAAACCGGTACCACTAACAGTTCg         <  1:561442/36‑1 (MQ=255)
                            atatCGGC.GGGAAAAGCGGTACCACTAACAGTTCGa        <  1:4389707/36‑1 (MQ=255)
                            atatCGGC.GGGAAAACCGGGACCACTAACAGTTCgt        >  1:1158358/1‑35 (MQ=255)
                             tagcggc.ggGAAAACCGGTACCACTAACAGTTCGaa       <  1:7719420/33‑1 (MQ=255)
                               tCGGC.GGGAAAACCGGTACCACTAACAGTTCGAAAg     >  1:3680526/1‑36 (MQ=255)
                               tCGGC.GGGAAAACCGGTACCACTAACAGTTCGAAAg     >  1:957362/1‑36 (MQ=255)
                                cggc.ggGAAAACCGGTACCACTAACAGTTCGAAAGa    <  1:7289805/36‑1 (MQ=255)
                                 ggc.ggGAAAACCGGTACCACTAACAGTTCGAAAGAt   >  1:543395/1‑36 (MQ=255)
                                 ggc.ggGAAAACCGGTACCACTAACAGTTCGAAAGAt   <  1:4065396/36‑1 (MQ=255)
                                 ggc.ggGAAAACCGGTACCACTAACAGTTCGAAAGAt   <  1:1756027/36‑1 (MQ=255)
                                 ggc.ggGAAAACCGGTACCACTAACAGTTCGAAAGAt   >  1:8441554/1‑36 (MQ=255)
                                 ggc.ggGAAAACCGGTACCACCAACAGTTCGAAAGAt   >  1:4312749/1‑36 (MQ=255)
                                  gc.ggGAAAACCGGTACCACTAACAGTTCGAAAGATg  <  1:6772383/36‑1 (MQ=255)
                                  gc.ggGAAAACCGGTACCACTAACAGTTCGAAAGATg  >  1:4938494/1‑36 (MQ=255)
                                  gc.ggGAAAACCGGTACCACTAACAGTTCGAAAGATg  >  1:7345659/1‑36 (MQ=255)
                                  gc.ggGAAAACCGGTACCACTAACAGTTCGAAAGATg  <  1:3956419/36‑1 (MQ=255)
                                  gc.ggGAAAACCGGTACCACTAACAGTTCGAAAGATg  <  1:1800639/36‑1 (MQ=255)
                                  gc.ggGAAAACCGGTACCACTAACAGTTCGAAAGATg  >  1:867785/1‑36 (MQ=255)
                                    |                                  
CGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC.GGGAAAACCGGTACCACTAACAGTTCGAAAGATG  >  REL606/3453659‑3453728

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG