Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
GATK DepthofCoverage Error bwubb Bioinformatics 5 04-10-2012 10:29 AM
GATK DepthOfCoverage at high depth lletourn Bioinformatics 2 03-28-2012 07:51 AM
GATK / DepthofCoverage nguyendofx Bioinformatics 0 11-07-2011 10:21 AM
GATK - DepthOfCoverage giverny Bioinformatics 2 09-14-2011 01:48 PM
GATK depthofcoverage foxyg Bioinformatics 1 08-21-2010 09:22 AM

Thread Tools
Old 07-01-2012, 05:11 PM   #1
Location: Shanghai

Join Date: Dec 2009
Posts: 16
Question GATK depthofcoverage problem


I am using GATK depthofcoverage to calculate my exome-seq data:

 nohup java -jar ~/GATK/GenomeAnalysisTK-1.3-23-g13905c0/GenomeAnalysisTK.jar \
-T DepthOfCoverage -o depthofcoverage_sample1 -I Sample_1.bam \
-ct 10 -ct 20 -ct 30 -R hg19.fa -geneList Genome_UCSC.refseq -omitBaseOutput &
But the output file depthofcoverage_sample1.sample_gene_summary is empty.

I thought it was because I didn't sort my refseq file, whose first few lines are as follow:

#bin    name    chrom   strand  txStart txEnd   cdsStart        cdsEnd  exonCount       exonStarts      exonEnds        score   name2   cdsStartStat    cdsEndStat      exonFrames
0       NM_032291       chr1    +       66999824        67210768        67000041        67208778        25      66999824,67091529,67098752,67101626,67105459,67108492,67109226,67126195,67133212,67136677,67137626,67138963,67142686,67145360,67147551,67154830,67155872,67161116,67184976,67194946,67199430,67205017,67206340,67206954,67208755,       67000051,67091593,67098777,67101698,67105516,67108547,67109402,67126207,67133224,67136702,67137678,67139049,67142779,67145435,67148052,67154958,67155999,67161176,67185088,67195102,67199563,67205220,67206405,67207119,67210768,       0       SGIP1   cmpl    cmpl    0,1,2,0,0,0,1,0,0,0,1,2,1,1,1,1,0,1,1,2,2,0,2,1,1,
1       NM_001080397    chr1    +       8384389 8404227 8384389 8404073 8       8384389,8385357,8385877,8390268,8395496,8397875,8399552,8403806,        8384786,8385450,8386102,8390996,8395650,8398052,8399758,8404227,        0       SLC45A1 cmpl    cmpl    0,1,1,1,0,1,1,0,
So I downloaded the GATK file and typed the following command:

 perl Genome_UCSC.refseq ~/ref/UCSC/hg19/hg19.fa.fai
then I got an error:
Can not open temporary file 1045: Too many open files at line 95, <$INPUT> line 1021.
What am I doing wrong?

Last edited by gary; 07-01-2012 at 05:17 PM.
gary is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 07:13 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO