Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000964	557,870	(T)5→6	intergenic (+421/‑3)	yddT → / → ydzN	putative exported phage protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000964	557,870	1	.	T	100.0%	51.8 / NA	14	intergenic (+421/‑3)	yddT/ydzN	putative exported phage protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (10/4);  total (10/4)

AGCCGGACGGTTGATGCCAAGCTCAAAAGTTTATATGGAAGTCAAAGCATAAAAAAGGGTGTTACAACGTGTAGTCTTCTTTTTATTTTGCCAACGTAACAAGTGTGAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTT‑AATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTAATTGTGTTGTGGAATATTGTTGGACTTTTTAGCAAAGAATCAGTAAAACAAGCACAGCAAGCAAATTAAAAGTCATAAAGTTGG  >  NC_000964/557624‑558073                                                                                                                                                                                                                                                        |                                                                                                                                                                                                            aGCCGGACGGTTGATGCCAAGCTCAAAAGTTTATATGGAAGTCAAAGCATAAAAAAGGGTGTTACAACGTGTAGTCTTCTTTTTATTTTGCCAACGTAACAAGTGTGAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGttttttaat                                                                                                                                                                                                          <  2:325870/251‑5 (MQ=255)             gATGCCAAGCTCAAAAGTTTATATGGAAGTCAAAGCATAAAAAAGGGTGTTACAACGTGTAGTCTTCTTTTTATTTTGCCAACGTAACAAGTGTGAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGt                                                                                                                                                                                               >  1:268038/1‑250 (MQ=255)                             gTTTATATGGAAGTCAAAGCATAAAAAAGGGTGTTACAACGTGTAGTCTTCTTTTTATTTTGCCAACGTAACAAGTGTGAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGtttt                                                                                                                                                                              >  1:200842/1‑251 (MQ=255)                                                                                        ttGCCAACGTAACAAGTGTGAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTc                                                                                                                    >  2:116533/1‑250 (MQ=255)                                                                                         tGCCAACGTAACAAGTGTGAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTc                                                                                                                    >  1:95718/1‑249 (MQ=255)                                                                                                               gTAATGATTTGTAGGATAAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTAAATAAAATAATGGTACTATAATAAATGTATACAGATTGTATGGTGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAAACAAGAAGTCAGTACTACTTGTTTAATAATAgtt                                                                                            >  2:51128/1‑249 (MQ=255)                                                                                                                                   ttAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTAATTGTGTTGTGGa                                                                         >  1:224521/1‑250 (MQ=255)                                                                                                                                          aGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGGCAGTACTACTCGTTTAATAATAGCTTCTGTAATTGTGTTGTGGAATAttgtt                                                                 <  2:249689/251‑1 (MQ=255)                                                                                                                                                      aGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTAATTGTGTTGTGGAATATTGTTGGACTTTTTAg                                                      <  1:169412/250‑1 (MQ=255)                                                                                                                                                            tGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTAATTGTGTTGTGGAATATTGTTGGACTTTTTAGCAAAGa                                                >  2:220203/1‑250 (MQ=255)                                                                                                                                                                        tatatTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTAATTGTGTTGTGGAATATTGTTGGACTTTTTAGCAAAGAATCAGTAAAACa                                    >  2:192893/1‑250 (MQ=255)                                                                                                                                                                         atatTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTAATTGTGTTGTGGAATATTGTTGGACTTTTTAGCAAAGAATCAGTAAAACaa                                   <  2:95718/250‑1 (MQ=255)                                                                                                                                                                                                       aaaTGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTAAttgtg                                                                                <  1:166072/175‑1 (MQ=255)                                                                                                                                                                                                       aaaTGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTAAttgtg                                                                                >  2:166072/1‑175 (MQ=255)                                                                                                                                                                                                         aTGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTAATTGTGTTGTGGAATATTGTTGGACTTTTTAGCAAAGAATCAGTAAAACAAGCACAGCAAGCAAATTAAAAGTCATAAAGTTgg  >  1:53229/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                            AGCCGGACGGTTGATGCCAAGCTCAAAAGTTTATATGGAAGTCAAAGCATAAAAAAGGGTGTTACAACGTGTAGTCTTCTTTTTATTTTGCCAACGTAACAAGTGTGAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTT‑AATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTAATTGTGTTGTGGAATATTGTTGGACTTTTTAGCAAAGAATCAGTAAAACAAGCACAGCAAGCAAATTAAAAGTCATAAAGTTGG  >  NC_000964/557624‑558073

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