Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,866,696 | G→T | 12.3% | 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,696 | 0 | G | T | 12.3% | 154.2 / 10.8 | 57 | intergenic (+19/+38) | PP_5140/thyA | MerR family transcriptional regulator/thymidylate synthase |
Reads supporting (aligned to +/- strand): ref base G (24/26); new base T (4/3); total (28/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.06e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01 |
GGAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCG > NC_002947/5866551‑5866842 | ggAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAAc < 5:217703/149‑1 (MQ=255) gcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAaga < 1:197763/149‑1 (MQ=255) gcgcTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGTGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGt < 7:212245/149‑1 (MQ=255) gcTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCa < 1:28376/149‑1 (MQ=255) tGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCcgcg < 5:178413/138‑1 (MQ=255) tGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCcgcg > 6:178413/1‑138 (MQ=255) tGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGt < 8:245417/145‑1 (MQ=255) tGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGt > 7:245417/1‑145 (MQ=255) aGGTGAGTGAATTCGCGTTCTTGCTGCAACTGGAGACACAGGGGGCGGTTTCGGTACCGGATACCGAACATTCCGATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGTTGAACCCGCGAAGA‑GGTCAGTTCTAACAATa < 6:11475/149‑1 (MQ=255) gcgATCTTGCAGCAACTGGAGCCAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGGTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGAc < 3:35037/149‑1 (MQ=255) aTCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGc > 7:67311/1‑149 (MQ=255) gcagcaACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCg > 2:58059/1‑149 (MQ=255) cagcaACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCgg < 2:196181/149‑1 (MQ=255) agcaACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGt > 3:208619/1‑111 (MQ=255) agcaACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGt < 4:208619/111‑1 (MQ=255) aaCTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACTGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCt < 3:160639/149‑1 (MQ=255) aaCTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCt > 7:46203/1‑149 (MQ=255) tGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCCGCGaa < 7:236445/100‑1 (MQ=255) tGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCCGCGaa > 8:236445/1‑100 (MQ=255) tGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGt > 5:37907/1‑148 (MQ=255) acaCGGGGCTGTCTCTGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCAGTGAGCCAAAAGCTTAGCGGGTGACCCCGCGAAGA‑GGTCAGTTCTAACAATATATATCTATCAGACCGCCATCGGCGCTGTCATCGgcg < 6:149451/147‑1 (MQ=255) aaaCGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAAGAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATAGgcg > 8:96996/1‑149 (MQ=255) aaaCGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACAGCCATCGGCGCTGTCATCGgcg > 5:227818/1‑149 (MQ=255) gCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTg < 6:37907/149‑1 (MQ=255) tATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCa < 8:46203/149‑1 (MQ=255) tACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACTGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCt < 3:235935/149‑1 (MQ=255) ggATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTcca < 8:233054/149‑1 (MQ=255) aCCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGcc < 5:198224/140‑1 (MQ=255) aCCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGcc > 6:198224/1‑140 (MQ=255) cGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTccagccag > 5:136548/1‑148 (MQ=255) cGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTccagccag > 6:49558/1‑148 (MQ=255) gAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGTCCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGaa < 1:70997/149‑1 (MQ=255) gAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCAAACAATa > 6:205999/1‑84 (MQ=255) gAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCAAACAATa < 5:205999/84‑1 (MQ=255) aaCATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGaaa < 6:227818/149‑1 (MQ=255) aTTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCa > 7:241913/1‑149 (MQ=255) tCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACt > 6:159922/1‑149 (MQ=255) cccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTc > 7:30685/1‑149 (MQ=255) ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTc < 8:241913/148‑1 (MQ=255) ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCg > 6:268046/1‑149 (MQ=255) ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCg > 8:203856/1‑149 (MQ=255) ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTCCTCATAGCCTTCTAGCCAGAAATCACTCg < 4:179712/149‑1 (MQ=255) ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTACTCATAGCCTTCTAGCCAGAAATCACTCg > 3:179712/1‑149 (MQ=255) cATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCgg > 2:84889/1‑149 (MQ=255) aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGc > 6:61562/1‑149 (MQ=255) cATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCATCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTc > 7:142493/1‑149 (MQ=255) ccAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCa > 6:8374/1‑149 (MQ=255) ccAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCa < 8:30685/149‑1 (MQ=255) aaaaGCTTCGCGGGTGAACCCGCGAAGAGGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCa > 1:152875/1‑148 (MQ=255) ctcTTCGCGGGTTCACCCGCGAAGC‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACtt > 1:157266/3‑149 (MQ=255) cTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGc < 3:148897/149‑1 (MQ=255) cTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGc < 7:203856/149‑1 (MQ=255) ttCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGc > 8:146764/1‑148 (MQ=255) tCGCGGGTTCACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGc > 2:196323/1‑147 (MQ=255) cgcgGGTTCACCCGCGAAGC‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCa < 2:157266/148‑1 (MQ=255) cgcgGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCa < 5:268046/148‑1 (MQ=255) gTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCg < 5:49558/149‑1 (MQ=255) | GGAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCG > NC_002947/5866551‑5866842 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |