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

GCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCG  >  REL606/649358‑649423
                                 |                                
gCGGAACCTTCACGATACCGTCATTAATCAGGaaca                                >  1:467963/1‑36 (MQ=255)
gCGGAACCTTCACGATACCGTCATTAATCAGGaaca                                >  1:6159708/1‑36 (MQ=255)
 cGGAACCTTCACGATACCGTCATTAATCAGGaacat                               >  1:7766543/1‑36 (MQ=255)
 cGGAACCTTCACGATACCGTCATTAATCAGGaacat                               >  1:2351772/1‑36 (MQ=255)
   gAACCTTCACGATACCGTCATTAATCAGGAacatca                             >  1:322602/1‑36 (MQ=255)
   gAACCTTCACGATACCGTCATTAATCAGGAacatca                             >  1:5499082/1‑36 (MQ=255)
    aaCCTTCACGATACCGTCATTAATCAGGAACATCAg                            >  1:1981371/1‑36 (MQ=255)
    aaCCTTCACGATACCGTCATTAATCAGGAACATCAg                            >  1:5065740/1‑36 (MQ=255)
     aCCTTCACGATACCGTCATTAATCAGGAACATCAGt                           <  1:903144/36‑1 (MQ=255)
     aCCTTCACGATACCGTCATTAATCAGGAACATCAGt                           >  1:7482982/1‑36 (MQ=255)
     aCCTTCACGATACCGTCATTAATCAGGAACATCAGt                           <  1:2382673/36‑1 (MQ=255)
     aCCTTCACGATACCGTCATTAATCAGGAACATCAGt                           <  1:4637259/36‑1 (MQ=255)
      ccTTCACGATACCGTCATTAATCAGGAACATCAGTg                          <  1:1802085/36‑1 (MQ=255)
      ccTTCACGATACCGTCATTAATCAGGAACATCAGTg                          >  1:3882507/1‑36 (MQ=255)
        ttCACGATACCGTCATTAATCAGGAACATCAGTGcc                        >  1:2036333/1‑36 (MQ=255)
        ttCACGATACCGTCATTAATCAGGAACATCAGTGcc                        <  1:2241591/36‑1 (MQ=255)
         tCACGATACCGTCATTAATCAGGAACATCAGTGCCt                       >  1:6996577/1‑36 (MQ=255)
          cACGATACCGTCATTAATCAGGAACATCAGTGccgt                      >  1:2524654/1‑34 (MQ=25)
          cACGATACCGTCATTAATCAGGAACATCAGTGCCtt                      <  1:5151219/36‑1 (MQ=255)
          cACGATACCGTCATTAATCAGGAACATCAGTGCCtt                      >  1:1645120/1‑36 (MQ=255)
          cACGATACCGTCATTAATCAGGAACATCAGTGCCtt                      <  1:3588751/36‑1 (MQ=255)
           aCGATACCGTCATTAATCAGGAACATCAGTGCCTTa                     <  1:3459700/36‑1 (MQ=255)
            cGATACCGTCATTAATAAGGAACATCAGAGCATTAc                    >  1:5217354/1‑36 (MQ=19)
             gATACCGTCATTAATCAGGAACATCAGTGCCTTACt                   >  1:451736/1‑36 (MQ=255)
             gATACCGTCATTAATCAGGAACATCAGTGCCTTACt                   >  1:3979373/1‑36 (MQ=255)
             gATACCGTCATTAATCAGGAACATCAGTGCCTTACt                   <  1:3423105/36‑1 (MQ=255)
             gATACCGTCATTAATCAGGAACATCAGTGCCTTACt                   <  1:1198987/36‑1 (MQ=255)
             gATACCGTCATTAATCAGGAACATCAGTGCCTTACt                   >  1:1927364/1‑36 (MQ=255)
              aTACCGTCATTAATCAGGAACATCAGTGCCTTACTc                  <  1:321868/36‑1 (MQ=255)
              aTACCGTCATTAATCAGGAACATCAGTGCCTTACTc                  >  1:733658/1‑36 (MQ=255)
              aTACCGTCATTAATCAGGAACATCAGGGCCTTACTc                  >  1:6701764/1‑36 (MQ=25)
               tACCGTCATTAATCAGGAACATCAGTGCCTTACTCa                 <  1:1878028/36‑1 (MQ=255)
                 ccGTCATTAATCAGGAACATCAGTGCCTTACTCATc               >  1:6285975/1‑36 (MQ=255)
                  cGTCATTAATCAGGAACATCAGTGCCTTACTCATCt              >  1:1884686/1‑36 (MQ=255)
                  cGTCATTAATCAGGAACATCAGTGCCTTACTCATCt              <  1:404953/36‑1 (MQ=255)
                    tCATTAATCAGGAACATCAGTGCCTTACTCATCTgg            <  1:1185125/36‑1 (MQ=255)
                    tCATTAATCAGGAACATCAGTGCCTTACTCATCTgg            <  1:7318028/36‑1 (MQ=255)
                    tCATTAATCAGGAACATCAGTGCCTTACTCATCTgg            >  1:2500852/1‑36 (MQ=255)
                    tCATTAATCAGGAACATCAGTGCCTTACTCATCTgg            >  1:5439272/1‑36 (MQ=255)
                    tCATTAATCAGGAACATCAGTGCCTTACTCATCTgg            >  1:409371/1‑36 (MQ=255)
                     cATTAATCAGGAACATCAGTGCCTTACTCATCTGGa           >  1:7540985/1‑36 (MQ=255)
                     cATTAATCAGGAACATCAGTGCCTTACTCATCTGGa           <  1:3145157/36‑1 (MQ=255)
                     cATTAATCAGGAACATCAGTGCCTTACTCATATGGa           >  1:4386467/1‑36 (MQ=37)
                      aTTAATCAGGAACATCAGTGCCTTACTCATCTGGAt          >  1:7879308/1‑36 (MQ=255)
                      aTTAATCAGGAACATCAGTGCCTTACTCATCTGGAt          <  1:5892533/36‑1 (MQ=255)
                       ttAATCAGGAACATCAGTGCCTTACTCATCTGGAtt         <  1:3584191/36‑1 (MQ=255)
                        tAATCAGGAACATCAGTGCCTTACTCATCTGGATTg        <  1:4958676/36‑1 (MQ=255)
                        tAATCAGGAACATCAGTGCCTTACTCATCTGGATTg        >  1:7448974/1‑36 (MQ=255)
                        tAATCAGGAACATCAGTGCCTTACTCATCTGGATTg        >  1:5542208/1‑36 (MQ=255)
                         aaTCAGGAACATCAGTGCCTTACTCATCTGGATTgg       <  1:3070568/36‑1 (MQ=255)
                         aaTCAGGAACATCAGTGCCTTACTCATCTGGATTgg       >  1:2931396/1‑36 (MQ=255)
                         aaTCAGGAACATCAGTGCCTTACTCATCTGGATTgg       >  1:517449/1‑36 (MQ=255)
                         aaTCAGGAACATCAGTGCCTTACTCATCTGGATTgg       <  1:3965601/36‑1 (MQ=255)
                          aTCAGGAACATCAGTGCCTTACTCATCTGGATTggg      <  1:3470305/36‑1 (MQ=255)
                           tCAGGAACATCAGTGCCTTACTCATCTGTATTgggg     <  1:3471565/36‑1 (MQ=38)
                            cAGGAACATCAGTGCCTTACTCATCTGGATTGGGGt    >  1:6352267/1‑36 (MQ=255)
                            cAGGAACATCAGTGCCTTACTCATCTGGATTGGGGt    >  1:5469927/1‑36 (MQ=255)
                            cAGGAACATCAGTGCCTTACTCATCTGGATTGGGGt    <  1:7718115/36‑1 (MQ=255)
                            cAGGAACATCAGTGCCTTACTCATCTGGATTGGGGt    >  1:2830385/1‑36 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   >  1:6546892/1‑36 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   >  1:3109343/1‑36 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   <  1:3192948/36‑1 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   >  1:2773528/1‑36 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   >  1:894421/1‑36 (MQ=255)
                             aGGAACATCAGTGCCTTACTCATATGGATTGGGGTc   >  1:3101472/1‑36 (MQ=37)
                             aGGAACATCAGGGCCTTACTCATCTGGATTGGGGTc   <  1:381307/36‑1 (MQ=37)
                              ggAACATCAGTGCCTTACTCATCTGGATTGGGGTCg  >  1:2798819/1‑36 (MQ=255)
                              ggAACATCAGTGCCTTACTCATCTGGATTGGGGTCg  >  1:836677/1‑36 (MQ=255)
                                 |                                
GCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCG  >  REL606/649358‑649423

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG