| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 3,217,770 | G→A | 12.5% | R187H (CGC→CAC) | oprJ → | outer membrane protein OprJ |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 3,217,770 | 0 | G | A | 12.5% | 81.2 / 3.1 | 32 | R187H (CGC→CAC) | oprJ | outer membrane protein OprJ |
| Reads supporting (aligned to +/- strand): ref base G (15/13); new base A (2/2); total (17/15) | |||||||||||
| 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.46e-01 | |||||||||||
AAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAA > NC_002947/3217646‑3217912 | aaaaGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCt < 6:97200/149‑1 (MQ=255) aaaGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTg > 3:26905/1‑149 (MQ=255) gCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGt < 1:214045/149‑1 (MQ=255) cAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCTTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGtt > 1:129380/1‑149 (MQ=255) ccTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTc < 3:173599/148‑1 (MQ=255) ccTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCg < 8:223284/149‑1 (MQ=255) cTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCgg > 6:52/1‑149 (MQ=255) gAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAgcggcg > 1:14654/1‑149 (MQ=255) ggAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGgcgcg > 4:228892/1‑149 (MQ=255) aTCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCt > 6:277887/1‑149 (MQ=255) ccTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCACCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGac > 8:11321/1‑148 (MQ=255) aTCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGaccac < 5:52/147‑1 (MQ=255) ccGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAagag > 7:146599/1‑149 (MQ=255) gAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCGTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGgc < 4:26905/149‑1 (MQ=255) aGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTgg < 3:242075/148‑1 (MQ=255) ccAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGgcgc < 4:282022/103‑1 (MQ=255) ccAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGgcgc > 3:282022/1‑103 (MQ=255) caaGGCTTACCTGGCCTTGGACGCCGCCGCGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGc < 2:129380/147‑1 (MQ=255) gCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTgg > 2:34529/1‑124 (MQ=255) gCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTgg < 1:34529/124‑1 (MQ=255) gCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGt > 4:26205/1‑149 (MQ=255) gCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTCGGGCTGGt > 4:48105/1‑149 (MQ=255) cTACCTGGCCTTGGACGGCGCCGAGCGGCACCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTgg < 2:275134/122‑1 (MQ=255) cTACCTGGCCTTGGACGGCGCCGAGCGGCACCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTgg > 1:275134/1‑122 (MQ=255) ggccTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGa < 1:251173/147‑1 (MQ=255) ccTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAg > 6:4195/1‑148 (MQ=255) tGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGGCCACCAAGAGGCGCTGGGGCTGGTCGAGCAg > 6:114576/1‑149 (MQ=255) cTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTccc > 7:115380/1‑149 (MQ=255) ttGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCg > 8:63336/1‑149 (MQ=255) aGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGaa < 8:115380/149‑1 (MQ=255) aGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGaa < 2:14654/149‑1 (MQ=255) aGCGGCACCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGaa < 7:11321/149‑1 (MQ=255) | AAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAA > NC_002947/3217646‑3217912 |
| 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: ‑ |