This function reads the output files from the somatic
and copynumber
programs, which are part of VarScan version 2. The resulting object can be used, in lieu of a “seqz” object, within Sequenza.
VarScan2seqz(varscan.somatic, varscan.copynumber = NULL, normal_var_freq = 0.25)
a data frame from the output file of somatic
.
a data frame from the output file of copynumber
(optional).
defines the variant allele frequency threshold on the normal genotype to filter out false positive heterozygous position.
A data frame with the seqz format, as returned by read.seqz
.
It is possible to run sequenza with the output of VarScan2 version 2.3.6. Earlier versions may also work but have not been tested.