Mercurial > repos > ziru-zhou > bamedit
diff split.sh @ 3:e623e9c7bb23 draft
Uploaded
author | ziru-zhou |
---|---|
date | Tue, 18 Dec 2012 09:20:12 -0500 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/split.sh Tue Dec 18 09:20:12 2012 -0500 @@ -0,0 +1,25 @@ +#!bin/bash + +#purpose: shell script which performs the split option +#author: Ziru Zhou +#date: October, 2012 + +main() +{ + output1="${1}" + output2="${2}" + input="${3}" + + $(samtools view -H ${input} > tmp) + linecount=$(samtools view ${input} | wc -l) + half=$(((${linecount} + 1) / 2)) + samtools view ${input} | shuf | split -d -l ${half} + + $(cat tmp x00 > x00h) + $(cat tmp x01 > x01h) + samtools view -bSo ${output1} x00h + samtools view -bSo ${output2} x01h + + rm tmp x00h x01h x00 x01 +} +main "${@}"