#!/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、不定长参数 不定长参数也叫可变参数。用于不确定调用的时候会传递多少个参数（不传参也可以）的场景。此时，可用包裹...
