Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,7800AG53.8% 1.3 / 10.5 13E244E (GAA→GAG) stfRputative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base G (5/2);  total (8/5)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
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.

ATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTC  >  NC_000913/1429576‑1429962
                                                                                                                                                                                                            |                                                                                                                                                                                      
aTCAACAAAGGCCACTGAAGCATCAAACAGTGCTGCAGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAA‑GAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGc                                                                                                                                         >  1:102070/1‑251 (MQ=255)
    gcaaAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAACACGCGGCGGCCACCAGTGCCGGTGCGGCGAAACCGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCATCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGc                                                                                                                                      >  1:185046/2‑251 (MQ=255)
                                  gcCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACTAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCcgc                                                                                                        <  1:131465/251‑1 (MQ=255)
                                      ncgCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATTAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAaacgaa                                                                                                                                                                 >  1:96683/3‑190 (MQ=255)
                                        gCAGAGTCCTCAAAAAACGCGGCGGCCACCAGGGCCGGTGCGGCGAACACATCAGAAACTAATGCTGCAGCGTCACAACAATCAGCCGCCGCGTCTGCCGCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGTCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAgn                                                                                                                                      <  1:205924/215‑2 (MQ=255)
                                             gTCCTCAAAAAACGCGGCGGCCACCATGGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGGGCAGCTTCCTCGGCAACGGCGGCAGAAAATTctgccagg                                                                                                    <  1:211310/244‑8 (MQ=255)
                                                                     caGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTAACTACAATCAGCAGCCACATCGGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCAACCTCGGCCCGGGATGCGGCGGCCTAAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAggn                                                                     <  1:203582/251‑2 (MQ=255)
                                                                            ggTGCGGCGAAAACGTCAGCAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTctgc                                                                                                        >  1:249672/1‑206 (MQ=255)
                                                                                                  cGAATGCTTCAGAGTAACTACAATCAGCAGCCACCTCTGCATCCACCGAGACCACGAAGGCATCAGAAGCTGAGACATCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGcc                                        >  1:272348/1‑251 (MQ=255)
                                                                                                                                        gcctCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCCAAAGAGGCAGCAAAATCATCAGAAACGAACGAATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGCCGCAAAAACAgcggcggcggggtg  >  1:303128/4‑246 (MQ=255)
                                                                                                                                        gcctCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGAAGCAAAATCATCAGAAACGCACGCATCATCACGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAgcggcggcgggg    >  1:48142/4‑246 (MQ=255)
                                                                                                                                                  cgACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGAAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGaa                                                                                                     >  1:161847/1‑142 (MQ=255)
                                                                                                                                                                                                    tCAAAGGAAGCGGCAAAATCAGCAGAAACGAACAGTTCATCAAGCTCCATTGGGTCATCTTCCGCTCTAAAGGCGGCTCAAAAGTCCGCGAAGCTGGCAAAAACGACCGAGACGAACGCCAGGTCTTCTGaa                                                             <  1:164152/132‑1 (MQ=255)
                                                                                                                                                                                                            |                                                                                                                                                                                      
ATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTC  >  NC_000913/1429576‑1429962

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