Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 649,391 T→A 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% 167.4 541411mrdAtranspeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  new base (25/29):  ref base (0/0):  total (25/29)

GCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCG  >  REL606/649358‑649423
                                 |                                
gCGGAACCTTCACGATACCGTCATTAATCAGGaaca                                <  1:3535213/36‑1 (MQ=255)
gCGGAACCTTCACGATACCGTCATTAATCAGGaaca                                >  2:3209541/1‑36 (MQ=255)
gCGGAACCTTCACGATACCGTCATTAATCAGGaaca                                >  2:118820/1‑36 (MQ=255)
 cGGAACCTTCACGATACCGTCATTAATCAGGaacat                               >  2:3664592/1‑36 (MQ=255)
  ggAACCTTCACGATACCGTCATTAATCAGGaacatc                              <  1:3672572/36‑1 (MQ=255)
  ggAACCTTCACGATACCGTCATTAATCAGGaacatc                              >  1:1824899/1‑36 (MQ=255)
  ggAACCTTCACGATACCGTCATTAATCAGGaacatc                              >  1:2903414/1‑36 (MQ=255)
   gAACCTTCACGATACCGTCATTAATCAGGAccctca                             >  2:1487787/1‑36 (MQ=25)
   gAACCTTCACGATACCGTCATTAATCAGGAccctca                             >  2:1485082/1‑36 (MQ=25)
   gAACCTTCACGATACCGTCATTAATCAGGAacatca                             >  2:2744146/1‑36 (MQ=255)
   aaaCCTTCACGATACCGTCATTATTCAGGaacataa                             <  1:1226132/35‑3 (MQ=25)
    aaCCTTCACGATACCGTCATTAATCAGGAACATCAg                            <  1:235356/36‑1 (MQ=255)
    aaCCTTCACGATACCGTCACTAATCAGGAacatcaa                            >  1:451453/1‑35 (MQ=25)
    aaCCCTCACGATACCGTCATTAATCAGGAACATCAg                            <  1:1428527/36‑1 (MQ=21)
     aCCTTCACGATACCGTCATTAATCAGGAACATCAGt                           <  1:3518439/36‑1 (MQ=255)
     aCCTTCACGATACCGTCATTAATCAGGAACATCAGt                           <  1:2132416/36‑1 (MQ=255)
      ccTTCACAATACCGTCATTAATCAGGAACATCAGTg                          <  1:2455308/36‑1 (MQ=25)
       cTTCACGATACCGTCATTAATCAGGAACATCAGTGc                         >  2:2985199/1‑36 (MQ=255)
       cTTCACGATACCGTCATTAATCAGGAACATCAGTGc                         >  2:3211024/1‑36 (MQ=255)
       cTTCACGATACCGTCATTAATCAGGAACATCAGTGc                         >  2:3348571/1‑36 (MQ=255)
       cTTCACGATACCGTCATTAATCAGGAACATCAGTGc                         <  2:543819/36‑1 (MQ=255)
        ttCACGATACCGTCATTAATCAGGAACATCAGTGcc                        >  1:900870/1‑36 (MQ=255)
        ttCACGATACCGTCATTAATCAGGAACATCAGTGcc                        <  1:570901/36‑1 (MQ=255)
         tCACGATACCGTCATTAATCAGGAACATCAGTGCCt                       <  2:3639763/36‑1 (MQ=255)
          cACGTTACCGTCATTAATCAGGAACATCAGTGCCtt                      <  1:919479/36‑1 (MQ=37)
          cACGATAGCGTCATTGATCAGGAACATCAGTGCCtt                      <  2:1029599/36‑1 (MQ=25)
          cACGATACCGTCATTAATCAGGAACATCAGTGCCtt                      <  1:138512/36‑1 (MQ=255)
          cACGATACCGTCATTAATCAGGAACATCAGCGCCtt                      >  2:824486/1‑36 (MQ=21)
          cACAATACCGTCATTAATCAGGAACATCAGTGCCtt                      <  1:1762378/36‑1 (MQ=25)
           aCGATACCGTCATTAATCAGGAACATCAGTGCCTTa                     >  1:3425091/1‑36 (MQ=255)
            cGATATCGTCATTAATCAGGAACATCAGTGCCTTAc                    <  2:727904/36‑1 (MQ=37)
            cGATACCGTCATTAGTCAGGAACATCAGTGCCTTAc                    >  1:465473/1‑36 (MQ=37)
            cGATACCGTCATTAATCCGGAACATCAGTGCCttcc                    >  1:647247/1‑34 (MQ=25)
            cGATACCGTCATTAATCAGGAACATCAGTGCCTTAc                    <  1:1849118/36‑1 (MQ=255)
             gATACCGTCATTAATCAGGAACATCAGTGCCTTTCt                   >  1:3661984/1‑36 (MQ=25)
             gATACCGTCATTAATCAGGAACATCAGTGCCTTACt                   <  1:269849/36‑1 (MQ=255)
             gATACCGTCATTAATCAGGAACATCAGTGCCTTACt                   >  2:3061105/1‑36 (MQ=255)
              aTACCGTCATTAATCAGGAACATCAGTGCCTTACTc                  >  2:3337936/1‑36 (MQ=255)
               tACCGTCATTAATCAGGAACATCAGTGCCTTACTCa                 >  2:2731711/1‑36 (MQ=255)
                aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt                >  2:1045735/1‑36 (MQ=255)
                aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt                >  1:1815763/1‑36 (MQ=255)
                aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt                <  2:3277647/36‑1 (MQ=255)
                 ccGTCATTAATCAGGAACATCAGTGCCTTACTCATc               >  2:806664/1‑36 (MQ=255)
                  cGTCATTAATCAGGAACATCAGTGCCTTACTCCTCt              >  2:2438432/1‑36 (MQ=37)
                  cGTCATTAATCAGGAACATCAGTGCCTTACTCATCt              <  1:1932502/36‑1 (MQ=255)
                  cGTCATTAATCAGGAACATCAGTGCCTTACTCATCt              <  2:2350696/36‑1 (MQ=255)
                   gTCATTAATCAGGAACATCAGTGCCTTACTCATCTg             >  1:1271576/1‑36 (MQ=255)
                   gTCATTAATCAGGAACATCAGTGCCTTACTCATCTg             <  1:2555996/36‑1 (MQ=255)
                   gTCATTAATCAGGAACATCAGTGCCTTACTCATCTg             <  1:3559105/36‑1 (MQ=255)
                   atCATAAATCAGGAACATCATTGCCTTACTCATCTg             <  1:3245120/35‑1 (MQ=25)
                     gaTTAATCAGTAACATCAGTGCCTTACTCATCTGGa           <  2:3671499/35‑1 (MQ=37)
                     cATTTATTAGGAACATCGGTGCCTTACTCATCTGGa           <  2:924247/36‑1 (MQ=19)
                     cATTAATCAGGAACATCAGTGCCTTACTCATCTGGa           >  1:2455306/1‑36 (MQ=255)
                      aTTAATCAGGAACATCAGTGCCTTACTCATCTGGAt          <  2:2149787/36‑1 (MQ=255)
                       ttAATCAGGAACATCAGTGCCTTACTCATCTGGAtt         >  2:2036261/1‑36 (MQ=255)
                       taaaTCAGGAACATCAGTGCCTTACTCATCTGCAtt         >  2:1216355/3‑36 (MQ=37)
                         aaTCAGGAACATCAGTGCCTTACTCATCTCGAttag       >  1:1938505/1‑34 (MQ=37)
                         aaTCAGGAACATCAGTGCCTTACTCATCTCGATTgg       >  2:737523/1‑36 (MQ=37)
                          aTCAGGAACATCAGTGCCTTACTCATCTCGCTTggg      >  1:3528447/1‑36 (MQ=25)
                           tCATGAACATCAGTGCCTTACTCATCTGGATTgggg     <  1:2809205/36‑1 (MQ=37)
                           tCAGGAACATCAGTGCCTTACTCATCTGGATTgggg     >  2:1336201/1‑36 (MQ=255)
                           tCAGGAACATCAGTGCCTTACTCATCTGGATTgggg     >  1:279810/1‑36 (MQ=255)
                            cAGGAACATCAGTGCCTTACTCATCTGGATTGGGGt    <  2:2712712/36‑1 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   <  1:1967186/36‑1 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   <  1:3588360/36‑1 (MQ=255)
                              ggAACATCAGTGCCTTACTCATCTGGATTGGGGTCg  <  2:3043241/36‑1 (MQ=255)
                                 |                                
GCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCG  >  REL606/649358‑649423

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 17 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG