powered by
Extract read counts from VCF file of a single sample.
extractCoverageFromVcf(vcfFileName, ADFieldIndex = 2)
A data.frame contains four columns: chromosomes, positions, reference allele count, alternative allele count.
Path of the VCF file.
Index of the AD field of the sample field. For example, if the format is "GT:AD:DP:GQ:PL", the AD index is 2 (by default).
vcfFile = system.file("extdata", "PG0390-C.test.vcf.gz", package = "DEploid") PG0390 = extractCoverageFromVcf(vcfFile)
Run the code above in your browser using DataLab