Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,328,536	G→A	A148V (GCA→GTA)	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,536	0	G	A	100.0%	120.5 / NA	33	A148V (GCA→GTA)	mreB	cell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/24);  total (9/24)

CAGTGGTACCACCACCGATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGC  >  REL606/3328465‑3328606                                                                        |                                                                       cAGTGGTACCACCACCGATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTACGGc                                                                    <  2:273411/76‑1 (MQ=255)       tACCACCACCGATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTACGGCAGccct                                                              >  2:995422/1‑74 (MQ=255)         cccccaccGATATCAACCACCATAGAGCCGGTCACTTCAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCATCg                                                            <  1:770697/73‑1 (MQ=255)          caccaccGATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCATCgg                                                           >  1:523266/1‑76 (MQ=255)               ccGATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCATCGGttctt                                                      <  2:1734009/76‑1 (MQ=255)                 gATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCATCGGTTCTTCa                                                    <  2:1802581/76‑1 (MQ=255)                    atCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATc                                                 <  1:1782124/76‑1 (MQ=255)                     tCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATCa                                                <  1:1502916/76‑1 (MQ=255)                     acAACCACCATAGACCCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATCa                                                <  1:1155923/75‑1 (MQ=255)                           accaTAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGa                                          <  1:1760855/76‑1 (MQ=255)                               tAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCCTCGGGTCTTCCACCACGAAGACTTc                                      >  2:377441/1‑76 (MQ=255)                                aGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCCTCGGTTCTTCAATCAGGAAGACTTcc                                     >  1:979166/1‑75 (MQ=255)                                aGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCa                                     <  1:319680/76‑1 (MQ=255)                                aGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCATCGGGTCTTCAATCAGGAAGACTTCa                                     >  2:809997/1‑76 (MQ=255)                                      ggTCGCTTCAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCa                               <  1:1761462/76‑1 (MQ=255)                                      ggTCGCTTCAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCa                               <  2:892736/76‑1 (MQ=255)                                        tCGCTTCAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCAcc                             <  1:1361729/76‑1 (MQ=255)                                             tCAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCATCGGGTCTTCAATCAGGAAGACTTCACGGGCACCAgcgc                        >  2:432082/1‑76 (MQ=255)                                              cAGAAACCGGCAGGCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGcc                       <  1:1251473/76‑1 (MQ=255)                                                 aaaCCGGCAGGCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCATGGGCACCAGCGCCCTg                    <  2:872053/76‑1 (MQ=255)                                                     cGGCAGGCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGcgcg                <  2:877137/76‑1 (MQ=255)                                                     cGGCAGGCCAGCACCAATTACGGCAGCCATAGGTTCTTCAATCAGGAAAACTTCACGGGCACCAGCGCCCTGcgcg                >  1:164667/1‑76 (MQ=255)                                                        cAGGCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGAt             <  2:623338/76‑1 (MQ=255)                                                         aGGCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGAtt            <  2:1289549/76‑1 (MQ=255)                                                          ggCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTc           <  1:88503/76‑1 (MQ=255)                                                          ggCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTc           <  2:644251/76‑1 (MQ=255)                                                          ggCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTc           <  2:833112/76‑1 (MQ=255)                                                           gCCAGCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCa          <  2:696061/76‑1 (MQ=255)                                                             cAGCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACg        <  1:660077/76‑1 (MQ=255)                                                              aGCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGa       <  1:1780359/76‑1 (MQ=255)                                                               gCACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGaa      <  2:887635/76‑1 (MQ=255)                                                                cACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAAt     <  2:262207/76‑1 (MQ=255)                                                                cACCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAAt     >  2:413822/1‑76 (MQ=255)                                                                 aCCAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGTGCGGATTCACGAAtt    <  2:1799265/76‑1 (MQ=255)                                                                   cAATTACGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCCGCGCGGGTTCACGAAATgg  >  1:15299/1‑75 (MQ=255)                                                                        |                                                                       CAGTGGTACCACCACCGATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGC  >  REL606/3328465‑3328606

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