Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	487,802	C→T	L423L (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	487,802	0	C	T	100.0%	39.8 / NA	15	L423L (CTG→TTG)	mscK	potassium dependent, small conductance mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/7);  total (8/7)

ACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATGCACTCTTCCAGAGCGATGCGTTCGTCAACAAACTGGAAGAAGGTCACACCAACGAAGTCAACAGCGAAGTTCACGATGCGTTATTGCAAGTGGTTGATATGCGTCGCGAATTGCTGGATCAACTCAACAAACAGTTGGGTAACCAGCTGATGATGGCCATTAACCTGCAAATCAACCAGCAGCAGTTAATGAGTGTGTCGAAAAACCTGAAATCCATCCTGACTCAGCAAATCTTTTGGGTGAACAGTAACCGTCCAATGGACTGGGACTGGATCAAAGCGTTCCCGCAAAGCCTGAAAGATGAATTTAAGTCGATGAAAATCACGGTGAACTGGCAAAAAGCCTGGCC  >  NC_000913/487594‑488037                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            acaaACGCTGCCCTCGGCGGATGAACTGGAAAACATGAGCATCCGCAGCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGTGTGATGCACTCTTCCAGAGCGATGCGTTCGTCAACAAACGGGAAGAATGCCACCCCGACGAAGGCAACAGCGAAGTTCACGATGCGTGATTGCAAGTGGGTGATATGGGGCGCGAATTGTGGGAGCAACTCAACAAACAGTTGGGTAACCAGCTGATGATgg                                                                                                                                                                                                   <  1:224805/251‑1 (MQ=255)       gCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATGCACTCTTCCAGAGCGATGCGTTCGTCAACAAACTGGAAGAAGGTCACACCAACGAAGTCAACAGCGAAGTTCACGATGCGTTATTGCAAGTGGTTGATATGCGTCGCGAATTGTTGGATCAACTCAACAAACAGTTGGGTAACCAGCTGATGATGGCCAtt                                                                                                                                                                                              >  1:354658/1‑250 (MQ=255)           cccacggcggATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATGCACTCTTCCAGAGCGATGCGTTCGTCAACAAACTGGAAGAAGGTCACACCAACGAAGTCAATAGCGAAGTTCACGATGCGTTATTGCAAGTGGTTGATATGCGTCGCGAATTGTTGGAtcaac                                                                                                                                                                                                                                    <  1:26270/204‑1 (MQ=255)                 gcggATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATGCACTCTTCCAGAGCGATGCGTTCGTCAACAAACTGGAAGAAGGTCACACCAACGAAGTCAACAGCGAAGTTCACGATGCGTTATTGCAAGTGGTTGATATGCGTCGCGAATTGTTGGATCAACTCAACAAACAGTTGGGTAACCAGCTGATGATGGCCATTAACCTGCAAAt                                                                                                                                                                                   >  1:141870/1‑251 (MQ=255)                    naTGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATGCACTCTTCCAGAGCGATGCGTTCGTCAACAAACTGGAAGAAGGTCACACCAACGAAGTCAACAGCGAAGTTCACGATGCGTTATTGCAAGTGGTTGATATGCGTCGCGAATTGTTGGATCAACTCAACAAACAGTTGGGTAACCAGCTGATGATGGCCAt                                                                                                                                                                                               >  1:130245/2‑236 (MQ=255)                                       ccAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATGCACTCTTCCAGAGCGCTGCGTTCGTCAACAAACTGGACGAAGGTCACACCAACGAAGTCAACAGCGAAGTTCACGATGCGTTATTGCAAGTGGTTGATATGCGTCGCGAATTGTTGGATCAACTCAACAAACAGTTGGGTAACCAGc                                                                                                                                                                                                            >  1:66722/1‑204 (MQ=255)                                                          gCGTCTCGAACAGTTTGAAGTTACCCAGCAGCGTGATGCACTCTTCCAGAGCGATGCGTTCGTCAACAAACTGGAAGAAGGTCACACCAACGAAGTCAACAGCGAAGTTCACGATGCGTTATTGCAAGTGGTTGATATGCGTCGCGAATTGTTGGATCAACTCAACAAACAGTTGGGTAACCAGCTGATGATGGCCATTAACCTGCAAATCAACCAGCAGCAGTTAATGAGTGTGTCGAAAAACCTGAAat                                                                                                                                          >  1:128281/1‑251 (MQ=255)                                                                                                                  cGTTCGTCAACAAACTGGAAGAAGGTCACACCAACGAAGGCAACAGCGAAGTTCACGATGCGTTATTGCAAGTGGTTGATATGCGTCGCGAATTGTTGGATCAACTCAACAAACAGTTGGGTAACCAGCTGATGATGGCCATTAACCTGCAAATCAACCAGCAGCAGTTAATGAGTGTGTCGAAAAACCTGAAATCCATCCTTACTCAGCAAATCTTTTGGGTGAACAGTAACCGTCCAATg                                                                                           <  1:317361/242‑1 (MQ=255)                                                                                                                         naacaaACTGGAAGAAGGTCACACCAACGAAGTCAACAGCGAAGTTCACGATGCGTTATTGCAAGTGGTTGATATGCGTCGCGAATTGTTGGATCAACTCAACAAACAGTTGGGTAACCAGCTGATGATGGCCATTAACCTGCAAATCAACCAGCAGCAGTTaa                                                                                                                                                                  >  1:264909/2‑164 (MQ=255)                                                                                                                                        agGTCACACCAACGAAGTCAACAGCGAAGTTCACGATGCGTTATTGCAAGTGGTTGATATGCGTCGCGAATTGTTGGATCAACTCAACAAACAGTTGGGTAACCAGCTGATGATGGCCATTAAcc                                                                                                                                                                                          >  1:14357/1‑125 (MQ=255)                                                                                                                                           nnacnnNNNCGANGTCAACAGCGAAGTTCACGATGCGTTATTGCAAGTGGTTGATATGCGTCGCGAATTGTTGGATCAACTCAACAAACAGTTGGGTAACCAGCTGATGATGGCCATTAACCTGCAAATCAACCAGCAGCAGTTaa                                                                                                                                                                  <  2:264909/144‑1 (MQ=255)                                                                                                                                                                    gatCACGATGCGTTATTGCAAGTGGTTGATATGCGTCGCGAATTGTTGGATCAACTCAACAAACAGTTGGGTAACCAGCTGATGATGGCCATTAACCTGCAAATCAACCAGCAGCAGTTAATGAGTGTGTCGAAAAACCTGAAATCCATCCTGACTCAGCAAATCTTTTGGGTGAACAGTAACCGTCCAATGGACTGGGACTGGATCAAAGCGTTCCCGCAAAGCCTGAAAGATGAATTTAAGTCGATGa                                 >  1:139850/3‑250 (MQ=255)                                                                                                                                                                      tCACGNNGCGTTATTGCAAGTGGTTGATATGCGTCGCGAATTGTTGGATCAACTCAACAAACAGTTGGGTAACCAGCTGATGATGGCCATTAAcc                                                                                                                                                                                          <  2:14357/95‑1 (MQ=255)                                                                                                                                                                          gATGCGTTATTGCAAGTGGTTGATATGCGTCGCGAATTGTTGGATCAACTCAACAAACAGTTGGGTAACCAGCTGATGATGGCCATTAACCTGCAAATCAACCAGCAGCAGTTAATGAGTGTGTCGAAAAACCTGAAATCCATCCTGACTCAGCAAATCttttn                                                                                                                 <  1:132960/164‑2 (MQ=255)                                                                                                                                                                                                  atGCGTCGTGAATTGTTGGATCAACGCAACAAACCGTTGGGTAACCAGGTGATGGTGGCCATTAACCTGCAAATCAACCGGCAGCAGTTAATGAGTGTGTCGAAAAACCTGAAATCGATCCCGACTCAGCAAATCTTTTGGGTGAACAGTAACCGTCCAATGGACTGGGACTGGATCAAAGCGTTCCCGCTAAGCCTGAAAGATGAATTTAAGTCGATGAAAATCACGGTGAACTGGCAAAAAGCCTGGcc  <  1:71038/251‑1 (MQ=255)                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            ACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATGCACTCTTCCAGAGCGATGCGTTCGTCAACAAACTGGAAGAAGGTCACACCAACGAAGTCAACAGCGAAGTTCACGATGCGTTATTGCAAGTGGTTGATATGCGTCGCGAATTGCTGGATCAACTCAACAAACAGTTGGGTAACCAGCTGATGATGGCCATTAACCTGCAAATCAACCAGCAGCAGTTAATGAGTGTGTCGAAAAACCTGAAATCCATCCTGACTCAGCAAATCTTTTGGGTGAACAGTAACCGTCCAATGGACTGGGACTGGATCAAAGCGTTCCCGCAAAGCCTGAAAGATGAATTTAAGTCGATGAAAATCACGGTGAACTGGCAAAAAGCCTGGCC  >  NC_000913/487594‑488037

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