Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,221,231 | G→T | 21.0% | intergenic (+46/+216) | PP_3699 → / ← PP_3700 | hypothetical protein/chromosome partitioning ATPase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,221,231 | 0 | G | T | 21.0% | 53.8 / 8.3 | 29 | intergenic (+46/+216) | PP_3699/PP_3700 | hypothetical protein/chromosome partitioning ATPase |
Reads supporting (aligned to +/- strand): ref base G (12/10); new base T (3/3); total (16/13) | |||||||||||
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.66e-01 |
ACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACAT > NC_002947/4221098‑4221371 | aCGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccccctga < 2:618813/150‑5 (MQ=255) gTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctc > 1:497148/1‑150 (MQ=255) cGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcag < 1:561088/150‑1 (MQ=255) gACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCg > 1:1083941/1‑124 (MQ=255) gACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCg < 2:1083941/124‑1 (MQ=255) tGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAATGAATTGTCAAATGCTGACCGTCGTGCATATCGCGCTCAGTGGGCATCTTACAGATGCCCCAAccccatcctcagcctca > 1:248611/1‑150 (MQ=255) gCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGTGGGCATCg < 1:365719/112‑1 (MQ=255) gCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGTGGGCATCg > 2:365719/1‑112 (MQ=255) gCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCaa < 1:1044250/150‑1 (MQ=255) gATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTa < 1:272962/150‑1 (MQ=255) cTGAAGCTTTGATCAGGGCTTATATGCTTGGGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGGTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTa > 2:486155/1‑150 (MQ=255) cAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACctcagcctcagcctcagc > 1:1031678/1‑122 (MQ=255) cAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACctcagcctcagcctcagc < 2:1031678/122‑1 (MQ=255) cAGGGCTGATATGCTTGAGGAGTTTCAACAAGACTGCAAGTGAAGTGGCAAATGCTAACCGTCGTACATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCCGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAgcgcgc > 2:867294/1‑150 (MQ=255) ggCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccccct < 1:714258/100‑3 (MQ=255) ggCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccccct > 2:714258/1‑98 (MQ=255) ggCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGa > 2:631057/1‑150 (MQ=255) ttGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACctcagcctcagcctcagc > 1:726474/1‑108 (MQ=255) ttGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACctcagcctcagcctcagc < 2:726474/108‑1 (MQ=255) ttGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCCAGCAGAGGCCCCAACCTCAGCCTCAGCGTCAGCATCAAAACCTAGCGATGGCCGAGCGCGCCGAAGCCACAAtct > 1:870391/1‑150 (MQ=255) gaggagTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCa < 1:631057/150‑1 (MQ=255) gTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCaa > 2:729855/1‑150 (MQ=255) aaGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACctcagcctcagcctcagc > 1:730244/1‑93 (MQ=255) aaGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACctcagcctcagcctcagc < 2:730244/93‑1 (MQ=255) aaGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGcc > 1:674245/1‑150 (MQ=255) gtgaagtgGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGc < 1:794733/96‑1 (MQ=255) gtgaagtgGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGc > 2:794733/1‑96 (MQ=255) tgaagtgGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCgg < 1:729855/150‑1 (MQ=255) gCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTa > 1:404463/1‑150 (MQ=255) tgaggttggggCATCTGTCCGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATg > 1:586478/8‑147 (MQ=255) tgaggttggggCATCTGTCCGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATg < 2:586478/140‑1 (MQ=255) cgcTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGccc < 2:497148/150‑1 (MQ=255) gggggCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACAt > 2:51147/1‑150 (MQ=255) | ACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACAT > NC_002947/4221098‑4221371 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |