Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 1,657,901 | A→T | 24.5% | intergenic (+31/‑168) | ynfD → / → ynfE | DUF1161 family periplasmic protein/putative selenate reductase, periplasmic |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,657,901 | 0 | A | T | 24.5% | 143.1 / 53.5 | 86 | intergenic (+31/‑168) | ynfD/ynfE | DUF1161 family periplasmic protein/putative selenate reductase, periplasmic |
Reads supporting (aligned to +/- strand): ref base A (38/27); new base T (10/11); total (48/38) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.53e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-02 |
CCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCC > NC_000913/1657803‑1657998 | ccaccaGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACaaa < 1:1981839/100‑1 (MQ=255) ccaGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGCAACACACTATCGACTGAACGCCGGATATGACAAATcc > 1:2194944/1‑100 (MQ=255) caGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCg < 1:1016038/100‑1 (MQ=255) gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGTTGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc < 1:1127528/100‑1 (MQ=255) gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGc < 1:105191/100‑1 (MQ=255) gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGc < 1:922805/100‑1 (MQ=255) gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGc < 2:2561282/100‑1 (MQ=255) gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc > 2:1094316/1‑100 (MQ=255) gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc < 1:2417611/100‑1 (MQ=255) cGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCg < 2:2074210/100‑1 (MQ=255) cGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCg < 2:2105591/100‑1 (MQ=255) ggTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgt < 1:1449607/100‑2 (MQ=255) ggTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgt < 1:1751393/100‑2 (MQ=255) ggTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgt < 1:903687/100‑2 (MQ=255) cGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTACCACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttcagt < 2:1073641/100‑7 (MQ=37) cGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttcagt < 2:1761173/100‑7 (MQ=38) gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttcagtc < 2:1045060/100‑8 (MQ=37) gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGaa > 2:749899/1‑100 (MQ=255) tctcTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAAc < 1:895122/100‑1 (MQ=255) tcgctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATCATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAAc < 1:1986014/97‑1 (MQ=255) ctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATa > 1:921602/1‑100 (MQ=255) ctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATa < 2:2495847/100‑1 (MQ=255) gCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATaca < 1:1932608/100‑1 (MQ=255) gCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATaca < 2:324441/100‑1 (MQ=255) gCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATaca > 2:769809/1‑100 (MQ=255) cTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaa > 2:136885/1‑100 (MQ=255) cTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaa > 2:1545490/1‑100 (MQ=255) ccGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACAt > 2:183234/1‑100 (MQ=255) gTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCa < 2:1094029/100‑1 (MQ=255) gTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCa > 2:367666/1‑100 (MQ=255) gTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCa > 2:770858/1‑100 (MQ=255) tCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATCATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAc < 1:626069/100‑1 (MQ=255) gCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTCAACATACAACATAATACCACCtta > 1:207042/1‑100 (MQ=255) cAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt > 1:347402/1‑100 (MQ=255) cAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt > 2:600472/1‑100 (MQ=255) cAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCAACCttagt > 1:1566596/1‑98 (MQ=255) ggATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTAc < 2:1850711/100‑1 (MQ=255) gATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACt < 1:1545481/100‑1 (MQ=255) gcgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCAt < 1:656282/99‑1 (MQ=255) gcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATAcc < 1:600469/100‑1 (MQ=255) gcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCt < 1:2347601/100‑1 (MQ=255) ccTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCtt > 2:185564/1‑100 (MQ=255) ccTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCtt > 1:2158787/1‑100 (MQ=255) gCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACCAATCCGGCGATTTCAACATACAACATAATCCCACCTTATTACTCATACCCTTCTa > 1:611210/1‑100 (MQ=255) cGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAtt > 2:94994/1‑100 (MQ=255) cGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAtt < 2:704208/100‑1 (MQ=255) cGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAtt > 1:1911811/1‑100 (MQ=255) aaCCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGa < 2:1112100/100‑1 (MQ=255) aCCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGat < 1:1888856/100‑1 (MQ=255) gCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATg < 1:1932318/100‑1 (MQ=255) gTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACCCATACCCTTCTATTGATATGGAt > 2:1225347/1‑100 (MQ=255) acacATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaa > 1:703048/1‑100 (MQ=255) cacaTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaat > 2:1971384/1‑100 (MQ=255) cacaTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaat < 2:2187283/100‑1 (MQ=255) tccgatctgACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaata < 1:582861/92‑1 (MQ=255) acaTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaata > 1:1563576/1‑100 (MQ=255) caTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaat > 1:1068428/1‑98 (MQ=255) aTCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCt < 2:1147966/100‑1 (MQ=255) aTCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCt < 2:1459353/100‑1 (MQ=255) aTCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCt > 1:2445232/1‑100 (MQ=255) tCGACTGAACGCCGGATATGACAAATCCGGCGCTTTGAACATACTACATAATCCCACCTTATTACTCATACCCTTCTATTGCTCTGGATTAATAATTCt < 1:2064472/99‑1 (MQ=255) ttcaaatcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTaa > 1:1274162/8‑100 (MQ=37) gACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaataagatcgga > 2:582863/1‑93 (MQ=255) aCTGAACGCCGGATATGACAAATCCGGCGATTTGAGCATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAc < 1:109281/100‑1 (MQ=255) cTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAcc < 2:1390109/100‑1 (MQ=255) gAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATAGCCTTATATTGATATGGCTTAATATTTCTTAACCCa > 1:2151422/1‑100 (MQ=255) gAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCa > 1:2432254/1‑100 (MQ=255) gAACGCCGGATATGACAAATCCGGCGATGTGAACATACAACATAATTCCACTTTATTACTCATACCCTTCTATAGCGATGGATTAATAATTCTTAAcaca > 1:1936160/1‑100 (MQ=255) aatcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCa > 1:1851047/4‑100 (MQ=38) atcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATACTTCTTAACCCaa > 1:472460/3‑100 (MQ=255) tcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCaaa > 2:414998/2‑100 (MQ=255) cGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATGCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCaaaa > 1:400101/1‑100 (MQ=255) cGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCaaaa > 2:460739/1‑100 (MQ=255) cGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCaaaa > 1:34191/1‑100 (MQ=255) gCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAAt > 1:939767/1‑100 (MQ=255) gCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAAt > 1:954254/1‑100 (MQ=255) gCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCAACTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAAt > 1:1438346/1‑100 (MQ=255) gCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAAt > 2:1711627/1‑100 (MQ=255) gATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATggct > 2:849198/1‑98 (MQ=255) tGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGAc > 2:882281/1‑100 (MQ=255) gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTa > 1:2023675/1‑96 (MQ=255) gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACt > 2:1963865/1‑100 (MQ=255) gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACt > 2:1589962/1‑100 (MQ=255) gACAAATCCGGCGATTTGAACATACAACATAATCCCACCATATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACt < 1:1539732/100‑1 (MQ=255) gACAAATCAGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACt > 1:2445986/1‑100 (MQ=255) cAAATCCGGCGATTTGAACATACATCATAATCCCACCGTATTACTCATACGCTTCTATTGATATGGATTCATGATTCTTAACCCAAAATGGGCAGAcgc > 1:1376030/1‑97 (MQ=255) cAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTcc > 2:1345306/1‑100 (MQ=255) | CCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCC > NC_000913/1657803‑1657998 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |