-
Notifications
You must be signed in to change notification settings - Fork 133
FindAVariation
##Motivation
Find a specific mutation in a list of VCF files
##Compilation
- java compiler SDK 1.8 http://www.oracle.com/technetwork/java/index.html (NOT the old java 1.7 or 1.6) . Please check that this java is in the
${PATH}
. Setting JAVA_HOME is not enough : (e.g: https://github.com/lindenb/jvarkit/issues/23 ) - GNU Make >= 3.81
- curl/wget
- git
- xsltproc http://xmlsoft.org/XSLT/xsltproc2.html (tested with "libxml 20706, libxslt 10126 and libexslt 815")
$ git clone "https://github.com/lindenb/jvarkit.git"
$ cd jvarkit
$ make findavariation
by default, the libraries are not included in the jar file, so you shouldn't move them (https://github.com/lindenb/jvarkit/issues/15#issuecomment-140099011 ). You can create a bigger but standalone executable jar by adding standalone=yes
on the command line:
$ git clone "https://github.com/lindenb/jvarkit.git"
$ cd jvarkit
$ make findavariation standalone=yes
The required libraries will be downloaded and installed in the dist
directory.
The a file local.mk can be created edited to override/add some definitions.
For example it can be used to set the HTTP proxy:
http.proxy.host=your.host.com
http.proxy.port=124567
##Synopsis
$ java -jar dist/findavariation.jar [options] (stdin|file)
- -o|--output (OUTPUT-FILE) Output file. Default:stdout.
- -p|--position (CHROM:POS) A list of 'chrom/position'Multiple calls to this option should end with double hyphen : --
- -f|--posfile (VALUE) Add this file containing chrom:positionMultiple calls to this option should end with double hyphen : --
- -h|--help print help
- -version|--version show version and exit
##Source Code
Main code is: https://github.com/lindenb/jvarkit/blob/master/src/main/java/com/github/lindenb/jvarkit/tools/misc/FindAVariation.java
The program will try to fix the chromosome names according to the dictionary in the VCF headers.
$ find ./ -name "*.vcf" -o -name "*.vcf.gz" |\
java -jar dist/findamutation.jar -p "chr1:1234"
htsjdk/testdata/htsjdk/samtools/intervallist/IntervalListFromVCFTestManual.vcf 1 8216713 8216713 yossi-1 NA12878 HET A G
htsjdk/testdata/htsjdk/samtools/intervallist/IntervalListFromVCFTestManual.vcf 1 8216713 8216713 yossi-1 NA12891 HET A G
htsjdk/testdata/htsjdk/samtools/intervallist/IntervalListFromVCFTestManual.vcf 1 8216713 8216713 yossi-1 NA12892 HET A G
htsjdk/testdata/htsjdk/samtools/intervallist/IntervalListFromVCFTestManual.vcf 2 2 2 . NA12878 HOM_REF C C
htsjdk/testdata/htsjdk/samtools/intervallist/IntervalListFromVCFTestManual.vcf 2 2 2 . NA12891 HET C T
htsjdk/testdata/htsjdk/samtools/intervallist/IntervalListFromVCFTestManual.vcf 2 2 2 . NA12892 HET C T
- Issue Tracker: http://github.com/lindenb/jvarkit/issues
- Source Code: http://github.com/lindenb/jvarkit
The project is licensed under the MIT license.
Should you cite findavariation ? https://github.com/mr-c/shouldacite/blob/master/should-I-cite-this-software.md
The current reference is:
http://dx.doi.org/10.6084/m9.figshare.1425030
Lindenbaum, Pierre (2015): JVarkit: java-based utilities for Bioinformatics. figshare. http://dx.doi.org/10.6084/m9.figshare.1425030