annotate test-data/convert.gvcf.vcf @ 10:7887b9d94c05 draft

"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bcftools commit 851f81495c875ac09d936537ffd2b32e6af2c8c5"

##fileformat=VCFv4.1
##FILTER=<ID=PASS,Description="All filters passed">
##FORMAT=<ID=GQX,Number=1,Type=Integer,Description="Minimum of {Genotype quality assuming variant position,Genotype quality assuming non-variant position}">
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
##FORMAT=<ID=GQ,Number=1,Type=Float,Description="Genotype Quality">
##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Filtered basecall depth used for site genotyping">
##FORMAT=<ID=DPF,Number=1,Type=Integer,Description="Basecalls filtered from input prior to site genotyping">
##FORMAT=<ID=AD,Number=.,Type=Integer,Description="Allelic depths for the ref and alt alleles in the order listed. For indels this value only includes reads which confidently support each allele (posterior prob 0.999 or higher that read contains indicated allele vs all other intersecting indel alleles)">
##FORMAT=<ID=DPI,Number=1,Type=Integer,Description="Read depth associated with indel, taken from the site preceding the indel.">
##INFO=<ID=END,Number=1,Type=Integer,Description="End position of the region described in this record">
##INFO=<ID=BLOCKAVG_min30p3a,Number=0,Type=Flag,Description="Non-variant site block. All sites in a block are constrained to be non-variant, have the same filter value, and have all sample values in range [x,y], y <= max(x+3,(x*1.3)). All printed site block sample values are the minimum observed in the region spanned by the block">
##INFO=<ID=SNVSB,Number=1,Type=Float,Description="SNV site strand bias">
##INFO=<ID=SNVHPOL,Number=1,Type=Integer,Description="SNV contextual homopolymer length">
##INFO=<ID=CIGAR,Number=A,Type=String,Description="CIGAR alignment for each alternate indel allele">
##INFO=<ID=RU,Number=A,Type=String,Description="Smallest repeating sequence unit extended or contracted in the indel allele relative to the reference. RUs are not reported if longer than 20 bases.">
##INFO=<ID=REFREP,Number=A,Type=Integer,Description="Number of times RU is repeated in reference.">
##INFO=<ID=IDREP,Number=A,Type=Integer,Description="Number of times RU is repeated in indel allele.">
##FILTER=<ID=IndelConflict,Description="Locus is in region with conflicting indel calls">
##FILTER=<ID=SiteConflict,Description="Site genotype conflicts with proximal indel call. This is typically a heterozygous SNV call made inside of a heterozygous deletion">
##FILTER=<ID=LowGQX,Description="Locus GQX is less than 30 or not present">
##FILTER=<ID=HighDPFRatio,Description="The fraction of basecalls filtered out at a site is greater than 0.3">
##FILTER=<ID=HighSNVSB,Description="SNV strand bias value (SNVSB) exceeds 10">
##FILTER=<ID=HighREFREP,Description="Locus contains an indel allele occurring in a homopolymer or dinucleotide track with a reference repeat greater than 8">
##FILTER=<ID=HighDepth,Description="Locus depth is greater than 3x the mean chromosome depth">
##contig=<ID=22,length=450>
##SnvTheta=0.001
##IndelTheta=0.0001
##MaxDepth_1=114.18
##MaxDepth_10=131.73
##MaxDepth_11=117.27
##MaxDepth_12=116.97
##MaxDepth_13=102.24
##MaxDepth_14=101.55
##MaxDepth_15=95.22
##MaxDepth_16=111.33
##MaxDepth_17=112.59
##MaxDepth_18=121.86
##MaxDepth_19=111.12
##MaxDepth_2=121.83
##MaxDepth_20=111.24
##MaxDepth_21=98.43
##MaxDepth_22=76.23
##MaxDepth_3=120.09
##MaxDepth_4=124.50
##MaxDepth_5=119.82
##MaxDepth_6=122.22
##MaxDepth_7=120.27
##MaxDepth_8=120.45
##MaxDepth_9=102.48
##MaxDepth_M=7005.66
##MaxDepth_X=61.05
##MaxDepth_Y=37.17
##FILTER=<ID=IndelSizeFilter,Description="Indel is outside reportable size range. Insertion range: [1,3], Deletion range: [1,11]">
##gvcftools_version="0.16"
##FILTER=<ID=HAPLOID_CONFLICT,Description="Locus has heterozygous genotype in a haploid region.">
##FORMAT=<ID=OPL,Number=.,Type=Integer,Description="Original PL value before ploidy correction">
##INFO=<ID=phastCons,Number=0,Type=Flag,Description="overlaps a phastCons element">
##INFO=<ID=AA,Number=1,Type=String,Description="Ancestral Allele, ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/pilot_data/technical/reference/ancestral_alignments/README">
##INFO=<ID=AF,Number=A,Type=String,Description="1000 Genomes Allele Frequency based on AC/AN; Format: Allele:AlleleFrequency">
##INFO=<ID=AMR_AF,Number=A,Type=String,Description="1000 Genomes Allele Frequency for samples from AMR population based on AC/AN; Format: Allele:AlleleFrequency">
##INFO=<ID=ASN_AF,Number=A,Type=String,Description="1000 Genomes Allele Frequency for samples from ASN population based on AC/AN; Format: Allele:AlleleFrequency">
##INFO=<ID=AFR_AF,Number=A,Type=String,Description="1000 Genomes Allele Frequency for samples from AFR population based on AC/AN; Format: Allele:AlleleFrequency">
##INFO=<ID=EUR_AF,Number=A,Type=String,Description="1000 Genomes Allele Frequency for samples from EUR population based on AC/AN; Format: Allele:AlleleFrequency">
##INFO=<ID=CLNACC,Number=.,Type=String,Description="Accession and version number assigned by ClinVar to the genotype/phenotype relationship. Multiple entries for an allele are pipe-delimited">
##INFO=<ID=CLNSIG,Number=.,Type=String,Description="String that describes the clinical significance. Possible values: unknown, untested, non-pathogenic, probable-non-pathogenic, probable-pathogenic, pathogenic, drug-response, histocompatibility, other. Multiple values are pipe-delimited">
##INFO=<ID=CLNORIGIN,Number=.,Type=String,Description="String that describes the origin of the variant allele. Possible values: unknown, germline, somatic, inherited, paternal, maternal, de-novo, biparental, uniparental, not-tested, test-inconclusive, other. Multiple values for an allele are pipe-delimited">
##INFO=<ID=CLNSRC,Number=.,Type=String,Description="Variant clinical source or channel. Multiple values for an allele are pipe-delimited">
##INFO=<ID=CLNSRCID,Number=.,Type=String,Description="Identifier used by source defined in CLNSRC. Multiple values are pipe-delimited">
##INFO=<ID=CLNGENEINFO,Number=.,Type=String,Description="Gene symbol(s) and NCBI GeneID(s). The gene symbol and ID are delimited by a colon and multiple such pairs are pipe-delimited, Example SYMBOL1:GeneID1|SYMBOL2:GeneID2">
##INFO=<ID=CLNDBN,Number=.,Type=String,Description="Disease name used by the database specified by CLNSRC. Values corresponding to each CLNACC entry are pipe-delimited">
##INFO=<ID=CLNDSDB,Number=.,Type=String,Description="Colon-delimited list of variant disease database name(s). Multiple values from a single database are pipe-delimited">
##INFO=<ID=CLNDSDBID,Number=.,Type=String,Description="Colon-delimited list of variant disease database identifier(s). Multiple values from a single database are pipe-delimited">
##INFO=<ID=CSQ,Number=A,Type=String,Description="Consequence type as predicted by VEP. Format: Allele|Gene|Feature|Feature_type|Consequence|cDNA_position|CDS_position|Protein_position|Amino_acids|Codons|Existing_variation|EXON|INTRON|HGNC|MOTIF_NAME|MOTIF_POS|HIGH_INF_POS|MOTIF_SCORE_CHANGE|DISTANCE|CANONICAL|SIFT|PolyPhen|GMAF|ENSP|DOMAINS|CCDS|HGVSc|HGVSp|CELL_TYPE">
#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	SAMPLE99
22	1	.	N	.	0	LowGQX	END=9;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	.:.:0:0
22	10	.	C	.	0	LowGQX	END=20;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:5:2:0
22	21	.	C	.	0	LowGQX	END=26;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:10:4:0
22	27	.	C	.	0	LowGQX	END=42;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:15:6:0
22	43	.	G	.	0	LowGQX	END=50;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:21:8:0
22	51	.	C	.	0	PASS	END=55;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:30:11:0
22	56	.	G	.	0	PASS	END=72;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:42:15:0
22	73	.	T	.	0	PASS	END=85;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:54:19:0
22	86	.	G	C	23	LowGQX	SNVSB=0;SNVHPOL=2	GT:GQ:GQX:DP:DPF:AD	0/1:56:23:22:0:16,6
22	87	.	T	.	0	PASS	END=101;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:69:24:0
22	102	.	A	.	0	PASS	END=140;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:84:29:0
22	141	.	G	.	0	PASS	END=185;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:90:31:0
22	186	.	T	.	0	PASS	END=187;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:93:32:2
22	188	.	T	.	0	PASS	END=204;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:102:35:0
22	205	.	T	.	0	PASS	.	GT:GQX:DP:DPF	0/0:72:36:0
22	206	.	T	.	0	PASS	END=231;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:87:30:0
22	232	.	A	.	0	PASS	END=249;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:66:23:0
22	250	.	G	.	0	PASS	END=257;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:87:30:0
22	258	.	A	.	0	PASS	END=269;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:114:39:0
22	270	.	G	.	0	PASS	END=279;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:150:51:0
22	280	.	A	.	0	PASS	END=314;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:166:63:0
22	315	.	C	.	0	PASS	END=316;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:220:74:0
22	317	.	T	.	0	HighDepth	END=342;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:229:77:0
22	343	.	T	.	0	HighDepth	END=377;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:241:81:0
22	378	.	T	.	0	HighDepth	END=384;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:200:75:0
22	385	.	G	.	0	PASS	END=388;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:217:73:0
22	389	.	C	.	0	HighDepth	END=390;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:223:75:0
22	391	.	T	.	0	PASS	.	GT:GQX:DP:DPF	0/0:223:75:1
22	392	.	T	.	0	HighDepth	END=397;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:223:75:0
22	398	.	T	.	0	PASS	END=420;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:178:60:0
22	421	.	C	.	0	PASS	END=450;BLOCKAVG_min30p3a	GT:GQX:DP:DPF	0/0:142:54:0