Entering edit mode
Hi all, I used a different reference genome to align instead of GRCh38 and got this error:
result <- dba.blacklist(result)
Genome detected: Hsapiens.1000genomes.hs37d5 Applying blacklist... No blacklist found for BSgenome.Hsapiens.1000genomes.hs37d5 No control reads specified, unable to generate greylist. No intervals removed.
result <- dba.blacklist(result, blacklist = BSgenome.Hsapiens.1000genomes.hs37d5)
Error in dba.blacklist(result, blacklist = BSgenome.Hsapiens.1000genomes.hs37d5) : object 'BSgenome.Hsapiens.1000genomes.hs37d5' not found
Would you have a suggestion? I have this file hg38-blacklist.v2.bed but don't know how to download blacklist file for BSgenome.Hsapiens.1000genomes.hs37d5. Thank you so much!