Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,168,644	T→G	K63N (AAA→AAC)	comR ←	DNA‑binding transcriptional repressor ComR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,168,644	0	T	G	100.0%	100.4 / NA	28	K63N (AAA→AAC)	comR	DNA‑binding transcriptional repressor ComR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (16/12);  total (16/12)

CAGTTGGCGATGGCAGCAAAATAGTCAGCCAGCGCAGACTCCACGCTTTTCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCTTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTCGACGAGGTCCGCAAGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCC  >  NC_000913/1168513‑1168770                                                                                                                                    |                                                                                                                               cAGTTGGCGATGGCAGCAAAATAGTCAGCCAGCGCAGACTCCACGCTTTTCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGt                                                                                                                        <  1:13222/140‑1 (MQ=255) cAGTTGGCGATGGCAGCAAAATAGTCAGCCAGCGCAGACTCCACGCTTTTCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGt                                                                                                                        <  2:7785/140‑1 (MQ=255)     tGGCGATGGCAGCAAAATAGTCAGCCAGCGCAGACTCCACGCTTTTCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAAt                                                                                                                    >  2:181870/1‑140 (MQ=255)                      tAGTCAGCCAGCGCAGACTCCACGCTTTTCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTggg                                                                                                   <  2:461162/140‑1 (MQ=255)                      tAGTCAGCCAGCGCAGACTCCACGCTTTTCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTggg                                                                                                   <  2:210218/140‑1 (MQ=255)                           agccagcGCAGACTCCACGCTTTTCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCtt                                                                                              >  2:28262/1‑140 (MQ=255)                           agccagcGCAGACTCCACGCTTTTCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCtt                                                                                              >  1:39709/1‑140 (MQ=255)                           agccagcGCAGACTCCACGCTTTTCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCtt                                                                                              >  2:168125/1‑140 (MQ=255)                              cagcGCAGACTCCACGCTTTTCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTgg                                                                                           >  1:291636/1‑140 (MQ=255)                                gcgcAGACTCCACGCTTTTCTATTCACAAAACAGCTGTGCTTCAGTCTTAGCGGAAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGCGGGCGCTTTGGCt                                                                                         >  2:100620/1‑140 (MQ=255)                                 cgcAGACTCCACGCTTTTCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTc                                                                                        >  1:297230/1‑140 (MQ=255)                                      aCTCCACGCTTTTCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTc                                                                                   >  2:18396/1‑140 (MQ=255)                                             gCTTTTCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTcga                                                                            >  2:141992/1‑140 (MQ=255)                                                tttCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACTGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTcgacga                                                                         >  2:44396/1‑140 (MQ=255)                                                                aGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTCGACGAGGTCCGCAAGAGAAGt                                                         <  1:198130/140‑1 (MQ=255)                                                                aGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTCGACGAGGTCCGCAAGAGAAGt                                                         >  2:55692/1‑140 (MQ=255)                                                                aGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTCGACGAGGTCCGCAAGAGAAGt                                                         <  1:141992/140‑1 (MQ=255)                                                                         tCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTCGACGAGGTCCGCAAGAGAAGTCGCTTCATa                                                <  1:18396/140‑1 (MQ=255)                                                                                 aGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTCGACGAGGTCCGCAAGAGAAGTCGCTTCATAACCGTGtt                                        >  2:77999/1‑140 (MQ=255)                                                                                 aGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTCGACGAGGTCCGCAAGAGAAGTCGCTTCATAACCGTGtt                                        >  2:160115/1‑140 (MQ=255)                                                                                 aGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTCGACGAGGTCCGCAAGAGAAGTCGCTTCATAACCGTGtt                                        >  1:229005/1‑140 (MQ=255)                                                                                   cGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTCGACGAGGTCCGCAAGAGAAGTCGCTTCATAACCGTGTTGc                                      >  1:384048/1‑140 (MQ=255)                                                                                          aCGATCGATATAGCGGTCGAGAACTGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTCGACGAGGTCCGCAAGAGAAGTCGCTTCATAACCGTGTTGCCAGAACa                               <  1:44396/140‑1 (MQ=255)                                                                                               cgatATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTCGACGAGGTCCGCAAGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTc                          <  1:100620/140‑1 (MQ=255)                                                                                                gatATAGCGGTCGAGAACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTCGACGAGGTCCGCAAGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCa                         <  1:359192/140‑1 (MQ=255)                                                                                                            gagaACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTCGACGAGGTCCGCAAGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCaa             >  1:357605/1‑140 (MQ=255)                                                                                                            gagaACGGCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTAGCTTCGACGAGGTCCGCAAGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCaa             >  1:41739/1‑140 (MQ=255)                                                                                                                  ggCACGAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTAGCTTCGACGAGGTCCGCAAGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCgg       <  2:41739/140‑1 (MQ=255)                                                                                                                       gAAATAACCCCTCGTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTCGACGAGGTCCGCAAGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATcc  <  1:28262/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                               CAGTTGGCGATGGCAGCAAAATAGTCAGCCAGCGCAGACTCCACGCTTTTCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAAATAACCCCTCTTTGTTGGTAAATTCCGCGTATAGCGTGGGCGCTTTGGCTCCGGTCGCTTCGACGAGGTCCGCAAGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCC  >  NC_000913/1168513‑1168770

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