Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 461,983 | C→T | 9.7% | G162D (GGC→GAC) | pqqF ← | coenzyme PQQ synthesis protein F |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 461,983 | 0 | C | T | 9.7% | 107.1 / 3.3 | 41 | G162D (GGC→GAC) | pqqF | coenzyme PQQ synthesis protein F |
Reads supporting (aligned to +/- strand): ref base C (16/21); new base T (2/2); total (18/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 = 4.68e-01 |
CCTGCTGCCGGGCAAGGCGTTCAAGTTCATCCAGCGGCTGCGGACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGATGCCCAGATCCGG > NC_002947/461834‑462131 | ccTGCTGCCGGGCAAGGCGTTCAAGTTCATCCAGCGGCTGCGGACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGc < 1:929629/150‑1 (MQ=255) ctgctgCCGGGCAAGGCGTTCAAGTTCATCCAGCGGCTGCGGACCGCAAAGGCTCAGGCAGATGTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGCCCGCCTCGTGCAAGGCCAGGGTGTATGGATTGcc > 2:549572/1‑150 (MQ=255) ttCATCCAGCGGCTGCGGACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATg > 2:716680/1‑150 (MQ=255) ccAGCGGCTGCGGACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCCTGAAAAGCGCCCAGTGGATGCCTgg > 2:104306/1‑150 (MQ=255) ccAGCGGCTGCGGACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTgg < 1:146797/150‑1 (MQ=255) ggCTGCGGACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGac > 1:313996/1‑150 (MQ=255) tGCGGACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGACAAGCGCCCAGTGGGTGCCTGGCTGACact > 1:602901/1‑150 (MQ=255) ggACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCa < 2:471867/114‑1 (MQ=255) ggACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCa > 1:471867/1‑114 (MQ=255) gACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACactgact < 1:529099/150‑1 (MQ=255) cGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCATGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTgca < 1:1069660/150‑1 (MQ=255) cGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATTAAAAGCGCCCAGGGGATGCCTGGCTGACACTGACGgca > 1:269925/1‑150 (MQ=255) cAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGTCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTgcagc > 2:1031203/1‑150 (MQ=255) cAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGTCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTgcagc > 2:1031124/1‑150 (MQ=255) aaaGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTgcagca > 1:318300/1‑150 (MQ=255) aaGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGACAGCTAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCaa < 1:996526/150‑1 (MQ=255) ggCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCGCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCCCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAgg < 1:1037496/150‑1 (MQ=255) tCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGa > 1:610424/1‑150 (MQ=255) aGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAAc < 1:549572/150‑1 (MQ=255) gCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCGGCGAAAGCCTGGCTGACACTGACTGCAGCAAGGCGAActg > 1:66977/1‑150 (MQ=255) ccTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCACGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGTCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGtt < 1:1031203/150‑1 (MQ=255) ccTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGATTGCGGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGTCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGtt < 1:1031124/150‑1 (MQ=255) gTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGa < 2:1201161/97‑1 (MQ=255) gTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGa > 1:1201161/1‑97 (MQ=255) aGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTgc < 2:792784/150‑1 (MQ=255) cGTTGATGAAAGACCGCAAGTGTTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGAc < 2:66977/150‑1 (MQ=255) atgaAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGa > 1:35233/1‑96 (MQ=255) atgaAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGa < 2:35233/96‑1 (MQ=255) tgaAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGTAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCg > 1:300684/1‑150 (MQ=255) cAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGACTGCGCGACCAGGCGATGAATTCgg < 2:610424/150‑1 (MQ=255) gCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCa < 1:192745/150‑1 (MQ=255) tGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACt < 2:313996/150‑1 (MQ=255) ccGCATCGTGCAAGGCAAGCGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCATGGCTGACCCTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGtt < 2:269925/150‑1 (MQ=255) gCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGTGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCa < 1:1068634/150‑1 (MQ=255) tGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACggcgt < 2:339437/150‑1 (MQ=255) aaGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACggcgttgg > 1:110281/1‑150 (MQ=255) cAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACggcgt < 2:1113750/142‑1 (MQ=255) cAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACggcgt > 1:1113750/1‑142 (MQ=255) gtgtATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGATg > 1:1089860/1‑150 (MQ=255) aTTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGATGCCCAGATc > 1:794690/1‑150 (MQ=255) ttGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGATGCCCAGATcc < 1:803950/150‑1 (MQ=255) gCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGATGCCCAGATCCgg < 2:1089860/150‑1 (MQ=255) gCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGATGCCCAGATCCgg < 1:93560/150‑1 (MQ=255) ccAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGAt < 2:1234499/137‑1 (MQ=255) ccAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGAt > 1:1234499/1‑137 (MQ=255) | CCTGCTGCCGGGCAAGGCGTTCAAGTTCATCCAGCGGCTGCGGACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGATGCCCAGATCCGG > NC_002947/461834‑462131 |
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: ‑ |