Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,328,582	C→A	A133S (GCT→TCT)	mreB ←	cell wall structural complex MreBCD, actin‑like component MreB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,328,582	0	C	A	100.0%	82.6 / NA	23	A133S (GCT→TCT)	mreB	cell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (8/15);  total (8/15)

AGAAACCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCG  >  REL606/3328511‑3328654                                                                        |                                                                         aGAAACCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGAGccc                                                                      <  2:365751/76‑1 (MQ=255)   aaaCCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTNAATCAGGAAGACTTCACGGGCACCAGAGCCCTg                                                                    <  2:1024061/76‑1 (MQ=255)       cGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGAGCCCTGcgcg                                                                >  2:605950/1‑76 (MQ=255)                    ccAATGGGGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGAGCCCTGCGCGGATTCACGAATTg                                                   <  1:239954/76‑1 (MQ=255)                      aaTTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGAGCCCTGCGCGGATTCACGAATTgag                                                 <  2:139091/76‑3 (MQ=255)                        ttGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGAGCCCTGCGCGGATTCACGAATTGcgcg                                               <  2:807169/76‑1 (MQ=255)                           cGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGAGCCCTGCGCGGATTCACGAATTGCgcggcg                                            <  1:692964/76‑1 (MQ=255)                                gCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGAGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCaa                                       <  2:549291/76‑1 (MQ=255)                                 ccATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGAGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAAc                                      >  2:206403/1‑76 (MQ=255)                                  cATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCNGAGCCCTGCGAGGATACACGAATTGCGCGGCGTTNAAcc                                     >  1:142755/1‑76 (MQ=255)                                     cGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGAGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTgg                                  <  2:1076487/76‑1 (MQ=255)                                      ggTTCTTCAATCAGGAAGACTTCACGGGCACCAGAGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTggg                                 <  1:474667/76‑1 (MQ=255)                                       gTTCTTCAATCAGGAAGACTTCACGGGCACCAGAGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGt                                >  1:488236/1‑76 (MQ=255)                                         tcttcAATCAGGAAGACTTCACGGGCACCAGAGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTcg                              >  1:290172/1‑76 (MQ=255)                                                 cAGGAAGACTTCACGGGCACCAGAGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAAcc                      <  2:818103/76‑1 (MQ=255)                                                 cAGGAAGACTTCACGGGCACCAGAGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAAcc                      <  2:942866/76‑1 (MQ=255)                                                       gACTTCACGGGCACCATAGCCCTGCGCCGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAAACGGcaca                >  1:60600/1‑76 (MQ=255)                                                       gACTTCACGGGCACCAGAGCCCTGCGCGGATTCACGAATTGTGCGGCGTTCAACCTGGGTCGCGCCAACCGGcaca                >  1:615194/1‑76 (MQ=255)                                                           tCACGGGCACCAGAGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACaaa            >  1:582438/1‑76 (MQ=255)                                                               gggCACCAGAGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACTACCCAg        >  1:1196559/1‑76 (MQ=255)                                                                ggCACCAGAGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGa       <  1:388216/76‑1 (MQ=255)                                                                 gCACCAGAGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCATCCTGGGTCGCGCCAACCGGCACATAAACCAGaa      <  2:1035379/76‑1 (MQ=255)                                                                   aCCAGAGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAAcg    <  1:564736/76‑1 (MQ=255)                                                                     cAGAGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAAcgcg  <  1:582384/76‑1 (MQ=255)                                                                        |                                                                         AGAAACCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCG  >  REL606/3328511‑3328654

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