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

TTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATA‑TCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGG  >  NC_000913/1652096‑1652568                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              ttCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGGACCAATCGGGGCATTCATATTATTTAAACAACAATCGACGCCCTCCGGAGTCAGATTACATTAAAGAGA‑CTTAGCAAAGAGCAAAACCATT‑ATGTTCACCTCAGAGAGTAAGGc                                                                                                                                                                                                                                 >  2:16749/1‑250 (MQ=255)          taCTGGAGTTCGCCCTGGTGAGCTGTGGGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAAAGGGTCATTCATCTTATTAAGCCCGCAATCGACGCTCTCCGGAGGCAGATGACATTAACGAGA‑CTGAGCAAAAAGCAAA‑TCATTGATGTTCACCTCAGAGAGTATGGcagaacagaa                                                                                                                                                                                                                        >  1:373794/1‑250 (MQ=255)                           gTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACCCAGGCGGGAACCAATCGGGTCATTCATCTGATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATA‑TCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACCAAAATGCACCttt                                                                                                                                                                                                      <  2:447653/251‑1 (MQ=255)                                    gtgCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATA‑TCATTGATGTTTAACTCAGAGAGTATGGCAGAACAGAAAAAAAAAAATGCACCTTTGtttttta                                                                                                                                                                                              >  1:456740/1‑248 (MQ=255)                                       cACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATA‑TCATTGATGTTCACCTCATAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCt                                                                                                                                                                                          <  2:456740/251‑1 (MQ=255)                                                  gaggaCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGATTCACGGTACCAAAAACACAGGCGGGAACCCATCGGGTCCTTCTTCTTATTTAGCCAGCCATCGACGCCCTCCGGAGTCAGATAACATTAAAGAGA‑CTGGGCAAAGAGGATA‑TCATTGATGTTCACCTCAGAGAGCATGGGAGAACAGAAAAACAAAAATGCAACTTTGTTTTTCAACCTGAAGTGACAgg                                                                                                                                                                               >  1:448529/1‑250 (MQ=255)                                                        aTAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGGACCAATCGGGTCATTCATATTATTAAGCCACCAAACAACGCCCTTCCGAGGCCGATGCAATTAACG‑GAGCAGAGCAAGGAGCAAC‑TCCTTGGTGTTTACCTCAAGGAGAATGGGAGAAAAGAAAAAAAAAAAAGAAACTTTGCTTTTCAAACTGAAATGTCAgaggg                                                                                                                                                                            >  1:216400/1‑244 (MQ=255)                                                                                                 aGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATA‑TCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGc                                                                                                                                <  1:273631/251‑1 (MQ=255)                                                                                                     aaaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCCATCGGGTCATTCATCTTATTAAGCCAGCAATAGACGCTATCCGGAGTCAGATGACATTAACGCGA‑CTGAGCAAAGAGCATA‑TCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAAAAAAAATGCACCTTTGTTTTTTAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAACATTTTACCGTTGTCTCAATAAGGc                                                                                                                                >  1:294045/1‑247 (MQ=255)                                                                                                         aCCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATA‑TCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTgg                                                                                                                        <  1:71275/251‑1 (MQ=255)                                                                                                                          aaaaCACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATA‑TCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACg                                                                                                       <  1:16749/251‑1 (MQ=255)                                                                                                                                                             tattaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATA‑TCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGAGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCg                                                                     >  2:22278/1‑250 (MQ=255)                                                                                                                                                                                  tctcCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATA‑TCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGt                                               >  1:22170/1‑251 (MQ=255)                                                                                                                                                                                            cAGATGACATTAACGAGA‑CTGATCAAAGAGCATA‑TCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAg                                      <  2:373794/250‑1 (MQ=255)                                                                                                                                                                                                                      aaGAGCATA‑TCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGAtgt                                                                                                                          >  1:437320/1‑141 (MQ=255)                                                                                                                                                                                                                      aaGAGCATA‑TCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGAtgt                                                                                                                          <  2:437320/141‑1 (MQ=255)                                                                                                                                                                                                                                 cATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCAGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAACAGGTGCTAACCCGGCATTTATAGCAAACCAGATgg  >  2:128400/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              TTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATA‑TCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGG  >  NC_000913/1652096‑1652568

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