aryarm / varCA

Use an ensemble of variant callers to call variants from ATAC-seq data
MIT License
23 stars 7 forks source link

skip rmdups and peak calling steps if input is BAM #17

Closed aryarm closed 4 years ago

aryarm commented 4 years ago

Instead of FASTQs, users can provide their own BAM files as input to the pipeline. However, the current behavior still calls samtools to remove duplicates from the BAM files (and MACS2 to call peaks from them).

I could have it skip the duplicate removal steps if a BAM file is provided, instead. And I could have it skip the peak calling step if a BED file is provided with the BAM file.