Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,429,558	T→C	S170S (TCT→TCC)	stfR →	putative prophage side tail fiber protein StfR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,429,558	0	T	C	100.0%	54.2 / NA	22	S170S (TCT→TCC)	stfR	putative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/9);  total (13/9)

CATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGC‑CCGTGAGGCGGCAACCCATGCGG‑CTGATGCTGCGGACTCAGCACGCGCAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACG  >  NC_000913/1429324‑1429801                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    cATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGCCGGCGATGCCAGTGCATCAGC‑TGCTCAGGTCGCGGCCCTTGTGA‑CTGATGCAACTGACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAgaagcgg                                                                                                                                                                                                                                     >  1:197314/1‑245 (MQ=255)     aCGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGCCGGCGATGCCAGTGCATCAGCTGCG‑CAGGTCGCGGCCCTTGTGA‑CTGATGCAACTGACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATc                                                                                                                                                                                                                                 <  1:157527/251‑1 (MQ=255)               atgCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGCCGGCGATGCCAGTGCATCAGC‑TGCTCAGGTCGCGGCCCTTGTGA‑CTGATGCAACTGACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGACGCGGCATCAGCAAAggac                                                                                                                                                                                                                       >  2:64283/1‑249 (MQ=255)                     ggcgggaggtgctgcgtcggcttGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGCCGGCGATGCCAGTGCATCAGC‑TGCTCAGGTCGCGGCCCGTGTGA‑CTGATGCAACTGACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGGCAGGTCAGGAAGCGTCCTCGGGCGCAGAAGCGGCAGCAGCAAAGGCCTCTGaa                                                                                                                                                                                                                 <  2:42057/230‑1 (MQ=255)                     ggccGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGCCGGCGATGCCAGTGCATCAGC‑TGCTCAGGTCGCGGCCCTTGTGA‑CTGATGCAACTGACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGa                                                                                                                                                                                                                  >  2:2884/3‑250 (MQ=255)                         ggatgtgcTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGCCGGCGCTGCCAGTGCATCAGC‑TGCTCAGGTCGCGGCCATTGTGA‑CTGATGCAACTGACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGAAAGCGTCCTCAGGAGCAGAAGCGGCCTCAGCAAAGGCCACTGCAGcgg                                                                                                                                                                                                             >  2:73486/8‑249 (MQ=255)                          gaggtgcTGCGTCGTCTTTAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGCCGGCGATGCCAGTGCATCAGC‑TGCTCAGGTCGCGGCCCTTGTGA‑CTGATGCAACGGACTCAGCACGCGCCGCCAGCACTTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGcgga                                                                                                                                                                                                            <  1:120666/245‑4 (MQ=255)                                 ggcgtcgtcttGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGCCGGCGATGCCAGTGCATCAGC‑TGCTCAGGTCGCGGCCCTTGTGA‑CTGATGCAACTGACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGcggaaaaaagt                                                                                                                                                                                                     <  2:104673/242‑11 (MQ=255)                                           tGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGTCACAGAGTACGGCAGACGCGAAGAAAGCAGCCGGCGATGCCAGTGCATCAGC‑TGCTCAGGTCGCGGCCC‑TGGGGTCTGATGCAACTGACTCAGCCCGCGCCGCCAGCACGTCCGCGGGACAGGCTGCAGCGTCAGCGCAGGAAGCGTCCTACGGCGCATAAGCGGCATCAGCAAAGGCCACTGAAGctgaaaatagtgccgcagcc                                                                                                                                                                                            <  2:238391/250‑20 (MQ=255)                                                                                                                                                                  gactgatgcaactgACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGaaa                                                                       <  1:99782/237‑2 (MQ=255)                                                                                                                                                                  gactgatgcaactgACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGaa                                                                        >  1:61930/14‑249 (MQ=255)                                                                                                                                                                      gatgcaactgACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGATCAGGAAGCGTCCTCCGGCGCAGAAGCGGAATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGaaagcgtc                                                                  >  1:138256/10‑245 (MQ=255)                                                                                                                                                                            actgACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCCCCGCGGCCACGaaagcggcagaggc                                                            >  2:225193/4‑239 (MQ=255)                                                                                                                                                                               gACTCAGAACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCCGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGCAAAAAGTGCCGCAGCCGCAGCGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGaggccgc                                                         >  1:35351/1‑245 (MQ=255)                                                                                                                                                                                  tCAGCACGCGCCGCCAGCACGTCCTCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCAGCAGCCGCAGAGTCCTCAACAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATAAGCCGCCCCGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGaggccgccac                                                      >  1:152939/1‑242 (MQ=255)                                                                                                                                                                                  tCAGCACGCGCCGCCAGCAAGTACGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGACATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGACCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCTCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGaggccgccac                                                      >  1:101583/1‑242 (MQ=255)                                                                                                                                                                                               ccagcACGTCCGCCGGACAGGCTGCAGCGTCAGCTCAGGAAGCGTCCTCCGGCGCAAAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGaggccgccacttcagcacgagat                                         <  1:69270/251‑23 (MQ=255)                                                                                                                                                                                                cagcACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGaggccgccacttcagcacgagatg                                        >  1:244063/1‑228 (MQ=255)                                                                                                                                                                                                cagcACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGACGCAACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGaggccgccacttcaccacgagt                                          >  2:62616/1‑228 (MQ=255)                                                                                                                                                                                                              gACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTAAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTAAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCaaaa                          <  2:136606/251‑1 (MQ=255)                                                                                                                                                                                                               aCAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAGGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAaag                         <  1:51890/251‑1 (MQ=255)                                                                                                                                                                                                                                      aagCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCCGCAACATCATCAGAaacg  >  2:132206/2‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    CATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGC‑CCGTGAGGCGGCAACCCATGCGG‑CTGATGCTGCGGACTCAGCACGCGCAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACG  >  NC_000913/1429324‑1429801

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