Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,797,421 | T→G | 10.5% | A191A (GCA→GCC) | vdh ← | vanillin dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,797,421 | 0 | T | G | 10.5% | 137.7 / 5.9 | 49 | A191A (GCA→GCC) | vdh | vanillin dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (23/19); new base G (2/3); total (26/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.54e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01 |
CCGCCAGTTCGCCGACGATGCGCCCGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGC > NC_002947/3797273‑3797563 | ccgccAGTTCGCCGACGATGCGCCCGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCt < 1:236314/149‑1 (MQ=255) ccgccAGTTCGCCGACGATGCGCCCGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCt > 2:236314/1‑149 (MQ=255) gCCGACGATGCGCCCGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTg > 2:294053/1‑147 (MQ=255) cgacgaTGCGCCCGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCAGGCATCggcgacg < 8:244471/149‑7 (MQ=255) gacgaTGCGCCCGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCCAGGAGTAGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTg > 1:290301/1‑147 (MQ=255) aTGCGCCCGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTg < 8:28713/149‑1 (MQ=255) cccGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCg > 3:103768/1‑148 (MQ=255) tCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATg > 1:4885/1‑148 (MQ=255) aaCCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGAcc > 4:154688/1‑149 (MQ=255) aaCCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGAcc > 7:334742/1‑149 (MQ=255) cGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGAccgcc < 5:144235/149‑1 (MQ=255) tGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCggg < 2:290301/149‑1 (MQ=255) gAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCC‑GCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCt > 7:203837/1‑149 (MQ=255) aGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCt > 1:174732/1‑148 (MQ=255) ttGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGAc > 1:185630/1‑137 (MQ=255) ttGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGAc < 2:185630/137‑1 (MQ=255) cGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGAc > 4:216543/1‑131 (MQ=255) cGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGAc < 3:216543/131‑1 (MQ=255) gTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGcc < 2:174732/149‑1 (MQ=255) cGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGccat > 6:27662/1‑145 (MQ=255) cGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCtt < 6:15024/147‑1 (MQ=255) cGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCtt > 5:15024/1‑147 (MQ=255) cGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGa < 8:334742/149‑1 (MQ=255) cGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGa < 1:230647/149‑1 (MQ=255) cGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCAGGCATCgg < 6:282308/95‑2 (MQ=255) cGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCAGGCATCgg > 5:282308/1‑94 (MQ=255) cAGGGTTGGCGATCAGCG‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCTTTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGc < 4:122740/149‑1 (MQ=255) aGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAg < 8:87880/139‑1 (MQ=255) aGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAg > 7:87880/1‑139 (MQ=255) ggTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGcacc < 3:154688/149‑1 (MQ=255) aGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGcc < 2:4885/149‑1 (MQ=255) c‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCC‑GCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAg < 8:203837/149‑1 (MQ=255) gggggCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGG < 7:310220‑M1/149‑1 (MQ=255) tCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGC > 5:13782‑M1/1‑149 (MQ=255) cTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCAc < 1:288567/149‑1 (MQ=255) ggcgCATTGCTGATGACATTGACCACGCCGTCGCAGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGt < 4:103768/149‑1 (MQ=255) tgatgaCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGc > 1:25624/1‑108 (MQ=255) tgatgaCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGc < 2:25624/108‑1 (MQ=255) atgaCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCAc > 2:323/1‑149 (MQ=255) tgaCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCAcc < 1:323/149‑1 (MQ=255) gaCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCg > 7:60687/1‑149 (MQ=255) cATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCggg > 5:331619/1‑149 (MQ=255) cACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCcctgcgcgcgccctgcggcgtggtgttgg > 6:267208‑M1/1‑120 (MQ=255) cACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTcc > 8:338934/1‑149 (MQ=255) cACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCA < 2:82174‑M1/115‑1 (MQ=255) cACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCA > 1:82174‑M1/1‑115 (MQ=255) gcacGATGCAGGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGtt > 3:168175/4‑139 (MQ=255) gcacGATGCAGGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGtt < 4:168175/136‑1 (MQ=255) ccGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCcctgcgcgcgccctgcggcgtggtgttgggcatcgcacc > 7:132058‑M1/1‑110 (MQ=255) cGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATg > 5:332649/1‑149 (MQ=255) gATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGc > 7:107131/1‑149 (MQ=255) tGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAAT > 7:108959‑M1/1‑99 (MQ=255) tGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAAT < 8:108959‑M1/99‑1 (MQ=255) | CCGCCAGTTCGCCGACGATGCGCCCGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGC > NC_002947/3797273‑3797563 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |