Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,199,503	(A)5→4	coding (809/894 nt)	lit →	e14 prophage; cell death peptidase Lit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,199,499	0	A	.	100.0%	64.9 / NA	15	coding (805/894 nt)	lit	e14 prophage; cell death peptidase Lit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (7/8);  total (7/8)

TATTTCATGCTACCCTGTCGGAAATGAAGAGTTGATTGAAGCTCTATGTACAGTGATGCTTCAATATCTTTTCCATGGCAAAAATATCAATGTGAATCTAGATGGGGAGTCCTTTTCATCGATTTTAGGTGATCTTCTCTGTGATATTTCA  >  NC_000913/1199420‑1199570                                                                                |                                                                        tatTTCATGCTACCCTGTCGGAAATGAAGAGTTGATTGAAGCTCTATGTACAGTGATGCTTCAATATCTTTTCCATGGCa                                                                         >  2:217524/1‑80 (MQ=255)       aTGCTACCCTGTCGGAAATGAAGAGTTGATTGAAGCTCTATGTACAGTGATGCTTCAATATCTTTTCCATGGCaaaatat                                                                   >  1:217541/1‑77 (MQ=255)       aTGCTACCCTGTCGGAAATGAAGAGTTGATTGAAGCTCTATGTACAGTGATGCTTCAATATCTTTTCCATGGCaaaatat                                                                   <  2:217511/80‑4 (MQ=255)        tGCTACCCTGTCGGAAATGAAGAGTTGATTGAAGCTCTATGTACAGTGATGCTTCAATATCTTTTCCATGGCaaaatatc                                                                  >  1:217531/1‑76 (MQ=255)               cTGTCGGAAATGAAGAGTTGATTGAAGCTCTATGTACAGTGATGCTTCAATATCTTTTCCATGGC‑AAAATATCAATGTGaa                                                         <  2:217512/81‑1 (MQ=255)                 gTCGGAAATGAAGAGTTGATTGAAGCTCTATGTACAGTGATGCTTCAATATCTTTTCCATGGC‑AAAATATCAATGTGAATc                                                       >  1:217529/1‑81 (MQ=255)                  tCGGAAATGAAGAGTTGATTGAAGCTCTATGTACAGTGATGCTTCAATATCTTTTCCATGGC‑AAAATATCAATGTGAATCt                                                      <  2:217513/81‑1 (MQ=255)                     gAAATGAAGAGTTGATTGAAGCTCTATGTACAGTGATGCTTCAATATCTTTTCCATGGC‑AAAATATCAATGTGAATCTAGa                                                   <  1:706936/81‑1 (MQ=255)                           aaGAGTTGATTGAAGCTCTATGTACAGTGATGCTTCAATATCTTTTCCATGGC‑AAAATATCAATGTGAATCTAGATGGGGa                                             <  1:217514/81‑1 (MQ=255)                             gagTTGATTGAAGCTCTATGTACAGTGATGCTTCAATATCTTTTCCATGGC‑AAAATATCAATGTGAATCTAGATGGGGAGt                                           <  2:217515/81‑1 (MQ=255)                              agTTGATTGAAGCTCTATGTACAGTGATGCTTCAATATCTTTTCCATGGC‑AAAATATCAATGTGAATCTAGATGGGGAGTc                                          <  1:217516/81‑1 (MQ=255)                               gTTGATTGAAGCTCTATGTACAGTGATGCTTCAATATCTTTTCCATGGC‑AAAATATCAATGTGAATCTAGATGGGGAGTcc                                         >  1:217535/1‑81 (MQ=255)                                               tGTACAGTGATGCTTCAATATCTTTTCCATGGC‑AAAATATCAATGTGAATCTAGATGGGGAGTCCTTTTCATCGATTTTAg                         >  1:217533/1‑81 (MQ=255)                                                   cAGTGATGCTTCAATATCTTTTCCATGGC‑AAAATATCAATGTGAATCTAGATGGGGAGTCCTTTTCATCGATTTTAGGTg                      <  2:217517/80‑1 (MQ=255)                                                    aGTGATGCTTCAATATCTTTTCCATGGC‑AAAATATCAATGTGAATCTAGATGGGGAGTCCTTTTCATCGATTTTAGGTGAt                    >  1:217537/1‑81 (MQ=255)                                                               aaTATCTTTTCCATGGC‑AAAATATCAATGTGAATCTAGATGGGGAGTCCTTTTCATCGATTTTAGGTGATCTTCTCtgtg          <  1:217518/80‑1 (MQ=255)                                                               aaTATCTTTTCCATGGC‑AAAATATCAATGTGAATCTAGATGGGGAGTCCTTTTCATCGATTTTAGGTGATCTTCTCTGTGa         >  2:217540/1‑81 (MQ=255)                                                                 tatCTTTTCCATGGC‑AAAATATCAATGTGAATCTAGATGGGGAGTCCTTTTCATCGATTTTAGGTGATCTTCTCTGTGata       >  2:217534/1‑81 (MQ=255)                                                                      tttCCATGGC‑AAAATATCAATGTGAATCTAGATGGGGAGTCCTTTTCATCGATTTTAGGTGATCTTCTCTGTGATATTTCa  >  2:217532/1‑81 (MQ=255)                                                                                |                                                                        TATTTCATGCTACCCTGTCGGAAATGAAGAGTTGATTGAAGCTCTATGTACAGTGATGCTTCAATATCTTTTCCATGGCAAAAATATCAATGTGAATCTAGATGGGGAGTCCTTTTCATCGATTTTAGGTGATCTTCTCTGTGATATTTCA  >  NC_000913/1199420‑1199570

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