Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	648,276	G→A	R208C (CGC→TGC)	mrdB ←	cell wall shape‑determining protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	648,276	0	G	A	100.0%	127.3 / NA	34	R208C (CGC→TGC)	mrdB	cell wall shape‑determining protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/19);  total (15/19)

TGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGTAATCATGCATCAGGAAGAACCACAGAATCGGAATGAACGCCGCTACCAGCACTACTGCGAC  >  REL606/648209‑648342                                                                    |                                                                   tGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAAt                                                            >  1:610116/1‑76 (MQ=255)  gAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATc                                                           <  2:1027170/76‑1 (MQ=255)   aataatGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCa                                                          >  2:905274/1‑76 (MQ=255)   aataatGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGCCCAGGAGCATCATTACGCGCTGGCACTGGTAATCa                                                          >  1:339089/1‑76 (MQ=255)    ataatGTGCTAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCAt                                                         >  1:1240744/1‑76 (MQ=255)    ataatGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGNTCTAGGAGCATCATTACGCGCTGGCACTGGTAATCAt                                                         <  1:1308933/76‑1 (MQ=255)      aatGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGc                                                       >  1:765628/1‑76 (MQ=255)         gtgATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCAtn                                                    <  1:229237/76‑2 (MQ=255)            aTAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAgg                                                 <  1:972511/76‑1 (MQ=255)            aTAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAgg                                                 <  1:972514/76‑1 (MQ=255)             tAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGga                                                >  1:358333/1‑76 (MQ=255)                 ccGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGgaagaa                                            >  1:421993/1‑76 (MQ=255)                  cgcgCCGANTGGTTCTTATTCCGGTTCTAGGAGTTTCATTACGCGCTGGCACTGGAAATCATACATCAGAAAGAAc                                           <  2:218375/76‑1 (MQ=255)                  cgcgCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAAc                                           <  1:419101/76‑1 (MQ=255)                  cgcgCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAAc                                           <  1:910582/76‑1 (MQ=255)                     gcCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGCAGAACcac                                        >  1:446212/1‑76 (MQ=255)                     gcCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAACcac                                        <  2:1041929/76‑1 (MQ=255)                      ccGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAACcaca                                       <  1:966067/76‑1 (MQ=255)                        gAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAACCACAGa                                     <  2:683715/76‑1 (MQ=255)                        gAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAACCACAGa                                     <  2:360746/76‑1 (MQ=255)                         aGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAACCACAGaa                                    <  1:180701/76‑1 (MQ=255)                             ggTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAACCACAGAATCgg                                <  1:1028039/76‑1 (MQ=255)                              gTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAACCACAGAATCGGa                               >  1:740344/1‑76 (MQ=255)                              gTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAACCACAGAATCGGa                               <  2:881769/76‑1 (MQ=255)                                   aTTCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAACCACCGAATCGGAATGaa                          >  2:564404/1‑76 (MQ=255)                                    ttCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTGATCATGCATCAGGAAGAACCACAGAATCGGAATGAAc                         <  1:994635/76‑1 (MQ=255)                                    ttCCGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAACCACAGAATCGGAATGAAc                         >  2:675052/1‑76 (MQ=255)                                      ncGGGTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAACCACAGAATCGGAATGAAcgc                       >  1:470560/2‑76 (MQ=255)                                          gTCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAACCACAGAATCGGAATGAACGCCGCt                   <  2:791466/76‑1 (MQ=255)                                           tCCAGGAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAACCACAGAATCGGAATGAACGCCGCTa                  <  1:1008913/76‑1 (MQ=255)                                               ggAGCATCATTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAACCACAGAATCGGAATGAACGCCGCTACCAg              <  2:1164142/76‑1 (MQ=255)                                                    atcntTACGCGCTGGCACTGGTAATCATGCATCAGGAAGAACCACAGAATCGGAATGAACGCCGCTACCAGCacta         >  1:1203548/1‑76 (MQ=255)                                                      catTACGCGCTGGCACTGGTAATCATGCATCAGGCAGAACCNCAGAATCGGAATGAACGCNGCNACCAGCactact       >  1:1096819/1‑76 (MQ=255)                                                           cgcgCTGGCACTGGTAATCATGCATCAGGAAGAACCACAGAATCGGAATGAACGCCGCTACCAGCACTACTGCGAc  >  1:1192935/1‑76 (MQ=255)                                                                    |                                                                   TGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGTAATCATGCATCAGGAAGAACCACAGAATCGGAATGAACGCCGCTACCAGCACTACTGCGAC  >  REL606/648209‑648342

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