Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,969,794 | 0 | A | C | 45.1% | 68.1 / 19.6 | 52 | G397G (GGT→GGG) | tap | methyl‑accepting protein IV |
Reads supporting (aligned to +/- strand): ref base A (24/4); new base C (0/23); total (24/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.27e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCT > NC_000913/1969647‑1969942 | ggTGACCGATCTGACAATATCCATCATGGTGGCGGCGGCGTTATTCGCCAGTTTCCAACCCTGCTGGACACGATTGACGCACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGTTACGGCTGGCAAGA‑TTGCGCACTTCACCt < 1:140728/151‑1 (MQ=255) aCCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCa < 1:140117/151‑1 (MQ=255) aCTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTAcc > 2:55073/1‑151 (MQ=255) gACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGc > 2:173728/1‑150 (MQ=255) cAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCaa > 2:478508/1‑150 (MQ=255) tatCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAAcc > 1:397450/1‑151 (MQ=255) ggccgggggggcgTTTTTCCCCAGTTTCAACCCCTGCTGGACACGTTTGACTGCCTCTTCGATCACCCTTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCt < 2:949892/140‑1 (MQ=255) ggTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTg > 2:573902/1‑150 (MQ=255) gTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGc > 1:765177/1‑150 (MQ=255) gcggcggcgTTATTCACCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGctct < 2:213753/150‑1 (MQ=255) cggcggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCtctc > 2:283287/1‑150 (MQ=255) gggtgctttATTCCCCATTTTCGACCCCTGCTGGACTCGATTGACTGACTCTTGGTTCAGCCTTTA‑GATCTCTTTGGCTGCCTGCGCTCTACGGCGGGCAAGA‑TTGCCCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTcc < 2:963020/143‑1 (MQ=255) ggcgTTATTCCCCAGTTTCGAACCCGGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTTTTTGGCTGCCGGCGCGCTACGGCTGGCAAGA‑TTGCGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgc < 2:597982/150‑1 (MQ=255) gcgTTATTCACCAGTTTCGAACCCTGCTGGACACGTTTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TGGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgcgc < 2:785706/151‑1 (MQ=255) aTTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGcc > 2:1008589/1‑151 (MQ=255) aGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTcc > 1:33839/1‑150 (MQ=255) gAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCg > 2:331582/1‑151 (MQ=255) aCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGc > 1:702018/1‑150 (MQ=255) cccTGTTGGACACGTTTGACTGCCTCTTCGATCAGCCTTTT‑GTTCTTTTTTGCTCCCTGCGCGCTCCGGCGGGAAAGA‑TGGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAt < 2:419398/151‑1 (MQ=255) cccTGCCGAACTGGATTGCCTGCCTCTTCACTCAGCCTTTT‑TTTCTTTTTTGCCGCGTGCCCGTTCCGGTGGGCAATC‑TTGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGc < 1:1058057/149‑1 (MQ=255) tgcggGACACGTTTACCTCCCTCTCCGACCGCCCTTTT‑GTCCTTTTTTGCGCCCTGCGCGCTGCGGGGGGAAAGT‑TGCCCCATTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt < 1:464256/149‑1 (MQ=255) gctgGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc > 1:120144/1‑149 (MQ=255) ggACACGTTTGCCTGCCTCTTCGATCACCCTTTT‑GATCTCTTTTGCTGCCTGCGCGCTCCGGCTGGCAAGA‑TGGCCCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc < 2:897243/146‑1 (MQ=255) gaacccgtttGACTGCCTCTCCGATCGCCCTTTT‑GATCTTTTTTGCTGCCGGCGCGCAAGGGCGGGAAAGT‑TGCCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc < 1:171598/138‑1 (MQ=255) acacTATTGACTGCCTCTCCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCGGGCAATT‑TTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc < 2:713273/144‑1 (MQ=255) cagtaTTGCCTGACTCTCCGATCACCCTTTT‑GTTCTCTTTTGCTCCCGGCGCGCTCCGGCTGGCAAGA‑TTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAg < 2:849283/147‑1 (MQ=255) gATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAat > 1:638368/1‑151 (MQ=255) ttGAGTGCCTCTAGGATCACCCCTTT‑GCTCTTTTTTGCGCCCTGCGCGCTCCGGCTGCAAAGA‑TTGCGCTCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGa < 2:721970/147‑1 (MQ=255) tCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGaa > 1:824316/1‑149 (MQ=255) gATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGc > 2:223674/1‑150 (MQ=255) aTCAGCCCTTT‑GATCTCTTTTGCTGCCCGCGCGCTACGGCGGGCAAGA‑TGGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAg < 2:12820/148‑1 (MQ=37) tCAGCCTTTT‑GATCTCTTTCCCTGCCTGCGCCCTAGGCATGGCAAGT‑TGGCCCACTTCGCCTGCCACTCCCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGa < 2:413615/144‑1 (MQ=25) accaccctttggctcttTTTTGCGCCCCGCGCGCTCCGGTGGGCAAGA‑TTGCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGa < 2:1104192/132‑1 (MQ=21) aGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTc > 2:865426/1‑151 (MQ=255) cccTTTGGATCTTTTTTCCTGCCTGCGCGCTCGGGCTGGCAAGA‑TTGCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAAt < 2:608172/148‑1 (MQ=25) ttt‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCg > 1:158737/1‑151 (MQ=255) gttcttTTTTGCGCCCGGCGCGCTCGGGCTGGAAAGA‑TTCCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa < 2:579209/149‑1 (MQ=25) aTCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa > 2:451890/1‑150 (MQ=255) ttttttGCTCCCTGCGCGCTCCGGCTGGAAAGA‑TTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa < 2:30282/145‑1 (MQ=37) tctTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGc > 1:923223/1‑151 (MQ=255) ctTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCt > 1:244336/1‑151 (MQ=255) ttttGCTGCCGCCGCGCTGCGGCGGACAAGA‑TGGCGCCTTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa < 1:738255/145‑1 (MQ=25) tttGCTGCCTGCGCGCTACGGCTGGCAGGA‑TTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCt < 1:808142/149‑1 (MQ=255) cctgcGCGCTCCGGCTGGCAAGT‑TTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATg < 2:951642/149‑1 (MQ=255) cctgcGCGCTCCGGC‑GGCAAAATTTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGt < 1:832537/150‑1 (MQ=25) tgcGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGc > 1:1103486/1‑151 (MQ=255) ccggcggacaaattggcgcACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAttttt < 1:895588/136‑1 (MQ=25) ttGCCCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCtggt < 1:9241/149‑1 (MQ=255) tgccccccTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCtgg < 1:429730/141‑1 (MQ=255) gcgcACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCg > 1:109339/1‑151 (MQ=255) aCTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGAtct > 1:282980/1‑151 (MQ=255) tCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCt > 1:764358/1‑151 (MQ=255) | GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCT > NC_000913/1969647‑1969942 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |