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%	79.6 / NA	20	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/10);  total (10/10)

GCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGA  >  NC_000913/281531‑281789                                                                                                                                  |                                                                                                                                  gCCCTTAAAATCCATCTGCCAGAGGCGGTTCGACGCTTCGTGTTCGAACCGGCCCGTGGCGGTAATGCCCGGTGAAGCGCCCGGCAGCACGCCATGGCGGGCCATCAGGATATGGTCGGTGCTGAAGGC‑GGCATGgtgt                                                                                                                         <  1:571699/139‑1 (MQ=255)  cccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGtgtg                                                                                                                        >  1:698262/1‑139 (MQ=255)  cccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCACGGtgtg                                                                                                                        >  2:692882/1‑139 (MQ=255)    cTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGcc                                                                                                                      <  1:310771/139‑1 (MQ=255)                  gCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGc                                                                                                        <  1:299420/139‑1 (MQ=255)                             ttCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAAt                                                                                             >  1:247786/1‑139 (MQ=255)                             ttCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAAt                                                                                             >  2:313424/1‑139 (MQ=255)                                          gtTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTccc                                                                                <  1:442171/139‑1 (MQ=255)                                                 ccGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGtt                                                                         <  1:735096/139‑1 (MQ=255)                                                                                      agcagGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGAt                                    <  2:48530/139‑1 (MQ=255)                                                                                           gCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATgtcgt                               <  2:698262/139‑1 (MQ=255)                                                                                                   gggCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg                       <  2:247786/139‑1 (MQ=255)                                                                                                   gggCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg                       <  2:387969/139‑1 (MQ=255)                                                                                                    ggCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGc                      >  1:265486/1‑139 (MQ=255)                                                                                                     gCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg                                      >  2:302970/1‑122 (MQ=255)                                                                                                     gCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg                                      <  1:302970/122‑1 (MQ=255)                                                                                                     gCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCg                     >  1:49115/1‑139 (MQ=255)                                                                                                         tCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTc                 >  2:628244/1‑139 (MQ=255)                                                                                                                       gTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCgg   >  2:189946/1‑139 (MQ=255)                                                                                                                       gTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCgg   >  1:497155/1‑139 (MQ=255)                                                                                                                        tGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGa  >  2:274793/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                  GCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGA  >  NC_000913/281531‑281789

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