Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	649,391	T→A	M471F (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
*	REL606	649,391	0	T→A	100.0%	160.0	53	M471F (ATC→TTC)	mrdA	transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  new base (25/28):  ref base (0/0):  total (25/28)

GCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCG  >  REL606/649358‑649423                                  |                                 gCGGAACCTTCACGATACCGTCATTAATCAGGaaca                                <  1:3535213/36‑1 (MQ=255) gCGGAACCTTCACGATACCGTCATTAATCAGGaaca                                >  1:3858030/1‑36 (MQ=255) gCGGAACCTTCACGATACCGTCATTAATCAGGaaca                                >  1:6948751/1‑36 (MQ=255)  cGGAACCTTCACGATACCGTCATTAATCAGGaacat                               >  1:7403802/1‑36 (MQ=255)   ggAACCTTCACGATACCGTCATTAATCAGGaacatc                              <  1:3672572/36‑1 (MQ=255)   ggAACCTTCACGATACCGTCATTAATCAGGaacatc                              >  1:2903414/1‑36 (MQ=255)   ggAACCTTCACGATACCGTCATTAATCAGGaacatc                              >  1:1824899/1‑36 (MQ=255)    gAACCTTCACGATACCGTCATTAATCAGGAccctca                             >  1:5224292/1‑36 (MQ=21)    gAACCTTCACGATACCGTCATTAATCAGGAccctca                             >  1:5226997/1‑36 (MQ=21)    gAACCTTCACGATACCGTCATTAATCAGGAacatca                             >  1:6483356/1‑36 (MQ=255)    aaaCCTTCACGATACCGTCATTATTCAGGaacataa                             <  1:1226132/35‑3 (MQ=21)     aaCCTTCACGATACCGTCATTAATCAGGAACATCAg                            <  1:235356/36‑1 (MQ=255)     aaCCTTCACGATACCGTCACTAATCAGGAacatcaa                            >  1:451453/1‑35 (MQ=25)     aaCCCTCACGATACCGTCATTAATCAGGAACATCAg                            <  1:1428527/36‑1 (MQ=18)      aCCTTCACGATACCGTCATTAATCAGGAACATCAGt                           <  1:2132416/36‑1 (MQ=255)      aCCTTCACGATACCGTCATTAATCAGGAACATCAGt                           <  1:3518439/36‑1 (MQ=255)       ccTTCACAATACCGTCATTAATCAGGAACATCAGTg                          <  1:2455308/36‑1 (MQ=21)        cTTCACGATACCGTCATTAATCAGGAACATCAGTGc                         >  1:6724409/1‑36 (MQ=255)        cTTCACGATACCGTCATTAATCAGGAACATCAGTGc                         >  1:6950234/1‑36 (MQ=255)        cTTCACGATACCGTCATTAATCAGGAACATCAGTGc                         >  1:7087781/1‑36 (MQ=255)        cTTCACGATACCGTCATTAATCAGGAACATCAGTGc                         <  1:4283029/36‑1 (MQ=255)         ttCACGATACCGTCATTAATCAGGAACATCAGTGcc                        <  1:570901/36‑1 (MQ=255)         ttCACGATACCGTCATTAATCAGGAACATCAGTGcc                        >  1:900870/1‑36 (MQ=255)          tCACGATACCGTCATTAATCAGGAACATCAGTGCCt                       <  1:7378973/36‑1 (MQ=255)           cACGTTACCGTCATTAATCAGGAACATCAGTGCCtt                      <  1:919479/36‑1 (MQ=21)           cACGATACCGTCATTAATCAGGAACATCAGTGCCtt                      <  1:138512/36‑1 (MQ=255)           cACGATACCGTCATTAATCAGGAACATCAGCGCCtt                      >  1:4563696/1‑36 (MQ=18)           cACAATACCGTCATTAATCAGGAACATCAGTGCCtt                      <  1:1762378/36‑1 (MQ=21)            aCGATACCGTCATTAATCAGGAACATCAGTGCCTTa                     >  1:3425091/1‑36 (MQ=255)             cGATATCGTCATTAATCAGGAACATCAGTGCCTTAc                    <  1:4467114/36‑1 (MQ=21)             cGATACCGTCATTAGTCAGGAACATCAGTGCCTTAc                    >  1:465473/1‑36 (MQ=25)             cGATACCGTCATTAATCCGGAACATCAGTGCCttcc                    >  1:647247/1‑34 (MQ=21)             cGATACCGTCATTAATCAGGAACATCAGTGCCTTAc                    <  1:1849118/36‑1 (MQ=255)              gATACCGTCATTAATCAGGAACATCAGTGCCTTTCt                   >  1:3661984/1‑36 (MQ=21)              gATACCGTCATTAATCAGGAACATCAGTGCCTTACt                   <  1:269849/36‑1 (MQ=255)              gATACCGTCATTAATCAGGAACATCAGTGCCTTACt                   >  1:6800315/1‑36 (MQ=255)               aTACCGTCATTAATCAGGAACATCAGTGCCTTACTc                  >  1:7077146/1‑36 (MQ=255)                tACCGTCATTAATCAGGAACATCAGTGCCTTACTCa                 >  1:6470921/1‑36 (MQ=255)                 aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt                >  1:1815763/1‑36 (MQ=255)                 aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt                <  1:7016857/36‑1 (MQ=255)                 aCCGTCATTAATCAGGAACATCAGTGCCTTACTCAt                >  1:4784945/1‑36 (MQ=255)                  ccGTCATTAATCAGGAACATCAGTGCCTTACTCATc               >  1:4545874/1‑36 (MQ=255)                   cGTCATTAATCAGGAACATCAGTGCCTTACTCCTCt              >  1:6177642/1‑36 (MQ=25)                   cGTCATTAATCAGGAACATCAGTGCCTTACTCATCt              <  1:6089906/36‑1 (MQ=255)                   cGTCATTAATCAGGAACATCAGTGCCTTACTCATCt              <  1:1932502/36‑1 (MQ=255)                    gTCATTAATCAGGAACATCAGTGCCTTACTCATCTg             >  1:1271576/1‑36 (MQ=255)                    gTCATTAATCAGGAACATCAGTGCCTTACTCATCTg             <  1:3559105/36‑1 (MQ=255)                    gTCATTAATCAGGAACATCAGTGCCTTACTCATCTg             <  1:2555996/36‑1 (MQ=255)                    atCATAAATCAGGAACATCATTGCCTTACTCATCTg             <  1:3245120/35‑1 (MQ=21)                      gaTTAATCAGTAACATCAGTGCCTTACTCATCTGGa           <  1:7410709/35‑1 (MQ=25)                      cATTTATTAGGAACATCGGTGCCTTACTCATCTGGa           <  1:4663457/36‑1 (MQ=18)                      cATTAATCAGGAACATCAGTGCCTTACTCATCTGGa           >  1:2455306/1‑36 (MQ=255)                       aTTAATCAGGAACATCAGTGCCTTACTCATCTGGAt          <  1:5888997/36‑1 (MQ=255)                        ttAATCAGGAACATCAGTGCCTTACTCATCTGGAtt         >  1:5775471/1‑36 (MQ=255)                        taaaTCAGGAACATCAGTGCCTTACTCATCTGCAtt         >  1:4955565/3‑36 (MQ=21)                          aaTCAGGAACATCAGTGCCTTACTCATCTCGAttag       >  1:1938505/1‑34 (MQ=21)                          aaTCAGGAACATCAGTGCCTTACTCATCTCGATTgg       >  1:4476733/1‑36 (MQ=25)                           aTCAGGAACATCAGTGCCTTACTCATCTCGCTTggg      >  1:3528447/1‑36 (MQ=21)                            tCATGAACATCAGTGCCTTACTCATCTGGATTgggg     <  1:2809205/36‑1 (MQ=25)                            tCAGGAACATCAGTGCCTTACTCATCTGGATTgggg     >  1:5075411/1‑36 (MQ=255)                            tCAGGAACATCAGTGCCTTACTCATCTGGATTgggg     >  1:279810/1‑36 (MQ=255)                             cAGGAACATCAGTGCCTTACTCATCTGGATTGGGGt    <  1:6451922/36‑1 (MQ=255)                              aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   <  1:1967186/36‑1 (MQ=255)                              aGGAACATCAGTGCCTTACTCATCTGGATTGGGGTc   <  1:3588360/36‑1 (MQ=255)                               ggAACATCAGTGCCTTACTCATCTGGATTGGGGTCg  <  1:6782451/36‑1 (MQ=255)                                  |                                 GCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCG  >  REL606/649358‑649423 Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG