| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,429,783 | 0 | G | A | 46.7% | ‑6.4 / 14.9 | 15 | A245A (GCG→GCA) | stfR | putative prophage side tail fiber protein StfR |
| Reads supporting (aligned to +/- strand): ref base G (6/2); new base A (5/2); total (11/4) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-02 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCA > NC_000913/1429555‑1430009 | gTCTTCCAGCGCAGGAACGGAATCAACAAAGGCCAATGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGACGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAaacgaacg > 1:198218/1‑251 (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) ggaaaaaaGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCAACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGc > 1:110498/4‑231 (MQ=255) cgcagccgCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGGGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTctn < 1:55440/244‑3 (MQ=255) agagTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTctgccagggcgg < 1:49061/251‑12 (MQ=255) ccaccaGTGCCGGTGAGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACACTCAGCCGCCACGTCTGCCTCAAACGCGGCCACAAACGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGACCGAGACGAATGcc > 1:84057/1‑251 (MQ=255) ccaccaGAGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAAc > 1:178236/1‑202 (MQ=255) aaCGAATGCTTCAGCGGCACGACAATCAGCAGCCACATCGGCATCCCCCGCGACCACTAAGGCATCAGAAGCGGCGACCGCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAGCGAACGCATCAGCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGTCGTCAGGAAATTCCACGAAGGCGTCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAgcn < 1:46761/251‑2 (MQ=255) aaTCAGCAGCCACATCTGCATCCACCGCGACCACGAAGACATCAGAAGCTGCGACCTAGGCCCGCGATACGGCGCCCTCAAAAGAAGCGACAAAATAATCCGAAACGAAAGCATCATCAAGCGCAAGTCGTGCAGCTTCCTCGGCAACGGAGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGGGGGCCAGAGCGACTCGGCTGCGGCAGGCTCaaa > 1:220550/1‑250 (MQ=255) cggcggCCTCAAAAGAAGCGGCNAAATCAGCANNAANGAACGNANNNTNNNGCNCCAGTAGTGCAGCTTCCTCNGCNNCGGCGGNAGGAAATTCCGCGAAGGCGGCaaaaa < 2:54605/111‑1 (MQ=255) cggcggCCTAAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCaaaaa > 1:54605/1‑111 (MQ=255) aaaaGAAGCGGCAAAATCATCAGAAACGAACGCATCATCACGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCGGGAAATTCCGCGAAGGCGGCCCAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGACAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCa > 1:235874/1‑236 (MQ=255) | GTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCA > NC_000913/1429555‑1430009 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |