Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 649,391 T→A 100% I471F (ATC→TTC)  mrdA ← transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606649,3910T→A100.0% 232.9 641411mrdAtranspeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  new base (31/31):  ref base (1/1):  total (32/32)

GCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCG  >  REL606/649358‑649423
                                 |                                
gCGGAACCTTCACGATACCGTCATTAATCAGGaaca                                >  1:621169/1‑36 (MQ=255)
gCGGAACCTTCACGATACCGTCATTAATCAGGaaca                                <  1:818597/36‑1 (MQ=255)
  ggAACCTTCACGATACCGTCATTAATCAGGaacatc                              >  1:4220679/1‑36 (MQ=255)
  ggAACCTTCACGATACCGTCATTAATCAGGaacatc                              >  1:2204278/1‑36 (MQ=255)
   gAACCTTCACGATACCGTCATTAATCACGAacatca                             >  1:1981267/1‑36 (MQ=37)
   aaaCCCTCACGATACCGTCATTAATCAGGAacatca                             <  1:5390291/35‑1 (MQ=21)
     aCCTTCACGATACCGTCATTAATCAGGAACATCAGt                           >  1:7487740/1‑36 (MQ=255)
      ccTTCACGATACCGTCATTAATCAGGAACATCAGTg                          <  1:3354115/36‑1 (MQ=255)
       cTTCACGATACCGTCATTAATCAGGAACATCAGTGc                         >  1:1128247/1‑36 (MQ=255)
       cTTCACGATACCGTCATTAATCAGGAACATCAGTGc                         >  1:4116473/1‑36 (MQ=255)
        ttCACGATACCGTCATTAATCAGGAACATCAGTGcc                        >  1:4540818/1‑36 (MQ=255)
         tCACGATACCGTCATTAATCAGGATCATCAGTGCCt                       >  1:5832384/1‑36 (MQ=255)
         tCACGATACCGTCATTAATCAGGAACATCAGTGCCt                       >  1:3656465/1‑36 (MQ=255)
         tCACGATACCGTCATTAATCAGGAACATCAGTGCCt                       >  1:8385368/1‑36 (MQ=255)
         tCACGATACCGTCATTAATCAGGAACATCAGTGCCt                       >  1:2065652/1‑36 (MQ=255)
          cACGATACCGTCATTAATCAGGAACATCAGTGCCtt                      <  1:5690051/36‑1 (MQ=255)
           aCGATACCGTCATTAATCAGGAACATCAGTGCCTTa                     <  1:2395260/36‑1 (MQ=255)
           aCGATACCGTCATTAATCAGGAACATCAGTGCCTTa                     >  1:7328012/1‑36 (MQ=255)
           aCGATACCGTCATTAATCAGGAACATCAGTGCCTTa                     <  1:7027942/36‑1 (MQ=255)
            cGATACCGTCATTAATCAGGAACATCAGTGCCTTAc                    <  1:2889699/36‑1 (MQ=255)
             gATACCGTCATTAATCAGGAACATCAGTGCCTTACt                   <  1:6562154/36‑1 (MQ=255)
             gATACCGTCATTAATCAGGAACATCAGTGCCTTACt                   >  1:7674984/1‑36 (MQ=255)
             gATACCGTCATTAATCAGGAACATCAGTGCCTTACt                   <  1:7665335/36‑1 (MQ=255)
               tACCGTCATTAATCAGGAACATCAGTGCCTTACTCa                 >  1:5665847/1‑36 (MQ=255)
               tACCGTCATTAATCAGGAACATCAGTGCCTTACTCa                 <  1:4511800/36‑1 (MQ=255)
                aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt                >  1:6032615/1‑36 (MQ=255)
                aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt                <  1:752264/36‑1 (MQ=255)
                aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt                >  1:5370255/1‑36 (MQ=255)
                aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt                <  1:5185365/36‑1 (MQ=255)
                aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt                <  1:3058475/36‑1 (MQ=255)
                aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt                >  1:3296672/1‑36 (MQ=255)
                 ccGTCATTAATCAGGAACATCAGTGCCTTACTCATc               <  1:7422279/36‑1 (MQ=255)
                 ccGTCATTAATCAGGAACATCAGTGCCTTACTCATc               <  1:6395862/36‑1 (MQ=255)
                 ccGTCATTAATCAGGAACATCAGTGCCTTACTCATc               <  1:2115405/36‑1 (MQ=255)
                  cGTCATTAATCAGTAACATCAGTGCCTTACTCATCt              <  1:7078640/36‑1 (MQ=37)
                  cGTCATTAATCAGGATCATCAGTGCCTTACTCATCt              <  1:2610138/36‑1 (MQ=255)
                  cGTCATTAATCAGGAACATCAGTGCCTTACTCATCt              >  1:5156534/1‑36 (MQ=255)
                  cGTCATTAATCAGGAACATCAGTGCCTTACTCATCt              >  1:8213977/1‑36 (MQ=255)
                  cGTCATTAATCAGGAACATCAGTGCCTTACTCATCt              <  1:2360434/36‑1 (MQ=255)
                   gTCATTAATCAGGAACATCAGTGCCTTACTCATCTg             >  1:4740249/1‑36 (MQ=255)
                   gTCATTAATCAGGAACATCAGTGCCTTACTCATCTg             >  1:764431/1‑36 (MQ=255)
                   gTCATTAATCAGGAACATCAGTGCCTTACTCATCTg             <  1:642490/36‑1 (MQ=255)
                   gTCATTAATCAGGAACATCAGTGCCTTACTCATCTg             >  1:1766427/1‑36 (MQ=255)
                   gTCATTAATCAGGAACATCAGTGCCTTACTCATCTg             >  1:1955955/1‑36 (MQ=255)
                    tCATTAATCAGGAACATCAGTGCCTTACTCATCTgg            >  1:1456268/1‑36 (MQ=255)
                     cATTAATCAGGAACATCAGTGCCTTACTCATCTGGa           <  1:2490340/36‑1 (MQ=255)
                      aTTAATCAGGAACATCAGTGCCTTACTCATCTGGAt          <  1:3455269/36‑1 (MQ=255)
                      aTTAATCAGGAACATCAGTGCCTTACTCATCTGGAt          <  1:587467/36‑1 (MQ=255)
                       ttAATCAGGAACATCAGTGCCTTACTCATCTGGAtt         >  1:3503929/1‑36 (MQ=255)
                       ttAATCAGGAACATCAGTGCCTTACTCATCTGGAtt         >  1:1799303/1‑36 (MQ=255)
                        tAATCAGGAACATCAGTGCCTTACTCATCTGGATTg        >  1:2878239/1‑36 (MQ=255)
                        tAATCAGGAACATCAGCGACTTACTCATATGGATTg        >  1:4376288/1‑36 (MQ=16)
                         aaTCAGGAACATCAGTGCCTTACTCATCTGGATTgg       >  1:5554222/1‑36 (MQ=255)
                         aaTCAGGAACATCAGTGCCTTACTCATCTGGATTgg       <  1:3919055/36‑1 (MQ=255)
                         aaTCAGGAACATCAGTGACTTACTCATCTGGATTgg       <  1:7389374/36‑1 (MQ=37)
                          aTCAGGAACATCAGTGCCTTACTCATCTGGATTggg      <  1:7888680/36‑1 (MQ=255)
                          aTCAGGAACATCAGTGCCTTACTCATCTGGATTggg      >  1:5265823/1‑36 (MQ=255)
                           tCAGGAACATCAGTGCCTTACTCATCTGGATTgggg     <  1:5406268/36‑1 (MQ=255)
                           tCAGGAACATCAGTGCCTTACTCATCTGGATTgggg     <  1:3588800/36‑1 (MQ=255)
                            cAGGAACATCAGTGCCTTACTCATCTGGATTGGGGt    >  1:5769799/1‑36 (MQ=255)
                            cAGGAACATCAGTGCCTTACTCATCTGGATTGGGGt    <  1:6899295/36‑1 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   >  1:7071110/1‑36 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   <  1:4742368/36‑1 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   <  1:3237676/36‑1 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   <  1:2879288/36‑1 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   <  1:1782939/36‑1 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   <  1:1104985/36‑1 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   >  1:1072461/1‑36 (MQ=255)
                              ggAACATCAGTGCCTTACTCATCTGGATTGGGGTCg  >  1:5474903/1‑36 (MQ=255)
                              ggAACATCAGTGCCTTACTCATCTGGATTGGGGTCg  >  1:3566818/1‑36 (MQ=255)
                                 |                                
GCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCG  >  REL606/649358‑649423

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG