Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	489,602	C→T	L1023L (CTG→TTG)	mscK →	potassium dependent, small conductance mechanosensitive channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	489,602	0	C	T	100.0%	32.1 / NA	12	L1023L (CTG→TTG)	mscK	potassium dependent, small conductance mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (5/7);  total (5/7)

CGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAACCAGCTTGAAGTGCATCTGCACA  >  NC_000913/489408‑489824                                                                                                                                                                                                   |                                                                                                                                                                                                                               cGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTATTGAAGGCGGCGACTGATCACCCAAGGGTGATGCACGAACCAATGCCGGAAGtc                                                                                                                                                                           <  1:86438/248‑1 (MQ=255)                                                  cGATCGCAAAGGAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGCCACTACTACGCGTCTGGTGATCCCTCTCGGCGTGGCCTATGGCTCCGATCTGGAACAAGTGCGTAAAGTGTTATTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACg                                                                                                                                                >  1:180195/1‑226 (MQ=255)                                                                 gATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACGGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGTCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTATTGAAGGCGGCGACTGAGCACCCAATGGTGAAGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTg                                                                                                                                             <  1:173251/214‑1 (MQ=255)                                                                   nnnnnnnnnncANNNNNTTTNTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGTCCTGTGGCTCCGATCTGGAAAAAGTGCGTAAAGTGGTATTGAAGGCGGCGACTGAGCCCCCaac                                                                                                                                                                                                       <  2:149572/144‑2 (MQ=255)                                                                                 gCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAACGTGTTATTGAAAGCGGCGAATGAGCACCCAAGGGTGATGCTCGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGt                                                                                        >  1:246922/1‑251 (MQ=255)                                                                                 gCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCCTCTAGGCGTGGCCTATGGATCCGATCTGGAAAAAGTGCGTAAAGTGTTATTGCAGGCGGCGACTGAGCACCCAAGGGTGATGCTCGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGt                                                                                        >  1:246924/1‑251 (MQ=255)                                                                                            cgggCGTCTGATCAACTTGTCGTTGACTGACACTTCTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTATTGAAggcgg                                                                                                                                                                                                                       <  1:101506/110‑1 (MQ=255)                                                                                                                    gactgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTATTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTg                                                                                                                                                        >  1:39898/1‑152 (MQ=255)                                                                                                                                         ggTGATCCGTCTCGGCGTGGCCTATGGCAGCGATCTGGAAAAAGTGCGTAAAGTGTTATGGAAGTTGGCGACTGAGCACCCAAGGGTGATGCACGAGCCAATGCCGGAAGTCTTCTTTACATCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCg                                                                                                          <  1:147077/177‑1 (MQ=255)                                                                                                                                                           ggCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTATTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCg                                                                                      <  1:21037/179‑1 (MQ=255)                                                                                                                                                                        naTCTGGAAAAAGTGCGTAAAGTGTTATTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGAATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCAGTGAAAACGACATCAACATTGCCTTTAACCAGCTTGAAGTGCATCTGcaca  >  1:50844/2‑250 (MQ=255)                                                                                                                                                                                          aaaGTGTTATTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCGTTACGGCATTTGGTTGGAGCACGTTGTAGCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGAc                                                                                                      <  1:229611/132‑1 (MQ=255)                                                                                                                                                                                                   |                                                                                                                                                                                                                               CGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAACCAGCTTGAAGTGCATCTGCACA  >  NC_000913/489408‑489824

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