Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,328,962	C→T	R6H (CGT→CAT)	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,962	0	C	T	100.0%	264.2 / NA	70	R6H (CGT→CAT)	mreB	cell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (44/26);  total (44/26)

CTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCGCTTACTTTACCAA  >  REL606/3328891‑3329034                                                                        |                                                                         cTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGaaa                                                                      >  2:1461506/1‑76 (MQ=255) cTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGaaa                                                                      >  2:819219/1‑76 (MQ=255) cTTGTCCCTTGACATAAATGAGGGGATGCGCAGCACCCANNTCAATGGNCAAGTCATTGGAAAACATGCCATGaaa                                                                      <  1:2011450/76‑1 (MQ=255)  ttGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAAt                                                                     <  1:1172371/76‑1 (MQ=255)    gTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAAttt                                                                   <  2:979588/76‑1 (MQ=255)     tCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAAtttt                                                                  >  2:1559899/1‑76 (MQ=255)      ccTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAAttttt                                                                 >  1:592156/1‑76 (MQ=255)          ttACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCaa                                                             <  2:1557614/76‑1 (MQ=255)          ttACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCaa                                                             >  2:1265841/1‑76 (MQ=255)          ttACATAAATGAGGAAATTCGCAATACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCaa                                                             >  1:185753/1‑76 (MQ=255)             cATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAAAAt                                                          >  2:1438887/1‑76 (MQ=255)              aTAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATa                                                         >  1:228513/1‑76 (MQ=255)               tAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATAc                                                        >  2:1743598/1‑76 (MQ=255)               tAAACGATGGTATCCTCAGTAACCAGGGCAATCGACAAGTGATTGCAAAGCATGGCATAAAATTTTTTCAACATAc                                                        <  1:2046705/76‑1 (MQ=255)                 aaTGAGGGTATTCGCAGTACCCTGGNCTATGGACAAGTCATTGGAATNCATGCCATGAAATTTTTTCAACATACTa                                                      <  1:299998/76‑1 (MQ=255)                 aaTGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTa                                                      >  2:711727/1‑76 (MQ=255)                  aTGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTaa                                                     >  2:1291783/1‑76 (MQ=255)                  aTGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTaa                                                     >  2:1291797/1‑76 (MQ=255)                  aTGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTaa                                                     >  2:587025/1‑76 (MQ=255)                   tcaGGGCATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGCACCCATGCCATGAAATTTTCGCAACATCCTAAg                                                    >  2:967828/3‑76 (MQ=255)                   tGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAg                                                    <  2:512342/76‑1 (MQ=255)                   tGAGGGTAGTCGCAGTACGCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAg                                                    <  2:512347/76‑1 (MQ=255)                    gAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAgg                                                   >  1:424472/1‑76 (MQ=255)                    gAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAgg                                                   >  1:1386953/1‑76 (MQ=255)                       ggTATTCGCATTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGAt                                                <  1:746481/76‑1 (MQ=255)                       ggTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGAt                                                <  2:1752977/76‑1 (MQ=255)                        gTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATa                                               <  2:1059636/76‑1 (MQ=255)                         tATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAACCATGCCATGAAATTTTTTCAACATACTAAGTGATaa                                              <  1:129127/76‑1 (MQ=255)                         tATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATaa                                              >  2:1468211/1‑76 (MQ=255)                         tATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATaa                                              >  2:678711/1‑76 (MQ=255)                         tATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATaa                                              <  2:1749257/76‑1 (MQ=255)                         tATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATaa                                              >  1:1543560/1‑76 (MQ=255)                         tATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATaa                                              <  2:773728/76‑1 (MQ=255)                          aTTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAAt                                             >  1:1192215/1‑76 (MQ=255)                          aTTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAAt                                             >  1:1396364/1‑76 (MQ=255)                            tCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATcc                                           >  2:1740698/1‑76 (MQ=255)                                aGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGaa                                       >  2:823093/1‑76 (MQ=255)                                aGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGaa                                       >  2:111478/1‑76 (MQ=255)                                 gTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGaaa                                      <  1:468656/76‑1 (MQ=255)                                  tACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAg                                     <  1:996534/76‑1 (MQ=255)                                  tACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAg                                     >  2:1436287/1‑76 (MQ=255)                                   aCCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAGATTTTTTCAGCATACTAAGGGATAATCCTGAAAGc                                    >  2:175924/1‑76 (MQ=255)                                   aCCCAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGc                                    >  1:1327286/1‑76 (MQ=255)                                    cccAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCt                                   >  2:1142707/1‑76 (MQ=255)                                     ccAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTg                                  >  1:828701/1‑76 (MQ=255)                                      cAGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTgg                                 <  1:1218883/76‑1 (MQ=255)                                       aGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTggg                                >  2:970664/1‑76 (MQ=255)                                       aGGTCAATGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTggg                                <  2:895622/76‑1 (MQ=255)                                             aTGGACCAGTCATTGGACTACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGNAANCTGGGGCGgaa                          >  2:322419/1‑76 (MQ=255)                                             aTGGACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGTGCGgaa                          <  1:1895908/76‑1 (MQ=255)                                             aTGGACAAGTCATTGGAAAACATACAATGAAATTTTTTCNACATACTAAGGGATANTCCTGAGAGCTCGGGCGgaa                          >  2:2018356/1‑76 (MQ=255)                                               ggACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaa                        >  1:1035063/1‑76 (MQ=255)                                                gACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaag                       >  1:489906/1‑76 (MQ=255)                                                gACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaag                       >  2:199124/1‑76 (MQ=255)                                                gACAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaag                       >  1:1734910/1‑76 (MQ=255)                                                 aCAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaaga                      >  1:1030600/1‑76 (MQ=255)                                                  cAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaa                     >  2:1940748/1‑76 (MQ=255)                                                  cAAGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaa                     >  2:834017/1‑76 (MQ=255)                                                   aaGTCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaaa                    >  1:1500769/1‑76 (MQ=255)                                                      tCATTGGAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATc                 <  2:1030599/76‑1 (MQ=255)                                                           ggAAAACGTGCCATGAAATTCTTTAAACATACTAAGGGATAATCCTGAAAGCTGGGGCTGATAAGAAAATCCGCtt            >  1:92552/1‑76 (MQ=255)                                                           ggAAAACATGCCATNAAATTTTTTCAACATACTAAGGGATAAGCCTGAAAGCTGGGGCGGTAATGAAAATCCGCtt            <  1:1753443/76‑1 (MQ=255)                                                           ggAAAACATACCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCGCtt            >  1:1327526/1‑76 (MQ=255)                                                            gAAAACATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCGCTTa           <  2:408534/76‑1 (MQ=255)                                                             naaaCATGCCATGAAATTTTTTCAACATACTAAGCGANAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCGCTTAc          >  1:1983279/2‑76 (MQ=255)                                                                 cATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCGCTTACTTTa      <  2:789922/76‑1 (MQ=255)                                                                 cATGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCGCTTACTTTa      >  1:1219931/1‑76 (MQ=255)                                                                  aTGCCATGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCGCTTACTTTac     <  2:190949/76‑1 (MQ=255)                                                                    gCCATGAAATTTTTTCAACATACTAAGGGATAATCCGGATAGCTGGGGNGGAAAAGAAAATCCGCTTACTTTacca   <  1:1799798/76‑1 (MQ=255)                                                                     ccATGAAATTTTNTCNACATACTAAGGGATAATCNTGAAAGTTGGGGCGGAAAAGAAAATCCGCTTACTTTaccaa  <  1:1666532/76‑1 (MQ=255)                                                                        |                                                                         CTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCGCTTACTTTACCAA  >  REL606/3328891‑3329034

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