Mercurial > repos > artbio > artbio_bam_cleaning

diff test-data/pipe.sh @ 0:65d6d2b554b3 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/artbio_bam_cleaning commit adfad19ff505ac7baa3688997bfa9f64243aaace"
author artbio
date Fri, 02 Oct 2020 00:17:33 +0000
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/pipe.sh	Fri Oct 02 00:17:33 2020 +0000
@@ -0,0 +1,9 @@
+input="match_chr21_DBA_974" && \
+genome="chr21.fa" && \
+samtools index $input".bam" && \
+sambamba view -h -t 8 --filter='mapping_quality >= 1 and not(unmapped) and not(mate_is_unmapped)' -f 'bam' $input".bam" \
+| samtools rmdup - - \
+tee $input".filt1.dedup.bam"| bamleftalign --fasta-reference $genome -c --max-iterations "5" - \
+| samtools calmd  -C 50 -b -@ 4 - $genome > $input".filt1.dedup.bamleft.calmd.bam" && \
+sambamba view -h -t 8 --filter='mapping_quality <= 254' -f 'bam' -o $input".filt1.dedup.bamleft.calmd.filt2.bam" $input".filt1.dedup.bamleft.calmd.bam"
+