Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	284,869	C→T	P557S (CCG→TCG)	yagF →	CP4‑6 prophage, dehydratase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	284,869	0	C	T	100.0%	59.7 / NA	20	P557S (CCG→TCG)	yagF	CP4‑6 prophage, dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/10);  total (10/10)

ACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGCCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGAACTTCATCGGCACCGCGGACAACCCGCTGACGCCGGAAGAGGGCGCG  >  NC_000913/284738‑284988                                                                                                                                    |                                                                                                                        aCCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTgg                                                                                                                  <  2:376999/139‑1 (MQ=255)  ccAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGc                                                                                                                 >  1:587996/1‑139 (MQ=255)  ccAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGc                                                                                                                 <  2:148155/139‑1 (MQ=255)  ccAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGc                                                                                                                 >  2:208067/1‑139 (MQ=255)                cgcTAAAGCATATCTCGTGGGGCATGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATaa                                                                                                   >  1:543228/1‑139 (MQ=255)                              tcGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCgaga                                                                                     <  2:662417/139‑1 (MQ=255)                              tcGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCgaga                                                                                     <  1:778566/139‑1 (MQ=255)                               cGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGAt                                                                                    >  2:587688/1‑139 (MQ=255)                                        aGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGa                                                                           >  2:628570/1‑139 (MQ=255)                                                        caCCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACt                                                           >  2:383329/1‑139 (MQ=255)                                                              tGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGcc                                                                                >  2:795126/1‑112 (MQ=255)                                                              tGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGcc                                                                                <  1:795126/112‑1 (MQ=255)                                                                    cTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGAAc                                               <  1:325592/139‑1 (MQ=255)                                                                                 tCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGa                                                                           >  1:383813/1‑98 (MQ=255)                                                                                 tCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGa                                                                           <  2:383813/98‑1 (MQ=255)                                                                                 tCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGAACTTCATCGGCACcg                                  >  2:582284/1‑139 (MQ=255)                                                                                              gCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGAACTTCATCGGCACCGCGGACAACCCGCt                     <  1:337428/139‑1 (MQ=255)                                                                                              gCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGAACTTCATCGGCACCGCGGACAACCCGCt                     <  2:742630/139‑1 (MQ=255)                                                                                              gCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGAACTTCATCGGCACCGCAGACAACCCGCt                     <  1:162458/139‑1 (MQ=255)                                                                                                                 ggAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGAACTTCATCGGCACCGCGGACAACCCGCTGACGCCGGAAGAGGgcgcg  >  1:373703/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                        ACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGCCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGAACTTCATCGGCACCGCGGACAACCCGCTGACGCCGGAAGAGGGCGCG  >  NC_000913/284738‑284988

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