Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,649,395 | G→C | 7.7% | intergenic (+872/+819) | PP_3214 → / ← PP_3216 | hypothetical protein/cysteine dioxygenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,649,395 | 0 | G | C | 7.7% | 163.4 / 4.7 | 52 | intergenic (+872/+819) | PP_3214/PP_3216 | hypothetical protein/cysteine dioxygenase |
Reads supporting (aligned to +/- strand): ref base G (22/26); new base C (2/2); total (24/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01 |
GACCACAGCGACCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGAGCCAGTTGCT > NC_002947/3649250‑3649535 | gACCACAGCGACCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGggc > 1:147809/1‑149 (MQ=255) cacaGCGACCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGggcggc > 6:36577/1‑149 (MQ=255) gACCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCg < 6:222081/149‑1 (MQ=255) aCCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCgg < 5:36577/149‑1 (MQ=255) cTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCa < 8:1372/149‑1 (MQ=255) cTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCa < 2:147809/149‑1 (MQ=255) ccAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGTCGGGCACCGGGCGGCTGCCCGGCAGAt > 3:295573/1‑149 (MQ=255) ccAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCAGCCGCCCGGTGCCCGGcgg < 2:123321/147‑3 (MQ=255) ccAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGAt > 3:94068/1‑149 (MQ=255) cccATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGAcc > 7:164382/1‑149 (MQ=255) tCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGAGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGgcg > 8:259232/1‑149 (MQ=255) agccaCTAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTc < 6:247678/149‑1 (MQ=255) ccaCGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGTCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCt < 4:295573/149‑1 (MQ=255) caCGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTc < 8:164382/149‑1 (MQ=255) aCGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGt > 3:205458/1‑141 (MQ=255) aCGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGt < 4:205458/141‑1 (MQ=255) aGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTc < 7:135427/149‑1 (MQ=255) gcTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCa > 7:213566/1‑119 (MQ=255) gcTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCa < 8:213566/119‑1 (MQ=255) acgatGGCGTGCGCGCCATCGTCGTCCACCAGCGCGATGCGCCCCCCGCGC‑ACACTGGCGACGAGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCCGCGCTGTATTCCTCGGTCCTACGCgg > 1:308120/5‑149 (MQ=255) ggCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGc > 5:334435/1‑149 (MQ=255) gcgcgcACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCt < 6:334435/149‑1 (MQ=255) cgcgcACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTa < 3:74271/149‑1 (MQ=255) aCGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGcc > 2:64733/1‑149 (MQ=255) gCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCg > 5:247546/1‑149 (MQ=255) aTCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGc < 6:247546/149‑1 (MQ=255) tcgtcCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTgg < 7:259232/149‑1 (MQ=255) caccaGCGCGATGCGCGCCCCGCGCAACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCg > 1:3337/1‑113 (MQ=255) caccaGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCg < 2:3337/112‑1 (MQ=255) caccaGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCcggc < 1:64733/149‑1 (MQ=255) caGCGCGATGCGCGCGCCGCGC‑ACTCTGGGGACGTGGTGGGTCTCTTGAACCAGTTGCCGGCCGGGGAGGGGGCGGCTGCCCGGGAGATGACCGGCGCTGTATTCCTCGGTCGTACGGGGTGTGCTACCGCCCGGGCGTGGCCggcggc < 2:308120/149‑1 (MQ=255) aGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGAc < 1:316161/91‑1 (MQ=255) aGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGAc > 2:316161/1‑91 (MQ=255) gcgcCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCggg < 3:80499/149‑1 (MQ=255) ccGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTAcc > 7:202439/1‑113 (MQ=255) ccGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTAcc < 8:202439/113‑1 (MQ=255) acacTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCAGCCGCCCGGTGCCCGGcg < 4:38634/64‑2 (MQ=21) acacTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCAGCCGCCCGGTGCCCGGcg > 3:38634/1‑63 (MQ=21) tGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCt > 8:43606/1‑149 (MQ=255) gCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTc < 3:307932/148‑1 (MQ=255) gACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTg > 6:1270/1‑110 (MQ=255) gACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTg < 5:1270/110‑1 (MQ=255) ctTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCt > 8:39791/1‑149 (MQ=255) ttGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGcc > 7:146950/1‑93 (MQ=255) ttGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGcc < 8:146950/93‑1 (MQ=255) ttGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAgcg < 4:206248/140‑1 (MQ=255) ttGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAgcg > 3:206248/1‑140 (MQ=255) gAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGa < 3:88391/148‑1 (MQ=255) aGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGcccct > 2:278286/1‑145 (MQ=255) aGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGa < 7:39791/149‑1 (MQ=255) cccGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGAGCCAg > 5:154282/1‑149 (MQ=255) gccgGGCAGCCGCCCGGTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAgcg > 4:329142/1‑125 (MQ=255) gggCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGAGCCAGTtgct < 7:171019/141‑1 (MQ=255) gggCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGAGCCAGTtgct > 8:171019/1‑141 (MQ=255) | GACCACAGCGACCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGAGCCAGTTGCT > NC_002947/3649250‑3649535 |
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: ‑ |