Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,013,339 | G→A | 25.0% | A38T (GCG→ACG) | PP_5678 → | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,013,339 | 0 | G | A | 25.0% | 21.2 / 5.2 | 16 | A38T (GCG→ACG) | PP_5678 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (6/6); new base A (2/2); total (8/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01 |
CGAGCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTACTCT > NC_002947/5013190‑5013483 | cGAGCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTCGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAgcg > 2:474024/1‑150 (MQ=255) ttGCAACAAGCCTCAGACTGGAGTTAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCACGGCGt < 1:32191/149‑1 (MQ=255) aCTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCgca < 1:535648/150‑1 (MQ=255) tCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGTGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAAc > 2:121190/1‑150 (MQ=255) aaGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTcc > 2:424541/1‑150 (MQ=255) gTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCtt > 1:335456/1‑150 (MQ=255) tgACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGAt < 1:121190/149‑1 (MQ=255) tGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGt < 1:559899/150‑1 (MQ=255) gTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTgg < 1:272070/146‑1 (MQ=255) gTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTgg > 2:272070/1‑146 (MQ=255) tcaGCACGGCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGcttct < 1:255138/88‑1 (MQ=255) tcaGCACGGCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGcttct > 2:255138/2‑89 (MQ=255) aGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCATTGCCGGTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCAAATTGATGTTCTCTGAGGAAGCTg > 1:195256/1‑150 (MQ=255) aCGCCGCGCTCAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACt > 1:135233/1‑150 (MQ=255) tGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTActct < 1:424541/150‑1 (MQ=255) gAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTc > 2:242525/1‑141 (MQ=255) gAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTc < 1:242525/141‑1 (MQ=255) | CGAGCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTACTCT > NC_002947/5013190‑5013483 |
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: ‑ |