| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 2,251,549 | G→T | 18.2% | V1014V (GTC→GTA) | PP_1983 ← | PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 2,251,549 | 0 | G | T | 18.2% | 59.6 / 3.1 | 22 | V1014V (GTC→GTA) | PP_1983 | PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase |
| Reads supporting (aligned to +/- strand): ref base G (9/9); new base T (2/2); total (11/11) | |||||||||||
| 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.83e-01 | |||||||||||
ACCGCCAGGTCCAGGCTGCGCGCCATGGCGACAATGGCGCGGACAATTTCTGCGTCGTTGGGGTCGACCGGAGCATCGCGTACGAAGCTCTGGTCGATCTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTCATCCATGGCAAAACTCACCCCGTAGCGCTTCAGCTCGCACATCTTGCTGATGGTGTCTTCC > NC_002947/2251439‑2251660 | aCCGCCAGGTCCAGGCTGCGCGCCATGGCGACAATGGCGCGGACAATTTCTGCGTCGTTGGGGTCGACCGGAGCATCGCGTACGAAGCTCTGGTCGATCTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCa < 2:45870/136‑1 (MQ=255) gcCATGGCGACAATGGCGCGGACAATTTCTGCGTCGTTGGGGTCGACCGGAGCATCGCGTACGAAGCTCTGGTCGATCTTCAGCGCATCTACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGa > 2:99773/1‑136 (MQ=255) gcCATGGCGACAATGGCGCGGACAATTTCTGCGTCGTTGGGGTCGACCGGAGCATCGCGTACGAAGCTCTGGTCGATCTTCAGCGCATCTACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGa > 2:499969/1‑136 (MQ=255) gcCATGGCGACAATGGCGCGGACAATTTCTGCGTCGTTGGGGTCGACCGGAGCATCGCGTACGAAGCTCTGGTCGATCTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGa > 2:5590/1‑136 (MQ=255) gcCATGGCGACAATGGCGCGGACAATTTCTGCGTCGTTGGGGTCGACCGGAGCATCGCGTACGAAGCTCTGGTCGATCTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGa > 1:322113/1‑136 (MQ=255) gcCATGGCGACAATGGCGCGGACAATTTCTGCGTCGTTGGGGTCGACCGGAGCATCGCGAACGAAGCTCTGGTCGATCTTCAGCGCATCGACCGGCAGCCCCTTCAGCTAAGTCAGCGAGGAATAGCCGATGCCGa > 1:309973/1‑136 (MQ=255) tGGCGACAATGGCGCGGACAATTTCTGCGTCGTTGGGGTCGACCGGAGCATCGCGTACGAAGCTCTGGTCGATCTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTc < 1:184639/136‑1 (MQ=255) tGGCGACAATGGCGCGGACAATTTCTGCGTCGTTGGGGTCGACCGGAGCATCGCGTACGAAGCTCTGGTCGATCTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTc < 1:308938/136‑1 (MQ=255) aaTGGCGCGGACAATTTCTGCGTCGTTGGGGTCGACCGGAGCATCGCGTACGAAGCTCTGGTCGATCTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTCATCCATg < 2:352405/136‑1 (MQ=255) gcgcGGACAATTTCTGCGTCGTTGGGGTCGACCGGAGCATCGCGTACGAAGCTCTGGTCGATCTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTCATCCATGGCaa > 2:7736/1‑136 (MQ=255) aaTTTCTGCGTCGTTGGGGTCGACCGGAGCATCGCGTACGAAGCTCTGGTCGATCTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTCATCCATGGCAAAACTCAcc > 2:251612/1‑136 (MQ=255) aTTTCTGCGTCGTTGGGGTCGACCGGAGCATCGCGTACGAAGCTCTGGTCGATCTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTCATCCATGGCAAAACTCAccc > 1:100310/1‑136 (MQ=255) tttCTGCGTCGTTGGGGTCGACCGGAGCATCGCGTACGAAGCTCTGGTCGATCTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTCATCCATGGCAAAACTCAcccc < 1:205120/136‑1 (MQ=255) cGGAGCATCGCGTACGAAGCTCTGGTCGATCTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTCATCCATGGCAAAACTCACCCCGTAGCGCTTCAGCTCGCACATc < 2:100310/136‑1 (MQ=255) aaGCTCTGGTCGATCTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTCATCCATGGCAAAACTCACCCCGTAGCGCTTCAGCTCGCACATCTTGCTGATGGTGTCtt > 2:177969/1‑136 (MQ=255) gCTCTGGTCGATCTTCAGCGCATCTACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTCATCCATGGCAAAACTCACCCCGTAGCGCTTCAGCTCGCACATCTTGCTGATGGTGTCTTcc < 1:499969/136‑1 (MQ=255) gCTCTGGTCGATCTTCAGCGCATCTACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTCATCCATGGCAAAACTCACCCCGTAGCGCTTCAGCTCGCACATCTTGCTGATGGTGTCTTcc < 1:99773/136‑1 (MQ=255) gCTCTGGTCGATCTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTCATCCATGGCAAAACTCACCCTGTAGCGCTTCAGCTCGCACATCTTGCTGATGGTGTCTTcc < 2:92004/136‑1 (MQ=255) cTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTCATCCATGGCaa > 2:409062/1‑74 (MQ=255) cTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTCATCCATGGCaa < 1:409062/74‑1 (MQ=255) tCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTCATCCATGGCAAAACTCACCCCGTAGCGCTTCAGCTCGCACATCTTGCt > 1:471858/1‑109 (MQ=255) tCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTCATCCATGGCAAAACTCACCCCGTAGCGCTTCAGCTCGCACATCTTGCt < 2:471858/109‑1 (MQ=255) | ACCGCCAGGTCCAGGCTGCGCGCCATGGCGACAATGGCGCGGACAATTTCTGCGTCGTTGGGGTCGACCGGAGCATCGCGTACGAAGCTCTGGTCGATCTTCAGCGCATCGACCGGCAGGCGCTTCAGGTAAGTCAGCGAGGAATAGCCGGTGCCGAAGTCATCCATGGCAAAACTCACCCCGTAGCGCTTCAGCTCGCACATCTTGCTGATGGTGTCTTCC > NC_002947/2251439‑2251660 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |