File: vcfmultiway

package info (click to toggle)
libvcflib 1.0.12%2Bdfsg-1
  • links: PTS, VCS
  • area: main
  • in suites: sid, trixie
  • size: 70,520 kB
  • sloc: cpp: 39,837; python: 532; perl: 474; ansic: 317; ruby: 295; sh: 254; lisp: 148; makefile: 123; javascript: 94
file content (21 lines) | stat: -rwxr-xr-x 467 bytes parent folder | download | duplicates (3) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
 
#!/bin/bash
# Multiway comparison
#
reference=$1
shift

echo comparing $@

for fileA in $@;
do
    for fileB in $@;
    do
	if [ "$fileA" = "$fileB" ]
	then
	    vcfstats $fileA >$(basename $fileA).stats.txt
	else
	    vcfintersect -r $reference -i $fileA $fileB | vcfstats >$(basename $fileA).common.$(basename $fileB).stats.txt
	    vcfintersect -r $reference -v -i $fileB $fileA | vcfstats >$(basename $fileA).unique.$(basename $fileB).stats.txt
	fi
    done
done