Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB122 | 6,038,658 | 0 | C | A | 31.1% | 74.0 / 16.5 | 61 | A37S (GCC→TCC) | puuD_3 | Gamma‑glutamyl‑gamma‑aminobutyrate hydrolasePuuD |
Reads supporting (aligned to +/- strand): ref base C (24/18); new base A (0/19); total (24/37) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.95e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACACTGGGGGCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTT > ACB122/6038522‑6038772 | acacTGGGGGCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGTTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGc > 1:27981/1‑137 (MQ=255) gggggCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAgc > 6:287564/1‑137 (MQ=255) gggggCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAgc > 8:225198/1‑137 (MQ=255) gcgcCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAgccggc > 5:112107/1‑137 (MQ=255) gcgcCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAgccggc > 2:177141/1‑137 (MQ=255) gcgcCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAgccggc > 7:281033/1‑137 (MQ=255) gcgcCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGGAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAgccggc > 1:245940/1‑137 (MQ=255) tttGTAATGAGGCGGCTCGATATTTGACGCCGAGCCGGTAAAAACAATGCCGTGCCGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCAATCACGGGCCGGCAAGATCAGGGGCAAGCCGGAAAAGCCGGCCGCCt < 6:50827/136‑1 (MQ=255) aTTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCACAGCCGgccgccg > 2:243994/1‑136 (MQ=255) ttGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTc > 5:154379/1‑137 (MQ=255) ttGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGGCGCGGCCGCCATCCCGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTc < 6:265638/137‑1 (MQ=255) atgCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACt > 1:22407/1‑137 (MQ=255) cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGt > 8:119936/1‑137 (MQ=255) cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGt > 8:182869/1‑137 (MQ=255) cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGt > 8:304637/1‑137 (MQ=255) cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGt > 6:267917/1‑137 (MQ=255) gCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCg > 6:182199/1‑137 (MQ=255) gCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCg > 6:17128/1‑137 (MQ=255) cTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGc > 1:58055/1‑137 (MQ=255) tCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGCGCCAGCAGCGCTTGGGTGTCGCGGCCGCCATCACGGGCCGGCAAGATCAGTGCCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGcc < 6:207809/137‑1 (MQ=255) ggCGAGCCGGTAAAAACAATGACGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTAGTCGCCGACCGTGTGCgc > 4:29473/1‑136 (MQ=255) gCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCgag < 5:228193/137‑1 (MQ=255) gCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGGCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCACGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCgag < 6:243956/137‑1 (MQ=255) gCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGt > 6:254444/1‑137 (MQ=255) gCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGt > 6:162130/1‑137 (MQ=255) gCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGt > 8:252856/1‑137 (MQ=255) aaaaCAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGAAATATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTccc < 6:191647/137‑1 (MQ=255) cGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGa > 3:245521/1‑137 (MQ=255) tGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACg < 1:177141/137‑1 (MQ=255) ggctggCGAGCCAGCAGGGCTTGCGGGGCGCGGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACg < 2:94275/133‑1 (MQ=255) aGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGcc > 5:36410/1‑136 (MQ=255) aGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGACGGCCGCCTCGCCATACATGTAGCCGACCGTGTGCGAGGAGTACTTCCCTACCTGCAGACGGcc > 1:288554/1‑136 (MQ=255) aGGCGACCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCATCCCTACCTGCTGACGGCa > 7:134788/1‑137 (MQ=255) ccAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTg < 7:98639/137‑1 (MQ=255) cagcagGGCTTGGGTGGCGCGGCCGCCATCACGGGCCGGCAATATCAGTGCCACGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGa < 6:259470/137‑1 (MQ=255) cagcagCGCTTGGGTGTCGGTGCCGCAATCCCGGGCCGGCAAGATCTGTGGCAAGCCGGAAAAGCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGa < 8:84233/137‑1 (MQ=255) cagcagCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGa < 2:253228/137‑1 (MQ=255) agcGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACAcc < 6:154379/137‑1 (MQ=255) gCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCg < 4:89839/137‑1 (MQ=255) ggccgccATCAGGGGCCGGCAAGATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCa < 4:245521/136‑1 (MQ=255) gccgccATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCAt < 5:254444/137‑1 (MQ=255) gccgcaATCCGGGGCCGTAAAGTTCAGTGGAACGCCGGAAAAGCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCAt < 2:271144/137‑1 (MQ=255) atacggggCCGGAAATATCTGGGCCACCCCGGAAATGCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCCAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 5:17128/132‑1 (MQ=255) agcccGGGGCGGCTAGATCATTGGAAAGCCGGAAAAGCGGCCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 4:253407/133‑1 (MQ=255) agaaggggcgggcAAGTTCTGTGGCACGCGGGCAATGCCGGCCGCCTCGACATACTTGTCGCCGACCGGGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 4:131609/132‑1 (MQ=255) accACGGGCCGGAAAGATCAGTGCCACGCGGCAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 3:8916/135‑1 (MQ=255) aTCAGGGGCGGCCAATATCTGTGGAAAGCCGCAAAAGCCGCCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 7:119936/137‑1 (MQ=255) aTCAGGGGCCGGCAAGATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 5:60767/137‑1 (MQ=255) aTCAGGGGCCGGAAATATCAGTGCCAACCCGCAAACGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 3:36543/137‑1 (MQ=255) aTCAGGGGCCGCCAATATCCGTGCCACGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 5:182199/137‑1 (MQ=255) aTCACGGGCCGGCAAGATCTGTGCCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 5:267917/137‑1 (MQ=255) aTCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 2:44900/137‑1 (MQ=255) aTCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 2:27981/137‑1 (MQ=255) aTCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 4:111467/137‑1 (MQ=255) aTCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 8:289118/137‑1 (MQ=255) aTCACGGGCCGGCAAGATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 7:182869/137‑1 (MQ=255) aTCACGGGCCGGCAAGATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 2:65559/137‑1 (MQ=255) aTCACGGGCCGGCAAGATCAGTGGAAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 7:225198/137‑1 (MQ=255) aTCACGGGCCGGCAAGATCAGGGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc < 8:35997/137‑1 (MQ=255) cACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTc < 2:245940/137‑1 (MQ=255) ggCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCtt < 8:281033/137‑1 (MQ=255) ggCAAGAGCAGTGGCACGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGAGCTCATGGGCtt < 8:225286/137‑1 (MQ=255) | ACACTGGGGGCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTT > ACB122/6038522‑6038772 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |