Masking Y chromosomes
1
0
Entering edit mode
6 months ago
lamia_203 ▴ 90

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:

  xyalign 
    --PREPARE_REFERENCE 
    --ref reference.fasta   
    --output_dir sample1_output 
    --sample_id sample1 
    --cpus 4 
    --reference_mask mask.bed      
    --x_chromosome chrX 
    --y_chromosome chrY

Any help understanding how to generate the mask.bed would be great! Many thanks.

genome whole generation sequencing wgs bash • 305 views
ADD COMMENT
3
Entering edit mode
6 months ago

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).

ADD COMMENT

Login before adding your answer.

Traffic: 1949 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6