搜尋
  • Prisma_J

HowTo : SNP array to vcf

生物資訊界已經廣泛使用vcf格式作為紀錄DNA變異的格式

這個格式包含SNV/indel/SV/CNV


而最近遇到一個資料轉換的問題跟大家分享解法


動機:

三級分析例如GWAS,起始資料格式是vcf格式

然而illumina SNP array的變異檔案格式是gtc

我們是否可以將gtc格式轉換成vcf格式,以利後續三級分析



GTCtoVCF

illumina github上提供一個小工具:GTCtoVCF

https://github.com/Illumina/GTCtoVCF

這是一個python 2.7的工具,請照網站說明安裝

由於幾個套件不能在Windows版本 conda安裝

這次環境在Linux環境執行


相依套件

bash Miniconda2-latest-Linux-x86_64.sh 
conda install -c miniconda numpy=1.11.2
conda install -c bioconda pyvcf=0.6.8
conda install -c bioconda pysam=0.9.0

下載

git clone https://github.com/Illumina/GTCtoVCF.git

其他準備資料

由於參考序列有一些規定,為了避免麻煩直接用內建腳本下載

## 下載參考序列
bash scripts/download_reference.sh genome.fa

GTC轉換VCF

不用擔心測試資料來源,在test資料夾中有測試用資料

tests/data/RegressionGTCDirectory/input 中有兩份gtc檔案

除了gtc檔案及參考序列外,還需要manifest(tests/data/small_manifest.csv)

python gtc_to_vcf.py \
--gtc-paths tests/data/RegressionGTCDirectory/input/ \
--manifest-file tests/data/small_manifest.csv \
--genome-fasta-file genome.fa \
--output-vcf-path tests/data/RegressionGTCDirectory/output/ 

檢查輸出

ls tests/data/RegressionGTCDirectory/output/








恭喜你完成VCF檔案轉換,接下來就可以進行GWAS 或 變異評估解讀等等三級分析


Reference

https://github.com/Illumina/GTCtoVCF












Ⓒ 2020 BY PRISMA BIOTECH. ALL RIGHTS RESERVED.

瑞林生物科技股份有限公司

​Prisma Biotech Corporation