Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,387,080	A→G	V764A (GTC→GCC)	mscM ←	miniconductance mechanosensitive channel MscM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,387,080	0	A	G	100.0%	67.3 / NA	19	V764A (GTC→GCC)	mscM	miniconductance mechanosensitive channel MscM
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/11);  total (8/11)

AAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGC  >  NC_000913/4386951‑4387215                                                                                                                                  |                                                                                                                                        aaGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGGCTTCAATTGc                                                                                                                               >  2:41563/1‑140 (MQ=255)  aGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGGCTTCAATTGCt                                                                                                                              <  1:41563/140‑1 (MQ=255)  aGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGGCTTCAATTGCt                                                                                                                              <  2:254515/140‑1 (MQ=255)          gCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGGCTTCAATTGCTCCTTCCgg                                                                                                                      <  1:360635/140‑1 (MQ=255)            ggAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGGCTTCAATTGCTCCTTCCGGGc                                                                                                                    <  1:88780/140‑1 (MQ=255)            ggAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGGCTTCAATTGCTCCTTCCGGGc                                                                                                                    <  1:152096/140‑1 (MQ=255)                              caGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGGCTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCtt                                                                                                  <  2:338842/140‑1 (MQ=255)                                                  gggCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGGCTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGt                                                                              <  2:152858/140‑1 (MQ=255)                                                   ggCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGGCTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGtt                                                                             >  2:295928/1‑140 (MQ=255)                                                     cGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGGCTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTgc                                                                           <  2:48105/140‑1 (MQ=255)                                                                         cgcACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGGCTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTa                                                                       >  1:172077/1‑124 (MQ=255)                                                                         cgcACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGGCTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTa                                                                       <  2:172077/124‑1 (MQ=255)                                                                         cgcACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGGCTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCg                                                       >  2:183103/1‑140 (MQ=255)                                                                                           tGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGGCTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCaaa                                     >  1:82716/1‑140 (MQ=255)                                                                                             cgcACTGATGGCATCGAGATCGACTTCGCTTTCATCGGCTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACg                                   <  2:159348/140‑1 (MQ=255)                                                                                                   gATGGCATCGAGATCGACTTCGCTTTCATCGGCTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGcc                             >  2:408792/1‑140 (MQ=255)                                                                                                            gagaTCGACTTCGCTTTCATCGGCTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCt                    >  2:12176/1‑140 (MQ=255)                                                                                                                          tttCATCGGCTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAg      <  2:82716/140‑1 (MQ=255)                                                                                                                              aTCGGCTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAgcggc  >  1:357086/1‑140 (MQ=255)                                                                                                                                  |                                                                                                                                        AAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGC  >  NC_000913/4386951‑4387215

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