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

CGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCG  >  REL606/649359‑649423
                                |                                
cGGAACCTTCACGATACCGTCATTAATCAGGaacat                               <  1:4662770/36‑1 (MQ=255)
   aaCCTTCACGATACCGTCATTAATCAGGAACATCAg                            >  1:498156/1‑36 (MQ=255)
    aCCTTCACGATACCGTCATTAATCAGGAACATCAGt                           >  1:1698975/1‑36 (MQ=255)
    aCCTTCACGATACCGTCATTAATCAGGAACATCAGt                           >  1:65222/1‑36 (MQ=255)
    aCCTTCACGATACCGTCATTAATCAGGAACATCAGt                           >  1:1846324/1‑36 (MQ=255)
       ttCACGATACCGTCATTAATCAGGAACATCAGTGcc                        >  1:1220144/1‑36 (MQ=255)
        tCACGATACCGTCATTAATCAGGAACATCAGTGCCt                       >  1:3049896/1‑36 (MQ=255)
         cACGAAACCTTCATTAATCAGTAACATCAGTGCCtt                      <  1:5443940/36‑1 (MQ=16)
          aCGATACCGTCATTAATCAGGAACATCAGTGCCTTa                     <  1:1107068/36‑1 (MQ=255)
          aCGATACCGTCATTAATCAGGAACATCAGTGCCTTa                     <  1:1899454/36‑1 (MQ=255)
           cGATACCGTCATTAATCAGGAACATCAGTGCCTTAc                    <  1:1260087/36‑1 (MQ=255)
           cGATAACGTCATTAATCAGGAACATCAGTGCCTTAc                    <  1:4044834/36‑1 (MQ=37)
            gATACCGTCATTAATCCGGAACATCAGTGCCTTCCt                   <  1:5559154/36‑1 (MQ=25)
            gATACCGTCATTAATCAGGAACATCAGTGCCTTACt                   >  1:3165427/1‑36 (MQ=255)
            gATACCGTCATTAATCAGGAACATCAGTGCCTTACt                   <  1:7520/36‑1 (MQ=255)
             aTACCGTCATTAATCAGGAACATCAGTGCCTTACTc                  >  1:1947712/1‑36 (MQ=255)
             aTACCGTCATTAATCAGGAACATCAGTGCCTTACTc                  >  1:1214550/1‑36 (MQ=255)
               aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt                >  1:2200685/1‑36 (MQ=255)
               aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt                <  1:5656656/36‑1 (MQ=255)
                ccGTCATTAATCAGGAACATCAGTGCCTTACTCATc               <  1:2447378/36‑1 (MQ=255)
                ccGTCATTAATCAGGAACATCAGTGCCTTACTCATc               <  1:5108638/36‑1 (MQ=255)
                 cGTCATTAATCAGGAACATCAGTGCCTTACTCATCt              <  1:4786816/36‑1 (MQ=255)
                 cGTCATTAATCAGGAACATCAGTGCCTTACTCATCt              >  1:4394116/1‑36 (MQ=255)
                 cGTCATTAATCAGGAACATCAGTGCCTTACTCATCt              <  1:3436077/36‑1 (MQ=255)
                 cGTCATTAATCAGGAACATCAGTGCCTTACTCATCt              <  1:1038759/36‑1 (MQ=255)
                  gTCATTAATCAGGAACATCAGTGCCTTACTCATCTg             >  1:5629629/1‑36 (MQ=255)
                  gTCATTAATCAGGAACATCAGTGCCTTACTCATCTg             >  1:4095524/1‑36 (MQ=255)
                    cATTAATCAGGAACATCAGTGCCTTACTCATCTGGa           <  1:5572936/36‑1 (MQ=255)
                    cATTAATCAGGAACATCAGTGCCTTACTCATCTGGa           >  1:3840644/1‑36 (MQ=255)
                     aTTAATCAGGAACATCAGTGCCTTACTCATCTGGAt          <  1:2058239/36‑1 (MQ=255)
                     aTTAAGCAGGAACATCAGTGCCTTAATCATCTGGAt          >  1:336962/1‑36 (MQ=25)
                      ttAATCAGGAACAACAGTGCCTTACTCATCTGGAtt         >  1:4495810/1‑36 (MQ=37)
                       tAATCAGGAACATCATTGCCTTACTCCTCTGGCTTg        <  1:5254827/36‑1 (MQ=16)
                       tAATCAGGAACATCAGTGCCTTACTCATCTGGATTg        <  1:945858/36‑1 (MQ=255)
                        aaTCAGGAACATCAGTGCCTTACTCATCTGGATTgg       >  1:4483268/1‑36 (MQ=255)
                        aaTCAGGAACATCAGTGCCTTACTAATCTGGATTgg       >  1:1942451/1‑36 (MQ=37)
                         aTCAGGAACATCAGTGCCTTACTCATCTGGATTggg      <  1:2802661/36‑1 (MQ=255)
                         aTCAGGAACATCAGTGCCTTACTCATCTGGATTggg      >  1:503620/1‑36 (MQ=255)
                         aTCAGGAACATCAGTGCCTTACTCATCTGGATTggg      >  1:4297355/1‑36 (MQ=255)
                          tCAGGAACATCAGTGCCTTACTCATCTGGATTgggg     <  1:2292086/36‑1 (MQ=255)
                           cAGGAACATCAGTGCCTTACTCATCTGGATTGGGGt    <  1:2133308/36‑1 (MQ=255)
                           cAGGAACATCAGTGCCTTACTCATCTGGATTGGGGt    <  1:4590146/36‑1 (MQ=255)
                            aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   <  1:2697801/36‑1 (MQ=255)
                            aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   <  1:5534586/36‑1 (MQ=255)
                            aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   >  1:3769570/1‑36 (MQ=255)
                            aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   <  1:1819472/36‑1 (MQ=255)
                            aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   >  1:181150/1‑36 (MQ=255)
                             ggAACATCAGTGCCTTACTCATCTGGATTGGGGTCg  >  1:5815014/1‑36 (MQ=255)
                                |                                
CGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCG  >  REL606/649359‑649423

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG