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%	47.8 / NA	16	intergenic (+2346/+397)	ydfD/ynfP	lysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/6);  total (10/6)

TTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTAT  >  NC_000913/1652085‑1652554                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                ttATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCGGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACTCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGATCATATCATTGATGTTCACCTCa                                                                                                                                                                                                                              <  2:120906/250‑1 (MQ=255)              cTGTTTATACTGGAGTTCGCCATGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCACTCGGGTCATTCATCTTATTAAGCCAACAATCGACGCTCTCCGGGGTCAGATGACATTAACGGGACTGAGCAAAGAGCCAATCATTTATGTTCACCCCAGAGAGTATGGccg                                                                                                                                                                                                                 >  1:240678/1‑248 (MQ=255)                      aCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCCGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaa                                                                                                                                                                                                         >  1:69236/1‑250 (MQ=255)                                           cTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGAAAACAATCGGGGCATTCCTATTATTAAACCAGCAATCGACGCCCTCCGGCGACAGATGACATTAACGAGACTGAGGAAAAAGGATATCCTTGTTGGTACCCCCAAAAAGTATGGTAGAAAAGAAAAACAAAAATGCACCTTTGtttt                                                                                                                                                                                   >  2:221570/1‑251 (MQ=255)                                                                     aGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaa                                                                                                                                                         <  2:126442/251‑1 (MQ=255)                                                                                                  ggAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGAc                                                                                                                            <  1:135534/251‑1 (MQ=255)                                                                                                                                 aCCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATaa                                                                                             <  1:253130/251‑1 (MQ=255)                                                                                                                                  ccAAAAACAAAGGCGGGCACCAATCGGGTCATTCATCTTATTAAGCCAACAATCGACGCCATCACGAGTCCGCTGACATTAAAGAGACTGACCCAAGAGCCTATCATTCATGTTCACCCCCGAGAGTATTGCCGAACCGAAAAACCAAAATGCACCTTTGTTTTTCAACCTGAAGTGGCCGCGAGTGTACAAAATTATGGTGACCATTTTACCGTTGAACCAATAAAGCAGATATGGGAAGCAGAGATaa                                                                                             >  2:323116/1‑250 (MQ=255)                                                                                                                                                        aTCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCg                                                                       >  2:86482/1‑250 (MQ=255)                                                                                                                                                            ggTCATTCATCTTATCAAGCCAGCACTCGACGCTCCCCGGAGTCAGCTGACATTAACGCGACTGAGCAAAAAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAACTGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGTTGACCATTTTACCGTTGACTCAATAAGGCAGATGTAGGACGCAGCGATAAAAAGTGCCGGACTCCGACATCGAAAAt                                                                  >  2:144941/1‑251 (MQ=255)                                                                                                                                                                                     aTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCGCCTTTGTTTTTCATCCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGc                                         <  2:6936/251‑1 (MQ=255)                                                                                                                                                                                          cGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctgct                                    >  1:59688/1‑251 (MQ=255)                                                                                                                                                                                                                  aaCGAGACTGAGCAAAGAGCATATCATTGATGGTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCGTTGTTGTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATGATGGTGACCATTTTACCGTTGACGCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGCCAGCTTGTGCTAAcc            <  2:240678/251‑1 (MQ=255)                                                                                                                                                                                                                    cGAGACTGAGCAAAGAGAATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCACCCTGAAGTGTAAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTCTGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCg          >  2:330642/1‑251 (MQ=255)                                                                                                                                                                                                                         cTGAGCAAAGCGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAACAACAAAAATGCACCTTTGTTTTTCCACCTGAAGTGTCAGCGCGAGTCAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCAtt     >  2:223302/1‑251 (MQ=255)                                                                                                                                                                                                                             gCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCGCCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTtat  >  2:217604/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                TTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTAT  >  NC_000913/1652085‑1652554

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