Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,866,752 | G→T | 12.7% | intergenic (+19/+38) | PP_5140 → / ← thyA | MerR family transcriptional regulator/thymidylate synthase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,866,752 | 0 | G | T | 12.7% | 145.1 / 10.1 | 55 | intergenic (+19/+38) | PP_5140/thyA | MerR family transcriptional regulator/thymidylate synthase |
Reads supporting (aligned to +/- strand): ref base G (24/24); new base T (4/3); total (28/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 = 5.18e-01 |
TGGTGGAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCA > NC_002947/5866603‑5866892 | tggtggAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTg > 1:230621/1‑150 (MQ=255) ggggAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGAGCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTGGCGGGTGaa > 2:1015538/3‑150 (MQ=255) gcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAagag > 1:591669/1‑150 (MQ=255) ggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGt < 2:230621/150‑1 (MQ=255) ggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGt < 2:85041/150‑1 (MQ=255) tGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGtt > 1:1089617/1‑150 (MQ=255) gCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTc < 1:741976/150‑1 (MQ=255) aaCGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTa > 2:979819/1‑150 (MQ=255) cgcgCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAAc < 2:646645/150‑1 (MQ=255) ggggAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATaga < 1:1078250/150‑1 (MQ=255) ggggAGGAAAGGCGGGATCTGGCAGAAACTGGAGACAAACGGGGCGGTATAGGTACCGGATATCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATaga < 2:107739/150‑1 (MQ=255) gAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGAt > 1:1061616/1‑150 (MQ=255) aGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATag > 2:1173034/1‑145 (MQ=255) aGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATag < 1:1173034/145‑1 (MQ=255) gcgcgcTCTTGCAGCAACTGGAGGCAAACGTGGCGGTATCGGTACCGGATACCGGACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCGTCGCGGGTGACCCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGa < 1:891446/150‑1 (MQ=255) cTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAt < 1:448430/150‑1 (MQ=255) gcaACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGgcgc < 2:1061616/150‑1 (MQ=255) caACTGGAGACAAACGGGGCGGGATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCt < 2:488422/150‑1 (MQ=255) gggCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTAAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGACGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGt < 1:428820/150‑1 (MQ=255) gTATCCGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCa < 1:1022944/150‑1 (MQ=255) tATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCAt > 2:1179270/1‑150 (MQ=255) ggTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCCGCGaa < 1:694499/78‑1 (MQ=255) ggTAACGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCCGCGaa > 2:694499/1‑78 (MQ=39) gTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCt < 2:1089617/150‑1 (MQ=255) cGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTcca < 2:195585/143‑1 (MQ=255) cGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTcca > 1:195585/1‑143 (MQ=255) aaCATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCcgc < 2:758816/61‑1 (MQ=39) aaCATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCcgc < 2:758793/61‑1 (MQ=39) aaCATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCcgc > 1:758816/1‑61 (MQ=255) aaCATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCcgc > 1:758793/1‑61 (MQ=255) aaCATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAaga < 1:796363/66‑1 (MQ=255) aaCATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAaga > 2:796363/1‑66 (MQ=255) ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTccagcca < 1:38532/138‑1 (MQ=255) ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTccagcca > 2:38532/1‑138 (MQ=255) ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCgg > 1:171555/1‑150 (MQ=255) tgtgGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTc > 2:329631/1‑150 (MQ=255) ggggCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAAt < 1:329631/150‑1 (MQ=255) gTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACTCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCGTCGGCGCTGTCATCGGCGCATGGTGCTCATCGCCTTCCAGCCAGAAATCACTCGGCTCAAGCTTCTCCa < 1:1015538/150‑1 (MQ=255) tCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAg > 1:253139/1‑150 (MQ=255) cACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGc < 1:979819/150‑1 (MQ=255) aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAGCCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAAtcttc < 1:67111/143‑1 (MQ=255) aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAAtcttc > 2:67111/1‑143 (MQ=255) aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGcc > 2:1271967/1‑150 (MQ=255) aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTAGAGcc > 1:983337/1‑150 (MQ=255) aTTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCg > 2:814057/1‑150 (MQ=255) ttGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCgg > 1:180333/1‑150 (MQ=255) aaaaGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATAc > 2:208873/1‑150 (MQ=255) aaaaGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATAc > 2:208786/1‑150 (MQ=255) aaGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACtt < 2:253139/150‑1 (MQ=255) ctcTTCGCGGGTTCACCCGCGAAGCGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTg > 1:232300/3‑150 (MQ=255) gCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGc > 1:528879/1‑150 (MQ=255) ttCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCg < 2:1028941/112‑1 (MQ=255) ttCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCg > 1:1028941/1‑112 (MQ=255) ttCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCACTCTTCTCCAGCCACTCCGGCTCATACTTGCCa > 2:25593/1‑150 (MQ=255) gTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTAt < 1:256408/41‑1 (MQ=255) gTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTAt > 2:256408/1‑41 (MQ=255) | TGGTGGAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCA > NC_002947/5866603‑5866892 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |