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%	90.9 / NA	28	T46I (ACT→ATT)	sieB	phage superinfection exclusion protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/17);  total (11/17)

TGTGCCCGGTTTTATTTCACCACCTCCGGGCTTCGGTGGTCTCGGCTATACCCCTACAGCGAGAGCTTGTGTTAACATTTCAATACCCTTACAGTTGAGAGTTATTGATATGTTGG‑ATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATACTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTTCTCCTTGCCGTAACCCGGTTATGAAAAA  >  NC_000913/1418562‑1419029                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                               tgtgCCCGGTTGTATTTCACCACCTCCGGGCTTCGGTGGTCTCGGCTATACCCCTACAGCGAGAGCTTGTGTTAACATTTCAATACCCTTACAGTTGAGAGTTATTGATATGTTGG‑ATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCGCGGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGaa                                                                                                                                                                                                                            <  1:537671/250‑1 (MQ=255)  gtgCCCGGTTTTATTTCACCACCGCCGGGCTTCGGTGGTCTCGGCTATACCCCTACAGCGAGAGCTTGTGTTAACATTTCAATACCCTTACAGTTGAGAGTTATTGATATGTTGG‑ATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGaaa                                                                                                                                                                                                                           <  2:77023/250‑1 (MQ=255)             tattcccccaccTCCGGGCTTCGGTGGTCTCGGCTATACCCCTACAGCAAGAGCTTGTGTTAACATTTCAATCCCCTTACAGTTGAGAGTTATTGATATGTTGG‑AGGTATTTACTCCATTGTTGAATCTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGct                                                                                                                                                                                                                 <  2:358188/242‑1 (MQ=255)             tATTTCACCACCTCCGGGCTTCGGTGGTCTCGGCTATACCCCTACAGCGAGAGCTTGTGTTAACATTTCAATACCCTTACAGTTGAGAGTTATTGATATGTTGG‑ATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGct                                                                                                                                                                                                                 >  1:358188/1‑249 (MQ=255)              aTTTCACCACCTCCGGGCTTCGGTGGTCTCGGCTATACCCCTACAGCGAGAGCTTGTGTTAACATTTCAATACCCTTACAGTTGAGAGTTATTGATATGTTGG‑ATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGctct                                                                                                                                                                                                               <  1:216665/250‑1 (MQ=255)                          ccGGGCTTCGGTGGTCTCGGCTATACCCCTACAGCGAGAGTTTGTGTTAACATTTCAATACCCTTACAGTTGAGAGTTATTGATATGTTGG‑ATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGt                                                                                                                                                                                                  <  1:473809/251‑1 (MQ=255)                                  cGGTGGTCGCGGCTATACCCCTACAGCGAGAGCTTGTGTTAACATTTCAATACCCTTACAGTTGAGAGTTATTGATATGTTGG‑ATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTGCAGATTATCGTTTTTGCCt                                                                                                                                                                                          <  1:517106/251‑1 (MQ=255)                                      ggtCTCGGCTATACCCCTACAGCGAGAGCTTGTGTTAACATTTCAATACCCTTACAGTTGAGAGTTATTGATATGTTGG‑ATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGAGTATCGTTTTTGCCTtttc                                                                                                                                                                                      <  1:158490/251‑1 (MQ=255)                                                ataCCCCAACAGCGAGAGCTTGTGTTAACATTTCAATACCCTTAAAGTGGAGATTTATTGATATGTGGG‑ATGTATTTACTCCATTTTTGAAATTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTGCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTc                                                                                                                                                                            <  2:109490/251‑1 (MQ=255)                                                                   ttGTGTTAACATTTCAATACCCTTACAGTTGAGAGTTATTGATATGTTGG‑ATGTATTTACTCCATTGTTGAAACTTTTTGCTAATGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGa                                                                                                                                                          <  1:546383/250‑1 (MQ=255)                                                                     gtgtTAACATGTTAAGCCCCTTACAGTTGATAGTTATTGATATGTTGG‑ATGTATGTACTTCATTGTTTAAGCTTTTTGCTAACGAGCCACTCGAAAGACGTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCGTGGCTCTTGCAGATTATCGGTTTTTCCTTTTCTTTCGTGTTCGCCATTTCCCTCGCAAGAtt                                                                                                                                                        <  2:238799/250‑1 (MQ=255)                                                                          aaCATTTCAATACCCTTACAGTTGAGAGTTATTGATA‑GTTGGGAGGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCGCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATCCCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGc                                                                                                                                                  <  2:58200/251‑1 (MQ=255)                                                                           aCATTTCAATACCCTTACAGTTGAGAGTTATTGATATGTTGG‑ATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATGGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCAGTTCCTTCTCAAGATTGCGAGca                                                                                                                                                 <  1:507312/251‑1 (MQ=255)                                                                               ttCAATACCCTTACAGTTGAGAGTTATTGATATGTTGG‑ATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCAATTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACat                                                                                                                                              >  1:356703/1‑250 (MQ=255)                                                                                     aCCCTTACAGTTGAGAGTTATTGATATGTTGG‑ATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCaa                                                                                                                                        >  1:453165/1‑250 (MQ=255)                                                                                       ccTTACAGTTGAGAGTTATTGATATGTTGG‑ATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCaaa                                                                                                                                       >  1:564141/1‑249 (MQ=255)                                                                                                     gTTATTGATATGTTGG‑ATGTATTTACTCCATTGTTGAAATTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATtaca                                                                                                                       <  2:108969/251‑2 (MQ=255)                                                                                                                              aCTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCACCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCg                                                                                               >  1:3651/1‑251 (MQ=255)                                                                                                                                     tgttgAAACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCt                                                                                        >  1:3094/1‑251 (MQ=255)                                                                                                                                             cTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAATTACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGagaa                                                                                   >  2:82650/1‑248 (MQ=255)                                                                                                                                                   gCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGc                                                                          <  1:82650/251‑1 (MQ=255)                                                                                                                                                                               ataCGTTTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGATTACTGTCGCATTCAATCATTATATTGAAATACCTTGGCTCTTTCAGCTTATCGTTTTTGCCTTTTCTTTAGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGAAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAgg                                               >  1:393070/1‑250 (MQ=255)                                                                                                                                                                               ataCGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGa                                              <  2:323272/251‑1 (MQ=255)                                                                                                                                                                                                     cACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCTATGCTTATATTGAAATCCCTTGGCTCTTTCAGATTATCGTTTTTGTCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGGTTTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTtctc                           <  1:557968/248‑1 (MQ=255)                                                                                                                                                                                                     cACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGAGTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTTCTCCtt                        <  2:380767/251‑1 (MQ=255)                                                                                                                                                                                                      aCTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTTCTCCTTg                       >  1:535041/1‑251 (MQ=255)                                                                                                                                                                                                                     ccGAAAGAGTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTTCTCCTTGCCGTAACCCGGTTAt        >  2:318005/1‑251 (MQ=255)                                                                                                                                                                                                                            agTTTACTGTCGCATTCAATGCTTATATTGAAATACCTTGGCTCTTTCCGATTATCGTTTTTGCCTTTTCTTTCGTAGTCGCCATTTCCTTCTCAAGATTGAGCGCACATATTCAAAAGCATTATTCATTACTACCCGAGAAACGAGTCTTGCTTCGTTTATCTGAGAACGAAATCGCTGTATTTAACGATTTCCTTAAAACAGGCAATCTTATTATCACTTCTCCTTGCCGTAACCCGGTTATGaaaac  >  1:45205/1‑249 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                               TGTGCCCGGTTTTATTTCACCACCTCCGGGCTTCGGTGGTCTCGGCTATACCCCTACAGCGAGAGCTTGTGTTAACATTTCAATACCCTTACAGTTGAGAGTTATTGATATGTTGG‑ATGTATTTACTCCATTGTTGAAACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATACTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTTCTCCTTGCCGTAACCCGGTTATGAAAAA  >  NC_000913/1418562‑1419029

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