Description Usage Arguments Value
View source: R/get_read_counts.R
Uses samtools pileup to get the read counts for each base in the genomic position specified
1 2 | get_read_counts(chr, pos, bam, tag = "", min_base_quality = 20,
max_depth = 1e+05, min_mapq = 30)
|
chr |
chromosome name |
pos |
genomic coordinate |
bam |
path to bam file |
tag |
the RG tag if the bam has more than one sample |
min_base_quality |
minimum base quality for a read to be counted |
max_depth |
maximum depth above which sampling will happen |
min_mapq |
the minimum mapping quality for a read to be counted |
a list, number of reads for each of the four basepairs
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.