使用snpflip校正基因组正负链
欢迎来到"bio生物信息"的世界 公众号有位小朋友留言希望我写一下 snpflip 的使用。 重新翻了一下之前写的推文,发现我是用意念写完的 snpflip ,总感觉我写过这个工具的使用了。 因此这里重新补充一下 snpflip 工具的使用。 关于 snpflip 的使用背景我就不多介绍了,之前的文章有写过 数据合并,踩不完的坑 下面直接讲一下怎么使用。 这个工具安装和使用都非常简单,老少皆懂。 第一步:下载、安装snpflip wget https://files.pythonhosted.org/packages/5d/58/c4e3427cd307c29c92631f0d2bc0ee599687d1d29019a390d9786ec6a44e/snpflip-0.0.6.tar.gz tar -zxvf snpflip-0.0.6.tar.gz 第二步:下载参考基因组 wget http://ftp-trace.ncbi.nih.gov/1000genomes/ftp/technical/reference/human_g1k_v37.fasta.gz gunzip human_g1k_v37.fasta.gz 第三步:校正正负链 查找混有正负链的SNP: /snpflip-0.0.6/bin/snpflip -b file.bim -f human_g1k_v37