Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,492,028 | 0 | C | T | 42.9% | 43.0 / 16.7 | 37 | Q90Q (CAG→CAA) | ppiA | peptidyl‑prolyl cis‑trans isomerase A (rotamase A) |
Reads supporting (aligned to +/- strand): ref base C (3/17); new base T (15/0); total (20/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.93e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CCATGGTCAAGGAAGGCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGACGC‑GGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCG > NC_000913/3491879‑3492178 | ccATGGTCAAGGAAGGCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTct < 2:2520866/151‑1 (MQ=255) ccATGGTCAAGGAAGGCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTc < 1:2214032/150‑1 (MQ=255) ggTCAGGAAAGGCGTTATCGGGGGCAGAGAGAGACACCTTGCTGGTGACGCTGGCGTTGTCGGCGGTACGTGAAAAAAAGATGGTACCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTGGTTTTTTttgt > 1:333058/1‑148 (MQ=255) ggTCAAGGAAGGCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTctgct < 1:3574690/150‑1 (MQ=255) cAAGGAAGGCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCa < 1:2614457/150‑1 (MQ=255) cGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTGGTTTTTTTTGCTGCATCTg > 2:2098828/1‑145 (MQ=255) tATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGa < 1:3905773/150‑1 (MQ=255) cGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAAcc < 1:1387627/151‑1 (MQ=255) gCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTTGTTTTTTTTTCTTCCTCTGCTC‑GGGTGAAccgc > 2:2953116/1‑151 (MQ=255) ttGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTTGGTTTTTT‑TTCTGCCTCTGCCC‑GGGGGAAACGCCGccc > 2:4061034/1‑150 (MQ=255) aTAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTTGTTTTTTTTGGTGCATCTGCTC‑GGGGAAACCGCCGCCCTg > 2:912164/1‑149 (MQ=255) cgcTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGAc < 2:2518856/151‑1 (MQ=255) gcTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTTTTCTGCATCTGCTC‑GGGGAAACCGCCGCCCCG‑AATCATAAAACCAGGaa > 2:1804032/1‑146 (MQ=255) gcTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGAcg < 2:3363309/151‑1 (MQ=255) gTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGAcgc < 1:3330944/149‑1 (MQ=255) tttGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTTGTTTTTTTTTCTTCCACTGCTC‑CGTGGAAACGCCGCCCCG‑TAAAATAAAAAAAAGAATGACGC‑gggg > 2:2338561/1‑148 (MQ=255) gCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTGGGTTTTTTTTCTTCATCTTCTC‑GGTGGAACCGCCGCCCCG‑AATCATAAAAACCGGGAAGACGC‑gggggaa > 1:625173/1‑141 (MQ=255) gTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTGGGTTTTTTTGCTGCATCTTCTC‑GGGGGAACCCCCGCCCCG‑GATCATAAAACCCGGGAAGACGC‑GGGGAAAGGGAggg > 2:3926618/1‑148 (MQ=255) tGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTTGTTTTTTTTTGTGGATCTGCTC‑GGGGAAAACCCCGCCCCG‑AATCAAAAAACCCGGAAAGACGC‑CGGGAAAAGGAGTgttgttt > 1:205914/1‑150 (MQ=255) gcgATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGACGC‑GGTGAAAGGTAGTGTTGTTATAAAAAc < 2:2845013/151‑1 (MQ=255) aTGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTGGGTTTTTTTTCTGCATCTTCTC‑GGTGGAAACGCCCCCCTG‑AATCATAAAACCAGGAAAGACGCGGGGGAAAGGAGGGGTTTTTTTAAAAAc > 2:1984923/1‑149 (MQ=255) tGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGACGC‑GGTGAAAGGTAGTGTTGTTATAAAAAcc < 2:1372281/148‑1 (MQ=255) tGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTGGGTTTTTTTTCTCTAAATGCTC‑GGGGGAAACGCCGCCCCGCATTCTTAAAG‑CAGGGATGACGG‑CGTGAAAGGAAGTGGTGTTTTtaaaa > 1:1372011/1‑146 (MQ=255) aCGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTGGTTTTTTTTGCTGCCTCTGCTC‑GGTGAAAACGCCGCCCTG‑AATCCTAAAACCAGGGATGAAGC‑GGGGGAAGGGGGGGGTGTTATTAAAACCGCTGtt > 2:508669/1‑147 (MQ=255) gtgtTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTGGGTTTTTCTGCTGCATCTGCTC‑GGTGAAAACCCCGCCCCG‑AATCATAAAACCAGGAATGACGC‑GGGGGAAGGGAGTGGTGTTTTAAAAACCGCTGTTCACATaa > 1:2000500/1‑150 (MQ=255) gtTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGACGC‑GGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATc < 2:3062356/150‑1 (MQ=255) gtTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGACGC‑GGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCg < 1:970298/151‑1 (MQ=255) ttGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTTGTTTTTTTTGATGCTTCTTCTC‑GGGGAAAACACCGCCCTT‑AATCATAAAAACCGGAAAGACGC‑GGTGGAAAGAAGGGTTGTTATAAAAACCCCCGTTTCCATaa > 2:1586414/1‑147 (MQ=255) ccGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGACGC‑GGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGtttt < 1:2331048/151‑1 (MQ=255) cGTTATCGGCTTCATTTGTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGTCGCCCTG‑AATCATAAAGCCAGGAATGACGC‑GGTGAAAGGTAGTGTTGTTATAAAACCCGCTGTTCACACAACCGTCAAAGtttta < 1:1323405/151‑2 (MQ=255) tATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTTTTCTTCATATGCTCTGGGGAAACC‑CCGACCCG‑AATCACAAAAACAGGGAAGACGC‑GGGGAAAGGTAGTGTTGTTATAAAAAACCCTGTTCACATAAACGACaaaaattt > 1:3447936/1‑142 (MQ=255) tATCGGCTTCATTTTTGATTGGC‑GGGGTTGGTTTTTTTTGTTGCATCTGCTC‑GGTGAAAACGCCGCCCTG‑CATCATAAAAACAGGGATGACGC‑GGGGGAAGGGAGTGTTGTTTTAAAAAACCCTGTTTAAAAAAACGACaaaaattttca > 1:3053106/1‑142 (MQ=255) aTCGGCTTCATTTTTGATTGGCGGGGTTTGGTTTTTT‑TGCTTCATCTTCTC‑CGGGGAAACGCCCCCCTG‑AATCATAAAACCAGGGATGAAGC‑GGGGGAAAGGAGGGTTGTTATAAAAAACGCTGTTTACATAATTGACAAAAATTTTcac > 1:3990526/1‑150 (MQ=255) cTTCATTTTTGATTGGC‑GGGGTTGGTTTTTTTTGCTCCATCTTCTG‑CGGGAAAACGCCGCCCCG‑AAACAAAAAAGCAGGAATGAAGC‑GGGGTAAGGGAGGGGTGTTATAAAAACCGCAGTTCAAAAAATTTACCAAGTTTTGcaccaaca > 2:1078579/1‑145 (MQ=255) aTTTTTGATTGGC‑GGGTTTTGATTTTTCTGCTGCGTCTGTCC‑GGGTAAAAAACCGCCCTC‑TATTATTAAATCACGGACGGAGC‑GGGGAAAAGGAGGGGTTCTATAAAAAAAGACGTTCACACAAACGACAAAATTTTGTACAGacaaccg > 1:603792/1‑146 (MQ=255) ttttGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGACGC‑GGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGgcg < 2:2683465/150‑1 (MQ=255) ttGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGACGC‑GGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCt < 2:2784758/150‑1 (MQ=255) tGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGACGC‑GGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTg < 2:2197623/150‑1 (MQ=255) ggC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGACGC‑GGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGtt < 1:1984695/151‑1 (MQ=255) ggtttggtttTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGACGC‑GGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGTTTATc < 1:3427618/151‑1 (MQ=255) ttttttCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGACGC‑GGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGTTTATCCAGCTcc > 1:1602856/1‑151 (MQ=255) ctgctgCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGACGC‑GGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCg > 2:430264/1‑151 (MQ=255) | CCATGGTCAAGGAAGGCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTG‑AATCATAAAGCCAGGAATGACGC‑GGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCG > NC_000913/3491879‑3492178 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |