Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	580,905	A→G	E24E (GAA→GAG)	nohD →	DLP12 prophage, DNA packaging protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	580,905	0	A	G	100.0%	126.6 / NA	37	E24E (GAA→GAG)	nohD	DLP12 prophage, DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/17);  total (20/17)

CTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAAC  >  NC_000913/580772‑581038                                                                                                                                      |                                                                                                                                      cTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAggg                                                                                                                                  >  2:31212/1‑139 (MQ=255) cTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAggg                                                                                                                                  <  1:850544/139‑1 (MQ=255)    aaaGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAAt                                                                                                                               >  2:992485/1‑139 (MQ=255)        cGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGc                                                                                                                             <  1:214341/137‑1 (MQ=255)        cGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGc                                                                                                                             >  2:214341/1‑137 (MQ=255)         gAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCTGGGAATGCCCg                                                                                                                          >  2:298837/1‑139 (MQ=255)         gAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCg                                                                                                                          >  1:410382/1‑139 (MQ=255)                        tGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTg                                                                                                           >  1:127532/1‑139 (MQ=255)                        tGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATACCCGTTCTGCGAGGCGGTg                                                                                                           >  1:110935/1‑139 (MQ=255)                         gTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTgg                                                                                                          <  2:413630/139‑1 (MQ=37)                          tCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGc                                                                                                         <  1:992485/139‑1 (MQ=255)                                 cTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTa                                                                                                  <  2:410382/139‑1 (MQ=255)                                    tctcTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGc                                                                                                               <  1:1198112/123‑1 (MQ=255)                                    tctcTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGc                                                                                                               >  2:1198112/1‑123 (MQ=255)                                    tctcTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCg                                                                                                              >  1:1194255/1‑124 (MQ=255)                                    tctcTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCg                                                                                                              <  2:1194255/124‑1 (MQ=255)                                      tctGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCGGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAg                                                                                             <  1:617256/139‑1 (MQ=255)                                      tctGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAg                                                                                             <  1:31212/139‑1 (MQ=255)                                                 ccGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGa                                                                                  <  1:453878/139‑1 (MQ=21)                                                           tgCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGc                                                                                                               <  2:1038681/100‑1 (MQ=17)                                                           tgCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGc                                                                                                               >  1:1038681/1‑100 (MQ=17)                                                            gCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGccgccg                                                                       <  2:127532/139‑1 (MQ=11)                                                            gCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATACCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGccgccg                                                                       <  2:110935/139‑1 (MQ=11)                                                                     gTCAACAAAAAGCAGCTGGATGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAAt                                                              >  1:825667/1‑139 (MQ=21)                                                                          caaAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCTGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTAt                                                         <  1:298837/139‑1 (MQ=255)                                                                            aaaaGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGc                                                       >  1:1159279/1‑139 (MQ=25)                                                                            aaaaGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGc                                                       >  1:715574/1‑139 (MQ=21)                                                                             aaaGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGcc                                                      >  2:611694/1‑139 (MQ=21)                                                                                       gctgACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATg                                            >  2:268470/1‑139 (MQ=25)                                                                                                                      ttCAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGAAGCTGAAATTGAGAACGAAAAGCTGCGCCGGGa             <  1:1079516/139‑1 (MQ=25)                                                                                                                        cAGAACTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAgg           >  1:709101/1‑139 (MQ=255)                                                                                                                           aaCTGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTg        <  1:268470/139‑1 (MQ=255)                                                                                                                              tGGCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTgaag     >  2:617992/1‑139 (MQ=38)                                                                                                                                gCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTgaagaa   >  2:607843/1‑139 (MQ=255)                                                                                                                                gCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTgaagaa   >  2:1015608/1‑139 (MQ=37)                                                                                                                                gCAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTgaagaa   >  1:99633/1‑139 (MQ=255)                                                                                                                                 cAGGAGCAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGGCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAAc  <  2:709101/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                      CTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAAC  >  NC_000913/580772‑581038

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