Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,328,101	A→G	L293P (CTC→CCC)	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,101	0	A	G	100.0%	157.6 / NA	42	L293P (CTC→CCC)	mreB	cell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/22);  total (20/22)

AACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTAC  >  REL606/3328028‑3328168                                                                          |                                                                    aacaacGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGc                                                                   >  1:1188042/1‑76 (MQ=255)   caacGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCAc                                                                 >  1:990597/1‑76 (MQ=255)    aacGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCAcc                                                                <  2:1271603/76‑1 (MQ=255)      cGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCAt                                                              <  2:337252/76‑1 (MQ=255)        aCTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGc                                                            <  1:1055647/76‑1 (MQ=255)         cTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGcc                                                           >  1:95980/1‑76 (MQ=255)              aTGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCt                                                      <  2:903428/76‑1 (MQ=255)               tGCCGGTTTCTTCCATTACCAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTc                                                     <  2:1116596/76‑1 (MQ=255)                gCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCg                                                    >  2:1402084/1‑76 (MQ=255)                gCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCg                                                    >  1:646341/1‑76 (MQ=255)                  cgttttCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGga                                                  <  2:1331898/73‑1 (MQ=255)                  cGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGga                                                  <  2:1571301/76‑1 (MQ=255)                   ggTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGgag                                                 >  2:564238/1‑76 (MQ=255)                    gTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGgaga                                                >  2:1444144/1‑76 (MQ=255)                      ttcttcCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATg                                              >  1:77734/1‑76 (MQ=255)                               tAACAAGCGGTCAAGGTTACGCAGCAGTGCGCCGCCGCCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGcc                                     <  1:1593386/76‑1 (MQ=255)                               tAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGcc                                     <  1:260997/76‑1 (MQ=255)                                  cAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGt                                  >  1:1201841/1‑76 (MQ=255)                                   aagcGGTCAAGGTTACGCAGCAGTGCGCCACCGCCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGtt                                 <  1:514850/73‑1 (MQ=255)                                    aaCGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGCAGCCAGTTc                                >  2:1592298/1‑76 (MQ=255)                                    aaCGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTc                                <  1:1404670/76‑1 (MQ=255)                                     aCGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTcc                               <  1:620289/76‑1 (MQ=255)                                     aCGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTcc                               <  1:418656/76‑1 (MQ=255)                                       ngTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCcgg                             >  1:672383/2‑76 (MQ=255)                                        gTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCcggc                            >  2:610896/1‑76 (MQ=255)                                        gTCAAGGTTACGCAGCAGTGCCCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCcggc                            <  2:586506/76‑1 (MQ=255)                                          cAAGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCcggcgg                          <  2:1403852/76‑1 (MQ=255)                                           aaGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCggg                         >  1:269378/1‑76 (MQ=255)                                            aGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGc                        >  2:1182347/1‑76 (MQ=255)                                            aGGTTACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGc                        <  1:1475307/76‑1 (MQ=255)                                               ttACGCAGCAGTGCGCCACCACCGGAGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCaac                     >  2:1520743/1‑74 (MQ=255)                                                tACGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCCCTg                    >  1:1513927/1‑76 (MQ=255)                                                 aCGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGGACTGt                   >  2:951809/1‑76 (MQ=255)                                                  cGCAGCAGTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGtt                  >  2:1467031/1‑76 (MQ=255)                                                    cagcagTGCGCCGCCACCGGTGGGCCCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTcc                <  1:636962/76‑1 (MQ=255)                                                    cagcagTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTcc                <  1:498130/76‑1 (MQ=255)                                                    cagcagTGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACCCTTcc                >  2:348860/1‑76 (MQ=255)                                                         gTGCGCCACCCCCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGc           <  2:710100/76‑1 (MQ=255)                                                          tGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCTGTGCa          >  2:776844/1‑76 (MQ=255)                                                          tGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCa          <  1:1397955/76‑1 (MQ=255)                                                          tGCGCCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACCGGTTCCGTGCa          >  2:1362740/1‑76 (MQ=255)                                                            cgcCACCACCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAAc        <  1:671093/76‑1 (MQ=255)                                                               caccacCGGTGGGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCAt     <  2:139238/76‑1 (MQ=255)                                                                  cgccGGTGGGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTAc  <  2:1332683/74‑1 (MQ=255)                                                                          |                                                                    AACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTAC  >  REL606/3328028‑3328168

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