| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 1,925,271 | C→A | 19.0% | V59V (GTG→GTT) | prc ← | tail‑specific protease Prc |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 1,925,271 | 0 | C | A | 19.0% | 49.7 / 2.5 | 21 | V59V (GTG→GTT) | prc | tail‑specific protease Prc |
| Reads supporting (aligned to +/- strand): ref base C (10/7); new base A (2/2); total (12/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01 | |||||||||||
TCGGCAGCGGTGAAGTAGCTGCGCGCCGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAA > NC_002947/1925152‑1925406 | tCGGCAGCGGTGAAGTAGCTGCGCGCCGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTgg < 1:64374/136‑1 (MQ=255) cGGCAGCGGTGAAGTAGCTGCGCGCCGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGc > 1:95365/1‑136 (MQ=255) cGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTgg < 2:255066/110‑1 (MQ=255) cGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTgg > 1:255066/1‑110 (MQ=255) gggTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCa < 1:356677/136‑1 (MQ=255) ccAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCAACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCt > 1:393748/1‑136 (MQ=255) ccAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCAACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCt > 1:419748/1‑136 (MQ=255) tGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACttgtt > 2:447718/1‑136 (MQ=255) tGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTgcc < 2:95365/136‑1 (MQ=255) gggCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCgg > 2:459855/1‑135 (MQ=255) gggCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGt > 1:220603/1‑136 (MQ=255) gggCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGt > 2:409171/1‑136 (MQ=255) cGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCt > 1:539493/1‑136 (MQ=255) tCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGc > 2:344360/1‑136 (MQ=255) ggcggcTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGc < 1:447718/136‑1 (MQ=255) ggcggcGTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGc < 1:406723/136‑1 (MQ=255) gcTTGCTGTAGTGGTGGCGCTTGAGCAATTCAACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGat < 2:393748/136‑1 (MQ=255) gcTTGCTGTAGTGGTGGCGCTTGAGCAATTCAACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGat < 2:419748/136‑1 (MQ=255) tGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCTGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGatcat > 1:133201/1‑136 (MQ=255) tGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGatcat > 1:36542/1‑136 (MQ=255) aTTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTggaag < 1:138426/136‑1 (MQ=255) caccacGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGaaa > 2:113638/1‑136 (MQ=255) caccacGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGaaa > 2:135481/1‑136 (MQ=255) caccacGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGaaa > 2:456219/1‑136 (MQ=255) caccacGTTGAGGCTGGCGACGATATATTAGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGaaa > 2:486757/1‑136 (MQ=255) | TCGGCAGCGGTGAAGTAGCTGCGCGCCGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAA > NC_002947/1925152‑1925406 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |