| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 282,112 | 0 | T | A | 26.7% | 44.6 / 18.2 | 30 | intergenic (‑129/‑166) | yagA/yagE | CP4‑6 prophage; putative DNA‑binding transcriptional regulator/2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage |
| Reads supporting (aligned to +/- strand): ref base T (11/11); new base A (4/4); total (15/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAG > NC_000913/281964‑282250 | cAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCtt > 5:37472/1‑149 (MQ=255) aaGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGgaga < 2:132734/149‑1 (MQ=255) aTAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTggg < 1:49613/135‑1 (MQ=255) aTAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTggg > 2:49613/1‑135 (MQ=255) aTGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATAtgt > 5:86704/1‑149 (MQ=255) aTGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATAtg < 6:37472/148‑1 (MQ=255) gtgtCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTggg < 1:31138/106‑1 (MQ=255) gtgtCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTggg > 2:31138/1‑106 (MQ=255) gAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTcttc < 4:43878/149‑1 (MQ=255) ttCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCTCCCCAAAGGGGCGA‑gggg < 6:169828/88‑2 (MQ=38) ttCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCTCCCCAAAGGGGCGA‑gggg > 5:169828/1‑87 (MQ=38) cGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTgg > 6:108066/1‑66 (MQ=255) cGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTgg < 5:108066/66‑1 (MQ=255) gACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGtt < 6:86704/149‑1 (MQ=255) ccAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAg > 4:120654/1‑149 (MQ=255) gggACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTtaaa > 7:66535/1‑149 (MQ=255) cATCTTTTCCCCCTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGtt > 1:46940/1‑149 (MQ=255) tCTTTTCCCCCTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTc > 1:50889/1‑148 (MQ=255) tCTTTTCCCCCTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCt > 3:22775/1‑149 (MQ=255) cccggcccTCTCCCCAAAGGGGCGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTtaaataaa < 7:101677/128‑1 (MQ=255) cccggcccTCTCCCCAAAGGGGCGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTtaaataaa > 8:101677/6‑133 (MQ=255) ttCCCCCTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCgct < 1:125785/125‑3 (MQ=255) ttCCCCCTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCgct > 2:125785/1‑123 (MQ=255) ccggcccTCTCCCCAAAGGGGCGAGGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGtt > 5:83086/5‑144 (MQ=255) ccggcccTCTCCCCAAAGGGGCGAGGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGtt < 6:83086/140‑1 (MQ=255) ccccTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGt < 2:46940/149‑1 (MQ=255) cccTCTCCCCAAAGGGGCGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTcttcccttc < 4:4233/81‑1 (MQ=38) cccTCTCCCCAAAGGGGCGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTcttcccttc > 3:4233/1‑81 (MQ=38) cccTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATg > 3:3793/1‑147 (MQ=255) cccTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTa < 4:3793/149‑1 (MQ=255) ccTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAg < 8:66535/149‑1 (MQ=255) | CAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGA‑GGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAG > NC_000913/281964‑282250 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |