This is a test version of Biostars. For the public version, visit https://www.biostars.org.
read counts per chromosome arm

I have DNAseq data and trying to get the read counts per chromosome arm. but I do not know the correct way. shall I align to every chromosome separately or the whole genome? how can I count the reads per chromosome arm? I mean what tool can I use for this goal?

genome alignment

1 answer

get a bed for each q/p for each chromosome, eg. see Programmatically convert cytoband location to genomic coordinates

then for each bed record, use samtools view -c to count the reads.

Log in to answer this question.