Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	281,662	Δ1 bp	coding (322/1155 nt)	yagA ←	CP4‑6 prophage, putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	281,660	0	G	.	100.0%	83.3 / NA	21	coding (324/1155 nt)	yagA	CP4‑6 prophage, putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (10/11);  total (10/11)

AAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAG  >  NC_000913/281525‑281768                                                                                                                                        |                                                                                                             aaaGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCTGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAGGGCggca                                                                                                           <  1:798211/139‑3 (MQ=255)       gCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGgtgt                                                                                                    >  2:807258/1‑139 (MQ=255)        cccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGtgtg                                                                                                   >  1:1045241/1‑139 (MQ=255)              aaaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTg                                                                                             <  1:1183288/139‑1 (MQ=255)              aaaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTg                                                                                             <  1:240699/139‑1 (MQ=255)                     tCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCg                                                                                      >  2:559678/1‑139 (MQ=255)                        gCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGc                                                                                   <  2:281898/139‑1 (MQ=255)                         ccAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGcc                                                                                  <  2:437808/139‑1 (MQ=255)                                cGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GACATGGTGTGCCCCTGGTCCTCGAGCCAGCGCtt                                                                           >  1:238515/1‑139 (MQ=255)                                 ggTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTa                                                                          >  1:105032/1‑139 (MQ=255)                                  gTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTaa                                                                         >  2:932617/1‑139 (MQ=255)                                                     aaCCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCg                                                      >  1:1180873/1‑139 (MQ=255)                                                                  cGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCa                                         >  2:642217/1‑139 (MQ=255)                                                                                       ccGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGcc                    <  1:932617/139‑1 (MQ=255)                                                                                          gcagcaGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg                 <  2:105032/139‑1 (MQ=255)                                                                                           cagcagGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCgg                                                                                                             >  2:967062/1‑47 (MQ=255)                                                                                           cagcagGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCgg                                                                                                             <  1:967062/47‑1 (MQ=255)                                                                                                ggCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATgtcg           >  1:295153/1‑139 (MQ=255)                                                                                                        cGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGa   <  1:190519/139‑1 (MQ=255)                                                                                                         gggCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg  <  2:1061355/139‑1 (MQ=255)                                                                                                         gggCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg  <  1:559678/139‑1 (MQ=255)                                                                                                         gggCCATCAGGTTATGGACGGTGCTGAAGGC‑GACATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg  <  2:238515/139‑1 (MQ=255)                                                                                                                             gTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg  >  2:453336/1‑119 (MQ=255)                                                                                                                             gTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg  <  1:453336/119‑1 (MQ=255)                                                                                                                                        |                                                                                                             AAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAG  >  NC_000913/281525‑281768

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