Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	648,692	C→T	R69H (CGC→CAC)	mrdB ←	cell wall shape‑determining protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	648,692	0	C	T	100.0%	60.7 / NA	16	R69H (CGC→CAC)	mrdB	cell wall shape‑determining protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/7);  total (9/7)

TCTACCGCCACCAGCAAAATAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCC  >  REL606/648629‑648760                                                                |                                                                     tCTACCGCCACCAGCAAAATAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGTGTGGAGGAAttt                                                          >  1:429374/1‑76 (MQ=255)            caGCAAAATAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGTGTGGAGGAATTTGCGCCATcacc                                               <  2:668814/76‑1 (MQ=255)              gCAAAATAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGTGTGGAGGAATTTGCGCCATcaccac                                             >  1:93984/1‑76 (MQ=255)              gCAAAATAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGTGTGGAGGAATTTGCGCCATcaccac                                             <  2:180689/76‑1 (MQ=255)               cAAAATAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGTGTGGAGGAATTTGCGCCATCaccacc                                            <  1:54650/76‑1 (MQ=255)               cAAAATAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGTGTGGAGGAATTTGCGCCATCaccacc                                            >  2:213922/1‑76 (MQ=255)                    taataCAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGTGTGGAGGAATTTGCGCCATCACCACCatgat                                       >  1:112527/1‑76 (MQ=255)                          aGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGTGTGGAGGAATTTGCGCCATCACCACCATGATGACCAg                                 <  1:489360/76‑1 (MQ=255)                                 ataGAGATAGGGGGCCCAGCCTTCATAAACGTGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATc                          >  2:104074/1‑76 (MQ=255)                                     agaTAGGGGGCCCAGCCTTCATAAACGTGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGa                      >  2:218313/1‑76 (MQ=255)                                        tAGGGGGCCCAGCCTTCATAAACGTGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGAttt                   <  1:76065/76‑1 (MQ=255)                                             ggCCCAGCCTTCATAAACGTGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCg              <  1:321168/76‑1 (MQ=255)                                                 cAGCCTTCATAAACGTGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGAttt          <  1:416418/76‑1 (MQ=255)                                                     cTTCATAAACGTGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACg      >  2:74642/1‑76 (MQ=255)                                                       tCATAAACGTGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCt    >  1:314669/1‑76 (MQ=255)                                                         aTAAACGTGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTcc  >  1:568904/1‑76 (MQ=255)                                                                |                                                                     TCTACCGCCACCAGCAAAATAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCC  >  REL606/648629‑648760

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