Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,647,969 | G→A | 74.5% | Y265Y (TAC→TAT) | oprB‑II ← | carbohydrate‑selective porin |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,647,969 | 0 | G | A | 74.5% | 87.0 / 30.4 | 55 | Y265Y (TAC→TAT) | oprB‑II | carbohydrate‑selective porin |
Reads supporting (aligned to +/- strand): ref base G (7/7); new base A (21/20); total (28/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.37e-01 |
AGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCGAAG > NC_002947/1647835‑1648113 | aGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCAc > 1:833846/1‑148 (MQ=255) aGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCAc < 2:833846/148‑1 (MQ=255) ctgctgctgGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGc > 2:1371654/1‑131 (MQ=255) ctgctgctgGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGc < 1:1371654/131‑1 (MQ=255) tgctgGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAg > 1:368857/1‑150 (MQ=255) ctgGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGcc > 1:1071550/1‑150 (MQ=255) tgGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCt > 1:1428262/1‑150 (MQ=255) tgGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCt > 1:540219/1‑150 (MQ=255) gggCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCtt > 1:1323649/1‑150 (MQ=255) gcgcCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTgg < 2:1323649/150‑1 (MQ=255) gcCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGAt < 1:554265/150‑1 (MQ=255) cAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATAc > 1:52326/1‑150 (MQ=255) aTCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGt > 2:163425/1‑150 (MQ=255) gCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCa > 1:390326/1‑150 (MQ=255) ttgctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGt < 2:934905/150‑1 (MQ=255) ctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCg > 2:1278926/1‑150 (MQ=255) tCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGa > 2:552877/1‑150 (MQ=255) tCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGa > 2:1531623/1‑150 (MQ=255) tgctgcGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTc < 2:1014662/142‑1 (MQ=255) tgctgcGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTc > 1:1014662/1‑142 (MQ=255) gcGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTAc > 1:1164535/1‑150 (MQ=255) gcggcgCCGCTGAGGGCGGCCTGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATATCAGTTCGAACGGCATTAc > 1:1325326/1‑143 (MQ=255) gcggcgCCGCTGAGGGCGGCCTGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATATCAGTTCGAACGGCATTAc < 2:1325326/143‑1 (MQ=255) cggcgCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCt < 1:552877/150‑1 (MQ=255) gcgcCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTgg < 1:954852/150‑1 (MQ=255) ggcggcCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGTTATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACGAGTTCGGACGGCATTACCGCACCCTGGGtgcagctgcc > 1:932874/1‑150 (MQ=255) ggcggcCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTggg < 1:1199133/140‑1 (MQ=255) ggcggcCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTggg > 2:1199133/1‑140 (MQ=255) ggccggcTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTgccgctgccgct < 2:1458677/150‑1 (MQ=255) ccGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGaaa < 1:1247255/150‑1 (MQ=255) tGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGc < 2:1071550/150‑1 (MQ=255) gTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGcc < 2:634806/150‑1 (MQ=255) tCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATAc < 2:415505/89‑1 (MQ=255) tCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATAc > 1:415505/1‑89 (MQ=255) ttGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGAt > 2:791129/1‑150 (MQ=255) gaACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGa < 2:52326/150‑1 (MQ=255) aaCATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAAc > 2:981723/1‑91 (MQ=255) aaCATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAAc < 1:981723/91‑1 (MQ=255) aaCATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAg < 1:791129/150‑1 (MQ=255) aCATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAgg < 2:540219/150‑1 (MQ=255) aCATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAgg < 2:1599258/150‑1 (MQ=255) cTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTgg < 2:1164535/150‑1 (MQ=255) cTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTgg > 2:28778/1‑150 (MQ=255) ccTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGa < 2:987245/120‑1 (MQ=255) ccTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGa > 1:987245/1‑120 (MQ=255) ccTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGCGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGt > 2:1360679/1‑150 (MQ=255) ccTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGt < 1:1346105/150‑1 (MQ=255) cTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGGGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCGGCTGGGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGtt > 1:310962/1‑150 (MQ=255) ttGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATCCCCGTCCGCACGGCATTACCCCACCCTGGGTGCGGCTGCCGGTGAGGTTGAAACCAATGCCGGATTCGGGGCTGGAAtgggac > 1:71189/1‑147 (MQ=255) aTTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGc < 1:1207559/122‑1 (MQ=255) aTTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGc > 2:1207559/1‑122 (MQ=255) aCTTTAATAGTATCCGGCACGATAGTCCCCTTGCAAGCCTTGGATACGTGGTGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCACTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCGaa < 1:746110/150‑1 (MQ=38) aCTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCGaa < 2:390326/150‑1 (MQ=255) aCTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCGaa < 2:1428262/150‑1 (MQ=255) gtGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCGAAg > 2:1125418/2‑150 (MQ=255) gtaatagtaGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGt > 2:668605/1‑128 (MQ=255) | AGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCGAAG > NC_002947/1647835‑1648113 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |