annotate scripts/prepare_statistics.sh @ 1:4c5de3a2c54e

toolshed2
author biomonika <biomonika@psu.edu>
date Sun, 07 Sep 2014 23:12:26 -0400
parents 05c27700e5ca
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
1 #!/bin/bash
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
2 set -e;
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
3
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
4 #variables
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
5 mkdir -p scripts/statistics; #create directory if it doesn't exist
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
6
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
7 echo "Calculation of statistics started."
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
8
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
9 #get statistics about mapped and unmapped reads for every contig
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
10
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
11 idxstats_bam_file ()
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
12 {
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
13 for ARG in "$@";
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
14 do
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
15 samtools idxstats bam/aln_cleaned_sorted_deduplicated_${ARG}_male.bam >scripts/statistics/stat_${ARG} &
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
16 done;
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
17 }
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
18
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
19 idxstats_bam_file mother father daughter son;
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
20
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
21 wait;
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
22 echo "Idxstats successfully created."
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
23
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
24
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
25
05c27700e5ca initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
26