The Perl Toolchain Summit needs more sponsors. If your company depends on Perl, please support this very important event.

NAME

Bio::Pipeline::Comparison::Generate::VCFWriter - Create a VCF with the differences between a reference and a single evolved genome. Outputs a gzipped VCF file and a tabix file.

VERSION

version 1.123050

SYNOPSIS

Create a VCF with the differences between a reference and a single evolved genome

use Bio::Pipeline::Comparison::Generate::VCFWriter; my $obj = Bio::Pipeline::Comparison::Generate::VCFWriter->new(output_filename => 'my_snps.vcf.gz'); $obj->add_snp(1234, 'C', 'A'); $obj->add_snp(1234, 'T', 'A'); $obj->create_file();

METHODS

add_snp

Stage a base position in the reference genome, the reference base and a new base for writing to the VCF file.

create_file

Create the VCF file and gzip it. Create an index for the file using tabix.

evolved_name

Optional input parameter. This is the name that goes in the column of the VCF file. The output filename is used as the base for this name by default.

SEE ALSO

AUTHOR

Andrew J. Page <ap13@sanger.ac.uk>

COPYRIGHT AND LICENSE

This software is Copyright (c) 2012 by Wellcome Trust Sanger Institute.

This is free software, licensed under:

  The GNU General Public License, Version 3, June 2007