Description Usage Arguments Value Examples
counts number of variants in certain affect categories use with plyr::ldply() on VCFList objects to create a table.
1 | variantCounts(data, unique = TRUE)
|
data |
tidyVCF with 'EFF' field parsed |
unique |
logical, if true only unique variants will be counted |
a single row dataframe with variant counts
1 2 3 4 5 6 7 8 9 10 | geneInfo <- getGeneInfo(genes = c("AT3G62980", "AT3G26810"))
## download and compile a list of VCF files
myVCFList <- VCFList(geneInfo)
## tabultate unique variants
plyr::ldply(myVCFList, variantCounts, unique=TRUE, .id="transcript_ID")
# tabulate total counts of variants relative to the reference sequence.
plyr::ldply(myVCFList, variantCounts, unique=TRUE, .id="transcript_ID")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.