Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,660,295 | C→G | 8.7% | intergenic (+97/‑132) | metK → / → ligB | methionine adenosyltransferase/DNA ligase B |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,660,295 | 0 | C | G | 8.7% | 141.5 / 3.4 | 46 | intergenic (+97/‑132) | metK/ligB | methionine adenosyltransferase/DNA ligase B |
Reads supporting (aligned to +/- strand): ref base C (21/21); new base G (2/2); total (23/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01 |
CCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGCCTACGCTAAGCGCTTCCCCGCCATGCAAGGATGCCAGGCCATGCCGTAC > NC_002947/5660146‑5660435 | ccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGcgcc < 1:558401/150‑1 (MQ=255) tGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGcgccgc < 2:602408/150‑1 (MQ=255) cACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCGGCGCTCgatc < 1:544380/150‑4 (MQ=255) accgaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCCGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTAAGATCGAGCGCCGCGCTCGttt > 1:339506/1‑149 (MQ=255) accgaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTc > 1:729956/1‑150 (MQ=255) accgaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTAAGATCGAGCGCCGCGCTCGttt < 2:339506/149‑1 (MQ=255) cgaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCGGCGCTCgatctga < 1:106497/150‑7 (MQ=255) cgcgCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTc > 1:301307/1‑150 (MQ=255) cAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGtctc < 2:660942/150‑1 (MQ=255) gACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGAc < 1:1108871/150‑1 (MQ=255) cGCTGCCGGCCTGTAAGCTTCCTGCAGCCCCAGACGCAAGCCTCTGCCGAGCGATCGGCAGGGGCTTTTTTCCGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGa < 1:870069/150‑1 (MQ=255) cGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCa < 2:840016/150‑1 (MQ=255) gCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGc > 1:1061620/1‑150 (MQ=255) tGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTg > 1:1061458/1‑116 (MQ=255) tGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTg < 2:1061458/116‑1 (MQ=255) tAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGaaa < 2:562162/149‑1 (MQ=255) tAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGaaa > 1:562162/1‑149 (MQ=255) tAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAat < 2:729956/150‑1 (MQ=255) ttCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCAc > 1:842092/1‑150 (MQ=255) ttCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCAc < 2:842092/150‑1 (MQ=255) tCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATCTCAcc < 2:301307/150‑1 (MQ=255) tCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCAcc < 1:18102/150‑1 (MQ=255) tGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGa < 2:626489/150‑1 (MQ=255) gCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCg < 2:915555/114‑1 (MQ=255) gCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCg > 1:915555/1‑114 (MQ=255) cAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACa > 2:479455/1‑150 (MQ=255) gAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTgcacgcac > 1:998451/1‑150 (MQ=255) gAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTgcacgcac > 1:1145970/1‑150 (MQ=255) gCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGcccgcccc < 1:413268/150‑1 (MQ=255) gATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCAccc > 1:431378/1‑150 (MQ=255) tCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGc < 2:652236/150‑1 (MQ=255) gggCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGCCTACGCt > 2:109306/1‑150 (MQ=255) ggCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGCCTACGCTa > 1:959032/1‑150 (MQ=255) cTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGCCTACGCTAAg > 1:164354/1‑150 (MQ=255) ttttttGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGCCTACGCTAAgc > 1:659786/1‑150 (MQ=255) gCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGCCTACGCTAAGCGCTTcccc > 1:883598/1‑150 (MQ=255) gCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGCCTACGCTAAGCGCTTcccc < 1:479455/150‑1 (MQ=255) aCGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGCCTACGCTAAGCGCTTCCCCGCCATGc < 1:109306/150‑1 (MQ=255) ttCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGCCTACGCTAAGCGCTTCCCCGCCATGCAAGGATGc < 1:1078707/150‑1 (MQ=255) cATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGCCTACGCTAAGCGCTTCCCCGCCATGCAAGGATGCCa < 2:959032/150‑1 (MQ=255) aTCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCgcccgccc < 2:609033/103‑1 (MQ=255) aTCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCgcccgccc > 1:609033/1‑103 (MQ=255) aTCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGCCTACGCTAAGCGCTTCCCCGCCATGCAAGGATGCCAg > 1:297596/1‑150 (MQ=255) cccGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGCCTACGCTAAGCGCTTCCCCGCCATGCAAGGATGCCAGGCCa > 2:813434/1‑150 (MQ=255) cGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACCACTGCACGCCCCGCCCGCCCCGCCACCCGCCTACGACAAGAGCATCCACGCCCGGCAAGGAAGCCAGGCCATg < 2:980096/150‑1 (MQ=255) gaaacGAGCGCGGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATTACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGCCTACGCTAAGCGCTTCCCCGCCATGCAAGGATGCCAGGCCATGCCGt > 2:421800/5‑150 (MQ=255) aaacGAGCGCGGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGCCTACGCTAAGCGCTTCCCCGCCATGCAAGGATGCCAGGCCATGCCGTa > 1:299060/4‑150 (MQ=255) aTCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGCCTACGCTAACCGCTTCCCCGCCATGCAAGGATGCCAGGCCATGCCGTAc > 1:969066/1‑150 (MQ=255) | CCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCGCCTACGCTAAGCGCTTCCCCGCCATGCAAGGATGCCAGGCCATGCCGTAC > NC_002947/5660146‑5660435 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |