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.