Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	648,086	T→C	E271G (GAG→GGG)	mrdB ←	cell wall shape‑determining protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	648,086	0	T	C	100.0%	84.1 / NA	22	E271G (GAG→GGG)	mrdB	cell wall shape‑determining protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/11);  total (11/11)

GGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGACTGAGT  >  REL606/648013‑648153                                                                          |                                                                    ggcTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGccct                                                                   <  1:82718/76‑1 (MQ=255)     aTCNACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCANAATCAGAATGCCCACTAATCCCAGCCCTTCCg                                                               <  2:1091530/76‑1 (MQ=255)       ccACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCCCTTCCGcc                                                             <  1:338278/76‑1 (MQ=255)        cacaGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCCCTTCCGCCa                                                            <  1:763316/76‑1 (MQ=255)         acaGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCCCTTCCGCCAg                                                           >  2:101450/1‑76 (MQ=255)          caGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCCCTTCCGCCAGt                                                          <  1:1140596/76‑1 (MQ=255)             cccGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCCCTTCCGCCAGTAcc                                                       <  2:53421/76‑1 (MQ=255)             cccGCGCATCATCAGCAAAATGTCGAGAGCGAGCAGAATCCCAATGCCCCCTAATCCCAGCCCTTCCGCCAGTAcc                                                       <  2:53415/76‑1 (MQ=255)                   cATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCCCTTCCGCCAGTACCGCGAAg                                                 >  1:412812/1‑76 (MQ=255)                   cATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCCCTTCCGCCAGTACCGCGAAg                                                 >  1:671632/1‑76 (MQ=255)                           gcagAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCCCTTCCGCCAGTACCGCGAAGATAAAGTc                                         <  1:425900/76‑1 (MQ=255)                               aaTGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCCCTTCCGCCAGTACCGCGAAGATAAAGTCGGTa                                     <  2:429326/76‑1 (MQ=255)                                    agagagCGAGCAGAATCAGAATGCCCACTAATCCCAGCCCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCg                                >  1:183457/1‑76 (MQ=255)                                      agagCGAGCAGAATCAGANTGCCCACTAATCCCAGCCCTCCCGCCAATACCGCGAAGATAAAGTCGGTATGGCGtt                              >  1:1112471/1‑76 (MQ=255)                                       gagcgagcAGAATCAGAATGCCCACTAATCCCAGCCCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTc                             >  2:198862/1‑76 (MQ=255)                                           gagcAGAATCAGAATGCCCACTAATCCCAGCCCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGt                         <  2:119473/76‑1 (MQ=255)                                              cagaatcagaatGCCCACTAATCCCAGCCCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGa                      <  2:543348/76‑1 (MQ=255)                                                         tGCCCACTAATCCCAGCCCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCtg           >  2:920359/1‑76 (MQ=255)                                                           cccACTAACCCCAGCCCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCtgtc         <  1:775290/76‑2 (MQ=255)                                                                tAATCCCAGCCCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGACTGa    >  2:315921/1‑76 (MQ=255)                                                                 aaTCCCAGCCCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGACTGAg   >  2:509510/1‑76 (MQ=255)                                                                 aaTCCCAGCCCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGACTGAg   >  2:900584/1‑76 (MQ=255)                                                                  aTCCCAGCCCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGACTGAGt  >  1:952011/1‑76 (MQ=255)                                                                          |                                                                    GGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGACTGAGT  >  REL606/648013‑648153

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