Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,580,095 | 0 | A | C | 22.9% | 103.7 / 12.5 | 48 | G456G (GGT→GGG) | hsdS | type I restriction enzyme EcoKI specificity protein |
Reads supporting (aligned to +/- strand): ref base A (30/7); new base C (0/11); total (30/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.41e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAGGATGGA > NC_000913/4579958‑4580229 | aaaTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAcc < 1:86331/140‑1 (MQ=255) aaTGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTg > 1:841682/1‑140 (MQ=255) gCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGt > 2:634092/1‑140 (MQ=255) tttCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAg > 2:559089/1‑140 (MQ=255) tttCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAg > 2:589342/1‑140 (MQ=255) aTTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGAttttt > 1:758459/1‑140 (MQ=255) tttATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAAcgc > 1:705980/1‑140 (MQ=255) tgtgGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACgcggcg > 2:1025080/1‑140 (MQ=255) aTGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGgcgc > 1:1129110/1‑136 (MQ=255) aTGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGgcgc < 2:1129110/136‑1 (MQ=255) ggCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTcc > 1:1009370/1‑140 (MQ=255) gCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCg > 1:454430/1‑140 (MQ=255) gCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCg > 2:1235192/1‑140 (MQ=255) gCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCg > 2:623796/1‑140 (MQ=255) cGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGc > 2:1227751/1‑140 (MQ=255) cgcACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGttt > 2:163018/1‑140 (MQ=255) cgcACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGC‑TTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGtttt < 2:595313/140‑1 (MQ=255) cccGAAAAGGTGCGCAAAAAAATAATGTTCAGGATTTTTTACGTGAGGTTTTTTTCCCCCCGCTAGCGGCGCGTTCAGTTTTGTTTTTTCCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGTTTTTCg < 2:1009370/139‑1 (MQ=255) ggAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGcc < 2:705980/140‑1 (MQ=255) ggTGCGCAAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCa < 1:457900/140‑1 (MQ=255) tGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACt > 2:503180/1‑140 (MQ=255) tGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACt > 1:418125/1‑140 (MQ=255) cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg > 2:137504/1‑140 (MQ=255) cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg > 2:69388/1‑140 (MQ=255) ccAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACt < 1:198523/136‑1 (MQ=255) ccAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACt > 2:198523/1‑136 (MQ=255) cAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGTTTTTCGGCCCGCCACTGGGCg < 1:1227751/140‑1 (MQ=255) cAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCg < 1:137504/140‑1 (MQ=255) cAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCg > 1:247869/1‑140 (MQ=255) aaaaaTAATGTTCAGGATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGt < 1:69388/139‑1 (MQ=255) aaaTAATGTTCAGGATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg < 1:865136/133‑1 (MQ=255) aaaTAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg > 2:865136/1‑133 (MQ=255) aaatGTTCAGGATTTTTTACGTGAGGTTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGTTTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGtt < 2:370964/139‑1 (MQ=255) tGTTCAGGATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATAAAACCCGGGTTTCCGGCCCGCCACTGGGCGGTAAGTTcac < 1:1235192/140‑1 (MQ=255) ttCAGGATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTTTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGt < 2:319537/130‑1 (MQ=255) ttCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGt > 1:319537/1‑130 (MQ=255) aGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACgg > 2:625462/1‑140 (MQ=255) gATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGaa > 1:290363/1‑140 (MQ=255) tttACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCtt > 2:868908/1‑140 (MQ=255) aCGTGAGGTTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGTTTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGTTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTg < 2:591925/140‑1 (MQ=255) gTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGcc < 1:1025080/140‑1 (MQ=255) aGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAgg > 1:822915/1‑140 (MQ=255) ggCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGc < 1:1220928/129‑1 (MQ=255) ggCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGc > 2:1220928/1‑129 (MQ=255) gCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAggat > 1:941177/1‑140 (MQ=255) gCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAggat > 1:1148302/1‑140 (MQ=255) ttttttACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAggatgg < 1:868908/140‑1 (MQ=255) tttttCCCCCCGTTAGCGGCGCTTTCAGTTTTGATTTTTCCCAGCAGCGGGGCGGCGCTGTTTTCCCCGTTGATAAAATCGGGTTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAggatgga < 2:180105/140‑1 (MQ=255) tttttCCCCCCGCTAGCTGCGCGTTCAGCTTTGTTTTTTCCCAGCACCGCGGCGGCGCTGTTTTCCCCGCTGATCAAACCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAggatgga < 2:290363/140‑1 (MQ=255) ttttACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTg > 2:1173006/1‑129 (MQ=255) ttttACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTg < 1:1173006/129‑1 (MQ=255) | AAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAGGATGGA > NC_000913/4579958‑4580229 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |