#!/usr/bin/env python
##coding=utf-8
import os
import sys
import argparse
import vcf
parser = argparse.ArgumentParser(description="")
parser.add_argument('-i', '--invcf', help = 'the input vcf file', required = True)
parser.add_argument('-o', '--outvcf', help = 'the output vcf file', required = True)
argv = vars(parser.parse_args())
inF = os.path.abspath(argv['invcf'].strip())
ouF = os.path.abspath(argv['outvcf'].strip())
def parser_vcf(infile,outfile):
samplename=str(str(str(infile).strip().split('/')[-1]).strip().split("_")[0])
vcf_r=vcf.Reader(filename=infile)
for record in vcf_r:
if len(record.FILTER) == 0 :
print("samplename: "+str(samplename))
sample1=samplename+".Mutect2"
sample2=samplename+".vardict"
print("sample1: "+str(record.genotype(sample1)['GT']))
break
# if record.samples[sample2]['DP'] == "." and record.samples[sample2]['GT'] == ".":
parser_vcf(inF,ouF)
Python41
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。
推荐阅读更多精彩内容
- 坦白来说,这一次的小练习,总体还是有点不明所以,看来的对函数以及类的章节,还需要好好的复习一下才行。 不多说了,看...
- 4、不定长参数 不定长参数也叫可变参数。用于不确定调用的时候会传递多少个参数(不传参也可以)的场景。此时,可用包裹...
- 大家好,我是杨数Tos,这是《从零基础到大神》系列课程的第41篇文章,第二阶段的课程:Python基础知识:Pyt...