If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We encourage our fourm members to be more involved, jump in and help out your fellow researchers with their questions. GATK forum is a community forum and helping each other with using GATK tools and research is the cornerstone of our success as a genomics research community.We appreciate your help!
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.
We will be out of the office on October 14, 2019, due to the U.S. holiday. We will return to monitoring the forum on October 15.
DepthofCoverage: More rows in output that in target file
I have used the following commands using DepthOfCoverage tool with two different bed files:
java -jar GenomeAnalysisTK.jar -T DepthOfCoverage -R ucsc_hg19.fa -I WT_recalibrated.bam -L coverage_summary.bed -ct 1 -ct 10 -ct 20 -ct 30 -ct 50 -ct 100 -o WT_cov
The line count for the input and output:
$wc -l WT_cov.sample_interval_summary 4988 WT_cov.sample_interval_summary $ wc -l coverage_summary.bed 10585 coverage_summary.bed
In the other case:
java -jar GenomeAnalysisTK.jar -T DepthOfCoverage -R ucsc_hg19.fa -I WT_recalibrated.bam -L exon.bed -ct 1 -ct 10 -ct 20 -ct 30 -ct 50 -ct 100 -o WT_exon
Line count for the input and output:
$ wc -l WT_exon.sample_interval_summary 5065 WT_exon.sample_interval_summary $ wc -l exon.bed 5065 exon.bed
The input in both the cases is of the standard format as shown below:
chr1 6529578 6529755 chr1 6530273 6530442 chr1 6530543 6530721 chr1 6530773 6530980 chr1 6531028 6531730 chr1 6531768 6531914 chr1 6532563 6532713 chr1 6533023 6533273
Could anyone help to interpret the discrepancy between number of target regions in bed file and _interval_summary file in the above two cases?