annotate scripts/prepare_statistics.sh @ 32:cf7cd22fc5d3 default tip

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