Skip to content

Commit

Permalink
No need to update reference fasta location
Browse files Browse the repository at this point in the history
  • Loading branch information
yueyaog committed Jun 12, 2024
1 parent 3c00997 commit ed6a906
Showing 1 changed file with 1 addition and 7 deletions.
8 changes: 1 addition & 7 deletions CoverageProfiler/CoverageProfiler.wdl
Original file line number Diff line number Diff line change
Expand Up @@ -47,14 +47,8 @@ workflow coverageProfile {
# Create directories for input, output, reference localization
mkdir input
mkdir output
mkdir reference

# Localize Reference Files and BAM
mv ~{referenceFasta} reference/reference.fasta
mv ~{referenceDict} reference/reference.dict
mv ~{referenceFai} reference/reference.fai
mv ~{alignedBam} input/~{sampleName}.bam

# Index BAM file
gatk BuildBamIndex \
--INPUT input/~{sampleName}.bam
Expand All @@ -63,7 +57,7 @@ workflow coverageProfile {
gatk --java-options "-Xmx~{command_mem_mb}m" DepthOfCoverage \
-L ~{intervals} \
--input input/~{sampleName}.bam \
--reference reference/reference.fasta \
--reference ~{referenceFasta} \
--output output/~{sampleName}

cat output/~{sampleName}.sample_interval_summary | awk 'BEGIN {FS = ","}{print $3}' | tail -n 1 > output/mean_coverage.txt
Expand Down

0 comments on commit ed6a906

Please sign in to comment.