Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913_3_hsa_pgi | 225,467 | 0 | G | T | 32.0% | 18.6 / 14.5 | 25 | intergenic (+10/‑33) | ileV/alaV | tRNA‑Ile/tRNA‑Ala |
Reads supporting (aligned to +/- strand): ref base G (10/7); new base T (7/1); total (17/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.05e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTAC > NC_000913_3_hsa_pgi/225443‑225614 | cACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTtaa > 1:356336/1‑149 (MQ=32) aCTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACTCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTtaaa > 1:907946/1‑149 (MQ=17) aCTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTtaaa > 1:1062807/1‑149 (MQ=32) aCTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTtaaa < 1:90772/149‑1 (MQ=32) aGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTtaaataaa > 1:910439/1‑149 (MQ=32) gCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCATGAGTACTGCGGTTAGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATaaaaa < 1:201210/149‑1 (MQ=255) gCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATAGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATaaaaa > 2:1589135/1‑149 (MQ=17) cTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGGAGTGGTTAAATAAAAAAt > 1:437735/1‑149 (MQ=17) tACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATATAAAATa > 1:913412/1‑149 (MQ=17) cttcAAATTTGCCGTGCAAATTTGAAGAGGTTTTATCTACATGTTATGGGGCTATAGCTCTGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGGTCGATGCCGCATAGCTCCACCTTCTCTGTAGTGGTTATATATAAAAtt > 1:1638481/4‑148 (MQ=11) ttcAAATTTGCCGTGCAAATTTGGAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATAc > 2:247435/3‑149 (MQ=255) ttcAAATTTGCCGTGCAAATTTGGAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATAc > 2:135664/3‑149 (MQ=255) ttcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTCGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACAATCTCTGTAGTGGTTAAATAAAAAATAc > 2:1657644/3‑149 (MQ=11) ccAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACt > 2:1033479/1‑149 (MQ=32) ccAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGATGGGAGAGCGCCTGCTTTTCACGCAGGAGGTCTGCGTTTCGATCCCGAATAGCTCCCACATCTCTGTAGTGGTTAAATAAAAAATACt > 2:95853/1‑149 (MQ=11) cAAATTTGCCGTGCAAATTTGGAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACtt > 1:10446/1‑149 (MQ=255) cAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACtt < 2:590129/149‑1 (MQ=32) cAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACtt > 1:1010357/1‑149 (MQ=32) aaaTTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTc > 1:782307/1‑149 (MQ=12) aaaTTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTc > 2:1547295/1‑149 (MQ=12) aaTTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTGCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCa < 2:204292/149‑1 (MQ=17) tttGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCGCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCaga < 1:1323162/149‑1 (MQ=17) gCACTGCAAATTTTAAGAGGTTTTTAGTAAATTTTATGTTGCTATATCTCAGATGGGAGAGCGCATGATTTGCACGCAGGAGGTCTGAGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAgtg < 1:1153177/149‑1 (MQ=11) cGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTAc < 2:221764/149‑1 (MQ=32) cGGCAAATTTGAAGAGGTATTAACTACATGTTATGGGGCTATAGCTGAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTAc < 2:693586/149‑1 (MQ=17) | CACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTAC > NC_000913_3_hsa_pgi/225443‑225614 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |