Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,652,855:1 | (G)6→7 | 100% | coding (362/846 nt) | sseA → | 3‑mercaptopyruvate sulfurtransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,652,849 | 1 | . | G | 93.9% | 208.5 / 8.4 | 66 | coding (356/846 nt) | sseA | 3‑mercaptopyruvate sulfurtransferase |
Reads supporting (aligned to +/- strand): ref base . (2/2); new base G (33/28); total (36/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01 |
GGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑G‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATC > NC_000913/2652751‑2652944 | ggATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTATAGAGAAAGTGTCGATTCT‑g > 1:500206/1‑100 (MQ=255) gaaTAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑g > 1:1762958/3‑100 (MQ=255) tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑g‑ggg > 2:197487/1‑100 (MQ=255) tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑g‑ggg > 1:762597/1‑100 (MQ=255) gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑g‑gggg‑gg < 1:2462019/100‑2 (MQ=255) cACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑g‑gggg‑ggt > 2:1997750/1‑98 (MQ=255) ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑g‑gggg‑ggtgg < 1:631393/100‑5 (MQ=255) tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGAc > 1:2088841/1‑100 (MQ=255) gTCTATGACGACGGTAAGCTTTTCTCGTCCCCGCGAGCATGGTGGATGCTGCGCACCTTTTGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCa < 2:927829/100‑1 (MQ=255) aTGACGAAGGTAATCTTTCCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCt > 2:1044454/1‑100 (MQ=255) tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCt > 1:2114986/1‑99 (MQ=255) tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTg < 2:1343844/100‑1 (MQ=255) tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTg > 2:2441667/1‑100 (MQ=255) gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTgg > 1:1194734/1‑100 (MQ=255) aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGc < 1:2528029/100‑1 (MQ=255) cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCa < 1:2361739/100‑1 (MQ=255) cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCa < 1:1979385/100‑1 (MQ=255) aaGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCGCCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAg < 2:163347/99‑1 (MQ=255) aaGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTGGCTTCTGG‑GGGG‑GTGGACTTGAAGGCTGGc > 2:2532503/1‑97 (MQ=255) aaGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCa > 2:2532076/1‑98 (MQ=255) aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAgataggaa > 2:1330722/1‑90 (MQ=255) aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAgcg < 2:1799859/100‑1 (MQ=255) gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGcgcg > 2:1340534/1‑100 (MQ=255) tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑G‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCgat > 1:432168/1‑100 (MQ=255) aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCgat > 1:1413816/1‑100 (MQ=255) aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCgatg > 1:2125029/1‑100 (MQ=255) tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCgatga > 2:1851257/1‑100 (MQ=255) tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATctgc < 2:1615082/100‑1 (MQ=255) ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATctgct > 1:136415/1‑100 (MQ=255) gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa > 1:125641/1‑100 (MQ=255) aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg < 2:2573835/100‑1 (MQ=255) gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTTGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt > 1:166310/1‑100 (MQ=255) gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt > 2:845377/1‑100 (MQ=255) aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTTGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc < 2:166308/100‑1 (MQ=255) aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc < 1:748094/100‑1 (MQ=255) ggtggATGCTGCGCACATTTGGTGTAGAGCAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc > 1:2054589/1‑98 (MQ=255) gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt > 2:448577/1‑100 (MQ=255) tggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑G‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag < 2:1139521/100‑1 (MQ=255) ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag < 2:336279/100‑1 (MQ=255) ntgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑G‑GGGGCGTTGACTTGCGGGCTGGCAGCGCTATTATGTTGGGTTAGTAGAAGGTGCATTag > 2:283864/2‑98 (MQ=255) aTGCTGCGCGCCTTTTGTGTAGAGAAAGTGGCGATTCT‑GTGGGG‑TTGGACTTGAAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag < 2:2263609/100‑1 (MQ=255) aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag < 2:1320832/100‑1 (MQ=255) tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc < 2:1906165/100‑1 (MQ=255) gctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCt < 2:1689177/100‑1 (MQ=255) gcgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc > 1:2533427/1‑100 (MQ=255) ggcACCTTTGGTGTAGAGAAAGTGTAGATTTTAG‑GGGG‑GAGGCCTTGGAGGCTGGCAGCGAAAGGATATGATGTTAGACGAAGGTGGAGTATAGTTg > 2:958593/2‑97 (MQ=255) aCCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GGGGGGTTGCAGGATGGCAGCGGGATGATCGGCGGTTAGAAGAAGGTGCAGTAGAGCTGCCGGaa > 2:1221413/1‑100 (MQ=255) ccTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg < 2:2543229/100‑1 (MQ=255) ccTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg > 1:94110/1‑100 (MQ=255) ttGGTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag < 1:1223922/100‑1 (MQ=255) ggTGTAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag < 2:1220871/100‑1 (MQ=255) tgtAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt > 1:2123423/1‑100 (MQ=255) gtAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt < 1:497476/100‑1 (MQ=255) gtAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt < 1:2001068/100‑1 (MQ=255) tAGAGAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa < 2:1213393/100‑1 (MQ=255) agagAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa < 1:1020154/100‑1 (MQ=255) agagAAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa < 1:1789618/100‑1 (MQ=255) gagaAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTATAGCTGCCGGAAGGAGAGTTTAAc < 2:2438177/100‑1 (MQ=255) gagaAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc > 1:1680762/1‑100 (MQ=255) gagaAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc > 1:2292486/1‑100 (MQ=255) gagaAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc > 2:1585818/1‑100 (MQ=255) agaAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg < 1:2427668/100‑1 (MQ=255) gaAAGTGTCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc < 2:237842/100‑1 (MQ=255) gtgtCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccgtg > 2:158085/1‑98 (MQ=255) tgtCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGt > 2:11731/1‑100 (MQ=255) tgtAGATTCTGG‑GGGG‑GTGCACTTGCAGTCTGGGAGCGCGATGATCTGCTGGTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGGTTAAcgccgc > 2:463683/1‑98 (MQ=255) tCGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt < 2:1619938/100‑1 (MQ=255) cGATTCTGG‑GGGG‑GTGGACTTGCGGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa < 2:156699/100‑1 (MQ=255) cGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa > 2:1725329/1‑100 (MQ=255) cGATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa > 1:1416320/1‑100 (MQ=255) gATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa > 1:2360401/1‑100 (MQ=255) gATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa > 1:556161/1‑100 (MQ=255) gATTCTGG‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa > 1:887898/1‑100 (MQ=255) ctt‑g‑gtgg‑GTTGACTTGCGGGCTGTCAGCTTGATGATTTGTTGTGAGGAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTCAACGCCTCGTTTAATc > 2:2338858/3‑98 (MQ=255) | GGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑G‑GGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATC > NC_000913/2652751‑2652944 |
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: ‑ |