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.4% | 108.3 / 12.9 | 49 | G456G (GGT→GGG) | hsdS | type I restriction enzyme EcoKI specificity protein |
Reads supporting (aligned to +/- strand): ref base A (29/9); new base C (0/11); total (29/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.76e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAGGATGGA > NC_000913/4579958‑4580229 | aaaTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAcc < 1:709933/140‑1 (MQ=255) aaTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAccc > 1:537596/1‑140 (MQ=255) aTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAcccc < 1:158353/140‑1 (MQ=255) aaTGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGCCCTTTTTACCCCCGCTAGCTg > 1:583751/1‑140 (MQ=255) gCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGt > 1:662948/1‑140 (MQ=255) aaTTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTc > 1:145186/1‑140 (MQ=255) aaTTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTc > 1:585852/1‑140 (MQ=255) cTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCtt < 1:688949/140‑1 (MQ=255) aTTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGAttttt > 1:512423/1‑140 (MQ=255) atTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAAcg > 2:329011/1‑140 (MQ=255) atTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAAcg > 2:588336/1‑140 (MQ=255) atTTATGGGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGGGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAAcg > 2:368286/1‑140 (MQ=255) tATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACgcgg < 1:99516/140‑1 (MQ=255) tgtgGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACgcggcg < 2:145186/140‑1 (MQ=255) tgtgGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACgcggcg < 2:585852/140‑1 (MQ=255) gTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCt > 1:84193/1‑140 (MQ=255) gTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCt > 1:741847/1‑140 (MQ=255) gTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCt > 2:653793/1‑140 (MQ=255) gCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCg > 1:3696/1‑140 (MQ=255) gCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCg < 2:203998/140‑1 (MQ=255) gCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCg > 1:405197/1‑140 (MQ=255) gCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCg > 1:376290/1‑140 (MQ=255) cgcACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGttt > 1:207070/1‑140 (MQ=255) tGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACt < 2:405197/140‑1 (MQ=255) cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg > 2:565735/1‑140 (MQ=255) cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg > 2:461313/1‑140 (MQ=255) cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg > 1:613698/1‑140 (MQ=255) cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg > 1:614600/1‑140 (MQ=255) cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg > 1:18910/1‑140 (MQ=255) cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg > 1:463353/1‑140 (MQ=255) cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg > 1:733262/1‑140 (MQ=255) cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg > 2:154975/1‑140 (MQ=255) cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTAATCAAATCCGGGTTTTCGGCCCGCCACTgg > 2:605974/1‑140 (MQ=255) ccAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGc > 1:400483/1‑140 (MQ=255) caaaaaaaTAATGTTCAGGATTTTTTAGGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGTTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCCCCGCTGATCAAATCCGGTTTTTCGGCCCGCCACTGGGc < 2:741847/136‑1 (MQ=255) aGAAAATAATGTTCAGGATTTTTTACGTGAGGTTTTTTTCCCCCCGCTAGCTGCGCTTTCAGCTTTGATTTTTCCCAGCAACGCGGCGGCGCTTTTTTCTCCGCTGATAAAATCCGGGTTTTCGGCCCGCCACTGGGCgg < 2:752005/140‑1 (MQ=255) taatGTCCAGGTTTTTTTACGGGGGGCTTTTTTCCCCCCGCTAGGTGCGGTTTTTTTTTTGATTTTTTCAGGAACCGCGGGGGCGCTGTTTCCTCCGCTGATCAAACCCGGTTTTTCGGCCCGCCACTGGGCGGTAAGtt < 2:376290/140‑1 (MQ=255) ttCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTcacca > 1:464996/1‑140 (MQ=255) aGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACgg > 1:234389/1‑140 (MQ=255) aGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACgg > 2:447784/1‑140 (MQ=255) aGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACgg > 1:27194/1‑140 (MQ=255) ggtttttttACGTGAGGTTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGa < 1:723819/137‑1 (MQ=255) ggATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGa < 1:447784/140‑1 (MQ=255) ttttACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCt < 1:274149/140‑1 (MQ=255) ttttACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGGGGCGGCGCTGTTTTCTCCGCTGATCAAATCTGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCt < 1:263008/140‑1 (MQ=255) aCGTGAGGTTTTTTTCCCCCCGCTAGCTGCGGGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTg < 2:613698/140‑1 (MQ=255) aCGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCCCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTg < 1:54221/140‑1 (MQ=255) cGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGc < 1:565735/140‑1 (MQ=255) gTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGcc < 1:374307/140‑1 (MQ=255) gTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGcc < 1:653793/140‑1 (MQ=255) gTGAGCCCTTTTTCCCCCCGCTAGCTGCGCGTTCAGTTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGTTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGcc < 2:583751/140‑1 (MQ=255) tttttCCCCCCGCTGCCTGCGCGTTCAGCTTTGTTTTTTTCCAGCAACGCGGCGGCGCTGTTTTCCCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAggatgga < 2:27194/140‑1 (MQ=255) tttttCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAggatgga < 2:234389/140‑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 < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |