Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 535,925 | 0 | T | G | 59.6% | 53.7 / 65.3 | 47 | G72G (GGT→GGG) | hyi | hydroxypyruvate isomerase |
Reads supporting (aligned to +/- strand): ref base T (3/16); new base G (28/0); total (31/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.45e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-09 |
TGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACG > NC_000913/535781‑536072 | tGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGCGGCGGGGGGGGGCGGGAtttgc > 1:351165/1‑148 (MQ=255) gcgcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGGGGGGGGGGGGGGCGGGATTTcc > 1:225166/1‑150 (MQ=255) cgcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTg < 2:397072/151‑1 (MQ=255) cAGTGTGGTTTTCGC‑GCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTGGGCGGGGGGGGGGGGGGGGATTGCCTGt > 1:569263/1‑149 (MQ=255) tgtgGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGATGGGCGGCGGGGGGGGGCGGGATTGCCCGTa > 2:29214/1‑148 (MQ=255) tgGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGGGGGGGGGGGGGGGGTTGTCTGTA‑tttc > 1:552684/1‑148 (MQ=255) ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTA‑TTCCTg < 2:72436/151‑1 (MQ=255) ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTTGGGGGCGGGGGGGGGCGGGTATGACTGTA‑tttc > 1:258146/1‑147 (MQ=255) gTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTA‑TTCCTgg < 2:705141/151‑1 (MQ=255) tttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGAGCGCGGGATTGCCTTTA‑TTCCTCGCCGTGGAGAAGAGgttt > 2:499484/1‑150 (MQ=255) ttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGAGGGCGGGATTTCCTGTA‑TTTCTGGGCGGGAAGAAGAGgtt > 1:678511/1‑149 (MQ=255) ttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGGGGCGGGGGGGGGCGGGTTTGCCCGTA‑TTTCTGGGCGTgaa > 2:544662/1‑140 (MQ=255) aTGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGAGCGGGGGATTGCCTGTA‑TTTCCGGCCGTGCAGAAGAGGTTCggg > 2:655283/1‑151 (MQ=255) tGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGAATGGGCGGCGGGGGGGCGCGGGGATGCCTGTA‑TTCCTTGGCGTGGAGAagagggttgggg > 1:727220/1‑143 (MQ=255) gTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATTGGGGGGGGGGGGGCGCGGGATTGTCCGTATTTCCTGGCC‑TGGAGAAGAGGTTTggg > 2:568652/1‑149 (MQ=255) tttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGCGGTGGGGGGGCGCGGGATTTCCTTTA‑TTTCTGGGCGCGGAGAaaagggtttg > 1:201334/1‑141 (MQ=255) cGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGGGCGCGGGATTGCCTGTA‑TTTCTTGGCGGGAAGAAGAGGTTCGGGGGGGGGGAGGAGCAGCg > 2:96520/1‑151 (MQ=255) aTTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTTGGCGGGGGGGGGGCGCCGGATTGCATGTA‑TTCCTGGCCGGGAAGAAGAGTTTCGGGGTGGGGTAGCAGCAGCGAtt > 2:208551/1‑150 (MQ=255) tGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGAGCGCGGGATTGCCTGTA‑TTCCTGGGCGTGAAGAAGAGGTTCGGGGTGGCGTGGCAGCAGCGGTTCg > 1:729044/1‑150 (MQ=255) agaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTAt < 1:447486/151‑1 (MQ=255) agaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGGATTGCCTGTA‑TTCCTGGCCGTGGAGAAGAGTTTTGGGGTGGCGTAGCAGCAGCGATTTGTTAt > 2:105916/1‑151 (MQ=255) gaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGGGGGGGGGCGCGGGATTGCCTGTA‑TTCCTGGCCGGGGAGAAGAGGTTTGGGGTGGCGGAGGAGCAGCGATTTGTTAt > 2:468036/1‑150 (MQ=255) tAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGccc < 2:2583/150‑1 (MQ=255) tAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGGGGGGGAGCGCGGTATTGCCTGTA‑TTTCTGGGCGGGGAGAAGAGGTTTGGGGTGGGGGAGCAGCAGCGATTTGTTTTGgcc > 1:75607/1‑150 (MQ=255) aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGAATGGGCGGCGGGGGGGCGCGGGATTGCCTGAA‑TTCCTGGCCGCGGAGAAGAGTTTCGGGATGGCGGGGCAGCAGCGAGTCGGTATGAACGTgcg > 2:600516/1‑151 (MQ=255) cATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGGGGGGGGGGGCGGGATTGCCTGTA‑TTTCTGGCCGGGGAGAAGGGTTTTGGGGTGGCGGAGCAGCAGCGATTCGTTATGCCCGTgcgc > 2:638192/1‑151 (MQ=255) tgtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGGGCGGGATTGCCTGTA‑TTCCTGGGCGTGAAGAAAAGGTTCGGGATGGCGGAGCAGAACCGGTTCGTTATTGCCCTGCGGTg > 1:737241/1‑151 (MQ=255) gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTgg < 2:476216/151‑1 (MQ=255) gCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGt < 1:679096/151‑1 (MQ=255) cTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCGGGGGGGCGGGGGGGGCCGGGATTGCCTTTA‑TTTCGTGCCGTGGAGAAGAGGTTTGGGGCGGGGGCGCCACAGCGATTTGTTATGTCCgctcg > 1:641664/1‑140 (MQ=255) tGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGGGGCGGGGGGGGGCGGGGTTGCCTGTA‑TTTCTGGCCGTGGAAAAGAGGTTCGGGGTGGGGGAGGAGCAGCGGTTCGTTATGTCCGTGGGCTggg > 1:215510/1‑148 (MQ=255) ggCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaat < 1:506770/151‑1 (MQ=255) ggCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaa < 1:723332/150‑1 (MQ=255) aataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATaaaaaa < 2:388866/150‑1 (MQ=255) aataaACTCGAACATACGCTGCACAATTTACCGGCGGGGGACCGGGGGGGGGGGGGGGGCGGGATTTGCTGTA‑TTTCTGGCCGGGGAGAAGAGATTTGGGGTGGCGTAGCAGCAGCGATTCGTTAATCCCCTGGGCTGGGGAATaaaaaa > 2:502255/1‑150 (MQ=255) ataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGGGCGCGGGGTTTCCTTTA‑TTCCTTGCCCGGAAGAAGAGGTTTGGGATGGGGTAGCAGCAGCGATTTGTTTTGCCCGTGCGGTGGGGAATaaaaaaa > 2:600314/1‑150 (MQ=255) aaaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACCGGGCGGCGGGGGGGCGCGGGATTTCCTGGA‑TTCCTTGGCGGGGAGAAGGGTTTTGGGGTGGCGGAGCAGCAGCGATTTGTTTTTCCCGGGCGGTGGGGAATaaaaaaaat > 2:8891/1‑148 (MQ=255) cTCGAACTTACGCTGCACAATTTACCGGCGGGGGTACGGGCGGCCGGGGGGCGCGGGTCCGCCTGTA‑TTCCTCGCCGTGGAAAAGAGGTTCGGGGTGGTCTCGCAGCACAGGTCCGTTTTGCACCGGCGCTGGGGGATagaaaaaa > 1:467325/1‑145 (MQ=255) cGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACt < 2:592544/149‑1 (MQ=255) cGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACt < 2:670233/149‑1 (MQ=255) cGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCgg < 2:421215/151‑1 (MQ=255) tGCACAATTTACCGGCGGGTGAATGGGCGGCGGGGGGGCGCGGGATTTCCTTTA‑TTCCTTGGCGTGAAGAAGAGTTTTGGGGGGGGGTAGCAGCAGCGCTTCGTTATGTCCGTGCGCTGGGGGATAAAAAAATTAACTGTCTGGTggg > 1:730776/1‑148 (MQ=255) gCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCTCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTa < 2:106976/149‑1 (MQ=255) cggGGGAGCGCGGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAAc > 1:146348/1‑147 (MQ=255) cgcgGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCa > 1:238007/1‑150 (MQ=255) gcgGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAAc < 2:602302/151‑1 (MQ=255) cgGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACg < 1:296424/151‑1 (MQ=255) | TGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACG > NC_000913/535781‑536072 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |