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%	133.4 / NA	38	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 (19/19);  total (19/19)

CGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCC  >  NC_000913/3400646‑3400899                                                                                                                                |                                                                                                                               cGGCAGCCATCGGTTCTTCATTCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCt                                                                                                                     >  2:472633/1‑139 (MQ=255)  ggCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTg                                                                                                                    <  1:682423/139‑1 (MQ=255)                  tcAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGtt                                                                                                    <  2:308024/139‑1 (MQ=255)                          gAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGATtgatgacg                                                                                            >  2:671470/1‑137 (MQ=255)                           aaGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCtgt                                                                                                             <  1:1022527/121‑1 (MQ=255)                           aaGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCtgt                                                                                                             >  2:1022527/1‑121 (MQ=255)                                     gggCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCAt                                                                                 >  2:263083/1‑139 (MQ=255)                                              cgcCCTGCGCGGACTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCAAACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGt                                                                        <  1:472633/139‑1 (MQ=255)                                                 ccTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACtt                                                                                                       <  1:759066/105‑1 (MQ=255)                                                 ccTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACtt                                                                                                       >  2:759066/1‑105 (MQ=255)                                                 ccTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCAc                                                                     >  1:508755/1‑139 (MQ=255)                                                                    tGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGtt                                                                                                    <  2:755811/89‑1 (MQ=255)                                                                    tGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGtt                                                                                                    >  1:755811/1‑89 (MQ=255)                                                                      cgcgGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACtt                                                                                                       >  1:565475/1‑84 (MQ=255)                                                                      cgcgGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACtt                                                                                                       <  2:565475/84‑1 (MQ=255)                                                                                   cTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTggg                                   <  1:261698/139‑1 (MQ=255)                                                                                     gggTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGcc                                                >  2:690800/1‑124 (MQ=255)                                                                                     gggTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGcc                                                <  1:690800/124‑1 (MQ=255)                                                                                       gTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTtgatg                                                                                               <  1:336450/75‑1 (MQ=255)                                                                                       gTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTtgatg                                                                                               >  2:336450/1‑75 (MQ=255)                                                                                           cgcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgc                           >  1:340326/1‑139 (MQ=255)                                                                                           cgcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgc                           >  1:578721/1‑139 (MQ=255)                                                                                            gcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTtgatg                                                                                               <  2:993847/70‑1 (MQ=255)                                                                                            gcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTtgatg                                                                                               >  1:993847/1‑70 (MQ=255)                                                                                             ccAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgcag                         <  2:340326/139‑1 (MQ=255)                                                                                                    gCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgca                          >  2:625597/1‑131 (MQ=255)                                                                                                    gCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgca                          <  1:625597/131‑1 (MQ=255)                                                                                                            aCCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCtgtg                                                                                                            >  2:399454/1‑41 (MQ=255)                                                                                                            aCCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCtgtg                                                                                                            <  1:399454/41‑1 (MQ=255)                                                                                                               aGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTAc       <  2:508755/139‑1 (MQ=255)                                                                                                                     cgcgGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                  >  1:727099/1‑90 (MQ=255)                                                                                                                     cgcgGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                  <  2:727099/90‑1 (MQ=255)                                                                                                                      gcgGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTcggt                                                                        <  2:1029453/67‑4 (MQ=255)                                                                                                                      gcgGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTcggt                                                                        >  1:1029453/1‑64 (MQ=255)                                                                                                                      gcgGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGccc  <  2:279079/137‑1 (MQ=255)                                                                                                                      gcgGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGccc  >  1:279079/1‑137 (MQ=255)                                                                                                                        gggCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                  >  1:840270/1‑87 (MQ=255)                                                                                                                        gggCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                  <  2:840270/87‑1 (MQ=255)                                                                                                                                |                                                                                                                               CGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCC  >  NC_000913/3400646‑3400899

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: ‑