Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | AE015451 | 4,022,292 | C→G | 12.5% | C5S (TGC→TCC) | PP_3547 ← | putative NAD(P)‑binding oxidoreductase with NAD(P)‑binding Rossmann‑fold domain |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AE015451 | 4,022,292 | 0 | C | G | 12.5% | 116.4 / 11.0 | 40 | C5S (TGC→TCC) | PP_3547 | putative NAD(P)‑binding oxidoreductase with NAD(P)‑binding Rossmann‑fold domain |
Reads supporting (aligned to +/- strand): ref base C (11/24); new base G (2/3); total (13/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01 |
TGGCGGCGCCGATACCGCGGCTGCCGCCGGTGATGACAATGACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTG > AE015451/4022059‑4022521 | tGGCGGCGCCGATACCGCGGCTGCCGCCGGTGATGACAATGACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAg < 1:537827/236‑1 (MQ=255) gcgcCGATACCGCGGCTGCCGCCGGTGATGACAATGACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGa < 1:709534/235‑1 (MQ=255) gATACCGCGGCTGCCGCCGGTGATGACAATGACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAg < 2:713967/226‑1 (MQ=255) cTGCCGCCGGTGATGACAATGACGTGCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATgg < 1:207133/229‑1 (MQ=255) caccgaggatgcgttgtctgaccggtatcacacccattgtgatccgcggatgaatgccgatcagtcgctggagttggcgttcctgattgccgaaacccacaagcaggtacggcgctaagcgcttttcgcttttgatcgtgcttctaagtgcgcggtagttcaggcaacgcagattgcgacttcaggaggccgagcgcagggcttgcggAGGGAGGTGACGGGCATGGATGCCCg > 2:384363‑M2/209‑234 (MQ=255) cgGTGATGACAATGACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgcg < 1:392788/236‑1 (MQ=255) ggTGATGACAATGACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCATCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgc < 2:310122/234‑1 (MQ=255) gTGATGACAATGACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgcgc < 1:306577/235‑1 (MQ=255) gTGATGACAATGACGTTATCCAGGGGTTTGCCAATGGGTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgcgc < 2:230615/235‑1 (MQ=255) tgatgaCAATGACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgcgc < 1:771363/234‑1 (MQ=255) cAATGACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAgg < 2:480860/233‑1 (MQ=255) ggtatcacacccattgtgatccgcggatgaatgccgatcagtcgctggagttggcgttcctgattgccgaaaccctcaagcaggtacggcgctaagcgcttttcgcttttgatcgtgcttctaagtgcgcggtagttcaggcaacgcagattgcgacttcaggaggccgagtggaggcattgcggAGGGAGATGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCccatg < 1:362235‑M2/46‑3 (MQ=255) gtatcacacccattgtgatccgcggatgaatgccgatcagtcgctggagttggcgttcctgattgccgaaacccacaagcaggtacggcgctaagcgcttttcgcttttgatcgtgcttctaagtgcgcggtagttcaggcaacgcagattgcgacttcaggaggccgagcgcagggcttgcggAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATGGGGc < 1:384365‑M2/52‑1 (MQ=255) cAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACgg > 2:789408/1‑183 (MQ=255) cAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGTGGCCTTCGGcgcg > 2:699713/1‑233 (MQ=255) cAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGAACTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCAGAGGCAGGTGACGGGCATCCATGCCCGTCacctccgt > 1:526676/1‑199 (MQ=255) aTGGTTTGGGGAAGGGAACATCGACACTAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCATCCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACgg < 1:225636/181‑1 (MQ=255) tGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgcgc < 2:771649/201‑1 (MQ=255) tGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGa < 2:714721/213‑1 (MQ=255) cgcggatgaatgccgatcagtcgctggagttggcgttcctgattgccgaaaccctcaagcaggtacggcgctaagcgcttttcgcttttgatcgtgcttctaagtgcgcggtagttcaggcaacgcagattgcgacttcaggaggccgagtggaggcattgcggAGGGAGATGACggg > 1:479806‑M2/165‑178 (MQ=255) cgcggatgaatgccgatcagtcgctggagttggcgttcctgattgccgaaaccctcaagcaggtacggcgctaagcgcttttcgcttttgatcgtgcttctaagtgcgcggtagttcaggcaacgcagattgcgacttcaggaggccgagtggaggcattgcggAGGGAGATGACggg < 2:479804‑M2/14‑1 (MQ=255) agggaaCATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCaa < 1:643599/235‑1 (MQ=255) gggaaCAACGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTTATCTTGATCTTGATCTTCGCGACTTCATGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGAt < 2:689590/176‑1 (MQ=255) cAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGccc > 1:257980/1‑186 (MQ=255) tAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGa > 1:300348/1‑235 (MQ=255) cGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGACTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGGGGCAGGTGACGGGCATGGAt > 2:141335/1‑158 (MQ=255) tGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGagcgcagc > 1:799055/1‑235 (MQ=255) gTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCGTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATgg < 2:256958/142‑1 (MQ=255) ttGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAgg < 2:673397/129‑1 (MQ=255) cacaGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAgc < 2:300346/234‑1 (MQ=255) aCCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGgatg > 1:594550/1‑218 (MQ=255) ctctGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCgggt > 2:139472/1‑212 (MQ=255) gatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCg < 1:468548/228‑1 (MQ=31) gatcttCGCGACTTCAGGAGGCCGAACGCAGGCCGTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAggg < 1:16776/228‑1 (MQ=14) tgcttCGCGACTCCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCAGGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTGCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGagcgcagcg < 2:28146/147‑1 (MQ=12) ttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGagcgcagc < 2:595310/145‑1 (MQ=32) aaCGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGttt > 1:256036/1‑156 (MQ=31) ccTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGtt < 1:139474/236‑1 (MQ=38) gAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTg < 2:676321‑M2/233‑1 (MQ=255) gAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTg > 1:676323‑M2/2‑234 (MQ=255) | TGGCGGCGCCGATACCGCGGCTGCCGCCGGTGATGACAATGACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTG > AE015451/4022059‑4022521 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |