Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,130,138 | G→T | 9.8% | intergenic (‑193/+111) | gcvH‑I ← / ← PP_0990 | glycine cleavage system protein H/DoxX family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,130,138 | 0 | G | T | 9.8% | 152.5 / 6.4 | 51 | intergenic (‑193/+111) | gcvH‑I/PP_0990 | glycine cleavage system protein H/DoxX family protein |
Reads supporting (aligned to +/- strand): ref base G (23/23); new base T (3/2); total (26/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 = 7.19e-01 |
GCTCTTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACG > NC_002947/1129991‑1130283 | gcTCTTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAg > 1:118951/1‑148 (MQ=255) ctctTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCa < 1:293496/149‑1 (MQ=255) aGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCaa < 1:329403/148‑1 (MQ=255) tAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGATGCTCAAACCCt < 6:135413/149‑1 (MQ=255) ccTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTtgt < 4:144973/144‑1 (MQ=255) ccTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTtgt > 3:144973/1‑144 (MQ=255) ccTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAc‑‑c > 8:277432/1‑149 (MQ=255) ccTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAc‑‑c > 8:20606/1‑149 (MQ=255) ttAGCAATTATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCa > 1:213743/1‑133 (MQ=255) ttAGCAATTATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCa < 2:213743/133‑1 (MQ=255) aataatGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATc > 7:132094/1‑127 (MQ=255) aataatGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATc < 8:132094/127‑1 (MQ=255) tGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCgagg‑ag < 7:20606/149‑1 (MQ=255) cGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGagg‑agg < 2:118951/148‑1 (MQ=255) ggCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCt < 6:214966/149‑1 (MQ=255) ccAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCtgt < 7:277432/149‑1 (MQ=255) aaCCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTgtgt < 4:167306/149‑1 (MQ=255) aCCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTa > 4:117933/1‑149 (MQ=255) cTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCATCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTATATc > 5:321761/1‑149 (MQ=255) aaTTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTc < 2:210096/149‑1 (MQ=255) aaTAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCg > 2:320856/1‑148 (MQ=255) ttACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGa < 3:103916/149‑1 (MQ=255) tACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGa > 1:61706/1‑148 (MQ=255) gTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCTCGCCAGTCCAGGCACCACATGTATCGTTTTCGCTACGTCGTAGCGAAAACGa > 5:117829/1‑149 (MQ=18) aaTTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACAtg < 2:61706/148‑1 (MQ=255) gCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTg > 3:218890/1‑149 (MQ=255) cGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGc < 5:1369/149‑1 (MQ=255) tatttattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGa < 3:117933/149‑1 (MQ=255) atttattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGtttt > 1:298094/1‑102 (MQ=255) atttattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGtttt < 2:298094/102‑1 (MQ=255) tattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGGGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTg > 1:104572/1‑149 (MQ=255) tattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTg > 3:292164/1‑149 (MQ=255) tattGTGTGGAATGATATCGGTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGt < 5:26168/148‑1 (MQ=255) atCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGt > 4:232683/1‑148 (MQ=255) tCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCtt > 7:97090/1‑143 (MQ=255) tCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCtt < 8:97090/143‑1 (MQ=255) cGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAt > 2:281762/1‑149 (MQ=255) gTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtc > 4:156355/1‑149 (MQ=255) gTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtc > 4:151348/1‑149 (MQ=255) gTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACttt < 4:218890/149‑1 (MQ=255) ctctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTcag > 2:89497/1‑148 (MQ=255) agagGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAgat < 1:89497/148‑1 (MQ=255) gAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATc > 2:173303/1‑148 (MQ=255) acaaGGGTTTGATGCTCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATc > 3:314056/3‑147 (MQ=255) acaaGGGTTTGATGCTCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATc > 8:31967/3‑147 (MQ=255) gggTTTGATGCTCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAt < 7:236429/123‑1 (MQ=255) gggTTTGATGCTCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAt > 8:236429/1‑123 (MQ=255) gTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGaaa < 1:281762/149‑1 (MQ=255) gTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGaa > 5:8036/1‑148 (MQ=255) gAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCAc < 5:69239/122‑1 (MQ=255) gAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCAc > 6:69239/1‑122 (MQ=255) gAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAAcg < 6:8036/148‑1 (MQ=255) | GCTCTTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACG > NC_002947/1129991‑1130283 |
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: ‑ |