Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,652,331	T→C	intergenic (+2346/+397)	ydfD → / ← ynfP	lysis protein/protein YnfP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,652,331	0	T	C	100.0%	19.5 / NA	9	intergenic (+2346/+397)	ydfD/ynfP	lysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (6/3);  total (6/3)

CTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGG‑TCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAA‑AATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACC  >  NC_000913/1652136‑1652562                                                                                                                                                                                                     |                                                                                                                                                                                                                                         cTTGGATGGGAGGACATAGATCTGAAAACTGGAACAATGCTGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGCGACCAATCGGGGTC‑TTCATCTTATTAACCCACCCCTCGACGCCCTCCGCAGTCACACGACACTAATGCGACCTAGCAAATAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGGACAGAAAAACAAA‑AATGCACCTTTGTTTTTCAACCTg                                                                                                                                                                                   >  1:9873/1‑250 (MQ=255)                                             cAGGAAAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGG‑TCATTCATCTTATTAAGCCAGCAATCGACGCTCCCCGGAGTCAGATGACATGAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAA‑AATGCACCTTCGTTTTGCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGttn                                                                                                                                      <  1:135798/251‑2 (MQ=255)                                                               aGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGG‑TCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCGCAGAGAGTATGGcagn                                                                                                                                                                                                                         <  1:133813/151‑2 (MQ=255)                                                                                                                                  aaTCGACGCTCTCCGGAGTCAGATGACATTAACGAGAATGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGACCCGAAAAACAAC‑AATGCACCTTTGTTTTTCAACCTGAAGAGTCAGCGAGAGTAAACAATTATGTTGACCATTTTACCGTTGACTCAATAAGGCAGAtgt                                                                                                                    >  1:114989/1‑185 (MQ=255)                                                                                                                                   aTCGACGCTATCCGGAGTCAGATGACATTAACGAGACTGCGCAAAGAGCATATAATTGATGTTCACCTCCGAGAGTCTGGCAGACCAGCAAA‑CAAACAATGCACCTTTGTTTTTCACCCTGAAGGGTCAGCGAGAGCAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCAATCGAAAATCATATCAGTCGCGACATACTTATGc                                                 >  1:222682/1‑251 (MQ=255)                                                                                                                                                       gATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAA‑AATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACNATTNNNNNGtnn                                                                                                                                      >  2:133296/1‑144 (MQ=255)                                                                                                                                                       gATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAA‑AATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGcn                                                                                                         <  1:133296/175‑2 (MQ=255)                                                                                                                                                                  cGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAA‑AATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAcc                                                                                                                                                  >  1:58198/1‑123 (MQ=255)                                                                                                                                                                                  gCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAA‑AATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAAcc  >  1:103726/1‑251 (MQ=255)                                                                                                                                                                                                     |                                                                                                                                                                                                                                         CTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGG‑TCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAA‑AATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACC  >  NC_000913/1652136‑1652562

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