| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 3,304,580 | C→A | 9.7% | P602H (CCT→CAT) | pvdQ → | acyl‑homoserine lactone acylase PvdQ |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 3,304,580 | 0 | C | A | 9.7% | 137.9 / 2.6 | 41 | P602H (CCT→CAT) | pvdQ | acyl‑homoserine lactone acylase PvdQ |
| Reads supporting (aligned to +/- strand): ref base C (21/16); new base A (2/2); total (23/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01 | |||||||||||
GCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCAGGCTGGTTCATCAGGCAGCACTGGCGTCGCTCAAGCAG > NC_002947/3304463‑3304701 | gcgcAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGc > 2:509474/1‑136 (MQ=255)gcgcAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGc > 2:553892/1‑136 (MQ=255) cAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCgg < 2:30969/136‑1 (MQ=255) cAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCgg < 1:73612/136‑1 (MQ=255) gtgtGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGc > 1:54860/1‑112 (MQ=255) gtgtGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGc < 2:54860/112‑1 (MQ=255) gcagcaGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCgca > 1:87273/1‑136 (MQ=255) cagcagCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCATGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCgcag > 1:5243/1‑136 (MQ=255) cagcagCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCATGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCgcag > 1:670996/1‑136 (MQ=255) tGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACAcccc < 1:449514/136‑1 (MQ=255) tGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACAcccc < 2:372593/136‑1 (MQ=255) tGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACAcccc < 2:131619/136‑1 (MQ=255) ggCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCcg > 1:334785/1‑136 (MQ=255) gccTGGAACGGCAAGGCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCcgcgcg > 2:609205/1‑136 (MQ=255) ggCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAg > 2:590243/1‑136 (MQ=255) ggCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAg > 2:125246/1‑136 (MQ=255) gCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGc < 1:525000/136‑1 (MQ=255) aCCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAggc < 2:18912/136‑1 (MQ=255) gACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGcggcgg > 2:542043/1‑136 (MQ=255) ggTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCATGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCa < 2:670996/136‑1 (MQ=255) ggTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCATGAAAGCTGGCGAGTGGCATTCAACCCGCCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCa < 2:5243/136‑1 (MQ=255) tATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCAgg > 1:170212/1‑136 (MQ=255) gCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCAGGCTGGt > 1:589669/1‑136 (MQ=255) gCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCAGGCTGGt > 2:7429/1‑136 (MQ=255) gCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCAGGCTGGt < 1:650419/136‑1 (MQ=255) gTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCAtt < 1:604252/62‑1 (MQ=255) gTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCAtt > 2:604252/1‑62 (MQ=255) gTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCAGGCTGGTtcatc < 1:7429/136‑1 (MQ=255) gTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCAGGCTGGTtcatc < 2:226934/136‑1 (MQ=255) cACTTCCAGAACCTGTTCAATGCGCTGGCCGAGTACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCAGGCTGGTTCATCAgg > 2:673798/1‑136 (MQ=255) cACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTc < 2:163860/102‑1 (MQ=255) cACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTc > 1:163860/1‑102 (MQ=255) ccAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCAGGCTGGTTCATCAGgcagca > 2:21154/1‑136 (MQ=255) ccTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCATGCTGGTTCATCAGGCAGCACTGGCg > 1:523101/1‑136 (MQ=255) cTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCAGGCTGGTTCATCAGGCAGCACTGGCGt < 2:36739/136‑1 (MQ=255) tCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCAGGCTGGTtcatc < 2:459459/117‑1 (MQ=255) tCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCAGGCTGGTtcatc > 1:459459/1‑117 (MQ=255) tCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCAGGCTGGTTCATCAGGCAGCACTGGCGTCGCt < 2:578695/136‑1 (MQ=255) tCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCAGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCAGGCTGGTTCATCAGGCAGCACTGGCGTCGCt > 1:219030/1‑136 (MQ=255) aaTGCGCTGGCCGAGCACCCTGAAAGCTGGAGAGTGGCATTAAACCCGGCCGATCAGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCAGGCTGGTTCATCAGGCAGCACTCGCGTCGCTCa > 2:512202/1‑136 (MQ=255) gcTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGTCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGAAAGCAGGCTGGTTCATCAGGCAGCACTCGCGTCGCTCAAGCAg > 2:674509/1‑136 (MQ=255) | GCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGTATGGGCCTGGTGCACTTCCAGAACCTGTTCAATGCGCTGGCCGAGCACCCTGAAAGCTGGCGAGTGGCATTCAACCCGGCCGATCCGCAGCACACCCCGCGCGGCCTGGCTGTCGAGCAGGCGGCGGTAAGCAGGCTGGTTCATCAGGCAGCACTGGCGTCGCTCAAGCAG > NC_002947/3304463‑3304701 |
| 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: ‑ |