Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,431,901 | G→T | 15.1% | E167D (GAG→GAT) | PP_3047 → | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,431,901 | 0 | G | T | 15.1% | 60.6 / 2.4 | 26 | E167D (GAG→GAT) | PP_3047 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (9/12); new base T (2/2); total (12/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01 |
CGAAGGCTCGTCGCCAGTTGGACGCTACGCATGAGTTTCATCGCATGGGCGCGCTGAACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAG > NC_002947/3431762‑3432050 | cGAAGGCTCGTCGCCAGTTGGACGCTACGCATGAGTTTCATCGCATGGGCGCGCTGAACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTgg < 1:112668/150‑1 (MQ=255) aaGGCTCGTCGCCAGTTGGACGCTACGCATGAGTTTCATCGCATGGGCGCGCTGGACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGATATGGAGCTGGGt > 1:667629/1‑150 (MQ=255) ggCTCGGCGCCAGTTGGACGCTACGCATGAGTTTCATCGCATGGGCGCGCTGAACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTAc < 2:322678/150‑1 (MQ=255) gCATGAGTTTCATCGCATGGGCGCGCTGAACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGaa < 2:152773/150‑1 (MQ=255) tGAGTTTCATCGCATGGGCGCGCTGAACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAAATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAgtg > 2:158982/1‑150 (MQ=255) tttCATCGCATGGGCGCGCTGAACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATATTCGGGTGAAGAGTGtt > 2:600489/1‑150 (MQ=255) tCGCATGGGCGCGCTGAACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGc > 1:288929/1‑150 (MQ=255) cgcTGAACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGcac < 2:288929/150‑1 (MQ=255) gcTGAACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGcaca > 1:614404/1‑150 (MQ=255) tGAACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGcaca < 1:26218/148‑1 (MQ=255) tGAACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGcaca > 2:26218/1‑148 (MQ=255) ggCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGcaca < 1:150249/143‑1 (MQ=255) ggCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGcaca > 2:150249/1‑143 (MQ=255) gttgttCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGAt < 1:616587/105‑1 (MQ=255) gttgttCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGAt > 2:616587/1‑105 (MQ=255) gttgttCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGgcg > 1:356792/1‑150 (MQ=255) ttgttCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGgcgc < 2:744914/150‑1 (MQ=255) ttCTTGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTgg < 2:399479/150‑1 (MQ=255) tGGATGCGTACGGTAGCACCGTGCTGCTCAATAACTATGATCGTTTCGGGCTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTggggg > 2:103109/1‑150 (MQ=255) tGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATCCACGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTgggtg < 1:406686/150‑3 (MQ=255) aTGCAGACTGTAGCACAGTGCTGATCAATATCTATGATCGTTTCGGACTTCAAAAAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGa < 2:247459/150‑1 (MQ=255) ggACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGATATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTAc < 2:667629/150‑1 (MQ=255) tCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGATATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGCACaggaggag > 1:428933/1‑113 (MQ=255) tCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGATATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGCACaggaggag < 2:428933/113‑1 (MQ=255) aTCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGa > 2:720100/1‑150 (MQ=255) ggcTTCACACAATCGAGATCGCGAGCGAGCGGGGTACGGCAGGTACCGATGTTCGGGTGACGTGTGTTGGCGCACTGGATGAACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAGCTTCTggc < 2:306827/148‑1 (MQ=255) tcgagatGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATc < 2:356792/150‑1 (MQ=255) gatGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCagag > 1:54000/1‑150 (MQ=255) | CGAAGGCTCGTCGCCAGTTGGACGCTACGCATGAGTTTCATCGCATGGGCGCGCTGAACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGACTTCAGCCAATCGAGATCGAGATGGAGCTGGGTACCGCAGGTACCGATGTTCGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAG > NC_002947/3431762‑3432050 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |