Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | AE015451 | 2,087,428 | G→A | 9.6% | intergenic (+241/+677) | PP_1865 → / ← aroH | transposase/Phospho‑2‑dehydro‑3‑deoxyheptonate aldolase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AE015451 | 2,087,428 | 0 | G | A | 9.6% | 115.8 / 5.2 | 42 | intergenic (+241/+677) | PP_1865/aroH | transposase/Phospho‑2‑dehydro‑3‑deoxyheptonate aldolase |
Reads supporting (aligned to +/- strand): ref base G (15/22); new base A (2/2); total (17/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01 |
TCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAG > AE015451/2087210‑2087647 | ncccTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGa < 2:833055/235‑1 (MQ=255) gtggTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGa < 2:480875/227‑1 (MQ=255) aaaaCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgcg > 2:651302/1‑199 (MQ=255) ccccTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGa < 2:492120/182‑1 (MQ=255) tCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGCCTACCgcgcgcct > 2:12440/1‑163 (MQ=255) tCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCCCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCt < 1:84533/236‑1 (MQ=255) cGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGc < 1:425937/236‑1 (MQ=255) tCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGt < 1:685884/236‑1 (MQ=255) cccTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGt < 2:374294/234‑1 (MQ=255) gCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCA‑‑CCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATcggc < 2:404607/190‑1 (MQ=255) gAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCt < 2:778394/203‑1 (MQ=255) gcgcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAg < 1:1034973/131‑1 (MQ=255) cgcTGCAGGCCCCTGCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTc < 2:772363/204‑1 (MQ=255) cAGGCCCCATCC‑TGGGGCC‑CAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGcc > 1:62248/1‑236 (MQ=21) ccccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCg > 1:791358/1‑213 (MQ=255) catccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCAcc < 2:696189/161‑1 (MQ=255) catGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGa > 1:179417/1‑96 (MQ=255) catGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTAc > 1:367740/1‑235 (MQ=255) tGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGc < 2:479738/148‑1 (MQ=255) gACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTgg < 1:46440/221‑1 (MQ=255) cGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCCTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGTCCGGCACCGCCGATGCATTCGGCGAGGGTCTGCTCGGTCATCTGGATAAGGATATCGCCGGCATGGc > 1:384128/1‑196 (MQ=255) cccGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGt > 2:618701/1‑155 (MQ=255) cccGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGt < 1:618704/155‑1 (MQ=255) aTCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCtt > 1:604419/1‑235 (MQ=255) tctcCTTCCGCCATGCATCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGACTGGAATACCGCGCACTTAGAAGCACGa < 1:123018/80‑1 (MQ=11) tccGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAACCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTAGCTCTTGCTTTgg < 1:564484/231‑1 (MQ=255) cGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTgg < 2:543428/192‑1 (MQ=255) aaTGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGgac > 1:213463/1‑190 (MQ=255) tCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGt < 1:35994/131‑1 (MQ=255) tCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAAc > 2:819564/1‑189 (MQ=255) ccACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAgc < 2:604416/234‑1 (MQ=255) ccACTCGGCCGCCTGAAGTCGCAATCGGCGTCGCCTGTACTCCCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAgcgc < 1:12440/236‑1 (MQ=255) cTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTtgctcttgctcttgctct < 2:433771/208‑1 (MQ=255) cTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTtgctcttgctcttgctct > 1:433774/1‑208 (MQ=255) tctcctGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGc < 1:790496/181‑1 (MQ=255) tcctGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGCAGCAAGATCaa > 2:98349/1‑56 (MQ=255) ctGAAGTCTCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGt < 2:930259/118‑1 (MQ=255) cAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCa > 1:1016825/1‑235 (MQ=255) aTCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAg < 2:367737/234‑1 (MQ=255) cTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCAGCTCGGTCATCTGGATAAGGATATCGCCGGCATGGCt > 1:803764/1‑135 (MQ=255) cTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCt > 1:168779/1‑196 (MQ=255) tCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTATTGCTTCTAAGCGCGCGGTAgg > 1:911414/1‑213 (MQ=255) | TCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAG > AE015451/2087210‑2087647 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |