Hi, I am working on masking chromosome on the Y chromosome. I can see that I need a mask.bed file for input but do not understand how to generate the mask.bed file itself. Is the mask.bed file the coordinates of the chromosome? This is the code I would be using:
I have used bedtools maskfasta for similar purposes, but the rationale behind should be the same: if you want to mask a region, then you must define that region. And of course, if you want to mask the entire Y chromosome, then the bed file should contain the Y chromosome coordinates, which you may find in the reference.fasta.fai file if you have your fasta reference already indexed (and if not I would suggest you to do so in order to perform the masking properly).