Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 3,011,010 | A→G | E531G (GAA→GGA) | ygfK → | putative selenate reductase subunit YgfK |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,011,010 | 0 | A | G | 96.7% | 187.1 / ‑1.0 | 61 | E531G (GAA→GGA) | ygfK | putative selenate reductase subunit YgfK |
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (35/24); total (36/25) | |||||||||||
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.84e-01 |
TTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGC > NZ_CP009273/3010923‑3011091 | tttGCGACCCCCAGTGCAAATAAAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGgaa < 2:446829/90‑1 (MQ=255) cACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTgg > 2:596417/1‑90 (MQ=255) cACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTgg > 1:1097561/1‑90 (MQ=255) ccAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGa < 2:373366/90‑1 (MQ=255) ccAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGa > 1:593896/1‑90 (MQ=255) ccAGTGCAAATACAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGa < 2:223173/90‑1 (MQ=255) cAGTGCAAATACACCTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGAt < 2:122779/90‑1 (MQ=255) agccccAATACAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATg < 1:202872/86‑1 (MQ=255) aaTACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAAtata > 2:1027052/1‑90 (MQ=255) tACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAg > 2:834416/1‑90 (MQ=255) aaCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCaa > 1:366558/1‑90 (MQ=255) tGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGc < 2:366558/90‑1 (MQ=255) gTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGc < 2:282819/89‑1 (MQ=255) gTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGc > 1:282819/1‑89 (MQ=255) gTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCt > 2:769650/1‑90 (MQ=255) aCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTgg > 1:1355953/1‑90 (MQ=255) cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac < 2:765179/90‑1 (MQ=255) cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac > 2:895436/1‑90 (MQ=255) cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGcac > 2:1354056/1‑90 (MQ=255) cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAATCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac > 1:377987/1‑90 (MQ=255) gCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcaca < 2:568329/90‑1 (MQ=255) ggATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGc < 1:310073/74‑1 (MQ=255) ggATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGc > 2:310073/1‑74 (MQ=255) tACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGcc < 2:377574/90‑1 (MQ=255) tACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGcc < 1:895436/90‑1 (MQ=255) aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGcc < 1:843610/89‑1 (MQ=255) aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGcc > 2:843610/1‑89 (MQ=255) aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCg > 1:1420228/1‑90 (MQ=255) aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCg > 2:1275111/1‑90 (MQ=255) cAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTc < 2:1210600/90‑1 (MQ=255) aGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac > 1:1440294/1‑75 (MQ=255) aGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac < 2:1440294/75‑1 (MQ=255) aGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCt > 1:405660/1‑90 (MQ=255) aGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCt > 2:1088307/1‑90 (MQ=255) gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt < 2:348983/53‑1 (MQ=255) gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt < 1:129448/53‑1 (MQ=255) gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt > 1:348983/1‑53 (MQ=255) gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt > 2:129448/1‑53 (MQ=255) gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg > 2:320554/1‑90 (MQ=255) gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg > 1:1186552/1‑90 (MQ=255) gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg > 2:713560/1‑90 (MQ=255) gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg > 1:268264/1‑90 (MQ=255) gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg > 2:1027028/1‑90 (MQ=255) gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg > 1:944630/1‑90 (MQ=255) gTGCGCTGAATATCCGCGAACTG‑AAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTgg > 2:632071/1‑90 (MQ=255) gcgcTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTggg > 1:370873/1‑89 (MQ=255) gcgcTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGt < 1:989530/90‑1 (MQ=255) gcTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTc < 2:1033409/90‑1 (MQ=255) cTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCa < 2:235797/59‑1 (MQ=255) cTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCa > 1:235797/1‑59 (MQ=255) atatCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGcc > 2:314056/1‑90 (MQ=255) atCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCAt < 1:769650/90‑1 (MQ=255) tCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATc > 2:278564/1‑90 (MQ=255) tCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATc > 2:1260958/1‑90 (MQ=255) gcgAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCgg > 1:130506/1‑90 (MQ=255) gAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGtt > 1:1155891/1‑90 (MQ=255) gAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGtt > 1:1118305/1‑90 (MQ=255) tGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCg < 2:293355/90‑1 (MQ=255) tGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCg < 1:320554/90‑1 (MQ=255) tGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCg < 1:1185317/90‑1 (MQ=255) aaaGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGAtt < 1:1275111/90‑1 (MQ=255) cgcgCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGc < 1:871305/90‑1 (MQ=255) | TTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGC > NZ_CP009273/3010923‑3011091 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |