Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	649,520	T→C	N428D (AAC→GAC)	mrdA ←	transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	649,520	0	T	C	100.0%	227.9 / NA	54	N428D (AAC→GAC)	mrdA	transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/34);  total (20/34)

CAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAG  >  REL606/649448‑649592                                                                         |                                                                         cAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGcc                                                                       <  1:475941/76‑1 (MQ=255)  aaCCGTAATGGTGTCACCCTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCg                                                                      <  2:1210502/76‑1 (MQ=255)   aCCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCgg                                                                     <  1:1123128/76‑1 (MQ=255)   aCCGAAATGGTGTCA‑CCAGATTCCACCGTTTTTTAAAAGCTTTCTGTTTCCCATCGCGGGTAGGCATGTCGCCGGa                                                                    <  2:948640/76‑1 (MQ=255)    ccGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGa                                                                    >  1:778794/1‑76 (MQ=255)     cGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGaa                                                                   <  2:1084506/76‑1 (MQ=255)      ggAATGGTGTCACCCTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAAc                                                                  >  2:936865/1‑76 (MQ=255)        taTGGTGTCACCCTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGt                                                                >  2:1760512/2‑76 (MQ=255)            gtgtCACCCTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGttctt                                                            <  1:1836433/76‑1 (MQ=255)            gtgtCACCCTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGttctt                                                            <  2:1693533/76‑1 (MQ=255)                 ccccTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGcc                                                       <  2:741665/75‑1 (MQ=255)                 aCCCTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGGAGGCATGTCGCCGGAACGTTCTTCCGcc                                                       >  2:1160825/1‑76 (MQ=255)                  cccTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCa                                                      >  2:1092045/1‑76 (MQ=255)                   ccTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAg                                                     <  1:656326/76‑1 (MQ=255)                   ccTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAg                                                     <  1:128860/76‑1 (MQ=255)                    cTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAgg                                                    <  2:928170/76‑1 (MQ=255)                     tGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGt                                                   <  2:852690/76‑1 (MQ=255)                          ccACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATa                                              <  2:589473/76‑1 (MQ=255)                          ccACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATa                                              <  2:1651635/76‑1 (MQ=255)                           cACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATAc                                             <  1:898605/76‑1 (MQ=255)                                ttttttAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGgtg                                        <  2:963663/76‑1 (MQ=255)                                  ttttAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTa                                      >  1:1291158/1‑76 (MQ=255)                                  ttttAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTa                                      <  1:166294/76‑1 (MQ=255)                                  ttttAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTa                                      <  2:1444226/76‑1 (MQ=255)                                  ttttAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTa                                      <  1:1738935/76‑1 (MQ=255)                                   tttAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTaa                                     >  1:16698/1‑76 (MQ=255)                                    ttAAGGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAAt                                    >  2:472508/1‑76 (MQ=255)                                     tAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATg                                   >  1:895942/1‑76 (MQ=255)                                     tAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATg                                   <  1:1760512/76‑1 (MQ=255)                                           gTTTCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATa                             <  1:224400/76‑1 (MQ=255)                                            tttCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATaa                            <  1:1132287/76‑1 (MQ=255)                                            tttCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATaa                            >  2:864490/1‑76 (MQ=255)                                            tttCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATaa                            <  1:501526/76‑1 (MQ=255)                                             ttCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCAtngc                           <  2:1496200/76‑4 (MQ=255)                                             ttCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAAc                           >  1:1007110/1‑76 (MQ=255)                                              tCTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAAcc                          >  2:813785/1‑76 (MQ=255)                                               cTGTTTCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCg                         >  2:369319/1‑76 (MQ=255)                                               cTGTTTCCATTCGCGGGCAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGCGTAATGCCCATAACCg                         >  1:531626/1‑76 (MQ=255)                                                  tttCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAAt                      >  2:947650/1‑76 (MQ=255)                                                   ttCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAAtt                     >  2:1018884/1‑76 (MQ=255)                                                    tCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGaangt                    <  2:1497063/76‑4 (MQ=255)                                                    tCCATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAAttt                    <  1:463093/76‑1 (MQ=255)                                                     ccATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTa                   >  2:1569051/1‑76 (MQ=255)                                                      cATTCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTAc                  >  2:476011/1‑76 (MQ=255)                                                         tCGCGGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACcca               <  2:1728853/76‑1 (MQ=255)                                                          cgcgGGTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACccat              <  1:579392/76‑1 (MQ=255)                                                             gggTAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACccatcca           <  2:1153256/76‑1 (MQ=255)                                                               gTAGGCATGTCGCCGGAGCGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCAtt         <  1:1171378/76‑1 (MQ=255)                                                                tAGGCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTc        >  2:1573757/1‑76 (MQ=255)                                                                  ggCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCgg      <  2:660688/76‑1 (MQ=255)                                                                  ggCATGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGCGTAATGACCATAACCGAATTTACCCATCCATTCgg      <  2:704093/76‑1 (MQ=255)                                                                     aTGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGgaga   <  2:1348714/76‑1 (MQ=255)                                                                      tGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGagag  >  2:529685/1‑76 (MQ=255)                                                                      tGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGagag  <  2:2124879/76‑1 (MQ=255)                                                                      tGTCGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGagag  <  2:2051383/76‑1 (MQ=255)                                                                         |                                                                         CAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAG  >  REL606/649448‑649592

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑