Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 372,628 | T→C | 12.3% | intergenic (+48/‑131) | dgcA → / → dgcB | dimethylglycine dehydrogenase subunit/dimethylglycine dehydrogenase subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 372,628 | 0 | T | C | 12.3% | 100.6 / 5.9 | 41 | intergenic (+48/‑131) | dgcA/dgcB | dimethylglycine dehydrogenase subunit/dimethylglycine dehydrogenase subunit |
Reads supporting (aligned to +/- strand): ref base T (15/20); new base C (2/3); total (17/24) | |||||||||||
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.73e-01 |
GCGAAGGCTACCTGCTGTACCGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAAC > NC_002947/372483‑372767 | gCGAAGGCTACCTGCTGTACCGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTAcc < 3:205104/149‑1 (MQ=255) cTACCTGCTGTACCGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAAt < 8:206603/149‑1 (MQ=255) ccTGCTGTACCGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTAAACCcgctcctaca < 4:208707/149‑8 (MQ=255) ctgTACCGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAg > 7:221970/1‑149 (MQ=255) ccGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACt > 8:179749/1‑148 (MQ=255) ggCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGacaaca < 1:137172/148‑1 (MQ=255) ggCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGacaaca < 6:233908/148‑1 (MQ=255) gCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACt > 2:96346/1‑149 (MQ=255) cccAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTCACCcgc < 8:98641/113‑1 (MQ=255) cccAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTCACCcgc > 7:98641/1‑113 (MQ=255) ccAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCg > 2:167182/1‑148 (MQ=255) ccAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCg > 5:76903/1‑148 (MQ=255) aGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGAc < 7:179749/149‑1 (MQ=255) aaCGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAAGAGCGGGTTCACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCAtt < 1:110699/149‑1 (MQ=18) cGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTcg > 7:113695/1‑149 (MQ=255) gTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAAGAGCGGGTTCACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTcgc < 6:182367/149‑1 (MQ=18) cATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCgg < 6:242169/148‑1 (MQ=255) cATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCgg < 8:39926/148‑1 (MQ=255) gcgATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTc < 4:9850/102‑1 (MQ=255) gcgATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTc > 3:9850/1‑102 (MQ=255) cgacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTTACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGTGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCGGACCGGACGCATTCGCGGGTGAACCCGc < 5:232406/144‑1 (MQ=255) cgacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGctc < 8:221970/146‑1 (MQ=255) ttGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATgggg > 2:126508/1‑149 (MQ=255) gtgCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGtctc < 3:36853/149‑1 (MQ=255) tgCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGtctc < 6:76903/148‑1 (MQ=255) aCTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTg > 4:69405/1‑149 (MQ=255) ttCTGATCGTTAACAGTCCCCGCCATCGGGCGCGATCCCTCTAGCAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCCCGGGTGAACCCGCTCTTACAAGGGCCCAGTCTCGCCCTCTGGt > 2:118658/1‑149 (MQ=255) aTCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGgcg > 6:230665/1‑149 (MQ=255) ttAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGgcgc > 8:112191/1‑146 (MQ=255) tCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGc < 3:213219/50‑1 (MQ=255) tCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGc > 4:213219/1‑50 (MQ=255) ccATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGttt < 8:113695/149‑1 (MQ=255) cTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATg < 3:69405/149‑1 (MQ=255) tAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTg < 1:167182/149‑1 (MQ=255) tAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTg < 5:262598/149‑1 (MQ=255) cccccacGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGtt < 1:96346/142‑1 (MQ=255) aGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCggg > 6:35329/1‑68 (MQ=255) aGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCggg < 5:35329/68‑1 (MQ=255) gAGCGGGTTCACCCGCGAATGCGGGCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAac > 2:75457/1‑149 (MQ=18) aGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTccc < 7:33851/139‑1 (MQ=255) aGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTccc > 8:33851/1‑139 (MQ=255) | GCGAAGGCTACCTGCTGTACCGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAAC > NC_002947/372483‑372767 |
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: ‑ |