Test-drive the GATK tools and Best Practices pipelines on Terra
Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.
DepthOfCoverage multiple samples
I have about 300 bam files (whole-genome sequence) and I'm trying to get the DepthOfCoverage per chomosome.
I'm using a script to submit a job for each chromosome, like that:
java -jar /GenomeAnalysisTK-3.5/GenomeAnalysisTK.jar -T DepthOfCoverage -R 'ref_gen.fa' -I samplesBam.list -L $Chr -o $Chr.DepthCov -nt 30 -omitIntervals -omitSampleSummary -omitLocusTable
The analysis started without problems, but after a while a got an error about "many files opened".
Sorry for not having the output, but I deleted before realise I was going to need them here.
Is there a better way to do that? I've tried sample by sample, but them I'll have 300x34 chromosomes files to get the average manually.