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% 199.9 561411mrdAtranspeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  new base (30/24):  ref base (0/2):  total (30/26)

GCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTC  >  REL606/649358‑649422
                                 |                               
gCGGAACATTCACGATACCGTCATTAATCTGGaaca                               >  1:5554285/1‑36 (MQ=21)
 cGGAACCTTCACGATACCGTCATTAATCAGGaacat                              >  1:3607721/1‑36 (MQ=255)
 cGGAACCTTCACGATACCGTCATTAATCAGGaacat                              >  1:871112/1‑36 (MQ=255)
 cGGAACCTTCACGATACCGTCATTAATCAGGaacat                              >  1:2666408/1‑36 (MQ=255)
 cGGAACCTTCACGATACCGTCATTAATCAGGaacat                              >  1:2602106/1‑36 (MQ=255)
 cGGAACCTTCACGATACCGTCATTAATCAGGaacat                              <  1:2551737/36‑1 (MQ=255)
 cGGAACCTTCACGATACCGTCATTAATCAGGaacat                              >  1:5686062/1‑36 (MQ=255)
 cGGAACCTTCACGATAACGTCATTAATAAGGaacat                              >  1:1499694/1‑36 (MQ=25)
  ggAACCTTCACGATACCGTCATTAATCAGGaacatc                             >  1:5318011/1‑36 (MQ=255)
  ggAACCTTCACAATACCGTCATTAATCAGGaacatc                             <  1:3741282/36‑1 (MQ=25)
   gAACCTTCACGATACCGTCATTAATCAGGAacatca                            <  1:7630322/36‑1 (MQ=255)
   gAACCTTCACGATACCGTCATTAATCAGGAacatca                            >  1:7040610/1‑36 (MQ=255)
    aaCCTTCACGATACCGTCATTAATCAGGAACATCAg                           >  1:4446618/1‑36 (MQ=255)
     aCCTTCACGATACCGTCATTAATCAGGATCATCAGt                          <  1:1356618/36‑1 (MQ=255)
     aCCTTCACGATACCGTCATTAATCAGGAACATCAGt                          <  1:5117287/36‑1 (MQ=255)
      ccTTCACGATACCGTCATTAATCAGGAACATCAGTg                         >  1:7402267/1‑36 (MQ=255)
       cTTCACGATACCGTCATTAATCAGGAACATCAGTGc                        <  1:107125/36‑1 (MQ=255)
       cTTCACGATACCGTCATTAATCAGGAACATCAGTGc                        >  1:4924038/1‑36 (MQ=255)
         tCACGATACCGTCATTAATCAGGAACATCAGTGCCt                      <  1:3706743/36‑1 (MQ=255)
          cACGATACCGTCATTAATCAGGAACATCAGTGCCtt                     >  1:3677047/1‑36 (MQ=255)
          cACGATACCGTCATTAATCAGGAACATCAGTGCCtt                     <  1:3551801/36‑1 (MQ=255)
           aCGATACCGTCATTAATCAGGAACATCAGTGCCTTa                    >  1:2311946/1‑36 (MQ=255)
           aCGATACCGTCATTAATCAGGAACATCAGTGCCTTa                    >  1:5691313/1‑36 (MQ=255)
           aCGATACCGTCATTAATCAGGAACATCAGTGCCTTa                    >  1:958533/1‑36 (MQ=255)
             gATACCGTCATTAATCAGGAACATCAGTGCCTTACt                  <  1:1277752/36‑1 (MQ=255)
             gATACCGTCATTAATCAGGAACATCAGTGCCTTACt                  <  1:386366/36‑1 (MQ=255)
              aTACCGTCATTAATCAGGAACATCAGTGCCTTACTc                 >  1:2839237/1‑36 (MQ=255)
              aTACCGTCATTAATCAGGAACATCAGTGCCTTACTc                 >  1:4032182/1‑36 (MQ=255)
              aTACCGTCATTAATCAGGAACATCAGTGCCTTACTc                 >  1:7346302/1‑36 (MQ=255)
              aTACCGTCATTAATCAGGAACATCAGTGCCTTACTc                 >  1:6028676/1‑36 (MQ=255)
              aTACCGTCATTAATCAGGAACATCAGTGCCTTACTc                 <  1:1216347/36‑1 (MQ=255)
               tACCGTCATTAATCAGGAACATCAGTGCCTTACTCa                >  1:1793144/1‑36 (MQ=255)
               tACCGTCATTAATCAGGAACATCAGTGCCTTACTCa                <  1:195497/36‑1 (MQ=255)
                aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt               <  1:449381/36‑1 (MQ=255)
                aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt               <  1:4986404/36‑1 (MQ=255)
                aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt               >  1:434187/1‑36 (MQ=255)
                aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt               >  1:1278937/1‑36 (MQ=255)
                aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt               >  1:23077/1‑36 (MQ=255)
                  cGTCATTAATCAGGAACATCAGTGCCTTACTCATCt             <  1:1939498/36‑1 (MQ=255)
                  cGTCATTAATCAGGAACATCAGTGCCTTACTCATCt             >  1:2982571/1‑36 (MQ=255)
                    tCATTAATCAGGAACATCAGTGCCTTACTCATCTgg           <  1:7156848/36‑1 (MQ=255)
                    tCATTAATCAGGAACATCAGTGCCTTACTCATCTgg           >  1:7156007/1‑36 (MQ=255)
                     cATTAATCAGGAACATCAGTGCCTTACTCATCTGGa          <  1:5521728/36‑1 (MQ=255)
                     cATTAATCAGGAACATCAGTGCCTTACTCATCTGGa          >  1:878218/1‑36 (MQ=255)
                     cATTAATCAGGAACATCAGTGCCTTACTCATCTGGa          <  1:5750652/36‑1 (MQ=255)
                     cATTAATCAGGAACATCAGTGCCTTACTCATCTGGa          <  1:4643773/36‑1 (MQ=255)
                      aTTAATCAGGAACATCAGTGCCTTACTCATCTGGAt         >  1:3419600/1‑36 (MQ=255)
                      aTTAATCAGGAACATCAGTGCCTTACTCATCTGGAt         <  1:589336/36‑1 (MQ=255)
                      aTTAATCAGGAACATCAGTGCCTTACTCATCTGGAt         <  1:7644992/36‑1 (MQ=255)
                       ttAATCAGGAACATCAGTGCCTTACTCATCTGGAtt        <  1:7885801/36‑1 (MQ=255)
                        tAATCAGGAACATCAGTGCCTTACTCATCTGGATTg       <  1:7233129/36‑1 (MQ=255)
                         aaTCAGGAACATCAGTGCCTTACTCATCTGGATTgg      >  1:3168020/1‑36 (MQ=255)
                         aaTCAGGAACATCAGTGCCTTACTCATCTGGATTgg      >  1:7819684/1‑36 (MQ=255)
                         aaTCAGGAACATCAGGGCCTTACTCATCTGGATTgg      >  1:7795278/1‑36 (MQ=37)
                          aTCAGGATCATCAGTGCCTTACTCATCTGGATTggg     <  1:1873789/36‑1 (MQ=255)
                          aTCAGGAACATCAGTGCCTTACTCATCTGGATTggg     <  1:1018422/36‑1 (MQ=255)
                          aTCAGGAACATCAGTGCCTTACTCATCTGGATTggg     <  1:1719117/36‑1 (MQ=255)
                          aTCAGGAACATCAGTGCCTTACTCATCTGGATTggg     >  1:3316005/1‑36 (MQ=255)
                          aTCAGGAACATCAGTGCCTTACTCATCTGGATTggg     <  1:1959008/36‑1 (MQ=255)
                          aTCAGGAACATCAGTGCCTAACTCATCGGGATTggg     >  1:6459599/1‑36 (MQ=25)
                           tCAGGAACATCAGTGCCTTACTCATCTGGATTgggg    <  1:2817032/36‑1 (MQ=255)
                            cAGGAACATCAGTGCCTTACTCATCTGGATTGGGGt   >  1:7236863/1‑36 (MQ=255)
                            cAGGAACATCAGTGCCTTACTCATCTGGATTGGGGt   >  1:192600/1‑36 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc  >  1:7257835/1‑36 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc  <  1:162547/36‑1 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc  >  1:7731238/1‑36 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc  >  1:2858128/1‑36 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATGGGGGTc  >  1:7888313/1‑36 (MQ=37)
                                 |                               
GCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTC  >  REL606/649358‑649422

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG