| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_hsa_tpiA | 59,109 | 0 | G | A | 31.0% | 30.8 / 19.8 | 29 | pseudogene (58/228 nt) | yabP | pseudogene, pentapeptide repeats‑containing |
| Reads supporting (aligned to +/- strand): ref base G (14/6); new base A (2/7); total (16/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.05e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCA > NC_000913_3_hsa_tpiA/58968‑59257 | 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:645900/1‑144 (MQ=255) ttACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTAcctcctt < 1:69001/149‑6 (MQ=255) ttACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTAcctcctt < 1:655847/149‑6 (MQ=255) ttACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTAcctcctt > 2:386339/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:108647/1‑149 (MQ=255) cccggacggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGa > 1:945732/1‑149 (MQ=255) cccggacggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCCCCTTATATGATGAGGTAATACTTATTAATAAAATCCCACCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAGGTTGCTACTAAAAAAGTAGTGCTGa > 1:1057982/1‑149 (MQ=255) ccggacggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGaa < 2:173277/149‑1 (MQ=255) ccggacggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGaa > 1:473184/1‑149 (MQ=255) ccggacggtgctaaccTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGaa > 2:656321/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) tGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACg > 1:565407/1‑149 (MQ=255) ggACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACt < 2:49665/141‑1 (MQ=255) ggACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACt > 1:49665/1‑141 (MQ=255) gACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCa > 1:495969/1‑149 (MQ=255) | GAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCA > NC_000913_3_hsa_tpiA/58968‑59257 |
| 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: ‑ |