Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,429,615	T→C	A189A (GCT→GCC)	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,615	0	T	C	95.5%	64.9 / ‑6.1	22	A189A (GCT→GCC)	stfR	putative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (15/6);  minor base G (0/1);  total (15/7)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTGATGCTGCGGACTCAGCACGCGCAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGG  >  NC_000913/1429483‑1429838                                                                                                                                     |                                                                                                                                                                                                                                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)                                                                             cccggcgcAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGcatc                               <  1:224359/247‑1 (MQ=255)                                                                               cggcgcAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGcatcat                             >  1:57351/3‑251 (MQ=255)                                                                                 gcgcAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGCGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCCACATCATCAGAACCGAACGCatcatc                            >  2:71444/1‑250 (MQ=255)                                                                                        agcGGCATCAGCAAAGGCCACGGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGccg                    <  1:128612/249‑2 (MQ=255)                                                                                              aTCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGACCTCAAAAAACGCGGCGGCCACCCGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATAAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCACAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGccggacgtg              >  1:30911/1‑244 (MQ=255)                                                                                                   caaAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAACGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCt         >  1:125993/1‑251 (MQ=255)                                                                                                     aaGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTAAAAAAAAGCGGCGGCCACCAGTGCCGGTGCGGCGAAACAGTCAGAAACGAATGCTGCAGCGTAACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAGGCGTCAAAGGCCGCCCCTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAACTCATCAGCAAAGAACGAATCATCAATGGCCGGTCGTGCAGCTTc       >  2:176063/1‑251 (MQ=255)                                                                                                         ccaatgaagcggaaaaaaGTGCCGCAGCGGCAGATACCCCAAAAACCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTCCACCGGCCCAACAACCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCg   <  2:90489/238‑1 (MQ=255)                                                                                                          cactgaagcggaaaaaaGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCgg  >  2:73900/13‑251 (MQ=255)                                                                                                                                     |                                                                                                                                                                                                                                GGCTGATGCTGCGGACTCAGCACGCGCAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGG  >  NC_000913/1429483‑1429838

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: ‑