Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | AE015451 | 1,126,645:1 | +C | 100% | intergenic (‑99/+50) | tdcG‑II ← / ← gcvP‑I | L‑serine dehydratase/glycine dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AE015451 | 1,126,645 | 1 | . | C | 97.3% | 131.7 / ‑2.5 | 37 | intergenic (‑99/+50) | tdcG‑II/gcvP‑I | L‑serine dehydratase/glycine dehydrogenase |
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (18/18); total (18/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑A‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAG > AE015451/1126433‑1126862 | gCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑A‑c > 1:833035/1‑214 (MQ=255) gCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑A‑c < 2:833032/214‑1 (MQ=255) cTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGcc < 1:1041509/215‑1 (MQ=255) cAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCa < 1:123724/214‑1 (MQ=255) cGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGc < 2:945362/186‑1 (MQ=255) ggCCGACCGTGTGGGGGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑A‑cccg < 2:310214/170‑3 (MQ=255) tGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCCCCATCGCATTCGCGGGTAA‑ACCCGCTGCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgat < 1:543185/201‑1 (MQ=255) cAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGGTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑A‑cccgc < 2:79570/128‑4 (MQ=255) gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCg > 2:622196/1‑157 (MQ=255) gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGt > 1:382506/1‑235 (MQ=255) gACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGc < 2:785358/132‑1 (MQ=255) tCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACg < 1:477239/188‑1 (MQ=255) aGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCAACCAAGCCAGGTGGGGGGGAGGGCTa > 1:43446/1‑157 (MQ=255) aCATGTGGGCCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCCCCATCGCATTCGCGGGTAA‑ACCCGCTCCTac < 2:166250/114‑1 (MQ=37) cATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCa > 1:963880/1‑185 (MQ=255) tcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGCTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcg > 1:458547/1‑225 (MQ=255) gCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGCGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgcgcg > 2:1031456/1‑176 (MQ=255) ccAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGGATGCGGGGGGATAGCCACCATCGCATTCGCGGGGAA‑A‑cc > 2:135108/1‑89 (MQ=38) aGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGTGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGCGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCa > 1:367402/1‑179 (MQ=255) agtGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTc < 1:332662/229‑1 (MQ=255) ccTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGta > 1:707950/1‑235 (MQ=255) ggCGCGGGTTTCCCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAACC‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTc < 1:135110/218‑1 (MQ=255) gcgcGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAgg < 2:752931/223‑1 (MQ=255) ttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGATGGGCAGGCGCACACCCGGTTGCGGTCGCCGAACACGTTGTc > 2:17966/1‑166 (MQ=255) ttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTccacca > 1:610543/1‑204 (MQ=255) ttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTccacca < 2:610540/204‑1 (MQ=255) ttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGGGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAacac > 1:1086870/1‑160 (MQ=255) cgcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCggg > 1:201078/1‑84 (MQ=255) cgcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTcc < 2:963877/195‑1 (MQ=255) ggTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAg < 1:1073504/114‑1 (MQ=255) ccaccaGCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCtt < 2:123949/174‑1 (MQ=255) accaTCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc > 1:892064/1‑156 (MQ=255) accaTCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc < 2:892061/156‑1 (MQ=255) ccaTCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAAc < 1:718864/235‑1 (MQ=255) aTCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCa > 2:659676/1‑184 (MQ=255) aTCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCa < 1:659679/184‑1 (MQ=255) aTCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAg > 1:82204/1‑236 (MQ=255) aTCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGCGGCGGGGAGGGCTATCAGGCGTCCTGATAGCCCTCGACCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGCTGTCGACCCGGCCGACCGGCGGCCAgc > 2:264740/1‑159 (MQ=255) cGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgg < 2:34535/57‑1 (MQ=25) cGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgg > 1:34536/1‑57 (MQ=25) cGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCCCGCGGCTGTAGCCATGGGGCCATTcat > 2:718861/1‑225 (MQ=255) ggTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGc < 1:622199/189‑1 (MQ=255) | GCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑A‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAG > AE015451/1126433‑1126862 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |