VCF tools¶
- Intersect two VCF files
- uses vcfPytools’ intersect command to generate the intersection of two VCF files. Two input files are required as input and the intersection of these two files is generated and sent to the output. These files must be sorted by genomic coordinate to function correctly, although the reference sequence order is no important. The intersection can be calculated on two VCF files or a VCF and a BED file. If the priority file argument is set (this must be equal to one of the input VCF files), then the record written to the output will come from this file. If this argument is not set, the record with the highest quality is written out.
- Annotate a VCF file (dbSNP, hapmap)
- uses vcfPytools’ annotate command annotate a VCF file. Currently, either a hapmap or a dbsnp file should be provided, not both. dbSNP option will annotate the VCF file with dbSNP rsid values. The input dbSNP file must also be in VCF v4.0 format. Only dbSNP entries with VC=SNP are included. hapmap option will annotate the VCF file info string to include HM3 if the record is included hapmap. If the ref/alt values do not match the hapmap file, the info string will be populated with HM3A.
- Filter a VCF file
- uses vcfPytools’ filter command. Quality option will check the variant quality for each record and if it is below the defined value, the filter field will be populated with the filter entry Q[value]. Any value in the info string can be used for filtering by using the ‘Filter by info’ option. This option takes three values: the info string tag, the cutoff value and whether to filter out those records with less than (lt) or greater than (gt) this value.
- Extract reads from a specified region
- uses vcfPytools’ extract command to extract reads from a specified region of a VCF file.
- VCF to MAF Custom Track for display at UCSC
- converts a Variant Call Format (VCF) file into a Multiple Alignment Format (MAF) custom track file suitable for display at genome browsers.
- Extract consensus Sequence from VCF
- extracts Consensus sequence from VCF file.
- BCFtools view
- converts BCF format to VCF format.