diff --git a/alignment/bamqc.sh b/alignment/bamqc.sh index 566eea451e0332e5868311645fcdd198ac22acd3..a556ce1b874614a5307e7794cb0bfa1d5f6a50a8 100644 --- a/alignment/bamqc.sh +++ b/alignment/bamqc.sh @@ -54,6 +54,6 @@ if [[ $nuctype == 'dna' ]]; then java -Xmx4g -jar $PICARD/picard.jar CollectInsertSizeMetrics INPUT=${sbam} HISTOGRAM_FILE=${pair_id}.hist.ps REFERENCE_SEQUENCE=${index_path}/genome.fa OUTPUT=${pair_id}.hist.txt samtools view -b -q 1 ${pair_id}.ontarget.bam | bedtools coverage -sorted -hist -g ${index_path}/genomefile.txt -b stdin -a ${bed} > ${pair_id}.mapqualcov.txt bedtools coverage -sorted -g ${index_path}/genomefile.txt -a ${bed} -b ${sbam} -hist > ${pair_id}.covhist.txt - perl $baseDir/scripts/calculate_depthcov.pl ${pair_id}.covhist.txt + perl $baseDir/calculate_depthcov.pl ${pair_id}.covhist.txt grep ^all ${pair_id}.covhist.txt > ${pair_id}.genomecov.txt fi