Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,328,543	C→T	A146T (GCT→ACT)	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,543	0	C	T	100.0%	151.2 / NA	41	A146T (GCT→ACT)	mreB	cell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/20);  total (21/20)

GTACCACCACCGATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTC  >  REL606/3328470‑3328615                                                                          |                                                                         gTACCACCACCGATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCa                                                                        <  1:1344447/76‑1 (MQ=255)   accaccaccGATATCAACCACCATAGAACCGGTNGCTTCAGAAACCGGCAGGCCAGNACCAATTGCGNCAGTCATc                                                                      >  2:1360626/1‑76 (MQ=255)   accaccaccGATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATc                                                                      <  1:642123/76‑1 (MQ=255)   accaccaccGATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATc                                                                      >  2:179362/1‑76 (MQ=255)    ccaccaccGATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCg                                                                     <  2:732614/76‑1 (MQ=255)       ccaccGATNTCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGtt                                                                  <  1:1415297/76‑1 (MQ=255)       ccaccGATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGtt                                                                  <  2:594207/76‑1 (MQ=255)           cGATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACGAAGTGAGGCAGCCATCGGctcttc                                                              >  1:1533611/1‑76 (MQ=255)            gATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCa                                                             >  1:691245/1‑76 (MQ=255)            gATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCa                                                             <  1:1019925/76‑1 (MQ=255)            gATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCa                                                             <  2:1197542/76‑1 (MQ=255)             atatCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCaa                                                            <  1:1085418/76‑1 (MQ=255)             atatCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCaa                                                            <  1:1343745/76‑1 (MQ=255)             atatCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCAGCAGTCATCGGTTCTTCaa                                                            <  2:1284167/76‑1 (MQ=255)               atCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATc                                                          >  1:955771/1‑76 (MQ=255)                tCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCa                                                         >  2:1497988/1‑76 (MQ=255)                  aaCCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAgg                                                       >  1:551117/1‑76 (MQ=255)                       ccaTAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTTAATCAGGAAGAc                                                  >  1:741717/1‑76 (MQ=255)                       ccaTAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGAc                                                  >  2:396770/1‑76 (MQ=255)                        caTAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACt                                                 >  1:563239/1‑76 (MQ=255)                        caTAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACt                                                 >  1:636957/1‑76 (MQ=255)                        caTAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACt                                                 >  1:1209701/1‑76 (MQ=255)                         aTAGAACCCGTCGCTTCAGAAACCGGCAGGCCAGCACCCATTGCCGGGATCATCGGTTTCTCAACCACGAAGagct                                                >  1:1464937/1‑73 (MQ=255)                             aaCCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCACg                                            >  2:639011/1‑76 (MQ=255)                                 ggTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCa                                        >  1:9219/1‑76 (MQ=255)                                    cGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATNGGTTCTTCAATCAGGAAGACTTCACGGGCACCa                                     <  2:1376526/76‑1 (MQ=255)                                    cGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCa                                     <  1:488448/76‑1 (MQ=255)                                    cGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCa                                     <  2:1052611/76‑1 (MQ=255)                                     gtttCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAg                                    <  1:795402/74‑1 (MQ=255)                                     gCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCTCGGGCACCAg                                    <  1:1497988/76‑1 (MQ=255)                                          aGAAACCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGccc                               <  1:865646/76‑1 (MQ=255)                                            aaaCCGGCAGGCCAGCACCNATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTg                             >  1:1504476/1‑76 (MQ=255)                                             aaCCGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTgc                            >  1:791648/1‑76 (MQ=255)                                               ccGGCAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAAACAGGAAGACTTCACGGGCACCAGTGCCCTgcgc                          <  2:1381490/76‑1 (MQ=255)                                                   cAGGCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGAt                      >  1:159142/1‑76 (MQ=255)                                                     ggCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTc                    <  1:744575/76‑1 (MQ=255)                                                     ggCCAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTc                    <  1:1359248/76‑1 (MQ=255)                                                        cAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACg                 >  2:1271458/1‑76 (MQ=255)                                                        cAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACg                 >  1:1497390/1‑76 (MQ=255)                                                        cAGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACg                 >  1:1301178/1‑76 (MQ=255)                                                         aGCACCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGATTTCACGGGCACCAGCGCCCTGCGCGGATTCACGa                >  2:802661/1‑76 (MQ=255)                                                            aCCAATTGCGGCAGTCATCGGTTCTTCAATCAGGAANACTTCAGGGGCACCAGCGCCCTGCGCGGATTCACNAAtt             >  1:193455/1‑76 (MQ=255)                                                             ccAATTGCGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTg            <  2:523506/76‑1 (MQ=255)                                                                    cGGCAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCgcggcg     >  1:1397923/1‑76 (MQ=255)                                                                      gCAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGtt   >  2:762877/1‑76 (MQ=255)                                                                       cAGTCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTc  >  2:958871/1‑76 (MQ=255)                                                                          |                                                                         GTACCACCACCGATATCAACCACCATAGAACCGGTCGCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTC  >  REL606/3328470‑3328615

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