SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



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

Reply
 
Thread Tools
Old 07-01-2012, 06:11 PM   #1
gary
Member
 
Location: Shanghai

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

Hi

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

Code:
 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:

Code:
#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 sortByRef.pl file and typed the following command:

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

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

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 05:40 PM.


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