Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,5980TG53.8% 1.6 / 12.0 13S184A (TCA→GCA) stfRputative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base T (6/0);  new base G (6/1);  total (12/1)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGAGGCGGCAACCCATGCGGCTGATGCTGCGGACTCAGCACGCGCAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGG  >  NC_000913/1429372‑1429844
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
naTGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGAGGCGGCAACCCATGCGGCTGATGCTGCGGACTCAGCACGCGCAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCCTCAAAAAGTGCTGCCGCTTCAGAGt                                                                                                                                                                                                                                >  1:145834/2‑251 (MQ=255)
    gtggAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCAGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGAGGCGGCAACCCATGCGGCTGATGCTGCGGACTCAGCACGCGCAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACACAGGCCACTGAAGCATCAAAAAGTGCTGCAGCTGCAGAGTCCTc                                                                                                                                                                                                                            >  1:55/1‑251 (MQ=255)
                                                                                                      ggcccttgtgactgatgcaactgACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAAGCAGCCGCCACGTCTGCCTCCACcgcgg                                                                                                                          <  1:79235/229‑2 (MQ=255)
                                                                                                        cccttgtgactgatgcaactgACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGAGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACcgcggcc                                                                                                                        >  1:17508/21‑248 (MQ=255)
                                                                                                                                     ctcgccgccagcACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGAGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGCATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGaggccgccacttcagc                                                                                           >  1:224908/7‑236 (MQ=255)
                                                                                                                                            ccagcACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACCGAATCAACACAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCAGCGGCCTCCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGCCCTCGGCCCGGGAt                                                                                    >  1:229482/1‑251 (MQ=255)
                                                                                                                                                               ggCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGCAAGGCATCAACCAAGGCCACTCACGAATCCAAAAGTGTTGACGCTGCCGCGTCCTCAAAAAgcg                                                                                                                                                                                                                    >  1:124172/1‑104 (MQ=255)
                                                                                                                                                                                       gTCTTCCAGCGCAGGAACGGAATCAACAAAGGCCAATGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGACGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAaacgaacg                                         >  1:198218/1‑251 (MQ=255)
                                                                                                                                                                                                cgcagaagcGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCTGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAAc                                                                                                                                                                            >  1:15429/9‑111 (MQ=255)
                                                                                                                                                                                                    naagcGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCATCAAAAAACGCAGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGCGATGCGGTGGCCTCAAAAGAGACAGCAAAATCATCAGCAACGAACGCATCATCAAgtgc                            >  1:24836/5‑248 (MQ=255)
                                                                                                                                                                                                    gaagcGGCATCAGCAACGGCCACTGAAGCGGAAAACAGTGCCGCAGCCGCAGAATCCTCACAAAACGCGGCGGACACCAGCGCCGGTGCGGCGAAAACGTCAGACACGACTGCTGCAGCGTCACAACAATCAGCCGACACGTCTGCCTTCACCGCGGCCACGAAAGCGTCAGAGGCACCCACTTCAGCACGAGATGCGGTGGCCTAAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAgtgc                            >  1:233429/5‑248 (MQ=255)
                                                                                                                                                                                                    gaagcGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAACAACCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTTCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAgtgc                            >  1:233441/5‑248 (MQ=255)
                                                                                                                                                                                                                              aGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCAACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAACGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAAcgg  >  1:132140/1‑251 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
GATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGAGGCGGCAACCCATGCGGCTGATGCTGCGGACTCAGCACGCGCAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGG  >  NC_000913/1429372‑1429844

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: