Question: Read depth calling for X chromsome - script modification
gravatar for mbhargava2003
7 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 written 7 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: 477 users visited in the last hour