| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 2,747,717 | Δ1 bp | 19.0% | coding (561/618 nt) | PP_2401 ← | membrane protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 2,747,717 | 0 | G | . | 19.0% | 56.8 / 9.9 | 21 | coding (561/618 nt) | PP_2401 | membrane protein |
| Reads supporting (aligned to +/- strand): ref base G (9/8); new base . (2/2); total (11/10) | |||||||||||
| 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.17e-01 | |||||||||||
CGGCGCTCCAGCTGCCGAGGCGATCAACGGTGGTCAGGGTACGGGTCAAGGTCGTCATGTGCGGTGCTCTCTTCAGGCAAGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTTGACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTACCACCCCCCCATGGGGATGCGCTTGTAGAACGACAGACAGAACACGGCAGCCAGCACCCACATGCCGCTGCTGTGCAGCGC > NC_002947/2747588‑2747839 | cGGCGCTCCAGCTGCCGAGGCGATCAACGGTGGTCAGGGTACGGGTCAAGGTCGTCATGTGCGGTGCTCTCTTCAGGCAAGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTTGACGGCg > 2:465183/1‑136 (MQ=255)cGGCGCTCCAGCTGCCGAGGCGATCAACGGTGGTCAGGGTACGGGTCAAGGTCGTCATGTGCGGTGCTCTCTTCAGGCAAGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTTGACGGCg > 2:339481/1‑136 (MQ=255) ggCGCTCCAGCTGCCGAGGCGATCAACGGTGGTCAGGGTACGGGTCAAGGTCGTCATGTGCGGTGCTCTCTTCAGGCAAGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTTGACGGCGa < 1:166071/136‑1 (MQ=255) cgcTCCAGCTGCCGAGGCGATCAACGGTGGTCAGGGTACGGGTCAAGGTCGTCATGTGCGGTGCTCTCTTCAGGCAAGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTTGACGGCGATg < 1:339481/136‑1 (MQ=255) cgcTCCAGCTGCCGAGGCGATCAACGGTGGTCAGGGTACGGGTCAAGGTCGTCATGTGCGGTGCTCTCTTCAGGCAAGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTTGACGGCGATg < 1:465183/136‑1 (MQ=255) gtcgtcATGTGCGGTGCTCTCTTCAGGCAAGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTTGACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTACCACCCCCCCATggg > 1:293632/1‑136 (MQ=255) gtcgtcATGTGCGGTGCTCTCTTCAGGCAAGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTTGACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTACCACCCCCCCATggg > 2:310859/1‑136 (MQ=255) tgtgCGGTGCTCTCTTCAGGCAAGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTTGACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTACCACCCCCCCATGGGGATgcgc > 1:372082/1‑136 (MQ=255) tgtgCGGTGCTCTCTTCAGGCAAGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTTGACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTACCACCCCCCCATGGGGATgcgc > 2:17956/1‑136 (MQ=255) tGCTCTCTTCAGGCAAGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTTGACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTACCACCCCCCCATGGGGATGCGCTTGTAGa < 2:99120/136‑1 (MQ=255) aGGCAAGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTT‑ACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTac < 1:566932/97‑1 (MQ=255) aGGCAAGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTT‑ACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTac < 1:564583/97‑1 (MQ=255) aGGCAAGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTT‑ACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTac > 2:564583/1‑97 (MQ=255) aGGCAAGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTT‑ACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTac > 2:566932/1‑97 (MQ=255) aGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTTGACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTACCACCCCCCCATGGGGATGCGCTTGTAGAACGACAGACAGAaca > 2:113125/1‑136 (MQ=255) gcgGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTTGACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTACCACCCCCCCATGGGGATGCGCTTGTAGAACGACAGACAGAACACGGcc > 2:567389/1‑135 (MQ=255) cggcgTTGTTTACGGGCCACTTTGGCGGTGAGCTTGACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTACCACCCCCCCATGGGGATGCGCTTGTAGAACGACAGACAGAACACGGCAGCCAGCAccc < 1:310859/136‑1 (MQ=255) cgTTGTTTACGGGCCACTTTGGCGGTGAGCTTGACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTACCACCCCCCCATGGGGATGCGCTTGTAGAACGACAGACAGAACACGGCAGCCAGCACCcaca < 1:389909/136‑1 (MQ=255) tGGCGGTGAGCTTGACGGCGATGGCGCTGAACATCAGGGCGACCGCCAGCGGGTACCACCCCCCCATGGGGATGCGCTTGTAGAACGACAGACAGAACACGGCAGCCAGCACCCACATGCCGCTGCTGTGCAgcgc < 2:192750/136‑1 (MQ=255) ttGACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTACCACCCCCCCATGGGGATGCGCTTGTAGAACGACAGACAGAaca > 2:377275/1‑88 (MQ=255) ttGACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTACCACCCCCCCATGGGGATGCGCTTGTAGAACGACAGACAGAaca < 1:377275/88‑1 (MQ=255) | CGGCGCTCCAGCTGCCGAGGCGATCAACGGTGGTCAGGGTACGGGTCAAGGTCGTCATGTGCGGTGCTCTCTTCAGGCAAGGGCGCGGGTGCTGCGGCGTTGTTTACGGGCCACTTTGGCGGTGAGCTTGACGGCGATGGCGCTGAACATCAGGGCGAACGCCAGCGGGTACCACCCCCCCATGGGGATGCGCTTGTAGAACGACAGACAGAACACGGCAGCCAGCACCCACATGCCGCTGCTGTGCAGCGC > NC_002947/2747588‑2747839 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |