Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,306,136 | 0 | A | C | 52.2% | 32.1 / 37.3 | 23 | G575G (GGT→GGG) | deaD | ATP‑dependent RNA helicase |
Reads supporting (aligned to +/- strand): ref base A (11/0); new base C (0/12); total (11/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.40e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-05 |
GAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGATGGTGGAGTGAGA > NC_000913/3305986‑3306282 | gAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGa < 1:401220/151‑1 (MQ=255) acgacgacgACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGacc > 2:302982/1‑150 (MQ=255) caCGACGGGGGGCGCGGCCTCCACGACGTTCGCCGCTAAAACGACGACCATCCCCACGGCCCCCTTCACGGCGTTCCCCGCTGAAGTT‑ACGCCCGCCTTCACGACGTTCGCCACCAAAA‑CCCCGCCCACCGCCACGACGCTCACCGCCa < 2:109701/149‑1 (MQ=255) caCGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGt > 2:474356/1‑151 (MQ=255) cctccccgactTCCGCCGCAAAACCGCCGCCCACCCCCCCGCCCCCCTCCCCGGCGTCCCCCGCTAAAGTT‑ACGACCCCCTCCCCGCCGTCCCCCCCCGAAC‑CCCCGCCCACCGCCACGACGCTCACCGCCAGTATGCGGCtgtg < 1:711037/139‑1 (MQ=255) agttcccccgcAAAACCGCCTACCATCCCCCCGGCCACCTCCCGGGGGTTCCCCGGAGAATTT‑CCGCCCCCCTCCACGAC‑TTCCCCCCCAAAACCCCCGCCCGCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTaa < 1:543883/143‑1 (MQ=255) cgTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACt > 1:646103/1‑151 (MQ=255) gTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTg > 1:357946/1‑151 (MQ=255) ccgcTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATg > 2:74806/1‑150 (MQ=255) ccgaaaaaaCGCCCACCACCCCCCCGCCCCCCCTCGCGGCTTCCCCCGCTAAGGTT‑CCGCCCCCCTTCCCGCCGTCCCCCCCCAAAC‑CCCCGCCCCCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCa < 2:897004/143‑1 (MQ=255) aaCGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCgg > 1:364087/1‑151 (MQ=255) acgacCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCtt > 1:99355/1‑150 (MQ=255) tCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTgaga > 2:769235/1‑150 (MQ=255) cGGCCCCCCCCACGGCGTCCCCCGTTAAAGTT‑CCGCCCCCCTTAACGACGTCCCCCCCCAAAC‑CCGCCCCCCCCGCCACGACCCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCg < 2:604324/149‑1 (MQ=255) cccctccccGGCGTCCCCCGCT‑AAGTTGGCGCCCCCCTCCACGACG‑TCCCCCCCGAAACCCCCGCCCCCCGCCCCGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGt < 2:760428/140‑1 (MQ=255) aCGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTaa > 2:57093/1‑149 (MQ=255) cccgtgaaattt‑ACGCCCCCCTCCCCGCCTTTCGCCCCCAAAC‑CCCCGCCCCCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTgcagc < 1:132022/141‑1 (MQ=255) aCCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTgcag > 2:938620/1‑150 (MQ=255) ttccgccccccTCCACGCCGTCCCCCCCCAAAC‑CCCCGCCCCCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCAcc < 1:583482/140‑1 (MQ=255) cgTTCCCCCCCGAAACCCACGCCCCCCCCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGc < 1:851104/148‑1 (MQ=255) cccccGAAC‑CCCCGCCCCCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAg < 2:323043/140‑1 (MQ=255) ccGAAC‑CCACGCCCCCCGCCACGACCCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGAtggtg < 1:290139/150‑1 (MQ=255) gAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGAtggagg > 1:945491/1‑149 (MQ=255) ccACGCCCCCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGATGGTGGAGTgag < 1:474438/151‑1 (MQ=255) aCGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGATGGTGGAGTgaga > 2:357787/1‑150 (MQ=255) | GAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGATGGTGGAGTGAGA > NC_000913/3305986‑3306282 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |