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%	116.0 / NA	27	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 . (13/14);  total (13/14)

GGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGG  >  NC_000913/281530‑281788                                                                                                                                   |                                                                                                                                 ggCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGgtg                                                                                                                         <  2:498794/139‑1 (MQ=255)        aaaaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggca                                                                                                                               >  2:525286/1‑125 (MQ=255)        aaaaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggca                                                                                                                               <  1:525286/127‑3 (MQ=255)                        agGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAgcg                                                                                                  >  1:806588/1‑139 (MQ=255)                          tcGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCt                                                                                                <  1:185167/138‑1 (MQ=255)                           cGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCtt                                                                                               <  2:449460/139‑1 (MQ=255)                             gTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggc                                                                                                                                <  1:401847/105‑2 (MQ=255)                             gTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggc                                                                                                                                >  2:401847/1‑104 (MQ=255)                              ttCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAAt                                                                                            >  1:404315/1‑139 (MQ=255)                                     cgtcgtGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCgg                                                                                     <  1:521057/139‑1 (MQ=255)                                      gtcgtGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCggg                                                                                    <  1:629377/139‑1 (MQ=255)                                        cgtGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCt                                                                                  <  2:697105/139‑1 (MQ=255)                                        cgtGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCt                                                                                  >  1:697105/1‑139 (MQ=255)                                                   cGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCgg                                                                                     >  1:402928/1‑125 (MQ=255)                                                   cGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCgg                                                                                     <  2:402928/125‑1 (MQ=255)                                                     gcacGTGGCGGGAATGCCCGGTGAAGCGCCCGGCCGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAACGTTCAt                                                                     >  1:423533/4‑139 (MQ=255)                                                             cGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCtt                                                                                         >  1:348273/1‑111 (MQ=255)                                                             cGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCtt                                                                                         <  2:348273/111‑1 (MQ=255)                                                               ggAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCg                                                                          >  1:65260/1‑124 (MQ=255)                                                               ggAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCg                                                                          <  2:65260/124‑1 (MQ=255)                                                                   tGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGcc                                                       <  2:404315/139‑1 (MQ=255)                                                                     cccGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGcc                                                       >  1:6791/1‑137 (MQ=255)                                                                     cccGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGcc                                                       <  2:6791/137‑1 (MQ=255)                                                                      ccGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATa                                                    >  1:634461/1‑139 (MQ=255)                                                                                    ggCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggcat                                                                                                                              >  2:381/1‑49 (MQ=255)                                                                                    ggCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggcat                                                                                                                              <  1:381/52‑4 (MQ=255)                                                                                               aTGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                                   >  1:805042/1‑99 (MQ=255)                                                                                               aTGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                                   <  2:805042/99‑1 (MQ=255)                                                                                                tGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGa                          >  2:590959/1‑139 (MQ=255)                                                                                                               ttATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa           <  2:423533/139‑1 (MQ=255)                                                                                                                  tGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTc                                                              >  2:360802/1‑85 (MQ=255)                                                                                                                  tGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTc                                                              <  1:360802/85‑1 (MQ=255)                                                                                                                        gTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCgg  >  2:282793/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                 GGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGG  >  NC_000913/281530‑281788

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