Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 3,002,903 | T→C | S261P (TCA→CCA) | ygeY → | YgeY family selenium metabolism‑linked hydrolase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,002,903 | 0 | T | C | 92.7% | 190.0 / 5.9 | 68 | S261P (TCA→CCA) | ygeY | YgeY family selenium metabolism‑linked hydrolase |
Reads supporting (aligned to +/- strand): ref base T (2/3); new base C (30/33); total (32/36) | |||||||||||
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.43e-01 |
GAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAAC > NZ_CP009273/3002816‑3002989 | gAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCa > 2:212234/1‑90 (MQ=255) gAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCa > 1:264588/1‑90 (MQ=255) gAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCa > 1:1265666/1‑90 (MQ=255) aaTTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCaa > 2:97743/1‑90 (MQ=255) aaTTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCaa > 1:1268865/1‑90 (MQ=255) cGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACcg < 2:697830/90‑1 (MQ=255) cGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACcg < 1:541778/90‑1 (MQ=255) cGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACcg < 2:295379/90‑1 (MQ=255) ggCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACcgt < 2:59317/90‑1 (MQ=255) aaGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc < 2:390943/90‑1 (MQ=255) aaGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc < 2:1159967/90‑1 (MQ=255) aaGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc < 1:126165/90‑1 (MQ=255) aaGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc < 1:440912/90‑1 (MQ=255) aCCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGAcc > 1:1435309/1‑90 (MQ=255) ccTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTg < 2:1214834/90‑1 (MQ=255) ccTCACCGTTTCTGAAATCTTCTTCACACCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTg > 2:1001104/1‑90 (MQ=255) cACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTgggg > 2:881365/1‑90 (MQ=255) cGTTTCTGAAATCTTCTTCACTTCCCCAAGCCGTTGCGCTGCAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGa < 1:1350391/90‑1 (MQ=255) cGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGa < 1:212234/90‑1 (MQ=255) cGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGa < 1:1318857/90‑1 (MQ=255) tttCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGaaa < 2:1331399/90‑1 (MQ=255) tttCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGaaa < 1:595090/90‑1 (MQ=255) ttCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAAc > 1:559011/1‑90 (MQ=255) ttCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAAc > 1:480611/1‑90 (MQ=255) ttCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAAc > 2:1275997/1‑90 (MQ=255) ttCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAAc > 2:184341/1‑90 (MQ=255) ttCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGACGAAAc > 1:1526765/1‑90 (MQ=255) tcACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGgcgcgc > 1:989735/1‑90 (MQ=255) tcACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGgcgcgc > 1:975098/1‑90 (MQ=255) cacaTCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCt < 2:159438/90‑1 (MQ=255) cacaTCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCt < 2:264588/90‑1 (MQ=255) tCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGAc > 1:847192/1‑90 (MQ=255) ccAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGaa < 1:913059/90‑1 (MQ=255) ccGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCcg > 1:853051/1‑90 (MQ=255) gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCcgcg > 2:757533/1‑90 (MQ=255) gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCcgcg > 1:1483002/1‑90 (MQ=255) gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCcgcg > 2:1203887/1‑90 (MQ=255) ttGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCgcgc > 1:1452184/1‑90 (MQ=255) ttGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCgcgc > 1:629347/1‑90 (MQ=255) tGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGcc < 1:869123/90‑1 (MQ=255) tGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGcc < 1:206737/90‑1 (MQ=255) tGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGcc < 1:574991/90‑1 (MQ=255) gcgcTGTAGCAGACAGCTGCGCAGTTTCAAGCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGccc < 1:1046468/90‑1 (MQ=255) gcgcTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGccc < 2:393938/90‑1 (MQ=255) gcTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctg < 2:853051/90‑1 (MQ=255) gcTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctg < 2:1570242/90‑1 (MQ=255) tGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctgcc < 1:893438/90‑1 (MQ=255) tGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctgcc < 1:1004685/90‑1 (MQ=255) gTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctgcct < 2:680590/90‑1 (MQ=255) aGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc < 1:1203887/90‑1 (MQ=255) aGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc < 2:173774/90‑1 (MQ=255) aGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc < 2:1270660/90‑1 (MQ=255) aGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc < 2:477724/90‑1 (MQ=255) cAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAg < 2:1175837/90‑1 (MQ=255) gACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTa > 1:713240/1‑90 (MQ=255) cAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACa > 2:1338884/1‑90 (MQ=255) cAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACa > 1:1321741/1‑90 (MQ=255) aGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAg < 2:576068/90‑1 (MQ=255) gCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa > 1:789871/1‑90 (MQ=255) gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa > 1:1318398/1‑90 (MQ=255) gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa < 2:1265666/90‑1 (MQ=255) gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa > 1:374494/1‑90 (MQ=255) gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa > 1:569860/1‑90 (MQ=255) gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa > 2:1213476/1‑90 (MQ=255) gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa > 1:726088/1‑90 (MQ=255) gCTGCGCAGTTCCAAACGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa < 1:757163/90‑1 (MQ=255) gcgcAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAg > 2:584261/1‑90 (MQ=255) gTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAAc < 2:1067407/90‑1 (MQ=255) | GAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAAC > NZ_CP009273/3002816‑3002989 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |