Question: Read depth calling for X chromsome - script modification
gravatar for mbhargava2003
13 months ago by
mbhargava20030 wrote:

Hello Everyone,

I am trying to do read depth calling of chromosomes from bam files of the exome sequencing data. I am still in the early stages of learning bioinformatics. Initially I thought of doing read depth calling for chromosomes 1 through 22 only. I used the following code. I split the job into different individual arrays using the option "-t 1-22". This option doesn't allow me to include the X-chromosome in it.


$ -q long

$ -t 1-22

$ -cwd

$ -j y

$ -l m_mem_free=32g

source /software/scripts/useuse reuse -q Tabix reuse -q Samtools

samtools mpileup -ABQ0 -d 250 -b finnish.sample.list -f cvar/reddy/human_g1k_v37.fasta -l /cvar/reddy/RefSeq_bed/RefSeq_exons_trim_chr${SGE_TASK_ID}.bed -t DP -uv > finnish_rdepths.chr$SGE_TASK_ID.vcf

Then now I want to include X chromosome for read depth calling. I tried to change this code like this. 1) I removed the option "-t 1-22" and launched the same job using qsub. This produced an error - "/cvar/reddy/RefSeq_bed/RefSeq_exons_trim_chrundefined.bed": No such file or directory. 2) I changed the code as follows - since I need the read depths of only X-chromsomes (other chromosomes already done) samtools mpileup -ABQ0 -d 250 -b finnish.sample.list -f cvar/reddy/human_g1k_v37.fasta -l /cvar/reddy/RefSeq_bed/RefSeq_exons_trim_chrX.bed -t DP -uv > finnish_rdepths.chrX.vcf

This produced an output file with all the headers and individual IDs and no data in it.

I am not sure where I am doing anything wrong and would appreciate any suggestions or modification regarding the code I am using to perform read depth calling either for X chromosome alone or 1-22 & X chromosome combined.

ADD COMMENTlink modified 10 weeks ago by Biostar ♦♦ 20 • written 13 months ago by mbhargava20030
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 490 users visited in the last hour