文件名称:样品
文件大小:10.6MB
文件格式:ZIP
更新时间:2024-03-04 02:36:58
Shell
样品 培训数据工作流程 1000个基因组高覆盖度补习班 下载在high_coverage索引上列出的Cram索引(.crai)ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/data_collections/1000_genomes_project/1000genomes.high_coverage.GRCh38DH.alignment.index 与GRCh38参考基因组比对 获取SV呼叫集VCF(使用GRCh38参考) 提取DEL / Non-DEL区域用于训练集 bcftools view -i ' SVTYPE="DEL" ' $VCF \ | bcftools query -f ' %CHROM\t%POS\t%INFO/END[\t%SAMPLE,%GT]\n ' \ | python sample_del.py > $