Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,050,011 | T→G | 8.0% | L15R (CTG→CGG) ‡ | oprB‑III → | carbohydrate‑selective porin |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,050,011 | 0 | T | G | 8.0% | 139.5 / 5.1 | 50 | L15R (CTG→CGG) ‡ | oprB‑III | carbohydrate‑selective porin |
Reads supporting (aligned to +/- strand): ref base T (20/25); new base G (2/2); total (23/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01 |
AAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACGGT > NC_002947/4049867‑4050159 | aaaCCTCGTCGCGTGCTTCGGGGGCAAGCCGGCTCCTACAAGGCGGGTGTAGGATCGGGCTTGCGCGCTAATGCC‑CCTGAAAGACAACCCGAGATTCTGTAATGCGATCCGCAATTCGCATCACCCCCACGCTGCTGctggccctggc < 4:118075/148‑1 (MQ=255) aaCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCa < 5:132419/149‑1 (MQ=255) cgcgTCGCCTGCTTCGCGGGCAACCCCGCTCCAACAAGGCCCCTGTCGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGca < 5:195358/149‑1 (MQ=255) gcgTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGcac < 3:295019/149‑1 (MQ=255) gTCGCCGGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAACGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAgcagca < 4:189669/149‑1 (MQ=38) cTGCTTCGCGGGCAACCCCCCACCAACAAAGCCCCGGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCACAAAAACAAACCCAGAATATGTAAATCCAAACGCAATTATCATAACCCCCCCCCAGCAGCTGGCCCAGGCcaccacaccagcc > 8:151295/1‑142 (MQ=255) cTGCTACGCGGCCACCCCCGCTCAAACACGCCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGcc < 5:276875/149‑1 (MQ=255) gCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGCC‑CCTCAAAGCCAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCGGGCCCTGGCCAGCACCCCTGCCCt > 5:224505/1‑149 (MQ=255) ttCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAgcagcagg < 2:97554/142‑3 (MQ=255) cgcgGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc < 8:279865/148‑1 (MQ=255) aaGCCCGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgac > 2:148047/1‑148 (MQ=255) ccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct > 1:290964/1‑149 (MQ=255) gCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa > 3:307629/1‑149 (MQ=255) gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa < 8:290521/149‑1 (MQ=255) gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa > 1:221646/1‑149 (MQ=255) cTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGACCTGGCCAGCACCACTGCCCTTGGCAGCGACGATCTCATGAc > 2:182283/1‑149 (MQ=255) cAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCAATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCg > 8:187349/1‑149 (MQ=255) aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgc < 1:166950/149‑1 (MQ=255) ccTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAcc < 2:290964/149‑1 (MQ=255) ccTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATACGCAATTCGCATCACCCCCACCCTCCTGCTGGCCCTGGCCAGCACCACTGCCCTTGACAGCCACGAACTCATGACCCgcaccccc > 2:318953/1‑144 (MQ=255) tGTAGGAGCGGGCTTGCCCGCGAATGCCGCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAccc > 1:19837/1‑149 (MQ=255) gAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCg > 5:135059/1‑149 (MQ=255) aGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCgg < 7:97547/149‑1 (MQ=255) aGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCgg < 2:221646/149‑1 (MQ=255) ggCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑c < 8:75452/149‑1 (MQ=255) cccGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTTCCAGCGACGATCTCATGTCCCGCAGCACACTGACCGGTGA‑CTgggg < 1:182283/148‑1 (MQ=255) gAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCt < 7:139179/149‑1 (MQ=255) cc‑ccTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGCTCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTgcgc > 3:203779/1‑148 (MQ=255) c‑ccTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCa < 4:307629/149‑1 (MQ=255) cAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCa > 3:250322/1‑148 (MQ=255) cgACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTc < 7:187349/148‑1 (MQ=255) aGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTc > 2:104678/1‑149 (MQ=255) cAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCgaa < 7:151295/149‑1 (MQ=255) cAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCgaa < 5:80665/149‑1 (MQ=255) cccGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCgaag > 2:75976/1‑147 (MQ=255) ccGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCgaag > 6:275738/1‑146 (MQ=255) gaTTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACgg < 6:135059/149‑1 (MQ=255) tAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTc > 2:216066/1‑146 (MQ=255) tAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCa < 4:203779/147‑1 (MQ=255) tAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCATGGCCAGCACCACTGCCCTTGCAAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTc > 1:164202/1‑146 (MQ=255) tCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACCTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCg < 1:75976/149‑1 (MQ=255) tCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCgg > 4:288503/1‑149 (MQ=255) cGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGc < 2:164202/148‑1 (MQ=255) cAATTCGCAGCCCCCCCCCCCTTCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGAc < 6:224505/149‑1 (MQ=255) tCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGt < 1:148047/149‑1 (MQ=255) cccACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAAc > 8:80479/1‑149 (MQ=255) cccTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGcc < 7:80479/149‑1 (MQ=255) agtgctgctgGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAg > 4:180755/3‑118 (MQ=255) agtgctgctgGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAg < 3:180755/116‑1 (MQ=255) ggtgctgGCCAGGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTAc > 3:2206/3‑148 (MQ=255) gtgctgGCCAGGGCCAGCAGCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CT‑GGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAAc > 2:33114/2‑149 (MQ=255) ctggccctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGc < 3:288503/149‑1 (MQ=255) ctggccctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGc > 2:95556/1‑149 (MQ=255) gccctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGccc < 5:275738/148‑1 (MQ=255) cctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACg < 1:104678/149‑1 (MQ=255) tggccAGCACCACTGCCCATGCCAGCGACGATCTCATGACCCGCAGCACCCAGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGAGAAACCGCCAACAACGCCCACGGt > 4:128930/1‑149 (MQ=255) | AAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCC‑CCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGA‑CTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACGGT > NC_002947/4049867‑4050159 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |