| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 579,146 | 0 | A | G | 70.6% | 43.7 / 15.9 | 35 | intergenic (‑253/+38) | borD/ybcV | prophage lipoprotein BorD/DUF1398 domain‑containing protein YbcV |
| Reads supporting (aligned to +/- strand): ref base A (6/4); new base G (10/14); total (16/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.57e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTT > NC_000913/578909‑579391 | gCTTTTTGAAGTACGCAGAATCGTAGGTTTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAAATTCATATTGTGAATATTTATTAAGGTATGTCAGGGGCGATGAATCGTCATTGTATTCCCGGATTAACTATTTCCACAGCCCTTTGGGGGAACTTCTCTTCGGGAGTGTCCTGTAAAAATTAAAA‑CGATGCAAGCAGGGTTAAGCGCGTGCGCGTATTGCATTATGCCAACGCCCCGGTGCTg < 1:395353/248‑1 (MQ=255) gCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATTTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCgg > 1:405546/1‑250 (MQ=255) ngTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTAAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCt > 1:80389/2‑226 (MQ=255) gAAAATTAAACAAACCCTAAACACTGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACg > 1:398070/1‑222 (MQ=255) aaaaTTAAACAACCCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCGGACGGGGAACTTCTCTTCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCAGTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTgg < 1:62363/251‑1 (MQ=255) aaaCAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGGTATGATTTAGCGTGGAAAGATGTGTGTAgtg < 1:342826/250‑1 (MQ=255) gTTGAAATTTCATATTGTTAATCTTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAAt > 1:36478/1‑249 (MQ=255) gAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGctn < 1:87026/251‑2 (MQ=255) cGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGACCTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGgaag > 1:78212/1‑155 (MQ=255) cGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTtg > 1:355217/1‑182 (MQ=255) aTTGTATTCCAGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTATGCGGGAGTGTCAGGGAATAATTCAAACCGATGCAAAAAGGGGTTAGCGAGCACATGTATCGTATTATGCCAACACACCGGTGCTGACACGGCAGAAACCGGACGTTATGAtt > 1:406221/1‑162 (MQ=255) gTATTCCCGGATGAACTATTACCACAGCCCTGACGGGGAACTTCTCTGAGGGAGGGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTAGTATGCCAACACCCCGTTGCTGACACGGAAGAAACCGGACGTTATGATTTATCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGGAATGAATTATCAACGGTATAGTAATATCTTTTATGTTCGTgta < 1:19596/250‑3 (MQ=255) tATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGct < 1:165207/192‑1 (MQ=255) ggATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTACGGGAATACTTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACACCGGTGCTGACACGGAAGCAACCGGACGTTATGATTTAGCGTGGAAAGATTtg > 1:388970/1‑168 (MQ=255) acaGCCCTGACGGGGAACTTCTCTGCGGGAGTGTNCGGGAATAANNAANAACGANGNNNANNNGGNTTAGCGCGTACATGTATTGNATNNTGCCAANACCCCGGTGCTGACACGGAAGAAAcc < 2:102527/123‑1 (MQ=255) acaGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAAcc > 1:102527/1‑123 (MQ=255) ncccTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATCATTAAAAACGATGCACCAAGGGTTTAGCCCGTACATGTATTGTATTATGCCAACACCCCGGCGCTGACACGGAAGAAACCGGACGTTCTGATTTAGCGTGGAACGTTTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTg > 1:38657/2‑222 (MQ=38) gCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATACTTAAAA‑CGCTGCACACAGGGTTTAGCGCGTCCACGTATTGCATTATGCCAACGCCCCg > 1:295327/1‑99 (MQ=255) gACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCAGTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCttt < 1:385864/206‑1 (MQ=38) ggAACTTATCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGATGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCa > 1:362897/1‑249 (MQ=255) gCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTGGCAGTAAATAGTAATGAATTATCGAATGTATAGTAATATCTTTGATGTGCATTGATATTTGTAGCCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTg < 1:328578/249‑1 (MQ=255) ttaaaaCGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAAc > 2:28311/3‑77 (MQ=255) ttaaaaCGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAAc < 1:28311/75‑1 (MQ=255) nanaCNATGCANACAGNGTTTAGCGCGTACATNTATTGTATGATGCCAACACNNCNGTGCTNACACg < 2:191627/66‑1 (MQ=255) aaaCGATGCACACAGGGTTTAACGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCATGGAAAGATTTGTGTAGTGTTCGGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTACGTTCATGGAGATTGGCAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAAtgtg < 1:434888/241‑1 (MQ=25) acacaGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATTAATTATCAAAGGGATAGTAATATCTTTTATGTTCATTGATATTTGTAACCCATCGTAAAACTCCGGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAAcc < 1:66993/251‑1 (MQ=255) naGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGaaa > 1:109513/2‑187 (MQ=255) gggTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATg < 1:13056/157‑1 (MQ=255) ttAGCGCGTACACGTATTGCATGATGCCAACTCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGATCACAGTAAATAGTAATTAATGATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTAGTAACCCATCGGAAAACTCCTGCGTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCGATCATAg < 1:315674/250‑1 (MQ=25) aGCGCGTACTCGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGAc < 1:402202/251‑1 (MQ=255) aGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGAc < 1:141167/251‑1 (MQ=255) gcattatGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATc < 1:16642/226‑1 (MQ=255) tatGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGAGTTCTCTTGATTTCAACNt < 1:348450/221‑1 (MQ=255) ngCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGt > 1:275728/2‑185 (MQ=255) ccAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCtt > 1:94623/1‑251 (MQ=255) | GCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTT > NC_000913/578909‑579391 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |