Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,400,773	A→G	R105R (CGT→CGC)	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
*	NC_000913	3,400,773	0	A	G	100.0%	104.4 / NA	29	R105R (CGT→CGC)	mreB	cell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/16);  total (13/16)

ATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACG  >  NC_000913/3400642‑3400895                                                                                                                                    |                                                                                                                           aTTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATaa                                                                                                                     >  1:129811/1‑139 (MQ=255) aTTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATaa                                                                                                                     >  1:195830/1‑139 (MQ=255)         aGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTg                                                                                                             <  1:66939/139‑1 (MQ=255)           ccATCGGTTCTTCAATCGGGACGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCt                                                                                                           <  1:247419/139‑1 (MQ=255)                 gTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCAc                                                                                                     <  1:588804/139‑1 (MQ=255)                   tcttcAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACtt                                                                                                   <  1:341689/139‑1 (MQ=255)                   tcttcAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCAAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACtt                                                                                                   <  1:290147/139‑1 (MQ=255)                    cttcAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTg                                                                                                  >  2:463579/1‑139 (MQ=255)                           cAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCAc                                                                                                     <  2:31316/129‑1 (MQ=255)                           cAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCAc                                                                                                     >  1:31316/1‑129 (MQ=255)                           cAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTtgatg                                                                                           <  2:129811/139‑1 (MQ=255)                                    ttCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTgg                                                                                  <  1:563134/139‑1 (MQ=255)                                         gggCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCa                                                                                                      >  2:740745/1‑114 (MQ=255)                                         gggCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCa                                                                                                      <  1:740745/114‑1 (MQ=255)                                              ccAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCAtttttt                                                                        >  2:179724/1‑139 (MQ=255)                                                   gcCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTc                                                                   >  2:729853/1‑139 (MQ=255)                                                         cgcgGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACgaag                                                             >  2:549094/1‑139 (MQ=255)                                                                        tGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                              <  1:8002/139‑1 (MQ=255)                                                                        tGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                              >  2:309033/1‑139 (MQ=255)                                                                                  tCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCa                                    <  1:309033/139‑1 (MQ=255)                                                                                         gggTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCg                             >  2:394362/1‑139 (MQ=255)                                                                                         gggTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCg                             >  2:556302/1‑139 (MQ=255)                                                                                          ggTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGa                            <  1:549094/139‑1 (MQ=255)                                                                                           gTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGaa                           >  2:729519/1‑139 (MQ=255)                                                                                            tCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAAt                          <  1:257481/139‑1 (MQ=255)                                                                                               cgcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgc                       >  1:419005/1‑139 (MQ=255)                                                                                                 ccAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgcag                     <  1:179724/139‑1 (MQ=255)                                                                                                                  cAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTa    <  1:394362/139‑1 (MQ=255)                                                                                                                    gAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGAGAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACg  <  1:420469/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                           ATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACG  >  NC_000913/3400642‑3400895

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