| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 578,427 | 0 | G | C | 57.1% | ‑6.1 / 13.4 | 14 | T107T (ACG→ACC) R34P (CGG→CCG) | rzpD rzoD | putative prophage endopeptidase RzpD putative prophage lysis lipoprotein RzoD |
| Reads supporting (aligned to +/- strand): ref base G (4/2); new base C (1/7); total (5/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.09e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.34e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCC > NC_000913/578202‑578644 | cTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCCGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGCGTTAGCTGATGCGAAAGCTGAAAATGATGCTTTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCt > 1:247164/1‑251 (MQ=255) tCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCg > 1:134637/1‑250 (MQ=35) aaTTACTGACATGCAGATGCGGCAGCGTGATGTTGCTGCGCGCTATGCAAAATACACGAAGGATTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCGGTCTGTCAGTCAGTGCGTGAACCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGc < 1:97123/251‑1 (MQ=11) gataTGCAGATGCGTCAGCGTGAAGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGGTGAGGTTGGCGCTGGTCGTCGTCGGTTGCAAATCAAAGCAGTCTGTCAGTCAGGGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGCCACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCa < 1:2481/247‑1 (MQ=11) gACATGCAGATGCGTCAGCGTGATGATGCTGCGCTCGATGAAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGAAGATGTTGCAGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGCGGATAATGCAGCCTCCCCCAGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCAc > 1:34066/1‑251 (MQ=255) cggcgcTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGGGTGATGATGTTGCCGCTGGACGTCGTCGGTTGCACATCAACGCAGTCTGTCAGTCAGTGCGTGAAGCCCCCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATcn < 1:150590/220‑2 (MQ=11) gtcgtcgtcGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCGCCGCCTCCGGCGTGAATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGc < 1:248930/217‑1 (MQ=255) gtcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCtt < 1:90492/251‑1 (MQ=255) nttGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCCGACACCGCTGAACGGGATTATTTCACACTCagagagag > 1:23807/2‑123 (MQ=11) gTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGACACAACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGCACGGGATTATTTCACCCTCCGAGAGAGGCTGTTCACTATGCAAAAACACCTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTAAGCAATACAAACGCGCTTCCAGCgg > 1:46383/1‑251 (MQ=37) aGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGATCAGTGCAGATAGATTTGCCCAGATCGATGGGCAACTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGcc < 1:43451/251‑1 (MQ=37) gCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTTACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATa < 1:125813/220‑1 (MQ=255) tGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATACACAcgcg < 1:6228/222‑1 (MQ=255) gTGAAGCCCCCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAACGAACCCAGAAGTATAGTAATGAGCATTGCAGAGAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATACACAcgcg < 1:135835/209‑1 (MQ=25) | CTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCC > NC_000913/578202‑578644 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |