Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,418,086	(T)5→6	coding (144/222 nt)	kilR ←	Rac prophage; inhibitor of FtsZ, killing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,418,081	1	.	T	100.0%	55.2 / NA	17	T50K (ACA→AAA)	kilR	Rac prophage; inhibitor of FtsZ, killing protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (5/12);  total (5/12)

GTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATG‑TTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTGTTTGCCGAGGCAATATATGTCC  >  NC_000913/1418004‑1418150                                                                               |                                                                      gTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATG‑ttt                                                                    >  2:270855/1‑81 (MQ=255)       aCCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATG‑ttttttgtg                                                              >  1:270860/1‑77 (MQ=255)       aCCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATG‑ttttttgtg                                                              <  2:270836/81‑5 (MQ=255)             aCTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTTGTGGTCAgg                                                         <  1:270837/81‑1 (MQ=255)              cTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGTGGTCAGTAATGAATGTTTTTTGTGGTCAGGt                                                        >  2:270859/1‑81 (MQ=255)              cTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTTGTGGTCAGGt                                                        <  2:270838/81‑1 (MQ=255)                  gCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTTGTGGTCAGGTTAc                                                     >  2:270854/1‑80 (MQ=255)                      ttACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTTGTGGTCAGGTTACGAata                                                >  2:270871/1‑81 (MQ=255)                       tACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTTGTGGTCAGGTTACGAat                                                 <  2:270839/79‑1 (MQ=255)                           ggTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTTGTGGTCAGGTTACGAATATATAg                                           <  1:270840/81‑1 (MQ=255)                               aaCCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTTGTGGTCAGGTTACGAATATATAGtttt                                       <  2:270841/81‑1 (MQ=255)                                aCCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTTGTGGTCAGGTTACGAATATATAGTTTTc                                      <  2:270842/81‑1 (MQ=255)                                aCCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTTGTGGTCAGGTTACGAATATATAGTTTTc                                      <  1:270843/81‑1 (MQ=255)                                  cATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGc                                    <  2:270844/81‑1 (MQ=255)                                   aTCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCt                                   <  1:270845/81‑1 (MQ=255)                                     cACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCttt                                 >  1:270876/1‑81 (MQ=255)                                          ccGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAat                            <  2:270846/81‑1 (MQ=255)                                             aGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATAtt                         <  1:270847/81‑1 (MQ=255)                                               gAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTGt                       <  2:270848/81‑1 (MQ=255)                                                                    aGTAATGAATGTTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTGTTTGCCGAGGCAATATATGTcc  >  1:270856/1‑81 (MQ=255)                                                                               |                                                                      GTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATG‑TTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTGTTTGCCGAGGCAATATATGTCC  >  NC_000913/1418004‑1418150

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