Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,418,807	C→T	T46I (ACT→ATT)	sieB →	phage superinfection exclusion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,418,807	0	C	T	100.0%	34.5 / NA	13	T46I (ACT→ATT)	sieB	phage superinfection exclusion protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/6);  total (7/6)

CCGGTTTTATTTCACCACCTCCGGGCTTCGGTGGTCTCGGCTATACCCCTACAGCGAGAGCTTGTGTTAACATTTCAATACCCTTACAGTTGAGAGTTATTGATATGTTGGATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGAC‑TTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATACTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTTCTCCTTGCCGTAACCCGGTTATGAAAAAATTAGAACGGAAGGGCA  >  NC_000913/1418567‑1419046                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                cctgTTTTATTTCACCACCTCAGTGCTTCGGTGGTCTCGGCTATACCCCTACAGCGAGAGTTTTTTTTAACTTTTCAATACCCTAACAGTTGAGATTTATTGATATTTTGGATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGAC‑TTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCt                                                                                                                                                                                                                                       <  1:363624/248‑1 (MQ=255)  cgtttttATTCCACCACCTCCGGGCTTCGGGGGTCTCGGCTATACCCCTACAGCGAGGGCTTGTGTTAACAGTTTAATACCCTTACAGTTGTGTGTTATTGATATGGTGTATGTTTTTACATCATTGTTGAACCTTTTTGCTAGCGAGCCACTCGAAAGAC‑TTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGGTTACTGTCGCATTCAATGCTTATATTGAAATACCtt                                                                                                                                                                                                                                      <  2:489516/248‑1 (MQ=255)    gTTTTATTTCACCACCTCCGGGCTTCGGTGGTCTCGGCTATACCCCTACAGCGAGAGCTTGTGTTAACATTTCAATACCCTTACAGTTGAGAGTTATTGATATGTTGGATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGAC‑TTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTgg                                                                                                                                                                                                                                    >  1:406692/1‑251 (MQ=255)                                                      gCGAGAGCTTGTGTTAACATTTCAATACCCTTACAGTTGGGCGTTATTGATATGTTGGATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGAC‑TTATGTATACGATTATCATTTTTGGTCTCAATCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCtt                                                                                                                                                                                  >  2:520473/1‑251 (MQ=255)                                                                                                            ttGGATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGAC‑TTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCATTTTATTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACAAGAGCAACg                                                                                                                            >  1:354189/1‑251 (MQ=255)                                                                                                                        aCTCCATTGTTGAACCTTTTTGCTAACGAGCCACTCGAAAGAC‑TTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCGGAGCAACGAGTATTGCTTCg                                                                                                                <  1:358980/251‑1 (MQ=255)                                                                                                                                 ttgAAACTTTTAGCTACCGAGCCACTCGGA‑CACATTATGTATACGATAATCATTTTTGGTCTCACTCTCTAGCTGATAACGAAAGAGTTTAATGTCGCATTCTATGCTTATATTGAAATACCTTGGCTCTTTAAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCACGATTGCGAGAACATATTCAACAGCATTATTCATTACTTCCAGAGAAACGGGAATTACTTCGTTTATCTg                                                                                                        >  2:537338/1‑250 (MQ=255)                                                                                                                                     aaCTTTTTGCTAACGAGCCACTCGAAAGAA‑GTTTGTATACGATTATCATTTTTGGTCTCCCTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGagaaa                                                                                                   <  1:567612/251‑1 (MQ=255)                                                                                                                                     aaCTTTTTGCTAACGAGCCACTAGAAAGAC‑TTATGTATACGCTTATCATTTTTGGTCTCACTCTCAGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTCCCAAAGCAACGGGTATTGCTTCGTTTATCTGagaa                                                                                                    >  2:273093/1‑250 (MQ=255)                                                                                                                                                                                         ttGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACt                                                >  1:546278/1‑251 (MQ=255)                                                                                                                                                                                                                               gTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTTCTCCTTGCCGTAACCCGGTTATGAAAAAATTAGAAc           <  1:273093/250‑1 (MQ=255)                                                                                                                                                                                                                                       cAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTTCTCCTTGCCGTAACCCGGTTATGAAAAAATTAGAACGGAAGGGca  <  1:520473/251‑1 (MQ=255)                                                                                                                                                                                                                                       cAATGCTTATATTGAAAGACCTTGGCTCTTTCAGATCATCGTTTTTGCCTTTTCTTTCCTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATTTTCAAAGGCATTATTCATTACTACCGGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTGATTATCCCTTCTCCGTGCCGTAACCCGGATATGACAAAATTAGAACggcatgga   >  2:161499/1‑244 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                CCGGTTTTATTTCACCACCTCCGGGCTTCGGTGGTCTCGGCTATACCCCTACAGCGAGAGCTTGTGTTAACATTTCAATACCCTTACAGTTGAGAGTTATTGATATGTTGGATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGAC‑TTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATACTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTTCTCCTTGCCGTAACCCGGTTATGAAAAAATTAGAACGGAAGGGCA  >  NC_000913/1418567‑1419046

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