drprg 0.1.1

Drug resistance prediction with reference graphs
Documentation
##fileformat=VCFv4.3
##FILTER=<ID=PASS,Description="All filters passed">
##fileDate==07/12/22
##ALT=<ID=SNP,Description="SNP">
##ALT=<ID=PH_SNPs,Description="Phased SNPs">
##ALT=<ID=INDEL,Description="Insertion-deletion">
##ALT=<ID=COMPLEX,Description="Complex variant, collection of SNPs and indels">
##INFO=<ID=VC,Number=1,Type=String,Description="Type (class) of variant">
##ALT=<ID=SIMPLE,Description="Graph bubble is simple">
##ALT=<ID=NESTED,Description="Variation site was a nested feature in the graph">
##ALT=<ID=TOO_MANY_ALTS,Description="Variation site was a multinested feature with too many alts to include all in the VCF">
##INFO=<ID=GRAPHTYPE,Number=1,Type=String,Description="Type of graph feature">
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
##FORMAT=<ID=MEAN_FWD_COVG,Number=R,Type=Integer,Description="Mean forward coverage">
##FORMAT=<ID=MEAN_REV_COVG,Number=R,Type=Integer,Description="Mean reverse coverage">
##FORMAT=<ID=MED_FWD_COVG,Number=R,Type=Integer,Description="Med forward coverage">
##FORMAT=<ID=MED_REV_COVG,Number=R,Type=Integer,Description="Med reverse coverage">
##FORMAT=<ID=SUM_FWD_COVG,Number=R,Type=Integer,Description="Sum forward coverage">
##FORMAT=<ID=SUM_REV_COVG,Number=R,Type=Integer,Description="Sum reverse coverage">
##FORMAT=<ID=GAPS,Number=R,Type=Float,Description="Number of gap bases">
##FORMAT=<ID=LIKELIHOOD,Number=R,Type=Float,Description="Likelihood">
##FORMAT=<ID=GT_CONF,Number=1,Type=Float,Description="Genotype confidence">
##contig=<ID=ahpC>
##contig=<ID=ethA>
##FILTER=<ID=ld,Description="Kmer coverage on called allele less than 3">
##FILTER=<ID=sb,Description="A strand on the called allele has less than 1.00% of the coverage for that allele">
##FILTER=<ID=lgc,Description="Genotype confidence score less than 5.0">
##FILTER=<ID=frs,Description="Fraction of read support on called allele is less than 0.51">
##INFO=<ID=VARID,Number=.,Type=String,Description="The identifier for the panel variant(s) the record overlaps with">
##INFO=<ID=PREDICT,Number=.,Type=String,Description="The drug resistance prediction(s) for the corresponding VARID(s), where 'R' = resistant, 'S' = susceptible, 'F' = failed, and 'U' = unknown">
##INFO=<ID=OGT,Number=1,Type=String,Description="Original genotype after adjusting for minor allele depth proportions of 0.1">
##INFO=<ID=PDP,Number=R,Type=Float,Description="Proportion of the total position depth found on this allele">
##bcftools_viewVersion=1.15.1+htslib-1.15.1
##bcftools_viewCommand=view -o out3.vcf ERR4811286.drprg.bcf; Date=Thu Dec  8 09:24:07 2022
#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	sample
ahpC	11	83d878f0	GTG	ATG,GTA	.	PASS	VC=PH_SNPs;GRAPHTYPE=SIMPLE;PDP=1,0,0;VARID=ahpC_G-88A;PREDICT=.	GT:MEAN_FWD_COVG:MEAN_REV_COVG:MED_FWD_COVG:MED_REV_COVG:SUM_FWD_COVG:SUM_REV_COVG:GAPS:LIKELIHOOD:GT_CONF	0:44,0,0:50,0,0:44,0,0:50,0,0:44,0,0:50,0,0:0,1,1:-61.6903,-914.886,-914.886:853.196
ahpC	19	9f136482	TCACCTTTGCC	TCACCTTGCC,TCACCTTTGCT,TCC	.	PASS	VC=COMPLEX;GRAPHTYPE=SIMPLE;PDP=0.736842,0,0.263158,0;VARID=ahpC_C-72T;PREDICT=.	GT:MEAN_FWD_COVG:MEAN_REV_COVG:MED_FWD_COVG:MED_REV_COVG:SUM_FWD_COVG:SUM_REV_COVG:GAPS:LIKELIHOOD:GT_CONF	0:41,0,14,0:43,0,16,0:44,0,0,0:50,0,0,0:207,0,44,0:219,0,50,0:0.2,1,0.666667,1:-257.956,-1006.99,-698.615,-1006.99:440.659
ahpC	49	1684c937	CG	CA,TG	.	PASS	VC=PH_SNPs;GRAPHTYPE=SIMPLE;PDP=0.0542169,0.216867,0.728916;VARID=ahpC_C-52T,ahpC_G-51A;PREDICT=U,u	GT:MEAN_FWD_COVG:MEAN_REV_COVG:MED_FWD_COVG:MED_REV_COVG:SUM_FWD_COVG:SUM_REV_COVG:GAPS:LIKELIHOOD:GT_CONF	2:3,18,63:6,18,58:3,18,63:6,18,58:3,36,63:6,37,58:1,0,0:-1168.45,-737.939,-247.18:490.759
ahpC	326	9374804e	G	A	.	PASS	VC=SNP;GRAPHTYPE=SIMPLE;PDP=0.993789,0.00621118	GT:MEAN_FWD_COVG:MEAN_REV_COVG:MED_FWD_COVG:MED_REV_COVG:SUM_FWD_COVG:SUM_REV_COVG:GAPS:LIKELIHOOD:GT_CONF	0:74,0:86,1:109,0:124,0:223,2:258,4:0.333333,1:-103.856,-1213.34:1109.49