| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 578,574 | 0 | A | C | 68.0% | 27.9 / 10.6 | 27 | intergenic (+6/+26) | rzpD/borD | putative prophage endopeptidase RzpD/prophage lipoprotein BorD |
| Reads supporting (aligned to +/- strand): ref base A (5/3); new base C (7/10); total (12/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.11e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAA > NC_000913/578338‑578773 | gAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGc < 1:91269/251‑1 (MQ=255) tgatgTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAattattgtg < 1:10178/250‑4 (MQ=255) atgTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCCGTGCGTGAAGCCACCACGGCCTCAGACGTGGATCATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCATCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATAGATGGGCAACTCATGc > 1:195195/1‑238 (MQ=25) atgTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTAGATTAATGAGCAGTTCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAgn < 1:23223/251‑2 (MQ=255) nacaTCAAAGCAGTCTGTCAGTCAGTGCATGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCAATATCGATGGGCAACTCATGCAATTATTGTGAGCAATa > 1:137217/2‑229 (MQ=25) aaaGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCACAGAGAGGCTGATCACTATGCAAAAACATCTGGGAGGAACCCAGAAGTATATTAATGAGCAGTGCAGAGAGAGCTTTCCATATCGATGGGCAACTCATGCAATGACTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATg < 1:216971/251‑1 (MQ=37) tGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGAtn < 1:137047/185‑2 (MQ=255) ggcagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGAAGCCTCCCCCCGACTGGAAGACACCGCGGAACGGGATTCTTTCCCCCTCACAGAGAGGCTGATCGCTATTCAAAGACAACTGGAAGGAACAAAGACGAATATTAATGAGCAGTGCAGACAGAGCTGACCATATCGATGAGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGGGTATAAATGCCTAAAGtaataa > 1:115009/3‑250 (MQ=255) gCGTGAAGCCACCACGGACTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGCACCCAGAAGTATATTAATGAGCCGTGCAGATAGAGCTGACCATATCGATGGGCAAATCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCg > 1:67711/1‑249 (MQ=255) gTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTgg > 1:3900/1‑251 (MQ=255) tGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATggg < 1:117080/150‑1 (MQ=255) tgATAATGCAGTCTCCCCCCGACTGGTAGACACCGCTGAACGGGATTATTTCACCCTCAGCGAGAGGCAGATCACTATGCAAAACCCACTGTAAGGAACCCAGTCGTATCTTAATGCGCAGTGCAAATCGAGCTGACCATATCGATAGGCAACTCATCCATTTATTTAAAGCAATACACACGCGCTCTCAGCGGAGTATCATTGCCTACAGTAATAAAACCGAGCAATCCATTTACGATTGTTTGTTGGGt > 1:215331/2‑251 (MQ=255) ggATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGTCCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGTAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTggg < 1:210961/250‑1 (MQ=255) gATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGCAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGACAACTCATGCAATTATTGTAAGCAATACACACGCGCTTCCAGCg > 1:118773/1‑192 (MQ=38) gATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAAATGGAAGGAACCCAGAAGTATATTACTGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATac > 1:11326/1‑176 (MQ=38) ggCAGACACCGCTGAACGGGATTATTTCTCCCACAAAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATGAAGGAGCAGTGCAGAGAGAGCTGACCATACCTATGGGCAACTCATGCAATTATTTTGAGCAATACCCACGCGCTTCCAGCg < 1:100802/169‑1 (MQ=255) aCGNGATTATTTCACCCTCAGAGAGAGGCTGACCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATa < 2:137217/137‑1 (MQ=255) attattTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGAGGGGCAACTCATGCAAGTATTTTGAGCAATACACACGCGCTTCCAGCTGAGTATAAATGCCTAAAGTAATAAAACCGATCAATCCATTTACGAATGTTTTCTGGGTTTCTGTTTTAACACCATTTTCTGCGCCGCCACAAATTTTAGCTg < 1:160777/251‑1 (MQ=255) tttCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTCTTTGGAGCAATacactcnctcnnnnncc > 2:137712/1‑131 (MQ=255) ttCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAAGGAGCAGTGCATATAGAGTTGCCCATATCGATGGGCAACGCATGCAATTATGGTGAGCAATACTGACGCGCTGCCAGCGGGGTATAAATGCCTAAGGTAATGAAGCCGAGTAATGGTGTTGCTAATGTTTTGGGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATc < 1:97003/250‑1 (MQ=255) gagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTttcttc < 1:13273/251‑1 (MQ=255) cACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCATTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATACCCACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTcn < 1:137712/251‑2 (MQ=255) gNNAAAACANNTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATac < 2:11326/99‑1 (MQ=255) caacTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTCTTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTACTAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACgaagaa > 1:207223/1‑250 (MQ=255) ggaaCCCAGAAGTATCTTAATCAGCAGTGCAGATAGAGTTGCCCATATCGATAGGCTACTCATGCAATTATTGTGAGCACTACACACGCGCTTCCAGCGGCGTATGAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGac > 1:217624/1‑144 (MQ=255) ccAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATg > 1:114273/1‑59 (MQ=255) tAATGANNAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCt < 2:3543/96‑1 (MQ=255) | GAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAA > NC_000913/578338‑578773 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |