Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913_3_hsa_tpiA | 59,104 | A→C | pseudogene (53/228 nt) | yabP → | pseudogene, pentapeptide repeats‑containing |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913_3_hsa_tpiA | 59,104 | 0 | A | C | 81.8% | 26.3 / ‑2.9 | 11 | pseudogene (53/228 nt) | yabP | pseudogene, pentapeptide repeats‑containing |
Reads supporting (aligned to +/- strand): ref base A (1/1); new base C (2/7); total (3/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.91e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGAC > NC_000913_3_hsa_tpiA/58968‑59251 | gAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGc < 2:600091/149‑1 (MQ=255) gggATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTa < 1:1022203/149‑1 (MQ=255) gggATTGTTAACTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTa < 2:810747/149‑1 (MQ=255) tGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTAcctcc < 1:732002/149‑4 (MQ=255) ttACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTAcctcctt > 2:386339/1‑144 (MQ=255) ttACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTAcctcctt < 1:655847/149‑6 (MQ=255) ttACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTAcctcctt < 1:69001/149‑6 (MQ=255) ttACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTAcctcctt > 2:645900/1‑144 (MQ=255) tACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTACTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTAcctcctta < 1:1077697/149‑7 (MQ=255) gTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATtaata > 2:1139388/1‑149 (MQ=255) ggctacaccggacggggcgaGACCGGACGGTGATAACCTGGAAGGAGCTACCGTGAACGGCGCATCACCATTAGTCGTTGAGGTAATCTTTATTAATAAAAAAAAACCAAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTg < 2:720018/148‑1 (MQ=255) cggacggtgctaacccggacggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTaaa > 1:297578/1‑149 (MQ=255) agggtgctaacctggacggtgctaacctggacGGTGCTACCGTGAACGGTGCTACCTCGTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAATGGAGTTGCTACTGAAGAAGTTGCTACTaaaaaa < 2:622678/147‑1 (MQ=255) cgaacccggacggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTActg < 1:400747/147‑1 (MQ=255) cccggacggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGa > 1:945732/1‑149 (MQ=255) cccggacggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGa > 1:108647/1‑149 (MQ=255) cccggacggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCCCCTTATATGATGAGGTAATACTTATTAATAAAATCCCACCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAGGTTGCTACTAAAAAAGTAGTGCTGa > 1:1057982/1‑149 (MQ=255) ccggacggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGaa > 2:656321/1‑149 (MQ=255) ccggacggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGaa < 2:173277/149‑1 (MQ=255) ccggacggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGaa > 1:473184/1‑149 (MQ=255) ggacggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGCGGTAATTATTATTAATAAAATCCCCCCCAAAAAGATTAATCCGAAAGCAGTTGCTACTGAAGACGTTGTTACTATCAACGTACTGACGAACa > 2:58149/1‑149 (MQ=255) gacggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACaa > 2:1013219/1‑149 (MQ=255) ggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTaaaaaa < 2:340784/132‑1 (MQ=255) ggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTaaaaaa > 1:340784/1‑132 (MQ=255) tgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTa < 2:1089337/107‑1 (MQ=255) tgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTa > 1:1089337/1‑107 (MQ=255) taaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTAc > 1:850773/1‑105 (MQ=255) taaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTAc < 2:850773/105‑1 (MQ=255) aaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGAc > 2:829999/1‑149 (MQ=255) | GAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGAC > NC_000913_3_hsa_tpiA/58968‑59251 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |