From b4bb3fd33e2f5f56c10895363b05fc2945c381d3 Mon Sep 17 00:00:00 2001
From: Sander Bollen <a.h.b.bollen@lumc.nl>
Date: Thu, 22 Oct 2015 16:43:19 +0200
Subject: [PATCH] fix for #222

---
 .../src/test/resources/VEP_oneline.vcf        |   8 +--
 .../src/test/resources/VEP_oneline.vcf.gz     | Bin 6141 -> 3550 bytes
 .../src/test/resources/VEP_oneline.vcf.gz.tbi | Bin 101 -> 101 bytes
 .../biopet-tools/src/test/resources/chrQ.vcf  |   8 +--
 .../src/test/resources/chrQ.vcf.gz            | Bin 5359 -> 2786 bytes
 .../src/test/resources/chrQ.vcf.gz.tbi        | Bin 95 -> 95 bytes
 .../biopet-tools/src/test/resources/chrQ2.vcf |   8 +--
 .../src/test/resources/chrQ2.vcf.gz           | Bin 5357 -> 2785 bytes
 .../src/test/resources/chrQ2.vcf.gz.tbi       | Bin 95 -> 95 bytes
 .../src/test/resources/unvep_online.vcf.gz    | Bin 5727 -> 3132 bytes
 .../test/resources/unvep_online.vcf.gz.tbi    | Bin 101 -> 101 bytes
 .../src/test/resources/unvepped.vcf           |   8 +--
 .../src/test/resources/unvepped.vcf.gz        | Bin 5094 -> 2514 bytes
 .../src/test/resources/unvepped.vcf.gz.tbi    | Bin 95 -> 95 bytes
 .../biopet/tools/BastyGenerateFastaTest.scala |  18 +++---
 .../sasc/biopet/tools/VcfFilterTest.scala     |  54 +++++++++---------
 16 files changed, 40 insertions(+), 64 deletions(-)

diff --git a/public/biopet-tools/src/test/resources/VEP_oneline.vcf b/public/biopet-tools/src/test/resources/VEP_oneline.vcf
index 6d575651f..e5cd64bcb 100644
--- a/public/biopet-tools/src/test/resources/VEP_oneline.vcf
+++ b/public/biopet-tools/src/test/resources/VEP_oneline.vcf
@@ -1,11 +1,5 @@
 ##fileformat=VCFv4.1
 ##reference=file:///data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta
-##UnifiedGenotyperCommandLine=<ID=ApplyRecalibration,Version=3.1-1-g07a4bf8,Date="Sat Jun 14 16:58:07 CEST 2014",Epoch=1402757887567,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=1000 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false input=[(RodBinding name=input source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/Child_7006504.ug.chrom_merged.vcf)] recal_file=(RodBinding name=recal_file source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.snp.recal) tranches_file=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.snp.tranches out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub ts_filter_level=99.0 lodCutoff=null ignore_filter=null excludeFiltered=false mode=SNP filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
-##UnifiedGenotyperCommandLine=<ID=ApplyRecalibration,Version=3.1-1-g07a4bf8,Date="Sat Jun 14 17:01:08 CEST 2014",Epoch=1402758068552,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=1000 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false input=[(RodBinding name=input source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.snp.recalibrated.vcf)] recal_file=(RodBinding name=recal_file source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.indel.recal) tranches_file=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.indel.tranches out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub ts_filter_level=99.0 lodCutoff=null ignore_filter=null excludeFiltered=false mode=INDEL filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
-##UnifiedGenotyperCommandLine=<ID=UnifiedGenotyper,CommandLineOptions="analysis_type=UnifiedGenotyper input_file=[/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Child_7006504.ready.bam, /data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Mother_7006508.ready.bam, /data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Father_7006506.ready.bam] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=[chrM] excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=250 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false genotype_likelihoods_model=BOTH pcr_error_rate=1.0E-4 computeSLOD=false annotateNDA=false pair_hmm_implementation=LOGLESS_CACHING min_base_quality_score=17 max_deletion_fraction=0.05 allSitePLs=false min_indel_count_for_genotyping=5 min_indel_fraction_per_sample=0.25 indelGapContinuationPenalty=10 indelGapOpenPenalty=45 indelHaplotypeSize=80 indelDebug=false ignoreSNPAlleles=false allReadsSP=false ignoreLaneInfo=false reference_sample_calls=(RodBinding name= source=UNBOUND) reference_sample_name=null sample_ploidy=2 min_quality_score=1 max_quality_score=40 site_quality_prior=20 min_power_threshold_for_calling=0.95 min_reference_depth=100 exclude_filtered_reference_sites=false output_mode=EMIT_VARIANTS_ONLY heterozygosity=0.001 indel_heterozygosity=1.25E-4 genotyping_mode=DISCOVERY standard_min_confidence_threshold_for_calling=20.0 standard_min_confidence_threshold_for_emitting=20.0 alleles=(RodBinding name= source=UNBOUND) max_alternate_alleles=6 input_prior=[] contamination_fraction_to_filter=0.0 contamination_fraction_per_sample_file=null p_nonref_model=EXACT_INDEPENDENT exactcallslog=null dbsnp=(RodBinding name=dbsnp source=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/dbsnp_137.hg19_nohap.vcf) comp=[] out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub onlyEmitSamples=[] debug_file=null metrics_file=null annotation=[] excludeAnnotation=[] filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Sat Jun 14 15:20:24 CEST 2014",Epoch=1402752024377,Version=3.1-1-g07a4bf8>
-##HaplotypeCallerCommandLine=<ID=ApplyRecalibration,Version=3.1-1-g07a4bf8,Date="Sat Jun 14 22:28:02 CEST 2014",Epoch=1402777682364,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=1000 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false input=[(RodBinding name=input source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/Child_7006504.hc.chrom_merged.vcf)] recal_file=(RodBinding name=recal_file source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.snp.recal) tranches_file=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.snp.tranches out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub ts_filter_level=99.0 lodCutoff=null ignore_filter=null excludeFiltered=false mode=SNP filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
-##HaplotypeCallerCommandLine=<ID=ApplyRecalibration,Version=3.1-1-g07a4bf8,Date="Sat Jun 14 22:31:13 CEST 2014",Epoch=1402777873043,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=1000 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false input=[(RodBinding name=input source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.snp.recalibrated.vcf)] recal_file=(RodBinding name=recal_file source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.indel.recal) tranches_file=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.indel.tranches out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub ts_filter_level=99.0 lodCutoff=null ignore_filter=null excludeFiltered=false mode=INDEL filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
-##HaplotypeCallerCommandLine=<ID=HaplotypeCaller,CommandLineOptions="analysis_type=HaplotypeCaller input_file=[/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Child_7006504.ready.bam, /data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Mother_7006508.ready.bam, /data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Father_7006506.ready.bam] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=[chrM] excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=250 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub likelihoodCalculationEngine=PairHMM heterogeneousKmerSizeResolution=COMBO_MIN graphOutput=null bamOutput=null bam_compression=null disable_bam_indexing=null generate_md5=null simplifyBAM=null bamWriterType=CALLED_HAPLOTYPES dbsnp=(RodBinding name=dbsnp source=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/dbsnp_137.hg19_nohap.vcf) dontTrimActiveRegions=false maxDiscARExtension=25 maxGGAARExtension=300 paddingAroundIndels=150 paddingAroundSNPs=20 comp=[] annotation=[ClippingRankSumTest, DepthPerSampleHC] excludeAnnotation=[SpanningDeletions, TandemRepeatAnnotator] heterozygosity=0.001 indel_heterozygosity=1.25E-4 genotyping_mode=DISCOVERY standard_min_confidence_threshold_for_calling=20.0 standard_min_confidence_threshold_for_emitting=20.0 alleles=(RodBinding name= source=UNBOUND) max_alternate_alleles=6 input_prior=[] contamination_fraction_to_filter=0.0 contamination_fraction_per_sample_file=null p_nonref_model=EXACT_INDEPENDENT exactcallslog=null kmerSize=[10, 25] dontIncreaseKmerSizesForCycles=false numPruningSamples=1 recoverDanglingHeads=false dontRecoverDanglingTails=false consensus=false emitRefConfidence=NONE GVCFGQBands=[5, 20, 60] indelSizeToEliminateInRefModel=10 min_base_quality_score=10 minPruning=2 gcpHMM=10 includeUmappedReads=false useAllelesTrigger=false useFilteredReadsForAnnotations=false phredScaledGlobalReadMismappingRate=45 maxNumHaplotypesInPopulation=128 mergeVariantsViaLD=false pair_hmm_implementation=LOGLESS_CACHING keepRG=null justDetermineActiveRegions=false dontGenotype=false errorCorrectKmers=false debug=false debugGraphTransformations=false dontUseSoftClippedBases=false captureAssemblyFailureBAM=false allowCyclesInKmerGraphToGeneratePaths=false noFpga=false errorCorrectReads=false kmerLengthForReadErrorCorrection=25 minObservationsForKmerToBeSolid=20 pcr_indel_model=CONSERVATIVE activityProfileOut=null activeRegionOut=null activeRegionIn=null activeRegionExtension=null forceActive=false activeRegionMaxSize=null bandPassSigma=null min_mapping_quality_score=20 filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Sat Jun 14 15:26:18 CEST 2014",Epoch=1402752378803,Version=3.1-1-g07a4bf8>
 ##INFO=<ID=DN,Number=1,Type=Integer,Description="inDbSNP">
 ##INFO=<ID=DT,Number=0,Type=Flag,Description="in1000Genomes">
 ##INFO=<ID=DA,Number=1,Type=String,Description="allelesDBSNP">
@@ -164,5 +158,5 @@
 ##contig=<ID=chrUn_gl000248,length=39786>
 ##contig=<ID=chrUn_gl000249,length=38502>
 ##INFO=<ID=CSQ,Number=.,Type=String,Description="Consequence type as predicted by VEP. Format: Allele|Gene|Feature|Feature_type|Consequence|cDNA_position|CDS_position|Protein_position|Amino_acids|Codons|Existing_variation|AA_MAF|EA_MAF|ALLELE_NUM|DISTANCE|STRAND|CLIN_SIG|SYMBOL|SYMBOL_SOURCE|GMAF|HGVSc|HGVSp|AFR_MAF|AMR_MAF|ASN_MAF|EUR_MAF|PUBMED">
-#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	Child_7006504	Father_7006506	Mother_7006508
+#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	Sample_101	Sample_102	Sample_103
 chr1	871042	rs199537431	C	CA	1541.12	PASS	FG=intron;FD=unknown;GM=NM_152486.2;GL=SAMD11;CP=0.000;CG=-1.630;CN=2294,3274,30362,112930;DSP=107;AC=2;AF=0.333;AN=6;BaseQRankSum=4.068;DB;DP=124;FS=1.322;MLEAC=2;MLEAF=0.333;MQ=60.0;MQ0=0;MQRankSum=-0.197;QD=19.03;RPA=1,2;RU=A;ReadPosRankSum=-0.424;STR;VQSLOD=0.079;culprit=FS;GATKCaller=UG,HC;CSQ=A|ENSESTG00000013623|ENSESTT00000034081|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||||A:0.0078|ENSESTT00000034081.1:c.306-110_306-109insA||||||,A|CCDS2.2|CCDS2.2|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||||A:0.0078|CCDS2.2:c.306-110_306-109insA||||||,A|ENSESTG00000013623|ENSESTT00000034116|Transcript|upstream_gene_variant||||||rs199537431|||1|3610|1||||A:0.0078|||||||,A|ENSESTG00000013623|ENSESTT00000034091|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||||A:0.0078|ENSESTT00000034091.1:c.306-110_306-109insA||||||,A|ENSESTG00000013623|ENSESTT00000034102|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||||A:0.0078|ENSESTT00000034102.1:c.29-110_29-109insA||||||,A|148398|XM_005244723.1|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|XM_005244723.1:c.306-110_306-109insA||||||,A|148398|XM_005244724.1|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|XM_005244724.1:c.306-110_306-109insA||||||,A|148398|XM_005244725.1|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|XM_005244725.1:c.306-110_306-109insA||||||,A|148398|NM_152486.2|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|NM_152486.2:c.306-110_306-109insA||||||,A|148398|XM_005244727.1|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|XM_005244727.1:c.306-110_306-109insA||||||,A|148398|XM_005244726.1|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|XM_005244726.1:c.306-110_306-109insA||||||	GT:AD:DP:GQ:PL	0/1:24,21:45:99:838,0,889	0/1:17,19:36:99:744,0,603	0/0:42,0:43:99:0,126,1717
diff --git a/public/biopet-tools/src/test/resources/VEP_oneline.vcf.gz b/public/biopet-tools/src/test/resources/VEP_oneline.vcf.gz
index 3dccfacc66deb4dfc034ce927addceeb9a716535..8f0e4cb3292cbb6ac9918f7f13633e21b9f2ed1a 100644
GIT binary patch
literal 3550
zcmV<44I%O$iwFb&00000{{{d;LjnN74ZT|HbD~NX|LpxMy6)8O&eSG`=1sir7UL_a
ztHy}QeSEP9wn`0x7jQD!@{ixsAhc<Fkh$|nrGkX>>ht)W(~WQ5Ed78iqj=+{?Z;04
zquRtjy?GOpC5cJskv4sNeRXxUa8viHd;6$e{W`cRA17BUH~lbw+J%dN%;lzbwO(Os
z9!6_-d$sctugUJ4OE*bf5Pm!E-?e|a?Y6t)i}7yrMB+BSn7wRC`!-C;io_RPl6bMd
zP5mfr|Llj|CucnQ`9H5%X9cU2v-Sgb#q(iFk_IG<HYDM=ZBDV1#(v20yFow#l62pn
zDDMvn_U7l=m%GrTp9YT($KCw`Zh`i}s0P35L2wfI$ng$q@UBQmhJKR1X0`239laYr
zjBL+ea2}sbYA|o(C?$S4i4s3OF}>5RBf{&BZ887krj7+s@G@BwUiG9#x+Go{lgS!#
zyc31o>N|rP0z6=IsE|Q#I^lg^hkxyFwsw}Zoc3{zm3PLEueomO=(Lxy@44ZPyB)dP
zeoX%2Wb|st=<VWYODcqP>KKyvpO?@boCOPr8%2r~p7L1nSi^e@icpVk-~&vaWZtw+
z!Nvr9kdEBs1Bp5ISq=8o4HslH&4h5#x6dx#3;Ynm7}!Y?dA?u;U5Bp#AyO{@G+gKk
zv&h*7zQ?8Ez^Q{hu^Xmqcas&uYyQ!+2LHzQ;%GXyxgav^vMcP*iwAkpE=-XhqE(R*
z7YI5eBCfYa_Hc&s&Tv6V(ltSeyCEow7bHe*f&y>?S#aw6ocewcxv5Y&&E6ff)|ao=
zS*y9uqSOsgm=U6Ac~Fy3H6H!EU`}`~{=J)!`>7j#aCRF`kzjFuUS{x!mT|N}zxja|
zeTqT^vHpOP-3Fy3NvU=_NYP``L{qXvPoMz(wS$^U{gaw}e>SJHibNwqsgkw7<pk8q
zXL}1B>5~tm5kWC=7idAY>AHmym?*M;BxvKlptbvvpeI5?1jAAagVBPop0n2e23*Ju
zU#Jrf;oZ_OETZQy|GfAwCzJ%qdB68#Hyb0U+-YVoG^f;sat~L?*}7r&0b~(M#d*0m
z?$+^=lSW&#J@njS0X|U1AKic(eVp?qPuAc-NEeqT-}Vw@><HQ=)Xu)2tV#T@1d&g$
ztUPPd2%?e_xJ^>Xrv-ZQU2YxF7J@8=Gf+8i0qx(kpr8oLc=(V4&wjc_Zi+S$)%zTw
zgsfr$12DC1T`1v{--b^yAqybsL}a-<ubnXf5CMnp`2^qysvMS^<>DeQBWhPsoPYoy
zNZ_xd2tvFi@nxPg=ow}i5VP63;8f(UAoA39;Jk6;#mnVyFm8c4EAZ$Rd?<ka(Xdw!
zJV!263e+;Q<X9+r3xeKdkVdt;M@XBXx3CBl5<0QEu^+rJgo9Pn&vxhC*TNxPM-3GI
zJF8e0kX>UT{jXoh{rMc4jvNG&(3U3AeFmW~yb!5=ci<$UKS><rXu=Eq;(pLS-w!?>
z`auT3n}aEz=?mlDz@FVc_U5yxeLJ2zx3k_^=$)<o1bu{X!Qg?O!RfmMN(Er$Y3B!v
zEGS`C<N)%3FgbK4B(v`cOggag`<**b{_QQvVOC(tie{?VlzSe<1@U-1Xa(q7`U|?V
zzbc~r^aA|>+8#HQPHb)(6DH8O94r93pr7SqD$0jlv<rQhvy#ItSsazdp|sxmp%?5j
zFf7j;-f@;Z)Yu-(PW8T}14hCNgCv9P<M)_rUgdGwE>F9aGn<~xlp<W<I~OJEC{AB{
z_O-;vduMo8BMFAl3PlTOATb5IEWZi$ybIX|J-Eosk1$Z}lEbs58^D;h3n^gH1kLi2
zPL1v^p<6Pv@js<I@X#h&5N;aqcEL9G(>md!U>gFs<|zEU_8~={519HK#o%RVdh`^F
z+H`=@fzX7D9NKeehFBz-y!`gK_KHysMLT{&75EtleNV7s1Nj>IH1eoQuq8`bZdYZ`
ze(lg|iEHeRC~>*UY`%|6kfFE7Uf}S|jnlhj2PO$<*^3&gep?d|79ZRU>tAGUjbeZ0
zhoI^O!Az2^+by1ZY&7)_ie5un|86?6XZh}iKScZdhR*l$xjiZSHb4oYqoNBQKicrn
zM|tZ#>bplezh={b-sby-Y=#`&w`ENOUHR=(*4Q8okV~QDka+KS^#i`~3)L5=7_a*N
zv6NYjdMkwR%#-r?a8L#M^fWo}|3?A+)ii%Rq|#WFTT43YEXXoLjzb_GJAIVE;N`;r
z%*N9G?QqtcX0{BYXSzfl!k!Pa!3$3hC^(%#(RpifJMQ*|merJ`G=hsIacC3oASww7
zvtCPkV4CxxfP5rDW_l(fM-;jx0O4pDEw0feZmI@W(JyWE3H^dDrKY0DY7JW~umwvJ
z*!1sjWewG@HziB?BK6<U7n$<X6}e|BY=I5?*9C>{_yuODO%uy{VTaU|W&OMvUlk9L
zA=Cwmr5b9}l1=Rk^nWuyFzsbba@X+2%%bSVU#+v(C$MP}9aFJA57cW$NQaRPZFA`w
zdam}5&j=kv<KZ`xc?)yLr*V{K4RFDoBNQW+Pk({#aHv9M(rsFPyd{))n;Zr#erHal
zWPLIlyp9{?D-m{_QrPm{ptFXaFP`0&S$_lLU}bK}fIP|((G#W<xd}zv7(mna&WQ9&
zLTEv&0z7Gd^;h}mwT?#@0SQ;>x{Y<sz^Bh~@my9dS(9{GKqE661zVP4D5o?^smC%U
z83?6Q8kNyttf{i0pVDY0jVW1JlXZa{y`(WP2$FQ;l*TA&Gz%-Ht_y57OBz*IbQR1#
zrLh<dmazg35cp6ki4;v!p!@^|VSGd+NxFswtH&}UQcMk-mL?!#zQ!QMG7VEe#GEY|
zifU=9z>%1FTaq;eE0(|mxOBHADVl0%0#D-7-6qDCCS#GZ(%ptB8HTKitS;RRQIk!}
z6zGymZ|e$Jt;qsia_MXpkk1mjQ7(P0VMS6kT@a7oSr%!sA!(v){DYYcl_3;^=c@pK
z4$^*34=({LO1JB(BJyX!Aae}hf7KFbPzwfy*I-(PX;m<&Jc9+*AYd>G2JJ)x1WI6%
zSup4vgJOb;3VMU5H?SdypH<MC9K9t;W);1~FMyViz@P(zEJ-ITHlPwJ82B(NIu>R-
zE*SVQs}h!kWmYioVOAh5L>fv4KFlUoghg61oP^m>HM3F$A7>Ssy9h1Kf`JdTD(fQr
zDcyP!XjvDfHe)!6G@yzox8))bm=87dC{tk5Ls&2{RM&N3?K4t_?WS&O0#`CphVZ&7
zD*}rdDZ_ciFmR<*hW3W6se%M%q(?1QmMuYXFjC`)6swiCvWBZ;ENYUBl=WT1!lKkJ
zNm)MtsW5v~k>VpfSF9=n<X6$s1$oR!S+q>(BUOHt)uUzUmY~ELX&Eg@GAtUo7%7Vu
zOcs^CV!@J4MFe_`lm$zMUM#>0M#_R^!LU_9sg$u&72T?Gs=@|S1ID*1N|~Xiu2qF!
zDFdZymRhM@WimC<$j2zlFe$29=}<*yz-TFzlm@fY6b-kGvJ8>}b9q(mDP@dcu#~H`
ztJ1NZRnb)0lm*~Y<$0CGNYUg<nJh+%QW+yvJzBe5wRV-o2tc}0W*H;JgdIv1rNLyX
zqKyNiJc^M~naL`Pk*VcC#*U+%`x?iw^qqFTyN1m)Jsxxuw2jHa_vriNrx*0tn>0~B
z`yS{u$~WZu0SU>zPh5I>UR>vtdw(QuPrq_EFJ5TuJ6(s}O!6=Ev47a~>w$BZem@Xo
zEI`V>_sLJ_8}xa$ozAK3`N;0?d-=s44tvAieEcxlcW<4UJ?`}O&TMLryZg>Cd$%^&
zJAaJc-wh9!xpVh01)>4vz8O3^UiNppxBJsv;^=U7#<`-0{C@KAe$?yIy>jPfdN*oJ
z?wrPLw=wPY8}@M4xPP#RjeJAipzd$vC!-A~J6D)v3A2wf`>6axM~{YSU`dr5aRN}N
z0m3S{(P?z-24oX%V!1J~9jDPBv|;ZaM`5eqZSTSlVe}lf2BY?PG>1A;O}#0%2E(>v
zkGdGQI+Hd%Ad*_0LHiOnVb{~@jN7y=UnsHxe<VehFEGF~Jn1?Uy5novownSv`)ztM
zskH2ITW=kIm#3{ZCEaXw-?zFzFRQJ-)5c9jmRs330owVsU?1JLbua<0Qk(uQgj`BZ
z7}s0(T>$>3q_n0JyNxg8*7Tumw_d+wqqHjMfsnTHH&t{|hSfT}9%}cU7JJFjei&Tb
zbXt(~ZF}DvJ3VJMkg^{PYgG>an&tmeRLR79daRaL!an!y{BUsd&*dSdNf3prJoU1l
zqZq**&VKiM`<go2Fa<iA_}XhKl75L{PRy>7<%fx#i@UJ*9ViRADVNu8Zqq^V88gm1
z3}gM+xZO5MW8!Y+^gF?YjbejfD<g46eSu!d`T>Em&IhW7hotWizEd5UrLvV7N-vz9
zux(T<bN~BjE=d3ksv#@QZ{@)pU5>1Kt>ugdPc2h_h-K>EuuS_QmT7e@J3{&IW!903
zFZP!4L%e1DjX3H*#4^3wvc_O`ZFjG`lk36#^<>zPuJAQ%>tuYbYS)%^ZNf}0U6`ho
zJ%QaZwyqVOJ~3col`eEifhW>6jKuIqp%0`B3|lyCVB@F%1K<hkHK{WI03VA81ONa4
Y009360763o02=@U0000000000031BX7ytkO

literal 6141
zcmV<Z7Xs)XiwFb&00000{{{d;LjnNc7wueoZ{tRi|J(T#M1s4^;5;(lr0nG#XuYh9
zUL(rR&ISQNku8Zk6seHZ*b4mNud4YJMa#<06Id-4Ba2<#Rb9_tHCy@Vr-kd0MG&s7
zSbJ!7KZ*_Rho625$%2H$w~2<s*H>3pjul&1?f!#w^-J$6AC9h;R{X)7Z+*ujCf|^*
zR!gpF`oYTDTy5>hZWP}e3oD8(Sp4317cOyn#1G=#hJ>wPy|#R3;QB=STfeR8n~k@d
z65I0Jd1%FM;9owFFoI7_XmD@2w@X&B#Q8$KY+EtWel)C@`M<5ta3aIW*OGe8Dom?m
z%ov{K#2+s^o4{UaoXGNuq^PPQ$;xGk<#>Y}BJD@Zx4d2CMkWeXBLw66n{8~OLbTt0
zXQEZ`*xh>Gt+gJwQLHU2FCt7xEC<xF0R5Pe`zO(SAYnGEz$fM^Si`!%8VuFh2eL~B
zV{566^kgP19+ME4cYXNdBdnPaktf^RI;6k<G$S%25}V-`1iK&g$D`u;nuJT&Uuybb
zP(0jNj)M<bCbuIj@W+<g`?DFE@dLk2faTitK|iiN1%(G|gAgaxmE%U%+#_u_UXc(5
zN*=GRr|H&@T~O9ruJjGuCXTdWe{?_V7>pA<`jNHXcrY=s*zJ4M(1(*jrx+t<9GHvH
zveABt%J#?X;3EmGCDAyRWtqA4hc@nZ@v~=b#v9^K2slN&(?qtxG7I;~T)PpdBXEdg
z{;>t$h<9dW2O)`4;VsV#9?eY<6Ob9c>|UL*jmS83!Ah3*j*h$ekV6*M){Aegh!jJh
z6Nniav0#@$uuS%K`_7mqe*qWZ{U@R9BCv9*XVdi^;-ya$24wBVW{`5<8`*$>;4QF7
zBCP3=2qI%;`I^iuy<iT~8aOfu-5>;^KtJGHu0un`jN;I>V;|;dMt{&5&A`>xplfpz
z2FuV|n<0V8PR8v6A)6P5Qd@x&T5Dpit<44vPJ@u{b_>>Aq<?1{wYFed-0hQTLs(IX
z!1+usus`7J#})BSj5que201T(yM`^|K%{Z>o4wg)PfXCdiNTRRNGaSPpTaWh0DLqE
zXY@g5v2bk{Y`-flKZgyX@Fr*yA}Nhy8vDyI*lx_ZXW1W24>L-VI*8*gLqbzd_O}t!
za49``$W=eXP%viGeruLGW1dOOM~F0Y>r>5uJAy2!WJ}^h7{tU_Sr8w$t#11KqZPW~
zOsOtul(pZF^kKhc8lBEBmGO-Q=EI0bAN7JIW*+dHkHpjZqwW|V(Rj>lw<6vq<CG@&
z<QvVj+HXHk1LqdjhIJS;zeb-kQLqgUGSvq#E4b-8yb8f3u0TV9sj#dpu_9zf-(=an
zXIGmQL=}J4TDhK+J#1{34SN*^YnnHS)A(pF-uw>Mfhv%QeIQXW_F2UW%Z10>E-P|I
z{-!}m-Y{_pX<!8`sC-|GdB&9qz)nH9Y|O&|!wK?yyp2f%KK*3_OGd*D;M;Z{H69XY
zX$9bqPw}sz3!Y-Y2ows_UOSM;F-4gxhy`$=PhB*Ezalg6z1=4-HMbX^yha?QwLJ}p
zW)l*o7dS0QRf|Ph&6hqT)C>cZ%YK(W$_wWjs)=TdCQQ0ylAmaFk4)QLT473<jzb9o
zcxSIL9_L_UnheX?KBdF42x@o?`8&wj{rEp$7_%wYEO*VSFJU&7l~qaN&zS9u+0K~l
zjM>hZ?cW`<oeGX9$#;z6Xl#Ew9G9>^Be`#m<ocs_XYd_ix$5`{$g4VK4|jbQXNJTV
zRU8;qwQIn@CG+`M3<E?q>5%$mcyz5254qs+6#?pR(B}?+KSR_n3{iRMj8M-A^^8!z
z3xv8%s}tzq1M%Ed05&qQ0(;u+cy`BZ?9e1(7(fxmGYW9N&RdbOgEf>gVhqOZEE;`S
z2Pj5uJ@aeJ4b9aW?1pZ&hA5?4(+1<-pko+jOK;uvM?D5|R#LS0Vpii6%ts~&La0~e
zsRi<bg!(oBeM|<KL*Xh()Y$>#FcYHIO!tNclu88IO4Q8~r-wy(iJ`-gju79t{+5cH
zr6rEdhxf8oI_u8bcvO!DmLfGJZj<>oOSK97F-8;J^N3f@VpDW!W3op#uzV;$iy$v{
zdGSbs6hJH_kORz~A@%#w?f8DwesegHl2MmPAAsDhv(tF0po$5o1@@nbEXAJr(+y2Z
zJWHo<f=5s(1}W5R03~Pyl$2J$Ml&%&A(BHj@d{5YvXYxtWeAv(+=z5jH2@q(1VsU2
zXV{;a5BjvPk7kBB9u3}ulf&fT&)qVBNjr2DmP-s_R)#n*H~MkmA<3fczR?;#bf)ha
zsH(nWg^r0L+kwAu9V$-!=6Ebxr%k~^5a)BO)EQp4B$~oP?f4L7W;RbwL1khFNK?Q!
zwm<-DFKE(afC&S+zCcvOLn(7?0BHDNjx;HB{;s!XCZh962mXO_L!Mwd^+zwOkj@<X
z&BHy>$MwzRY2b%0U~+;|MnKrf(u9QDKR0Zrt$~7cz!ME>?g&-qV6rM19_kAqi>UM@
z&GDFof6HOAzW4C!G;B5J<<d3JUi0Eh=jA-hi-MxOaQ==BV4lHS7+at7ET8AEdG(s*
zUyw~vWR(|W@qCtlZj{bv`R7jQ+$nvZozl0uf3dQ^%lj8A`y1ZBsAYaWwg2X)_Fwn@
z#TSKa0(Z>`FG4m|5m-?;L$)(yJ43cJWIIE)e|N}sDmMDgu3sDo`gN~gyukhp<i0tO
z``XtpUV(Cz;S+FIWybmL#rf{V8K0i-UYz0S8J>PO@bnCCPFth6ngYe#-g?vnI{p&T
z!NhXIyWuc>_m371w$U$85AY(!lth8IrAe_h9^Q^k$f^uP-)1$YHwhB^fP+-Nn)o^y
z(m_<D(=_IAkh%eSCLPI~%3M2Ada(j8yts?qtv<|GOm;)PwgpA1(>Cw)$zVKtKfxUT
z)nAu8pyOHSu63vwA3;q^3Ze44)OI6VpLU)gy;BwN5<czq^zuW2Wtp_Z>JZbu)5pXM
zDMhM0H%5~PZ-`{?=*zE>Tb@f(9Td$E#&$g;QGCg?@m|OTOh}XBU90w8v$277U}ihL
zAQW9PGe|aMJtZ4RhAC+f{(im|^1WOO`H+UC_8Z4uGQ9LV^~1i8)hi-d&_>-LZ0+pg
ziWC%-Nr?G;ncY_4uyaO8)3*F2YWj{=rW8)Vn$!K^ndRmsA2ta{CEAw3!IUgoxnZ=?
zc+_Ef53O$R-7T07u&D&ng2XcWd*U}J<}B!VE>#wwRajt1^}w+&-e;pD=|(i3S=t-O
zTgkmL>fHBO^@%erZ3RG6dRq*Fb-5&A@jPqVC^=}X2#p-ufPXm#Sb#V^ie!*<h^<;0
zaWTQ`Yv4em?K*FlqP{-~HffI4I9_GwYo|0fM-Q$w$gcdIdchB2&$O4gz(2Q9oS|h>
zmwnKB*^T+kLU{eK6@*~(7@efRjm7Oc`gadA`z*BlDEV)9RF=Juh!HGensSJPa5MJ_
zYZGrnqDK)~&%Iq2+y%a3{w~4Y#3%bcijZ&$dTH*TK>w6`YtY>+ts2qGhQ&x25Puo3
zz;rm!DG_EF*Y(G9M5oC{fDk#)f?L=Hz#z=1csVmkWJwgZ#v`LMeb8t9hYo|$`w{$a
z5(XHy<E*|}C9T&!?E8n$inLFMK)-fszDz|W;-U3q<Xzz$LTF+|k>M`aR(i7tA|!*j
z`$599-}DBNe9gW1i(6h$RFxHu+#o{qM^7u-qe}>>IfN^BNlQMqvk<+PZ4%j`OAjox
zA6>sahtl}t{?u7El})C0J!@H=hc6j0M%E;%%+@Q669~OOKdlk%+XLp^UN*h)^7O^l
zr!LxiFe<a#e_?hF_ulY?`Ic=bu&gq7aKc<fYIumTKULRTr|<=z4RqUeDjJ_mPM8cW
zUrZqPBgi-qz12Qt3A;Vg>*r7IPLUw+b`xj^E379cxC;%;$;yflm{)e+>YcCv7F^mx
z_d3%_RrjaNUqO9I%3y{2=!BiOMi2Y5?oQ#UFQBove5i^;Yts$MAC)yaC#=!gh5;ba
zmltWBB1q)Eyh6M8YO(;kW!1ErD;`eJo`V^x=(_h;n4S`_Aw}auEBZjf%JkU@)8i7P
zuUKDyb@Mi0KQN;Qbv>}%njN%_Y6R3~YH0pqJcXJC1~hhd#SMGLDY*wwNxZVwG!XXZ
z52q*0U%Pe~Oh<afh(+$wkI0_{4?}+#0cbP-%?GqhP{<i#*(*jL%+jvo5=r6}J(D6#
z5IQ7etcbz08sb#gcPs3>USP$w%+Y)C<aPJSJWndFGXTRLb`C(u1dFVpLN?fw5T{VW
z{_^w`GVMkyjA}ObkGC{EGK(-+Grzi?9Xtg-1GfIaMB6nJLpzJybuVTfq6PzP0)v2r
z`6D~$toY;Yt3FL5k{J@DidOEXvcO6HtZ$$qJ-L9NlB0Zv{t)8^ddqe4AcC%7$E>W6
zgqaiKGuh8Fs?qbRw7<eE==XQ1gw(v)z?Dv?o&T*YO1$WGy?4{z6xXG=n!!^e%f;S`
z$xa5C{AggG1vIFqXvv6!4YNsWZnlB!l(7PRjMqh;U^@oxg8dk?c?)Z7v5`YLyKb~1
z;eSO*&zyKo9Dr3~=mwYwUxzx(+_ftED1LK<KKQx>Sido#AP94NxZnZ^vKauv>j3#a
z224bjAz4~@905*{HHQ1l`~Ey61P>-!0a+|w=Uz;JL`2d<;~1kHwo84{Nx}%_%7U|W
z^cUg#fVYX$FnF|0&yy@rDq}4zJXxeO@2t?-z5Ny77KpP1jc!1PHK~6%=$w{3p_y^l
z5X*%m<wcpFVWBrPr6G5#lxZX8XV?Ve6GW`gK6|?&!@;icZ@u;AQ_~^t<Af~yuk2!x
zrR)<n(*5*}ynEe;rbB}nCXFMLyeisN7hzAY4waRde~>UpvPpH(zg!Qxuh)Z*X+59`
zIPHdCstcn|kM{KDY^wK1rqQ2uUQNBTl^bPmE|^D9`Zj`4fwVH;x}HOg;?p39Bo7%T
zNu3d)^36e{bZWKV83XgrFG-SSYb05sr`C(|uA0Sb?2#T@F~#ZLRc7t^E;(nefYOOx
zL1VQJG|GVmzOjWEpfHt5z%H}(UBFpU`jt4P&`3k;=j6-}H<};izo|%ae?ogOJLLN_
zJ@Ke622KX;<CnlSdw!gk%fo7A%%-n;N@gz5orPBq!gyb^`zt=Y!`~mjzL`O=WC8~Y
z2)%Vg{hOfXE%3H&{`QZ)o7+a|SbEOtW5PlQrGAO1(b`7Xk}+%d36BivZJ7!hu(#eO
zbmLRZw<4vtAyCy^aFNF}r#=QDXgPfejxwdGv{x&_Ws=&H)QsopbiL$jIrcIrhc6#1
z4b}dum@ks-3P8SsI?dS0Q7py86?cJz_e*-&=k@LgiSu3B$%UUR$>7(r7ApD73ghvj
z1xNzlRGq-;*A>BQ9I2Y|6+TsKcDn!!0GNdC+sel60-D+>qy3cD9Z!e)Ea`5lhiET1
zbgrG?_9(C0kV+VYiWX?R)bQX?QhG~ux5V>Png-_Qq+g&Ia`;Z;B?)*XU+1*2fg3PZ
zNOH$}Q)S<)7XDiH*(t_e-7l{)J0ae(MXHl0`SD@zSklLr3ZDBXCiFctUk0gYi@dc&
ztmBXc%^ayA9$H-{0^sEW04`u@w?CM5rc{<e@Q7Qirm)8p8|<p|kOil+<iaGqw1;>9
z*f@ZXCR@oi;TEidFTdR6%|Ld+lcBd4sOZVPKbldnC4j}5LEv07Z@Gr3a3b?Jjd^1J
z#=K=4g20O>$kGH^I2PZz{`cpS2KmbkwkiB2_McG~lzFa-L^DTZfrkCljDj71jTmA>
z<#;)hgKhA<{JI#QWRE&QPz6F$RK!M;SEawe|Fh|#;GWwgQB75tSrFX0%N2cn0!>rh
zF;&s?K)g~y1S4xKb9gI^k^bX_MSH<$@Y&YHn-ou@Af^S-shlHZ5X<%H%jl={PzBe-
zHZ3{cs$H=kr2wnind2(7PEz1iE|hyq=wXa$%eA}~{vfd)yDeya1L0t2Rz88DKVb~;
zgy~@2h~D;ur0LqPEa~QiP=c0O^2GY;E)(dr3Wt{-eX*)>vZQc_<6Jh*i%ni)Wxj@n
zFK7g=*%TDvkVeS)G*y-dLiUhGEND1R5_v^Fq>*wOm2GkoFW0Oi=QIij3$e0tNTcL5
zQj-%@S+0>;&1pngkVO#tkfvGCa6Bh~0%~-~=0t)d2@rlY0&!f42yYg1HM{5df=Ey$
zPHjpxM1`(7Sf!~dY7J4LY*rD(rX<!Vk}I^$@{+&_%^D4Gxw@OIAc=}pqe(7Tx5{x%
ziRbE=<?2>cR#AAdPU>9UU^QNCsx`d$T-&k$QcHXdFFsec2$`>0t42Q8wZsXmD9JVU
z_<P|+60fjQoj1NOM25(yU4}E49wa(&do!fhp$oOkqEM&LYyw}I0Qwi3H5|ljf>7n4
zHWjsbWP(_o(1d8HnV@77qyr9Ms2Y*fY=T^wAgI9N2)$CJS2(4{e$9+tt<X1FRy{`F
ztZo29s79c40?)DsJ60ePj!dYUS&+F}Z|AZJRWpk$$JT~fHlb=}0o<aFLq4HuW|b3a
zgEXITU}i;>)T1n_b{3(y*JVpJn@}~h$jf#4Cs*r%p?SH^wTlS{mWHfS=i7W{2(D^s
z9{+1JhZFAoY(kOMWw|!?3(_LDtFkK9sInj}GQ2GELXE@)X_4mzMd6OFT4Z~Lm&6(e
zE=Wrymgk!_!BLPZB~ng2Dr-@2Wsa*Wk_Bl|ca<hr=h``GQ4e5Mz#hj)xl%qCn#Tm-
z_pBw?_;Ep6SWAUEa!jwhSlX7{tO@agG`AKw8CM6n3ev(_fGm#csxTI>3Ux`ZAT5l=
zLoKe!6$NQwtR{f1Ba}jJD^ZY}$CN4*z*GVFc8roQvZgE@Gry1<N|c)7QSM@~rc?*{
z3d-C}f+!wUs2~@Ku_+v-R0=uOI&fQ1=0*|#mmiBgA-53#Oa2&lkuA%4%$j0>vLLw}
z)4W*NNRarWYZf*VgrhbR#nRsSV|y118$psjx@K-8L4^+G7^PCIDb_U(1!c)b%27`i
z3md6YlE{kVXye@p$FcZMJN-!?Lyrfo2!GV&xHi5|p6{54&ZNO~>2=_1CTYl@@!IZl
zmsoguo_*4v{XCbJd&VnwX7)nkxz#p`FO%d#U-3Z4>w%_)_XA-K2UdCRJh}8oE+%cK
z6Dr*t>fPr~@_|1C8Fb9i{qVWnH)i^%)p<5%Q+?EaZVl+WwcfMwet0_`q#vd+zMle7
z4`<)?9t@lQyLr~T(`3bA`e}?3M)%40$^GrH)5c!8bvGRkFD7H-qTjxlcDfh(V0Q8D
zULRZ}4fzGC{~|dVz1aJ4>&5<WL@)M#mihytkPB7eSdqU7Bgl#pBwK;IXkE1Q3vev1
z!SNRp-7qe?Jq_CUFz|1>ZEfp+@PkMHrZ?0^LlZ(uROJSL(;H}pK5TQ`O>07bufg85
zdid)*Xn<~7BMmF`rNArj4=c$0B?lQBM%u;%o4*^qrSUg<SHnX};YJ^6@=f_#k0v%)
zS-olB-n4<97jL?T#x(?<zoC}`wB$3JK76OiAOd`{8vdIt@|JA?a=&@kh8oadg`4R_
z*SJgmW_qvbH~a6~kX8gfVBMSKr4?>UY2Ku7inOkAQ@roc?t7PatsC%vO@Ho;jE*tu
zvGji&>{Un~&5}og$g12kz4K0k;yKZ+nV#YN=OXoP;syRP@k08)l8xXSNB{r)tY4$D
z6}5&(gS)mH0xQ4e05Z}~w(0tjo-B8%Ker%4_y(VUes-DCg<lclb%k-9To!J-iQ<r0
zYy7VTlT5~XfuO5lD?<GRzHIX=7^?ZYp-#{c`vuH*$b(ACH>sfbQIQimNTI1d|9xn(
zEF=a|;f2QM(qNKZO499fc}0VVl8IkMGVx!KO!_L4NvBFy%JN@ItP+dAtS#lMXiNDQ
z>?nT~$>igbUG!$xdi%OPx$eEYo(wM7EAARPJf6E2rR!$%S_OE{UaD%7jzFW#HLnF3
zM-%|B?4`^KFv4B~I);A)9AGaw=;k<uQ-1h=U|4pQcX<E+ABzYC000000RIL6LPG)o
P8vp|U0000000000ht}Bs

diff --git a/public/biopet-tools/src/test/resources/VEP_oneline.vcf.gz.tbi b/public/biopet-tools/src/test/resources/VEP_oneline.vcf.gz.tbi
index 1e25b0dccdd45daf908188094e8a998f875094df..1d1bc2ce9351b5f9d267a75fd6f7df3fc34326df 100644
GIT binary patch
delta 30
mcmYdIouDta@>W1(oAVxpOFfAqJfaNU)(i_=CLEmTssI4DMhjK|

delta 30
lcmYdIouDsvtTiCAt^JY0rJh6)9#MvFYlbWa?&lL-6#%s|3a<bF

diff --git a/public/biopet-tools/src/test/resources/chrQ.vcf b/public/biopet-tools/src/test/resources/chrQ.vcf
index 4b97e4867..2ed17f44b 100644
--- a/public/biopet-tools/src/test/resources/chrQ.vcf
+++ b/public/biopet-tools/src/test/resources/chrQ.vcf
@@ -1,11 +1,5 @@
 ##fileformat=VCFv4.1
 ##reference=file:///data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta
-##UnifiedGenotyperCommandLine=<ID=ApplyRecalibration,Version=3.1-1-g07a4bf8,Date="Sat Jun 14 16:58:07 CEST 2014",Epoch=1402757887567,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=1000 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false input=[(RodBinding name=input source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/Child_7006504.ug.chrom_merged.vcf)] recal_file=(RodBinding name=recal_file source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.snp.recal) tranches_file=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.snp.tranches out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub ts_filter_level=99.0 lodCutoff=null ignore_filter=null excludeFiltered=false mode=SNP filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
-##UnifiedGenotyperCommandLine=<ID=ApplyRecalibration,Version=3.1-1-g07a4bf8,Date="Sat Jun 14 17:01:08 CEST 2014",Epoch=1402758068552,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=1000 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false input=[(RodBinding name=input source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.snp.recalibrated.vcf)] recal_file=(RodBinding name=recal_file source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.indel.recal) tranches_file=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.indel.tranches out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub ts_filter_level=99.0 lodCutoff=null ignore_filter=null excludeFiltered=false mode=INDEL filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
-##UnifiedGenotyperCommandLine=<ID=UnifiedGenotyper,CommandLineOptions="analysis_type=UnifiedGenotyper input_file=[/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Child_7006504.ready.bam, /data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Mother_7006508.ready.bam, /data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Father_7006506.ready.bam] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=[chrM] excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=250 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false genotype_likelihoods_model=BOTH pcr_error_rate=1.0E-4 computeSLOD=false annotateNDA=false pair_hmm_implementation=LOGLESS_CACHING min_base_quality_score=17 max_deletion_fraction=0.05 allSitePLs=false min_indel_count_for_genotyping=5 min_indel_fraction_per_sample=0.25 indelGapContinuationPenalty=10 indelGapOpenPenalty=45 indelHaplotypeSize=80 indelDebug=false ignoreSNPAlleles=false allReadsSP=false ignoreLaneInfo=false reference_sample_calls=(RodBinding name= source=UNBOUND) reference_sample_name=null sample_ploidy=2 min_quality_score=1 max_quality_score=40 site_quality_prior=20 min_power_threshold_for_calling=0.95 min_reference_depth=100 exclude_filtered_reference_sites=false output_mode=EMIT_VARIANTS_ONLY heterozygosity=0.001 indel_heterozygosity=1.25E-4 genotyping_mode=DISCOVERY standard_min_confidence_threshold_for_calling=20.0 standard_min_confidence_threshold_for_emitting=20.0 alleles=(RodBinding name= source=UNBOUND) max_alternate_alleles=6 input_prior=[] contamination_fraction_to_filter=0.0 contamination_fraction_per_sample_file=null p_nonref_model=EXACT_INDEPENDENT exactcallslog=null dbsnp=(RodBinding name=dbsnp source=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/dbsnp_137.hg19_nohap.vcf) comp=[] out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub onlyEmitSamples=[] debug_file=null metrics_file=null annotation=[] excludeAnnotation=[] filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Sat Jun 14 15:20:24 CEST 2014",Epoch=1402752024377,Version=3.1-1-g07a4bf8>
-##HaplotypeCallerCommandLine=<ID=ApplyRecalibration,Version=3.1-1-g07a4bf8,Date="Sat Jun 14 22:28:02 CEST 2014",Epoch=1402777682364,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=1000 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false input=[(RodBinding name=input source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/Child_7006504.hc.chrom_merged.vcf)] recal_file=(RodBinding name=recal_file source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.snp.recal) tranches_file=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.snp.tranches out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub ts_filter_level=99.0 lodCutoff=null ignore_filter=null excludeFiltered=false mode=SNP filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
-##HaplotypeCallerCommandLine=<ID=ApplyRecalibration,Version=3.1-1-g07a4bf8,Date="Sat Jun 14 22:31:13 CEST 2014",Epoch=1402777873043,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=1000 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false input=[(RodBinding name=input source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.snp.recalibrated.vcf)] recal_file=(RodBinding name=recal_file source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.indel.recal) tranches_file=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.indel.tranches out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub ts_filter_level=99.0 lodCutoff=null ignore_filter=null excludeFiltered=false mode=INDEL filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
-##HaplotypeCallerCommandLine=<ID=HaplotypeCaller,CommandLineOptions="analysis_type=HaplotypeCaller input_file=[/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Child_7006504.ready.bam, /data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Mother_7006508.ready.bam, /data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Father_7006506.ready.bam] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=[chrM] excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=250 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub likelihoodCalculationEngine=PairHMM heterogeneousKmerSizeResolution=COMBO_MIN graphOutput=null bamOutput=null bam_compression=null disable_bam_indexing=null generate_md5=null simplifyBAM=null bamWriterType=CALLED_HAPLOTYPES dbsnp=(RodBinding name=dbsnp source=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/dbsnp_137.hg19_nohap.vcf) dontTrimActiveRegions=false maxDiscARExtension=25 maxGGAARExtension=300 paddingAroundIndels=150 paddingAroundSNPs=20 comp=[] annotation=[ClippingRankSumTest, DepthPerSampleHC] excludeAnnotation=[SpanningDeletions, TandemRepeatAnnotator] heterozygosity=0.001 indel_heterozygosity=1.25E-4 genotyping_mode=DISCOVERY standard_min_confidence_threshold_for_calling=20.0 standard_min_confidence_threshold_for_emitting=20.0 alleles=(RodBinding name= source=UNBOUND) max_alternate_alleles=6 input_prior=[] contamination_fraction_to_filter=0.0 contamination_fraction_per_sample_file=null p_nonref_model=EXACT_INDEPENDENT exactcallslog=null kmerSize=[10, 25] dontIncreaseKmerSizesForCycles=false numPruningSamples=1 recoverDanglingHeads=false dontRecoverDanglingTails=false consensus=false emitRefConfidence=NONE GVCFGQBands=[5, 20, 60] indelSizeToEliminateInRefModel=10 min_base_quality_score=10 minPruning=2 gcpHMM=10 includeUmappedReads=false useAllelesTrigger=false useFilteredReadsForAnnotations=false phredScaledGlobalReadMismappingRate=45 maxNumHaplotypesInPopulation=128 mergeVariantsViaLD=false pair_hmm_implementation=LOGLESS_CACHING keepRG=null justDetermineActiveRegions=false dontGenotype=false errorCorrectKmers=false debug=false debugGraphTransformations=false dontUseSoftClippedBases=false captureAssemblyFailureBAM=false allowCyclesInKmerGraphToGeneratePaths=false noFpga=false errorCorrectReads=false kmerLengthForReadErrorCorrection=25 minObservationsForKmerToBeSolid=20 pcr_indel_model=CONSERVATIVE activityProfileOut=null activeRegionOut=null activeRegionIn=null activeRegionExtension=null forceActive=false activeRegionMaxSize=null bandPassSigma=null min_mapping_quality_score=20 filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Sat Jun 14 15:26:18 CEST 2014",Epoch=1402752378803,Version=3.1-1-g07a4bf8>
 ##INFO=<ID=DN,Number=1,Type=Integer,Description="inDbSNP">
 ##INFO=<ID=DT,Number=0,Type=Flag,Description="in1000Genomes">
 ##INFO=<ID=DA,Number=1,Type=String,Description="allelesDBSNP">
@@ -81,5 +75,5 @@
 ##FILTER=<ID=IsIntergenic,Description="Variant found in intergenic region">
 ##contig=<ID=chrQ,length=16571>
 ##INFO=<ID=CSQ,Number=.,Type=String,Description="Consequence type as predicted by VEP. Format: Allele|Gene|Feature|Feature_type|Consequence|cDNA_position|CDS_position|Protein_position|Amino_acids|Codons|Existing_variation|AA_MAF|EA_MAF|ALLELE_NUM|DISTANCE|STRAND|CLIN_SIG|SYMBOL|SYMBOL_SOURCE|GMAF|HGVSc|HGVSp|AFR_MAF|AMR_MAF|ASN_MAF|EUR_MAF|PUBMED">
-#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	Child_7006504	Father_7006506	Mother_7006508
+#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	Sample_101	Sample_102	Sample_103
 chrQ	1042	rs199537431	C	CA	1541.12	PASS	FG=intron;FD=unknown;GM=NM_152486.2;GL=SAMD11;CP=0.000;CG=-1.630;CN=2294,3274,30362,112930;DSP=107;AC=2;AF=0.333;AN=6;BaseQRankSum=4.068;DB;DP=124;FS=1.322;MLEAC=2;MLEAF=0.333;MQ=60.0;MQ0=0;MQRankSum=-0.197;QD=19.03;RPA=1,2;RU=A;ReadPosRankSum=-0.424;STR;VQSLOD=0.079;culprit=FS;GATKCaller=UG,HC;CSQ=A|ENSESTG00000013623|ENSESTT00000034081|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||||A:0.0078|ENSESTT00000034081.1:c.306-110_306-109insA||||||,A|CCDS2.2|CCDS2.2|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||||A:0.0078|CCDS2.2:c.306-110_306-109insA||||||,A|ENSESTG00000013623|ENSESTT00000034116|Transcript|upstream_gene_variant||||||rs199537431|||1|3610|1||||A:0.0078|||||||,A|ENSESTG00000013623|ENSESTT00000034091|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||||A:0.0078|ENSESTT00000034091.1:c.306-110_306-109insA||||||,A|ENSESTG00000013623|ENSESTT00000034102|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||||A:0.0078|ENSESTT00000034102.1:c.29-110_29-109insA||||||,A|148398|XM_005244723.1|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|XM_005244723.1:c.306-110_306-109insA||||||,A|148398|XM_005244724.1|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|XM_005244724.1:c.306-110_306-109insA||||||,A|148398|XM_005244725.1|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|XM_005244725.1:c.306-110_306-109insA||||||,A|148398|NM_152486.2|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|NM_152486.2:c.306-110_306-109insA||||||,A|148398|XM_005244727.1|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|XM_005244727.1:c.306-110_306-109insA||||||,A|148398|XM_005244726.1|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|XM_005244726.1:c.306-110_306-109insA||||||	GT:AD:DP:GQ:PL	0/1:24,21:45:99:838,0,889	0/1:17,19:36:99:744,0,603	0/0:42,0:43:99:0,126,1717
diff --git a/public/biopet-tools/src/test/resources/chrQ.vcf.gz b/public/biopet-tools/src/test/resources/chrQ.vcf.gz
index f4d4d695f057d29e4744ede187c8c1b155929366..4068353b8e36ee701149416a96adcd97dd23706b 100644
GIT binary patch
literal 2786
zcmV<83LW(yiwFb&00000{{{d;LjnNB3cXotbK*J@{mlFd)t%binOZ_@W8T-h<zR5A
z3xr_8?A|YwVhgN_v3>o>V~ZcZE!n&jW^iUQlS*M5tGkc7Tc=y%<HsfS$TEo5IML?a
z{-auvKYaWcktK<UZxfBZZeLzrIyk|X#$aB%{KLH5y^b$eIC)syr@rHng;J?qx+}T0
z@B<fzm#H1wmHl&NiQ@!=^1--&qy0WGG-F&Ir|WwXX>xh?5)y6TCuBvUvO!`yqG3V<
zU;B;v#=SM3{PwRS(b-m1%0&AfUWszDBuOUmgEff-X<aaEB@y)ndF*+_Be8My%6Q-0
zidViazD#|a-I;SskT%{ReZsmqDj<(-2u=*Qg6yz>?27ngNaN&KRM)#j?5?RF=r(l(
zk548Agu^IE2=ymHOp{m8yG9WewlUUE%TKP0APBsdiAzM&lLF}y+d)JoF64MO@P*lT
z%>oK+5VM#Pvp1cHzAqy0(siiwq!p}>3(UMbo*zqH7qO`?BWh#+8iyks_9OC_prcnn
zM=y<nkbFZ)w+N7!eqDoMo+S&28)J$D-ill?FW^0gG1MG;6ln5Q=1q$%*oc4+k`ay{
zNF<2Q3Wz7zcgT9ml?c(-&raU<s1ITE^f(S|dSV8KC9Z&zOzi+^V5kbSzzRKT3u$Ot
zMYt!ze&XUaUkJzYqiF&8HMOH)I@X0C?zv0<2L5dEAUjC?1W_NYwi!`IP$3b)wu|)P
z3}uy}%!VW`K`~wv6hsb*5RQ=tP9Tn8zAu>XdjU>P8fV#?LTP>ZVVi}D>nupHhkQN|
z1<Tw_%*@#A-q}0hSp5}_$?X*T4_3Ms44s(VpBJ1S(J~6w=nM7i;5qOS#QFin=^7;@
zPMCGwOVB*7peb3RdoY0hN?}bU^wpZYI-ApcMWPX5qhgnaf`UT%tcOsMo+-3O1Vscp
z$RT0kc2Eo*1^*&OYy5&-{7BF}AwGg;DS^hw5m(Pyv|m69VgH3WArEiJ!r%l?es=Ht
zTTn{8_`Kiy5$0_K(>vvWp*STNmb<?~R)~Fm1G)$sC3v|vHi~%3N`erDdClEUfWpY=
z5qm=GW1WXQalwI*E|?`>cmX<6g3_4b?9<pK(LZ8Do?%*9)+7N$B_Z$_Cy-ALx~Eu}
z2b2Ymr!Xg#^9We|W(x`y;VvE&Qs9XuF2V^~2h8qMfMT+W2sFS9o4Yk47$5le5g`r?
z=>}xEJa2t307L{5)TRXJ2v#{vm*=9BEhC0&yEuUX9*9TXAb=2uB)Z6w20cM918UYG
z2B#vtg2*%MAbE`==jGxHv|9jY1s)B-hbO3iH0%{4&k=@|0&97X<WMR41d2X$N@LdX
z0n^6l6HEg6gmtVqqTb7%;b7M6UhjVXQ95MvD8RzsGmGVxvI|V4|Kk(6J)c9<k%gcW
z3Rx1}ateJBl}N>_gP;lhNunTQ6H)1#>p}l~J$TIP0Vm*f4&^I#VcawI*<jvV%%=Ka
zys!qd-dXCMxim(P5H4su&=WX4jbW((t=y;7b9hjEzQ}>(fnhT0j0wkg50G+Ealdl|
z#@`=FMzbf7tXQU=hH@*i_=$MTbE`moOX{$l{pB{=_b*T%pzH}v=|mU0v6Bfj-Zd6L
zT~N=mHWg*fE=qk0eO8=566auP<V))l^=&WZWVkzX_$(l~E3iG7y|(*_H5jp|43Z4C
zk3XZYIm+W*xqMx%tl9KzrfkCnzQZVXgD5%n?6JoAtu?$UkOaeEg#rf(NW{pF=Qm-V
z$B=DMgSVOa2n|&l=kJ!-gElSo8DX&m<$1}vMx4fMONQ3;IZ44o>%bvGH(;k;7}2Ci
z_!!yxK&}}JKe-fA<SFOWryv3^L(yZWSPWB6rRGTp<&4@hYKB<ETwnHhxc0WK912o#
zL*-FULO+x2_6@lUbsE{sBw*nwE9|QD*&kiDTA~8GBQ|(9OWxl{JCb3Khh7l)3`fb$
zvJ0IAlx(MfsV^=8Wl`>Cn7_l_8bowOeK7TcpeNz;#@qMyzBRRTL&ucXznPBoS+=_o
z57EB9p`&)zx5vA>4OD_ytBAql2OA#tDl5GQb@#yLk8B#yr)<B#d&tqPrqpU+D|@`>
zg$>dGVMZl+;(Zp)n_}TVX?}Z(aa8vYwaf~vx2+OUc(OY_G`~gq>@?Y<|6_!H)Xg6@
zsVo+|ttIPq9J1uhkq2Vl?V}hPFA5FtzAfz!hO^$3V;KfdY>Fa<J>|W@i#Q%wa5{rU
zXS-xDHhRNWt0GBB0GE~|uuYUgRAS=uTFYu+lF1>DJQ9yX=PEKrVOs)Fj)sBLMi+8L
zZOAJ6T|>|4cXT0Ds#QfTfTaakWJv-{|9mer%znKhwW@D2|Bk-k#&@pBJo61KuwnnP
zps*eP2^eanDJ%6Y4ymFj_46=(C?2PPFc(x?YD2BGlxFP>_TSA9d+WQHWUdjx%!1&W
zu3Y}}32d599aHW1Jg{E*K-P?`S(!^*sJXSjeMPAmjEC>0a{ztE^EgO&0d$0Ogl&sO
z**}cS4^@~>woS{9w@zvt#Cd}y?#!7MUMG3ubyz5mROnH{Xp3U6>q5;x-McM$eFJbX
zGrSW(Hj@b03DYaRF=N{ZNR!%URMOo*C_$?&d9wOSS6S=jMz>{;_^ZU#<a(_kA9py`
zZGj5Uj!CkBL3~)n0TK^^?WxU<C}7(#?@cPG&wpBMqb$Z7s4QgDhZYqYvh8)jq&ElZ
zHtbiN#kNP@bPa3&Fv)%++JB+5Uz`@0{p18X9MH1qJwyElybB(oY?Qtj>HSSFyXeDV
zZ`fOm??xMAV9oS#x3_`5TOS*n?vQWS%#HQm(bdf`zbveqyD1DZne??ew`~3yZuI^%
z(>Thn);KeCmpxDJu0}nBxxafoy&08YUsxI#rD?BUg57`V_D&y`vKW_``%77|Em=HO
z7qTSpZ<YOR^#coFNtRTl6vc9@RjW4Cs$A-px_U{jsd7bDN)z3(N_|s<WQYR4(>Ju#
zfAE7RzhjQH@n|8}6t!8eC>?XCS^CJ3<xY2^F~*TPT~h<Tu2<o1tSL%MEmxHWd`Q*0
zQkG?<1#b*%qRCRDqjxo>qxUsd>Z%=mtkpY*hr^m$k?PHkan&(kyrOpcmL^xKiqhdn
zOImi_ibHEu2MBPLH1@ew0)<p=H9EJ3CPRs+c3}Sr4UE#6-f4R0c-O|ps$d5K+R3&)
zY*LL@C*N;qeXFzYnzcK#eBJFp%4_<jH@14#%#`?_40EmKU(K?ws;bnKH;kvVb+E~N
zyU3Z~@5?-)i5K{*Eb;i?L5$!@=Kq_G-e%5jG*8&5$Zfk)mFgF=EG_s|YEeJdGj(Ns
z(}hK$RFvKI-7)2g&%ikEFj=l2f=k0Vflhe6U_UZ%F`fogt;>=C>J4_K)-MQD>wKUJ
zcu4vQ;d^a^BUM@)D7y$aWwlvtH8=kmEhGtuL2W42%6oY*W0wPT$5zgG@HH~^mmpKW
z16l2tAgdKccEIvK3)X>&H+!q`OT5+ijyTqT39|aPk(JC@TQ}Oqq;1}|C&Q9-DYq51
ztjKM()^4@h&1$nOm7C2Le<L@_a;sgfvo{S@g*SDn3U8#gs+8fQ%3eriX!px<LvDQd
oKgbb)>G~xA03VA81ONa4009360763o02=@U00000000000J;%G{Qv*}

literal 5359
zcmV<L6cFnliwFb&00000{{{d;LjnNO6zyDlbK1z3|J(B^x-zx9NiD&Wc;nG-$>L#l
zco|{m=1!$lLTW(QLapc#n<;<zJ*Opsknun!j+3fOYRaH)_36{+`8(as{QB$McIce>
zD<hO|)y|i4QT*lCUwt~KK6OnhV|VlX{M<4^<Gj_ot(^bfJ>PT(=L;je(`S)uIaHU5
zmGk97Z0Igq8s2$i24->lT$~$0XaMmW*Ph$d>Qa}5YmfRWTdfS&>f0`rf9thm#q*r?
zgqnt9&wL}anOnG}egL0xsVKe|-!FukQJ&4~g_aRg`Hg0T<o_a<h-D&Ho0WP~s1dcT
zO^GCk<u`@4$IPWHmIbL+snzSXO0`zVOgZv!hCqH}xQ4S1>_A7R_Ao(g*NZ|O1tS0P
zi3Ce_--#UO(pdHFAe85Z6Hww)!vb;4K|Xq5|4C&RULoF+xl~`W70~^l5iYYkx{f<T
zV<8Wfcq9-HsgLAs7ruM}n(<5Y!E_>v_I6v-m-K*!x*q{!H^bg&xNToie_^`|S?Tw;
zjh<my*eH;V9)ZB;i0k{s6ir;_wkS-wvR#mmZB9VoY^70ZC9<;Yz?eC-Wrs`ZBSUfX
z%6OPWZfJwB-iw(hFq<%t2fg9Vw5<_~-MfLY@*Ef#YwYr)t|^0YzrF3FddT#-Z<wgR
zSY*3>Ci_BtV?kw65Ck$aKFgy{2S2-pH}a@EreGA|S}qGkfz0g(ePstA4rWnH{~UpB
zglj!8nNNd+dBbtoz3#D)0?+WUJ~$!@=*YK0ONMjJ`|VAaMdwE3gqKD@w_P9;@EPv0
zpqD;d#Ou1e)+X^cm<8tklk({T&~hSY-F7YNB(3oTbY+J+OJ?8c@qj_F79bKct2;CR
z&sZ9+TqO&K&461Cd&a)ad|(RX1GZ&b+*R}-^ldYAVT`Qx`t9KqOl<|S);*ssd}F2i
z6b4)AY<HA;&vYfa0wXk5R9_jMhYIIT$X6SIHs{IrbfHEBs>Rhl=q9)o=Ln2XcbIvH
zqwkm0)zRPZ!DpMf+-L<$#Ew7~`IG5IX^Re0*U>qWhU5*KZH5SB#lS|HKjjVW`P?>b
z(EU0?eg+Ff=5>%Hcv9lW-1iqgi#&bi80MYsU_^;S2es^lPkG2mziT0PmrTp;M)Xs3
z1#LR%snf(5H<3iY2T#)@mrDlB5qL>BTTmC=Af(#T0ROmDJ4yQ&!?(eh5?OLDYkeFl
zgPy8u?e_1ve$N2)p~s_*I&6WF2khnxb>!Z#Gr}hBj~m@BspF+HB>_Hu;*nPV<L!i5
zmp~in5LCa++eyG8|4F2}1GR#guKaT!OyV3Q#PpgVR4YOmBBQGd!cN<{w*;?}&ef&u
zSZQN1S`^Kt&sIEcQmgpIoWJ`7+JP((v;BmlZQqNs<!2L*9xZZw2Ci4+ue>9n4`E;l
zQ*i!XYw~nf!az^VUleCPLwABWA4VZ9!l%0^V#+Am44$G{P`r(eMP*=*58(&j220VP
z2NH#Dt}F=T7^3tg_yU;Fw-OEDTc9)7S${LBnK}Pv8ex#+_QW9?4G5SHvs4IG^LdiZ
z7cK<U6a!?_eup=066XrCiL4FBBq15cC+^)t-Lw~mpS+m#Lka?TXD-nnXP{yp471Tb
zc@O&n$l)QxZ?>84&HsF7%vNg(VpFI;huP{vwO*-6C(L%jY$wci!fYqZ_V13_js-{j
zm7f^Jao_%NI4++4gyeoWlIsmy?fy@M<?{VUAg}z89o+RI#teZk$Qdw-YF&hm0(rS9
z1`Lr+(pCRHb99W%9M$a{uLw~8fI4^Z=>$=~GengtCxm)Js3(N_Qy|nulAS;fchs?$
z40NPp26p7j(e#RVrms`qXOM((j{=OZ{k}|0wt`efwf?A;dZP<;;EiERN$uLOeSNtC
zy`fpHz)QK*<o>AJZ)=*aD(Y2l*d-9N;-tOfv$9yjc%%a(lv~v%v_O2Qpu9CeALD*%
zP)H^AI+H;R)4^-?WNoOxO6G->5;wEh=z*wIA-o&H5#k%$jX1k$SQ3Rz_fE7*MqL@6
z!{w-9Dyk<hx9BWNLv0LywBcBB9O`7F*aS^l8}Gc=H(W?SbGAu#o8%EYDS%iAAWtxR
ziqvn0m!q3u>)q3V{55V8Ndxe0TWeY30?Mg?YhbsnEO6|(Y4vzek_0~7WA`9ZbW+IK
z07_5^$SKc&#YU`zElU>l!X@rlq$xMa%HS{=z7gprVgNXf2#N#5_MkV_Z<R?;8BR5Q
zH0*x_BZtB4&-H@ApfwtbAjXQ&b6p~+8|`>&Au-WbPg6&??a4<1S=F@+-_ns~leu%-
z;_U2SoP<g1xFMJb!p#^XF^0pML{%6l9T&VzPsddgP#LQM!W7IK8o+?D<1|Syz<`0X
zKJ(JnA2N}{1EAr8I+CE!{)eJYbwuamHvEBfLmyx`w?`+<kk$<9&8KVP&HIbV<G>Fm
zpo^ti76D-;%L5WF|5UIYM*|sYgC%NQ-2sZu!eEtAJmeQZ7D1*ZiSZbO|JcA}WvB6V
z3bsA_<&~x+G^O%$`{j}#l}okSGyCsY0B$0<ioW%dd-;;ol<G}EdWN@Jty-5#)$(aC
z|5PZQ_VQ1a(y3DVIV+_fb^l^%{*?DGmgWz<f3bJ+)2{sw-?e|;`xoC8vX#W9Sb7$+
z)oUf8TslFv6J$F<wi9GKLAHN)$aX9?`pK?eJYn>8uU|Yf{RznZa3J^E*Dnr0xm@=V
zxGOi}bob(P_u_<4Pj@d)@bm;ve;Rmtf;Y#}*j`P6WNt<dw}7_0Ky)xRZ2xL7NY4II
z!z>DZhkSq+F(x!%PQ-(vIvQM#bcm`1yw6*X_)UUXKVTraM;(uoAshtTaGLlWc5*X7
z%j7+&QRyqIl3c973orJ3eW?sK6yw!!t5qRMwOjg?GVYJ2AIBKuU;Vh;0vS(zd!<0e
z_yS^Da0s=@OD#Jvl}Y;n!aEm%RKeD6SIIV(1c4+eRso-OtsVwm2q~3ZyEYsLcta#T
zqt6~AtB%b>9VE>=Em}=!5Ee)a?}dy(g*+%;se8|wH4o^($X0SeC@7FA1RJ`VP!ED(
z@+$K`ovwxaEZ0Kr61SB95QPGfDxbI=_FT+f0ZpAY=rCVho7*c=kWj`x#`8sbTS3Ii
z86i!}a2F`)E1sDWH~};#yWLa6-lTk3Bp{U_%7TLlovRzg$ivaFO}e*gr+a+~ssn6V
z0d9fgs_-ec8)S3J+K$bI1!xrr47eOb;o19ayeC<RERlujLEMV(m2u;~!K_cMNoFYk
znv&aM;H--U^|$S5(Z*i~iEUkDgKfaREDbb3tu9A0`04;lwZ!8>ir3e`fQHd(Q!WKP
zcg(yb#>%2pC;ZqciOs>SZS>PCf5%?%L)bIv#wPHuC<s%uO!q|})LwdHKGhIjKUA3y
zIuFrEwzzS7yN-YFVq~BCh8x8H?~WYVn}BL;9`cYwErgpJn=rgE@~IL8bTxC<9WWPo
z#Q2?oyRl97Tx21hiFK3MKZg2gW38;?EsQ;^XBCT{(5LPqT!QMbqn&v%jkvZunjt!k
z7XmNh%u{v=s{j~;5fv|I#(^yMLUlCM+LK#l+PiHN^xiLEe`BAa+m6!wW@Mzkx3T9w
zZQF)@-UaeC6ZNGc%Dg-<9<)tWI0F|N8$qDii<ObwECLTnA?~h|NWu@jK~!yu&;H_8
zD%I+Bq4XCwh!Fkp-HO(*04_BHcNGgf<zqPu-dkwV!1QgtVIjY<-PR0J<D1=~({!j1
z5A8U{B0mlf8PG>oG{}usaudf8dUt%1BU+bFChv69;l<~N&m)(cX!lmjjc$F#=smM{
z2S<!IOf!aMxv~8t#v)R~O^n^4ilQDfFW4+oOxwyyd^|p4FqC{Dh1icE;|c3(>lhMd
zYpCoWKfXGK1Lmy9P!8s%9v@*Y+%d;XBS2uDTfN#nf&eC5UPE`=lW|`5$BbV>eu>jy
zZua33D_4iNyQ8j-nNyiVVQ;vQ6$ggb@#*IrjrI{V+L6xyi9Sz6J%*9Mex5?BdoWml
z-qLKEj}^B^NYBCuwXM2$DNK%WSf8Wuff3wMKR0}O#PBc!>2umw4ldpV><4OeAg?pi
z-m`+1miK_Yks6A>5O<-bOoPJC%$Z?VJ0|u3DhZdyiaWyY_`&3e@hjW(*<`5Xl(>yu
z%3tswI}hPMi~zK``OOEEOpwS4HOwVZ`qQN9DBvsMlJ7}T!hDPR#0UuP)le%peJ3}4
z$6-decXG5|+<9Gp^PEQ|*C~Ku2P+47k+6AMP~kLKlTd3bgkAD{7c!|va}(`Z+&|y*
z@JQx9Tagd8W3mV4643P>38EDVp`69p6(=OOL6Jb2KoF3S&*?sA&K@rh+BEk_GNAZY
zu(Z7#furnM@gO5T*npqnt$f0Nh;ae6<tlCnpek6CrSXN58Ko{sf0mK=o`c+efKgEI
zuTcnzda;1Z>`q(%n<Gk{;IP~~No|UBIj$zSYh>71Tak2UfNZt~cBg;>^%yB>A@hir
z<lJ-t=`Ldq`e=t)9%DHM=7RMY@!kWC2n#vnvttKK>i<^|*UYiigbcJ2LN!1FJPx(U
z%r<iC*#71SZSXJ$n7`4XzzH|{u)zcnWD@|wD~9vlGZN5+PZtJmM}QHeh2iexJ$L3)
ziW?Iuqx1RU*|!5A5Yc4QI7Bap<#JoJVmCs*Qs>NU{RMa)@IE#gf?L~sKgj^05@TWD
z&LSUqW%$<m{Re<sz|I0B>VXXRg#JOleO&N_V#Z#9FK;=?5+!ee&^sQ|aCRdTX#?^W
z7C~H!h!x6bXT6PZuxk9SsPDct9O62Th{Ato74tM?AF+_mw=d-SunkQH8X3feBZ-fS
z)_Eer@!26qiTsoLERH65qF<d4I*0SYmn0wX0Gw3A&*g<-yUS~OeL7KkLtX1l+XqAM
zbZH0a$pyIwp+^C@3WSweWIGmDip!lGf;>c+ICKV-^EU&QlA-x}X9SbKJSTCO-NVTO
zEp<OD*ZC;EM<2<y71{3You}HKt>b;>94H+tIW)G{fQB1jfhV@$101H380?bBwE<@Z
z$s@Hgu93Lb+j!3p7rNQXf0yIr=7{oO`gGp6e8(eC42%rQ$M1n_cI<eQE}v#AZ8|w<
zDXF?Zb_QNO@WWln?ozzH#@`<wUQC}Y2(uu8@LNaRz9~xH0Begjr+@rxE((&qWS`Z=
zfQ1Ii?GjO=5d~P1k(K>`TZa6$Ob!j0k>mMxc#Qg%G;<mPS<MC$xldy1J@Y}z$ssuM
zl!VgmUKR>*XpciP?x%Bl@nhEaQYeRqkGX<s_bBGO1iLLDUqYTHX1o>4F>%gZAmII;
zT=qF!9pTHHRhsdIpESteV^#{ecxL$FXs!a1z>}&Yruwj?c#R_wGalh{v8J~RKmmYB
z__~d>-`qk|Gnr^Nq;*D<fijJ&oBSr)^97x4#<)G$<ZTEg1fik<63-Mo*cGSVOy14R
z`7K2Qc^lUY977JSWvNnux#GtePi$ZY#PI3H@ZRO8@8%Q#-s#g_jGernr7}Ij-%=v^
z!Q;*LVfU|sJ{~GK_MaHgchGzmq@pf1r6pn=i_UrENEC6ac1Qrg%LV|vg{7Tdf7+gK
zUi$1FmzWP>_X##w=ii4YIGM&5Ci$g3y!$7F3_cn{F537JXax_y?D%9Lz2Hgs?FG(y
zeD9A(6f6mVIO#L1N#2Xaa!o9g-(>PYek1RNVyPsRkKm;Nyodsxx&HS{af9<$i$bIH
zRnxzqEcoObQN)t@3tpgL|2Cyy#lMG*a<MK-)szpRC`r}BY<!bG_A!DYC^gEpa<L)R
zD_=4Hi^*X-`$i|RX!6WV8N0F<OMZL;Mbo}*s%^~!{>ooOFtWllhquB=<v*Sw+GWH3
zi_yfJ6c0ld@&st*_7PHuWxM<^`XSj=L7G^m#oJqZDfWT{VC5@wq{8zg2VU7kxkEzn
zLkwHC;i&k7#QoT9&hr~E2P-o+0|@^KBf}l0Cv*dT+Y^GOZ5|+LZ(f8Hv`B*|=2v?W
zL$9TOU2yopsw`G3HE|c=XxB$%c-$sQZiy0}aUh6$5bbRfw<w@&xNVP%q{B~`n<Uop
zBc235b^xLRAWJ{_O|r)fxko(3(Nm24sJ682!#F-8+BPV7;-njR<ODP<pz_#$u=%|y
zoogt5OVI~P=dm4s;LW0bTOZyG9$P(astnckBj9dj*m_j^yj<%(Y99xeqki(CYonV9
zywt_fSKV99<bS<Kr89{s4w6r87*BK)KaX!N2kjP$U%i@)24_$gp7mO1lXm9}s{ga=
z8>N31>-Y@Ce-<a(vz`0WXS>&q&vtJ*|AM+a6NR#L<_BV<Q7P5RCGku>Q<XEZQWlG%
zbT(Er?X1(4!7hB}UUXV=<ledL-o5A!<l#UUD^j^$ElL;NzN{&OmMC7RV}1chxKO)z
zowHhkZ$nv<8s$Pss=*(jRFw*%C^ev`rHy4#s9h+kEL|ua8M9vLLK(`{i|jtJTrLXL
z`bF#Vq6P0u<%^Cci^Y;8UGQBdIsQzC14^m_3-Bq(_%|g2IaO@bF0Na$2-%`^0i`Hl
z8R=qjBP$oX1suLxhIv4@7jYGYi>fs)lF~!&Xct=mE#GttSLy{=zpOmAhgw^kb_M=V
zgteBErfJ+%DhqY-kzb<VF7X)4R!@S$f6f!zrVeu#u?_Nn86Cls$p3qMRGKL4T7AzP
zMX_lXOG5R%C<;3N6dJZ0C^5N$@~DD;NJVM$`Qkh!gb%QBSYV=9%`z8xK?tCDh5!30
zn)0}xp;A>8a!h^2yh7s@6xBGas3RmKd<XSCor7~KH8`XAk>gV=*GrB1<39$vAV6R!
z*Q8SMr6d?fmkf8ivmB7%Q_jk-!ddxma8`L0&ML=pmWlG;i>=HQUoEZLt4OQ%H}qJ2
z70#-E&DmLZ+EiN2*0|ZdZjSqB!nxR#$^}VmmMhIhqggN23qqk@Z}1+mRuCJ_QWbk@
z<udeCg%b1#&9YR0KPBuC3V{3zVoj|5^8e9&MC8b8001A02m}BC000301^_}s0stET
N0{{R3000000051NLwf)K

diff --git a/public/biopet-tools/src/test/resources/chrQ.vcf.gz.tbi b/public/biopet-tools/src/test/resources/chrQ.vcf.gz.tbi
index 845dc40b2ef04fb67c5e54e06798e822c4c54bec..0cf22460546e4c2e1f9fab6448e4c95be79cb8b7 100644
GIT binary patch
delta 24
gcma!#pP(;bow9Hf&+#1+5&=gU-hPt%Hql-I0C-0VGXMYp

delta 24
fcma!#pP(<0dom!h?X(gP&x$05jl3qY6YUiMZO;fn

diff --git a/public/biopet-tools/src/test/resources/chrQ2.vcf b/public/biopet-tools/src/test/resources/chrQ2.vcf
index e49f468d7..5c98366c9 100644
--- a/public/biopet-tools/src/test/resources/chrQ2.vcf
+++ b/public/biopet-tools/src/test/resources/chrQ2.vcf
@@ -1,11 +1,5 @@
 ##fileformat=VCFv4.1
 ##reference=file:///data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta
-##UnifiedGenotyperCommandLine=<ID=ApplyRecalibration,Version=3.1-1-g07a4bf8,Date="Sat Jun 14 16:58:07 CEST 2014",Epoch=1402757887567,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=1000 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false input=[(RodBinding name=input source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/Child_7006504.ug.chrom_merged.vcf)] recal_file=(RodBinding name=recal_file source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.snp.recal) tranches_file=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.snp.tranches out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub ts_filter_level=99.0 lodCutoff=null ignore_filter=null excludeFiltered=false mode=SNP filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
-##UnifiedGenotyperCommandLine=<ID=ApplyRecalibration,Version=3.1-1-g07a4bf8,Date="Sat Jun 14 17:01:08 CEST 2014",Epoch=1402758068552,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=1000 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false input=[(RodBinding name=input source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.snp.recalibrated.vcf)] recal_file=(RodBinding name=recal_file source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.indel.recal) tranches_file=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.indel.tranches out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub ts_filter_level=99.0 lodCutoff=null ignore_filter=null excludeFiltered=false mode=INDEL filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
-##UnifiedGenotyperCommandLine=<ID=UnifiedGenotyper,CommandLineOptions="analysis_type=UnifiedGenotyper input_file=[/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Child_7006504.ready.bam, /data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Mother_7006508.ready.bam, /data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Father_7006506.ready.bam] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=[chrM] excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=250 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false genotype_likelihoods_model=BOTH pcr_error_rate=1.0E-4 computeSLOD=false annotateNDA=false pair_hmm_implementation=LOGLESS_CACHING min_base_quality_score=17 max_deletion_fraction=0.05 allSitePLs=false min_indel_count_for_genotyping=5 min_indel_fraction_per_sample=0.25 indelGapContinuationPenalty=10 indelGapOpenPenalty=45 indelHaplotypeSize=80 indelDebug=false ignoreSNPAlleles=false allReadsSP=false ignoreLaneInfo=false reference_sample_calls=(RodBinding name= source=UNBOUND) reference_sample_name=null sample_ploidy=2 min_quality_score=1 max_quality_score=40 site_quality_prior=20 min_power_threshold_for_calling=0.95 min_reference_depth=100 exclude_filtered_reference_sites=false output_mode=EMIT_VARIANTS_ONLY heterozygosity=0.001 indel_heterozygosity=1.25E-4 genotyping_mode=DISCOVERY standard_min_confidence_threshold_for_calling=20.0 standard_min_confidence_threshold_for_emitting=20.0 alleles=(RodBinding name= source=UNBOUND) max_alternate_alleles=6 input_prior=[] contamination_fraction_to_filter=0.0 contamination_fraction_per_sample_file=null p_nonref_model=EXACT_INDEPENDENT exactcallslog=null dbsnp=(RodBinding name=dbsnp source=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/dbsnp_137.hg19_nohap.vcf) comp=[] out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub onlyEmitSamples=[] debug_file=null metrics_file=null annotation=[] excludeAnnotation=[] filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Sat Jun 14 15:20:24 CEST 2014",Epoch=1402752024377,Version=3.1-1-g07a4bf8>
-##HaplotypeCallerCommandLine=<ID=ApplyRecalibration,Version=3.1-1-g07a4bf8,Date="Sat Jun 14 22:28:02 CEST 2014",Epoch=1402777682364,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=1000 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false input=[(RodBinding name=input source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/Child_7006504.hc.chrom_merged.vcf)] recal_file=(RodBinding name=recal_file source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.snp.recal) tranches_file=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.snp.tranches out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub ts_filter_level=99.0 lodCutoff=null ignore_filter=null excludeFiltered=false mode=SNP filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
-##HaplotypeCallerCommandLine=<ID=ApplyRecalibration,Version=3.1-1-g07a4bf8,Date="Sat Jun 14 22:31:13 CEST 2014",Epoch=1402777873043,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=1000 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false input=[(RodBinding name=input source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.snp.recalibrated.vcf)] recal_file=(RodBinding name=recal_file source=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.indel.recal) tranches_file=/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.indel.tranches out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub ts_filter_level=99.0 lodCutoff=null ignore_filter=null excludeFiltered=false mode=INDEL filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
-##HaplotypeCallerCommandLine=<ID=HaplotypeCaller,CommandLineOptions="analysis_type=HaplotypeCaller input_file=[/data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Child_7006504.ready.bam, /data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Mother_7006508.ready.bam, /data/DIV5/KG/kg_wes_mr/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Father_7006506.ready.bam] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=[chrM] excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=250 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub likelihoodCalculationEngine=PairHMM heterogeneousKmerSizeResolution=COMBO_MIN graphOutput=null bamOutput=null bam_compression=null disable_bam_indexing=null generate_md5=null simplifyBAM=null bamWriterType=CALLED_HAPLOTYPES dbsnp=(RodBinding name=dbsnp source=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/dbsnp_137.hg19_nohap.vcf) dontTrimActiveRegions=false maxDiscARExtension=25 maxGGAARExtension=300 paddingAroundIndels=150 paddingAroundSNPs=20 comp=[] annotation=[ClippingRankSumTest, DepthPerSampleHC] excludeAnnotation=[SpanningDeletions, TandemRepeatAnnotator] heterozygosity=0.001 indel_heterozygosity=1.25E-4 genotyping_mode=DISCOVERY standard_min_confidence_threshold_for_calling=20.0 standard_min_confidence_threshold_for_emitting=20.0 alleles=(RodBinding name= source=UNBOUND) max_alternate_alleles=6 input_prior=[] contamination_fraction_to_filter=0.0 contamination_fraction_per_sample_file=null p_nonref_model=EXACT_INDEPENDENT exactcallslog=null kmerSize=[10, 25] dontIncreaseKmerSizesForCycles=false numPruningSamples=1 recoverDanglingHeads=false dontRecoverDanglingTails=false consensus=false emitRefConfidence=NONE GVCFGQBands=[5, 20, 60] indelSizeToEliminateInRefModel=10 min_base_quality_score=10 minPruning=2 gcpHMM=10 includeUmappedReads=false useAllelesTrigger=false useFilteredReadsForAnnotations=false phredScaledGlobalReadMismappingRate=45 maxNumHaplotypesInPopulation=128 mergeVariantsViaLD=false pair_hmm_implementation=LOGLESS_CACHING keepRG=null justDetermineActiveRegions=false dontGenotype=false errorCorrectKmers=false debug=false debugGraphTransformations=false dontUseSoftClippedBases=false captureAssemblyFailureBAM=false allowCyclesInKmerGraphToGeneratePaths=false noFpga=false errorCorrectReads=false kmerLengthForReadErrorCorrection=25 minObservationsForKmerToBeSolid=20 pcr_indel_model=CONSERVATIVE activityProfileOut=null activeRegionOut=null activeRegionIn=null activeRegionExtension=null forceActive=false activeRegionMaxSize=null bandPassSigma=null min_mapping_quality_score=20 filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Sat Jun 14 15:26:18 CEST 2014",Epoch=1402752378803,Version=3.1-1-g07a4bf8>
 ##INFO=<ID=DN,Number=1,Type=Integer,Description="inDbSNP">
 ##INFO=<ID=DT,Number=0,Type=Flag,Description="in1000Genomes">
 ##INFO=<ID=DA,Number=1,Type=String,Description="allelesDBSNP">
@@ -81,5 +75,5 @@
 ##FILTER=<ID=IsIntergenic,Description="Variant found in intergenic region">
 ##contig=<ID=chrQ,length=16571>
 ##INFO=<ID=CSQ,Number=.,Type=String,Description="Consequence type as predicted by VEP. Format: Allele|Gene|Feature|Feature_type|Consequence|cDNA_position|CDS_position|Protein_position|Amino_acids|Codons|Existing_variation|AA_MAF|EA_MAF|ALLELE_NUM|DISTANCE|STRAND|CLIN_SIG|SYMBOL|SYMBOL_SOURCE|GMAF|HGVSc|HGVSp|AFR_MAF|AMR_MAF|ASN_MAF|EUR_MAF|PUBMED">
-#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	Child_7006504	Father_7006506	Mother_7006508
+#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	Sample_101	Sample_102	Sample_103
 chrQ	50	rs199537431	T	A	1541.12	PASS	FG=intron;FD=unknown;GM=NM_152486.2;GL=SAMD11;CP=0.000;CG=-1.630;CN=2294,3274,30362,112930;DSP=107;AC=2;AF=0.333;AN=6;BaseQRankSum=4.068;DB;DP=124;FS=1.322;MLEAC=2;MLEAF=0.333;MQ=60.0;MQ0=0;MQRankSum=-0.197;QD=19.03;RPA=1,2;RU=A;ReadPosRankSum=-0.424;STR;VQSLOD=0.079;culprit=FS;GATKCaller=UG,HC;CSQ=A|ENSESTG00000013623|ENSESTT00000034081|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||||A:0.0078|ENSESTT00000034081.1:c.306-110_306-109insA||||||,A|CCDS2.2|CCDS2.2|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||||A:0.0078|CCDS2.2:c.306-110_306-109insA||||||,A|ENSESTG00000013623|ENSESTT00000034116|Transcript|upstream_gene_variant||||||rs199537431|||1|3610|1||||A:0.0078|||||||,A|ENSESTG00000013623|ENSESTT00000034091|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||||A:0.0078|ENSESTT00000034091.1:c.306-110_306-109insA||||||,A|ENSESTG00000013623|ENSESTT00000034102|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||||A:0.0078|ENSESTT00000034102.1:c.29-110_29-109insA||||||,A|148398|XM_005244723.1|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|XM_005244723.1:c.306-110_306-109insA||||||,A|148398|XM_005244724.1|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|XM_005244724.1:c.306-110_306-109insA||||||,A|148398|XM_005244725.1|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|XM_005244725.1:c.306-110_306-109insA||||||,A|148398|NM_152486.2|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|NM_152486.2:c.306-110_306-109insA||||||,A|148398|XM_005244727.1|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|XM_005244727.1:c.306-110_306-109insA||||||,A|148398|XM_005244726.1|Transcript|intron_variant&feature_elongation||||||rs199537431|||1||1||SAMD11||A:0.0078|XM_005244726.1:c.306-110_306-109insA||||||	GT:AD:DP:GQ:PL	0/1:24,21:45:99:838,0,889	0/1:17,19:36:99:744,0,603	0/0:42,0:43:99:0,126,1717
diff --git a/public/biopet-tools/src/test/resources/chrQ2.vcf.gz b/public/biopet-tools/src/test/resources/chrQ2.vcf.gz
index 22435b2c513dc40a2f9632f1970395188292aa67..94a70075ea0258cf6132cf9a7a7dfedb18542fac 100644
GIT binary patch
literal 2785
zcmV<73Lf<ziwFb&00000{{{d;LjnNA3cXotbJ|K0{Y-wv)J^SfQj4)9UIwmqOF<x0
z7cfSc?A|ZbSQ^ljg+^BoKeqVs+cOfcF$qfI#N|Q~n(jW+-E+EUc>jJ$eX?ZH8YlX^
z+kdQ8)OYXSM`TGN5;#QXuiKZGmo855r8$__FaNMEcdz5i6;2)&_i5nzWT92+m)=Tk
zEdu7@@G^B`r?P*pEODG*P(B#<Z}i^>rf!bQ<8*ybB3&)dUP7V|f`qI{R5nTML^Mn&
z3-sS;VBXu~$#4HU5}j>Dl}xnn<CQF@DvDx}fUQX^NgI-3JBer@$z$IqK8ejwuZ;Jt
zt$5|j;>$E}_?<PkC28{;q)%A4Mg`=t1Hp;mR+1eSkX@003~8JkiyB6^h}|^}nBh=Y
z^7v#@Ksbz8LTE5yF-=}U@0vwaIOf<mEkC&~f`Iui6OYKICk4_ac34Cv9^`nJ1=8%h
zRsjVLh}m3;)tgRa-xrbh=sGk+(n{9H1!mqI&yS_9i`X=l5p{5Ijl&TR`w{s|($OoR
zqnAc3Bwta|EdnH_U)Eq+XUPKM#+f34w=!4E3wY1v3^m6-1)6-7dD9{bHX`7IWQ5}f
z5=r8-0^$h{T(X`DB~tW_vy*px8bBC*BaWFvPt3ry<rQ#}sa+rq3{_#q?9iu<l!lgF
zgnJ?!BpzOig>WoCnih~>Qzv56u^|O<&t1k>@Mnt$IV=qlL<6+iW<(i5g+v599x{eA
zlvRc@ACh<k#du8+i(C>R93vl`KwQauUozkK8BR_b=h>S=X?^)&n}v$&j3wAdff$I`
zGB*=*GY-FZ_f9xg{|U$Bc8Y@sJ6%hLPE7943qg-)8L>6`Onrwvvj9P?A5fgGQ9|N`
zTQ~d!&EpE1k|nwa1L&_5)>J}Yt;tVkb6TuOG$MRd?9otCP$-{`5Gv9$h1Q6mh+r4F
zBuu;xilL(rU&LsQUyz3%3A!gFK+r5D&=|S$>N$({Ge{vEyl^Mv;SG5hT=o=X_wK(X
zrNocV`@J7w(ME8+Qvn!?Q-Wc+2P<TUI1o3Wi||pBmwRKgh?neyg(%Ex?sftcMn;d=
zms%hDJmiT74uo{UJo&;4(2)|9#vEs##vY0O5hL;p)5@|YVGxysz+;?1KDp?gVrd>w
z7JQMyf>h2U;Psm?C|HEMcu+`zCz^N&Cuq&M-6w`(vWf^az#N;mH6j@w1osgkE)3~1
zvRt0Gz83%@0txC+0(1nc9HuLB(an|-$F*IYzyJ@#rygSv;*dlaS<;{<=w(37I>g{q
zgjW!GjvXYgapb;Se1>)l;H<!-A^7kF^^b<VV&pl(uu@<x?~xoTMIS)XM?q=aIzC|9
z7=3_Apn&j>6-U&6*)tr>n%^7Uk3UL>d>#c@_<Lrt+){RdiS&PbBDd#rXgabHbV4Cd
zqFX_sFR~J)cy*98p+8B)GB%NwesewOpRWgxc|8yWyw0I~p)QPjmN6U5dyCoB7>pP8
zVAeZJy)%!-=n=vNjR$%Hr>8M26`+;-l=`j+N+1?FkUTI<Mx8Mc`0fExE-LSLZov54
zBgts?1d<ic)YDLIWfngXk9lqtsBcMKzO%pFM*IE+>I0NLsVSWpQa5%ofyTSW0;mh>
zS=Oeatl33rK%viy^GD(yERB3=eV~Ekr-BT3XAU1FBzFb22ea39Kkx=4mX$%0!S?ZI
z^fgC$yepTltCc;Qp3RhPxWIQ9#U6{2W6vIIoZs5Rn*vENWGlp6C?F9hJCWandmcl!
zK@Hw!<|8yzX`H`XVjtSHG~k5A6IA3S?;3F$^DP-#)8`}w53QL?q;9}T{V<|Qk??V{
z4S-xT7Jl+5q{vgwsZT5dFGJDer&t_QPNmjK2<430Girud#6n;8c)0eqtsG*hyrJ@`
zAfca0cKe3hgF1~IZW6GFl$CZ>#_W$SUoBCA-4P$Wo2BUQqaDfc$3rhje1@asX4!>K
z0!p@9z|?1tfU+oeGtA!=Ze<Z&(Ev=nAm~ZNyz%zEvu{nE+|V(l^>3ylW0vi1<U_PC
zZ|JC#_3iPlZUdDd-YR17_`!yUzsgGQLESyD`6HVK^dZ|Xh#qout84W-*vcO7MPY+9
zK$ufWo_HT+^OjurPnzGJVjR`|LoKrc>usw<7M|>m53R3}K0i(N>Hj#PA9eGGO)8JY
zZfnVV9hWQxbL4@Tcl#)Y#*0D&yl+eUgW;?<6<CJsiBD0cu&2B?c#+2g3r=US=p2s>
z#%6EWYE={^VQ^_F65B*6L?tGHsI|NXCYc=a$s_RvbfF?+6uu<@<!H#<Ho8zNwWeA_
zzw783{f;h_O0}xh3Semg7FAII(?8z|4Y%K@D6Q%@nSVoH2;)0fWS;p77TB==SWx(m
z{{)O$<x15WTO3M7(;DYt{7^hj0pTvFwrb5<rKMffzrp^S`C)H;7n962GME`-*L3BH
zpHE=Zbn2LDzvqGVDhBdqWX;Q5;z7-=|LqG(EjAv$na%<99nWKyhyv(J=Lp*ti}HUM
zmmjKdoqU^?9dDh~IEeEGOWv7tEuv0}#_O<99;q;*gwqzqe%FJVf4X;DiuwlNU}ktH
zfE+Ht_zBZ1y)kFo2uPDUXH?SNKqx`0EqU_#N>^Fy<wdt;p9HJK)73`3sUCMY_HBU*
z&yPv6fI)nC#Q_ozf$gcok0@Z<Fz-z&s4sq6Y@;m38>lQ~(}xxn8nW$m!KF6`>NfmW
zoW-_B-gHfS|1imZBiesq@L!x3nE&JiI$Y4Q={-aJ2D}Rqp=^|~7#aOdFS{7SVQ<)5
zjPFJpb70SmaksaDzS|g^o9<9-*Q|~G-_fU=VSZWIH+NGQWO3<hYi>K@Gu#;cX{K?M
zU+r;b=q`Jn+<hALOz!^f_4H;`f_-6WV3wx6ehGH}rQ17WSju8t;_feH#kORNR9&cw
zy1&)-x7Bw%fTg-piek0Zs#lw}s#=<rjFMWfsTEZ#O$^&E^(`GTA!0$NZ|Z6A5U{79
zV~zCjXrb1%+Et^Xb*!Oo8zWOyJKc%S=|<^vEgi_ZQH8s)u4%1WxvDkcL#Z~jvZ`t=
zcw^cVT~(SLqpNEjqp$NiSM3;Mz0o;58P;nRrE%3UKXptPuhlwzTURSpP3wqrB|W=t
z#i1{100g)yI{(}%fhwxDnw?uySD`*sJFxkL?nUcN?{uScyk+BKYhVYW*~#`jd{WI;
zC*N%7eY>-7ne{uXeBJFpzU#)OH@18B%u>Xk3UjUIU(K?wsx{?G-Efl57QrU-?INdv
zzc2G_CO!*RS=NcagAl=!D*iVcqs^V&ygFf{qPCq%RcTzPs<IGQrA33-$kdgMO&3;#
zR?&9XH^-DKJ_F;t!&J3#2rdod1lr*Bg8#z2#dsP}wV^5!sBf^Vw0=RLTIT~*z(dMU
z2;XZP0;$##K>0<&sn)Klt*gy{Mhis&VyHE>YUQmwn6b+Nx??M6Jop+}?Ux{{eFw7o
zFF{r>itK>pe-^9*6W{Eu<}dM9^E=|$_$A01Uq@E5W^Kc4o0GP6+nx+d%B9-YYGqAr
z*Xr$7t9@0yDl6rytCo19Hp^<OU2X6;%~}oKG?Xg5QQ9@F3?EhgLMcP9UsjuH^WFae
nbCk<Q_ay)TABzYC000000RIL6LPG)o8vp|U0000000000<RUl{

literal 5357
zcmV<J6cXzniwFb&00000{{{d;LjnNM6zyDjbKAI*|J(B^P?_4@q*jqh>XwJQMeDM5
zbXk<0mp7GSQ6wbsnIae-u{Grnzun*=ik3&p#Boyfl$x?cfM_)O{51g0ufHxVn=Y8Q
z)&u#r)%{W~iog8&t49~qqmDsk>~5Z)pPPE1pSSzBmGj?w=V`ZpzSM&|Z5}$NO*N@l
zIbSWshUT!9?w*H+ZxpxB#f9z%IuO5ctc69*9(7o-ajDm0>$UEf1IwZEZ~eBcxURjK
zQbV_`xu*vfa|*ZA^Wjr26~*`B`=wCR%kxFO(AEPgzftvo{9ou0u}s8jvr=ygHPY&+
zGa?CM`AwnYGGir+WkIS{YV~@pQmqv-Q;uDn!I$6Yj&5&!%h!;pJxmbGal=4EfyjS+
zBL0fqcSGC0)Yk*c59Ebz`;>T8H$fZ=kdNkDe^S|jSBSe}4%Jp{4Rk-Kgv;!XZlcaW
zU&<pT8VSS$>LGc{fiGWxX7mz$FznEz{oU5&CC#US=7qr6&8R;fZQIw>TUyRiRtAG@
zqpO=HHVP!wBM|r;a(%y;qKU(tHiaqImILy!j4245tyN0RSXQRx>vNm7tzbnxWGHH0
z>kre=2`muSdol9_W)t}Gus^z)byQ-qd&k$;t_=etja`1!RAo3Bbhdp|3z)X>bOZGl
ziEOveU|*=GFR3gFf<Wf_XL;Q1;%85H$1Zgy6pSL+$Yr4@kcIW2tt}tK!AxpupF^;X
zV59j4^Qa#)uiG}e*IX7*;29n^2S<cH9eWmNNw=?gzm;~GbfJfKaH;!r+XXTKpWz-0
zdKs`~w64o*bsBwxSzzuzDWA><Eyr@!EXSmF+!{?l*H)mhc=nwh7Z?O<0U{BznoWK1
zjFs-lRkF0%9Jp1nXX07R1ExSeU|W{ST}AT)&oTlB#>i@a&>79Z)Yc$t&Gp#Q)7P3u
zVX&FZc1NlEOjoQcFhYGzwYBcLsBrFte6=BHa}j?}7OID!T3qdeW`J99j==adn;CaF
z`hG<n4gC!tJeJPoglkwLcKEW$pA0ulS~QTlhRzW;#BW$E9U_o50~=-Dj5l-^3(K%T
z_nQp)IV=#F*FcisNwFVu-(PwxbhWvy8+V$G5hV&8)U=i!<sm2et_Iv)GA*~M=x68(
z>TKF?&0=FrBZ+noo~DHkmkgLA@Dg*jqz<@2K-HBF{&Csr#_eBp&jMqLWy!s){c)rW
z`z=lFbbinEyE>>3Jsxe;W=o7bU^ic=E%!&=F*b33Om(-Swwufp2l(iTM_T!hw^L?b
z0&So}Q2jD*Cq4_kCz0w7)Cy+0_Rc*piF1$;(`tfHtq5g^jE*J<J8kFg3cN}>Z>=ob
zOd5;fvS_S4w&rn@n#C{1;@v0E4rGCd?I#><`(BhSFPnI@aGB%NcibX><sAt;2m>pa
zg7f!UlP9wh26|%NvN-n`x)a3tAPi^`KAmL|Q%2EZ@D$Gd;%#IsEe7`Z5Pa|~uoM+~
zAW>+>+Jr!kAxc|;FMtVsD^VZ5`5JTV%{P;p8;ftI5%_U#j~$}XfPiT;vjw4Qv52$z
z(t&`QV1R7e@A9TJajqeo$m(c9Vv<pO;@&;d3~Q--@r!Xkq#%HI#tQv$4l3rsFdOaT
z_pr~093DXYX6bZq{^vVmwpvpVn?n6L%vKkw^-4uLVYU-yJ7KmHW;<cFe|OAwEI8t?
z{KP1Z`}U8+anbZAB=^ISTz}N=41OXkm+wCUdF6-f;I0=jW(a(K&VW%=n<8`+$jePJ
zWQc6yuKM?xqpN4;sBY(YMS%JT)Vae?Cy4r;A*xh4A=DE>Jt5Se0--MB>;!VSqqenT
zpd$@4uq|JXXII2EJdJuDgCvZ56kvRv_hn+RHKa1C4#w@o8y%nnZ;aYXV%NImX{$Bp
z4b5r|Udp8=55~PgM^&|!(z@!8dIVxtl(ctzRu*d*k2GL}a;r*13&e*C%3B5WF&QKV
zg;XN1GZ@4$4ZK#1*M<tLWL`)paWjjI9*9a6!n+|HA-=Jkkh7bFB~eJbccN81>PmNQ
zE=LtpQ9XXSP3K_}Y9sigjwXt2Q#%{QrfAaYWaqts?mz-sur%4F$s=-70I?82o?!L_
zso#t)$2X()yQc&BYuqH_2H@K=H?qV9lv4rMz;0Vv;Mg;5b$L*d1U}ql_aIVqQpni=
zN>B;NDbIk#Mx=x-OD1)L74BFhDL2l_;4m4!5$VQa062~aiUY*Xus_pom1$oY%~WkX
z8hivJhr#U6&62^O4H}9dMvBmKT_UI(?RaY;5z%&EZH;d`(~ktQs$=S&sUgb-a~77#
z+1bB136s`wLog8p=@>mWhQpdfRp=-k2fR#6##Lib8L0un6wDjwz<|EvG;uJ%fPu3<
z^U~HIGLgdtpy7Zz;-Jv^htis9h|VV+_yg&NKEQBpk9LwF%{kPYPuIko_ZO4Lfgesl
z6HB!$0>Vm`2P9nnsbD*f1~SqCOH{eKeH5LE!78J8$S;5_{7g$6<1q;Tk-}tUr}1?P
zwmtjhm8K*#rSfz8<&q$kOSRfF`|nr)rV+e_zV(uO`I6L>>P<m<hPPU+T9-=I@@X&s
zR4ASH@=ul0sZ#nmE2STG|6*nQl=m-I#t*!Iv3K&*uKf?+wSV3F7vB}KmBgl4dKR+P
zYbBvvIzhG*WII8&6J$F<wtsiXb}Tmf$*x~KVf1ycUpzDY3CR6$AotqWFAhMtT=x;U
zD>veF_u_Q-;)G96cP~!x^aM|T8hCnwH^<S~UQK~yZiF_sfR3|7bTH8^?`k-V&;C)v
zEcAbee1I1*rqpM4$b({QJiHuh5LF3ypSv3Kn*@=5z(8`38XhM@IPkaOH1;{{<Ys`D
z$$Ju`($;1rzF2`5UaZCDQW>Tcqt$S$Z9$Uiw6!Z`G8oT3PB6y5`f<4lGM;(XT7iu5
z1;n)E5Gu_}ZOb>5Y3BjLI~Rde!PZ_+$u^b*fy60R0iSlvJ_cS0DV1EiI-2-+LnJw)
z&mJSUY>S6FNSb$QxSmlzD3CVZ3z>ilc~HD+?LBK&U7!Oa+wlb<zd&XXZ0LGQT?mHp
ztIYd!x)$=YTno92-BSKT6beMDeByT4cQAYTG;vzL%e>ac*j|x>gfj6ko-dQz3L;j{
z2x;27vqVW>@yry%37|RM?VjmYn(|?hfK>c23l64q(MlB~kH(`8>D{)vz3WR*9bnT6
za0?t)g-?;)Ae%GRu`MnvK&wDt$mJjk&)#R_J@G<hi7X8l;#PF8j2rh2W_@Z-GfM%`
z6yFvDXI(C-w{1^~HvT$DZ0i~+wgLMxRnP!6dmPE&t3xc+Vvh?bUS9(P8inh$T=M(Q
zgt>8yl|`ve__0$QoBdl$A0$`)j=kWAuxHwfOyFOkA0%j*?u$OCz2wGxq9MF~*kT^&
zJU}Da;>PXmI{v+fk$vXrjvxKMJ91=ie5$fVz(Wo-5pJe7p}RroQN{P^dTwvJU@q{8
z@jC-|Bb)3y$U-y|>&3Bu0`*gBt*q-V^*yX-6^ouQpw2Q_f$FfMlX)?TxRx`XBRY*1
z0x#msGj<8902qW36)$H-fh_XE)_A0Lrnkzhf7>DGy<fooCLTk#9Vhus&q#f5W8Zn&
zwhjBd3*>9W>PtkFd3mTmsA*L=2N#;?zOP!#wI1Ir0uM<b?yi$a!VkScRBejS{^C|D
z)#`Pj^cOdX5dHDpiuR}gE;R>t6$?D&V>t`nTWC|?@GQPzA-}Pl_8d~<o86(aWT+4g
z?b`Y>KMoHW&_~wP&y7}c6GsqwcYK^9+Luoz@Ai`6#pj1FLWi4Z?^ex?Zhyt-J+t?Q
zM~v4EBZ6hQv4bPVB2vRmjNPG%(mG~duvw-UmYI|IWOBq{DER^ku^&Ok6V_YpV@Me7
zk+OgM<mwm>n7x@mIhdPza)h~X$DFKmAAxyp^{w6!1Tf+98oJk+PV%xpX8a2BOOyt4
zvyYBgd24jLJL>9~Ih6$z_PPUEaj3gpkABY4=o~?#6M77g=<`Hc$1vhs&r@jk4h9R*
zTar!lvEud!>6sXzwpI5oh3PR4>v1$b)crf^<%Z9W7#?IGeNOwz!NnVZ{XmU2<aK6P
zdsfg^^B%A_QbX|<;4ajRsZiJ%IWz33$HX2$CBaHxb4S=6Kb#&her*{Zn~s#661TBS
z`3wFd=OO%u5r8(+-+VyH1c{tb-B=N2FpH~>0=^Qg_?{Fc%rmJ+bf4f}4K;JqcXQKs
zZKel%Cr9hWo!8Ab&v{gGodFoOv2uVH30ouu6;6XS2{pGu*d@<*A>(Q^H_@KO{qsE!
zk7VJoHThuK276!*0bSn_KU|Xl%2}LUu>*4J7YURJ1OW;8oa}Sv?D6uTO>>VVLyB+t
zE6dFhILe+C7c$a=1^6l2$|wAX7#C1muA>GYs)7w!>0c<BQ|gf9XBm0#Imqn?7zOqI
z289r-7Yn$|?zH*8Iil3|56it9*QQ99<7$GtM!JQy6-jmmNV+w!I|US|$4E&Hm`mI^
z=OznCb{TWfM?K8)1j{in7p%vK`yOb7SjZutUCUol@4x)0W{#{TV4#%%ssZBTai~e=
zmY!S3_BTgpgNHf5{EZF;PMGS$0uw-xO#le58P0pph)<UuUFx_U0Y;D%hP#vZow-LT
zZcMZoT`UgIz8wI8h{l`70eU$sm)oKlxe@Y}IA><-&&TtC_mR;M+}h^*NjeCX=t~`U
z7Wv34-7`1uKLFeUc9tMf7i736^bZG}<ANs?Gu9e>dCN(bD0vHn-tmxzv+J2i>yx*z
z2;xvgtWZANn{9-HRpWQ1_3m54A+F<yDExO;u}DJp5ew;l`$Db{+t75Vl3`RhlIW;t
zlP4k^pB-|P$UmvaqG*yQ`qlZMdpI9_iSq#uz;QMFTwWM;dc3CBW>cj<($xN}b1?MI
zR+gWfT#$PZdgz0zKv<cFmThvSINZr0$U}sQLZ?qTe{*0d9-6Or#xVKIa}tHwJ)A7j
zQuniRosZ&s^bv1ck?ro@d7|z4CfaAtfzpYRLt}dlXp{mAJh25I;4qa$V3&lB1vtx(
zAE}vfjo7u`Mtgp^&~z*RU5=BRBg%u>(|O<W9gjRQFfu3~zXz__vEwvdKFwC@Y<kdA
z5_N&>bi8`t1-p{nrFeUdzdt^_m;qZ7W<mntw~n}dQ<S_8))uCxfBbAN^y9vGpVh&D
zg$By)5>cZb`dE^Ywe^5ohWxfn4h<Ng?Rr*ljQUnIa~c9!%>om-k7MdR^FYe+Avp3B
zhtl3&779^lk3uu<r*nGIW7hW)D2IoSxq@o<DCWBayDcDJL7pZ?v=z%San4*I;Qbz7
z_BmV~;mheNjp)Kp5@hf(D}`J<)4gE4XaSPIld2=8`mmySjUyH_9^rGbCbtVf0f0&P
zy7i>r*g{hyo@h6ub;r}8GK;F4{3hD-1)XI?xZO|lHiQy_P*DeoX9^zdic)VT?`G!w
zmZE{Yjp_xCA&1wpRH?vR(c_#aHZTLCdo(q?clqgi`NY3>`eYYlC+}ye%#QH4gh+ny
zB;7vj{Z-J%Lj~LV69f7Vn$Ln%)J0laBGxhKf=7;65x1=_@d0>Q0D!l!wA&xdI#bTe
zfZgK~^C9d$#s-`G`w#`Ev*^MkzqE&U|Ac_SM?=U(8!rT{;Nh1QoeU%wJPE(Oz*&#(
z{qcx`B>@m817<eKd$CxqiDmMeOdiN@<h@WVm89|!yflCpQNT0T|9&ZMaQ<phXq3Kc
z`WKW1pF9;sB$>b91q$|WQwmo6d)O!!>!MUm_z;ScR6WebH|b*^BPfDWqg*Q&8&bXU
z74yHC9JaHkI*CM+XJ*FOm9<>);}a;F_H9#bYaZ}d{vv{rHKsYd6-Fxm@eI)(8x3BJ
zCf=lY7_opSKr^?GkU%WU;fK)=@umvW#4;_~-r7sC@5cZuUzsBno+ml*$|lMk5{egK
z*s^rHg+EB#kKGnLzX5ZwGCdtY_)i!a?l3)}>+{>55Hv0007)x-5mL}H37(i=tz`tg
zR^D~N<_D{?Sgq8=U4)}vACcj4n<Ty^N_fVBAnro6w+!5(fU@DXGbxfTKVfc?NXL(O
z68zW&hzfu#`Q$gr9y8<~@f1f(F!E!ot!^JC(HYUULBSIzO~)fApkV@)$IgSr?@ei3
zL(yA`HdMNgo#+E^77aSu=w|rX?yECp)apC}?p8+a$JT(CYrRMH<M47kh(9!Sd^3fY
zdN}&3cdHuwuluNUrxC?r{HczjiEg6j$<5`k(?;>PuBPMR8PtVm{r1_k(>;Ug|Lpok
z8JtBrK11=JMalMT=f3pW?zQ8y-J8z8pf1lU!kOobjYg$ZE0@Hx*_m=CR?1>gl+GrK
zs-AUwG8l!&oQrN-4xKxP-8&b(p*$LDVnr&~t3~OeH;`3j*cQc$)`Z{w5iVLiyvtcF
z!MBktNsV%$B-P-LP^wA=QIs0c(^e<4DAX>LmMmQ;T^Une=|UOF)r;&Vv0N?+)%r#I
z@}dpzOXZ8MDvQOEBwg?=COP^{h66;Z0t@gd$oMxQ0tr=Y)Gn^uvIwc7bOCiJ02%3G
zdLt_ryVV=MTZVZ+u@_MZgo~;*F5<dF?y47C^eo@>3RkTQFnw8h?2OcoI_nAip9pI$
z#Z9xQsZ<u~;v>I7!TsSelC2g;g#TQ`hD~keEF%Nt|1v6qCz1d6_^32d*tPneIf`P_
zD3*lkdr=fL{wXvp$5$e91?8~?&LI`0^z+4eiU}WJ<FLR)v6^Krbo~Gj@f!c%Q#9dm
zKSQOeDCC&>ig|^`D=4aQSW!nvNcaxwdpZZ_RBCWW@gv8lSgw~E^~Zk<H9>&DP_9X(
z;!8;|iY^)Mc4s*t!Ka*+Uxl;s-{7qBDx6i0<t!8BzZYAXDZW}-wO5f=?QiI@`YN1N
z|C+P2-mIy#o9#)ncio%}&V+NZDU}P7*eq9?jYhLxsuzSpz24wGVyz%Hnx!iC)XHV(
zsR||N5t?PG0Dnr@Art`i7sQ%a`{n-uhkYa;#cKcnABzYC000000RIL6LPG)o8vp|U
L0000000000a9c)|

diff --git a/public/biopet-tools/src/test/resources/chrQ2.vcf.gz.tbi b/public/biopet-tools/src/test/resources/chrQ2.vcf.gz.tbi
index d376218edbf3aeb9bcbf9a16275c36a6005c57b2..a1c3a9073ca4fab9c79fb6f8070905d35567de33 100644
GIT binary patch
delta 24
gcma!#pP(;bow9HfkNZvuiGZUF{rWRxCfX|i0BI@+GXMYp

delta 24
fcma!#pP(<0dom!hty_tQXGIc&QQCuuiS`NrY%~ah

diff --git a/public/biopet-tools/src/test/resources/unvep_online.vcf.gz b/public/biopet-tools/src/test/resources/unvep_online.vcf.gz
index f102295f99d0bc62e25de296b84dc6610930a683..210f83e016e5d83d303e755b4772cb8299089997 100644
GIT binary patch
literal 3132
zcmV-C48!vuiwFb&00000{{{d;LjnLF47FNqbJ|D}{>=W0sZCv7Qirid?=q{qQovwp
z0b?M}-F|V&(g3OmiIs#M-><(O329)M2<4Kk+OlEQZ_`ir(>-$Dy;}t#S;fiL&FshC
z;6?4=U*5e-$ciK+@`z0z-&|i``)=l5_eYP~_3y*$>T!C#cC#-}&-=&^$&=jCt~YCJ
zJw@@x-Cgg!)a$VCot2wrE*Kw82lw`GqrTmrwx|2;Gf8aRUcT*!J&H22CP}+bQZEU1
zSrA9|uR+v*UQB1d{_{k2S*l8f>L7I2ydIV$X-J}YOHxkT;T$irB#1bDHw;Ng(*DOY
z=YwIX-ubrrY9D#@)9`V@Y4^WFTi|^-X`t_V08R=YIo)vs-8G5GI7qWoRmbTy@w*M8
z*zp3N%lK^8KzWzM8403UoCeuh=)Hau6JCGni1lZ8O(KZHx7mj9u4fI(CH3Nj%r;Qt
zy*T1d-y1eC;DOqL8X4ZsXMFCP=x>7U&dH0G^FD2G^4|3ERO_yZPiK_`o*Uh{yNSCS
zB;-%d#%%)|xBDdCks2etCV`~Ew=MLC7s&!}qfC*(Q=Tgx8)VO?4E5-S0np@G<;|Nk
zY(gLh*~Cr1kc3lTHc-#p$S2!*ZiGYMxp;Ul3?hIrbka2T0>KIT3%&z{Ozi_{xG)u#
z@nRPS9#@9LMHB8x+$h_)+q@A@^(XTN`rE)u;`!9!KxEwI)Wlyj5Ax!Dl%XI(>#`!+
z2qq*VuD3zXc!`S1(552UhM?5l5)>ysNsybO5RyQA&i#OMKL}$t6FR5WyNA*G_QO6K
z4cBFyxgm;jK@_hJZc^&Tqn~}2gj4e$-IP4c-RR3=zvUbW4)@1ZPLF7n#9Q<^2)+1g
z93jB^3rhD}l#w)}-kmT*k7)<Z$qGG#1N7$}dMXRfdh+ANn$9~CO$e1rH^GiG(5RlB
z9ZaOJ0jx#@CB*fSPj=a+i&EGq^8b;dt^0;H?h8TBghU9Ir3@A$pYNWF(S8Ok<VJ5a
z2?uz06b3(jjf&6yzc`~LOfTnskhu98LG8}-fMGagF0^~JMvI*r<sZNnp;BCyZ>Rkx
zSz2WA4($#zw>&@qlF5r3a;wkca>%m{BoNBQrNy_q0vmgR_9=~X5TqND{5?hFD;%q+
znk)uXG6G-I4C={8&w<OG1I9v_mvBxhmm{F_n>G|QVTDHkCGZ+#8{}qa8&kipF-plg
zA+P|`*fymT&iN>MP6+Wqq!*Lb>ausH07L{Dffo><Bj|ECZeELi(MB|`vN?ePzK}53
z#4$j;Bgs`!H0Tv}88EZmxsX)kt^s)(J80gyiT`%>8P+WbXAK$MK@J7fKN;UPBhQHo
zodUhgA~`mS-h-hJIi*qW?h(_b=sg?)MTBmwZW4rVjN#za^t03Z@S}7{&(VN||DRK=
zO3H3<kim~1<l%A+%_j>4o6wFH(L+w5Z@dwy`E=kcp+86x7i_{C{qA`%xI7PD4)Y)<
z;O!xlZ_I`1?a*0{9&ewPb7wStT8x&r7pZr-2~zX|aKYk%ULooG6j}vn<#`{3ehx~M
zH#v|zFib(6Daqq|hL8@b{C(#hod5ht3Yrx}vZj?P4&{MoaRDBW2d_YVtH7r>`|A?z
z=Qo%SF!s2mbmnl|n6QAR)nWnE1@o*}Q&F+(qJ0#=o|PWH691?*4z2Y*h`ey0lVSDF
z;R6@Rrv}%9<+<PYbiqh@V^Cypef%r-nv*)N#^w36S}f-mE2V@Ba_6FS6DQeeW=}0X
zJ}k!f4T@kKuTkv707)p><@HTy=3S^Zn89UbzQ97YPY=&lZU}4IKB9z03pB4wx;47{
zlwQfuHu#$DAw%2PC)_sR?ZaIXWKH0sWE%mw7A*X_37|w?51jfMCy-?rdh{(8jp;z8
zLtzN*g4zpe1}xItUh(y~_sUfc#e4oj6$Uv8{gq_L7V-_uY2;Ct5KCUN+^x!4{@$bA
zk~Fv-QQ_*8`F@{NBtyR*XMy80H_7f-J=i2*Wcv+Vecli-mK@Ry=lAop#!0XaB5?JF
zU?<7X?Uv6ywwiheN2ioFxSvm)WpTUVU!r~cLKk?&-kw%-8>j@)RndivA6<Cpqhj<P
z&D|rPKXPe6?~D6FzC%tPY+2L5SMl{ZA8b$t$fZ<rD7+86`yoI0h3?C@7$<Z8*vhg&
zzm-OK<yrOma9BtB^lfq&{Ff5?$u@tyq*5%Zt0moae6q@!;{e2CZ-7!*yaHH&*;+an
zjhDCcJeG0%N{`4>*y~|8c;n@P1?Nj>I&VWp)Bf$)vO1EK#qi6LIJQalfJ#cDeAdz#
zm=$^`Brhb)L(feVj6$yjU>uEO{{~&*j%r{P{br-D=r?pFbreNb8^mHkELf5trhojb
zXsCa^BU#FKx&Mi|$emZIC^AzM3tZTL3>13f7lfg9Of2hV98yP?^~+)W5FVmHXbKcd
zHPntJo7#8y|7m$(-YZOs)bPP9<M=LEZ}LB%z@<s_nu^_dpkH%Ax{Pe;n9DXWbG5&H
z!{{)cj(?iX5$qjbr*W1KK%aY$P_9@3{SV|0FIA{bdQB@{ZwW1q(!+wq-^{6%d`{+z
z*YTh{G2tW`rLDjXdmEVf;@xeP&o>ATPUcn$$fFuDeZzESH>GTw0BHu^1(QJ~2qS1+
zk|&+7!Ma$zHp!$Nl4zZ6Y^-YrK7WqO=dx<anxxAD8JWo_*s>HuIhRo?Kb9%UASj*7
zs7wZ9O_dG(Tt=&8Ov%EUtPAYum5hPGkfa;uGDan%Sy(Z3T@bTb$*8)bs}T0NjKyTI
zj1@?LAcs;Vq-dG~?I#Eb<0BzS(lsnNJ(ih}Vrtm5G=UJyH3ln|X_x{bmTbvTR7+C@
ziNq}1lB_9Mu>={wRk|%n(NsedWD-~DHZis|8H=1%={8KsFl1E}b(L;FO*So4;7hKu
ztt$|<CJTJYRkBrJK1-NJxyrSM6-m`}0UrNmRiw#=q=~BWKP+VE451nPv<`vjpzNRM
z%S+&jD($+ei1Jw~$eaS?U$q1t)KWp=J(!kZS~UtPuV6tp2o#J`K|AvRKnX%JO9h=%
zP)u-9BX98X1~vrvS*5(m$y<_S*2!D^0cZ;e0y-$jl5}=r13IBbfrnYqv9Q{4sldam
zN>~!MS*gIotUy_aJX8ui%qCWZO<F0O!EC6SS?hwwS%u*)VoS4B;DJ_UUBo|CTF-!%
zbx~<Eg)^jqRYbKdmq1`1YM4=`Af``IslZrW*M+^$L>af6x~U0L$wV2$>#D2>B4(nD
z=M}@iwN@G18?vSf3Ydu=jaXT>1kJ%jjU!R4*2c;Pu8y&2NHS42cMS`RO1lze^8l>E
z?olU-kN8}%>H?5Y(9#8U%tRSlCd`q#yvpj)v2;t&;!Lzc3yKVjRxT#W(1OjPHdhQR
z*;GWN$3z)eGR$HDS1?fqmIcdJjigdxr7F5rmsEu<rUtBUb&@h;O<k*lU#WmnHA}7a
zt}>gNXys#)6-<h%)+SWZ88KQ)t)#)?G)2oTldM2eU@x!hJ*C137E8I#yDA;WS%;>|
zwk#l*y3DH#BSn*IZ8D4$r4}PqJvzHwcXpLw1SDN+v%*L*;f7KtX)v3r=;FX6k1#T7
zD_Lb2nOaU{y}SASq&2%=v_}2b{C3cC#>>{jCuiI$Zj3GZ=C)P5x@#@+i|-SbF#9O8
zkIFA}HnvOyORC&TQs`9;C_urjUaRM{;JCPh<<`tuELwx14VTL#j=F=sy^p>`@oUr_
zPVDLA34m5jy(4#rV|(FD`WSb6Gn>A_kh;C0eT6%)mUnwoo6e24A{+2SQgpeEfk)s;
ze=(z1lCIOU<*qZZ>8mrP>r8FEd;E`0Y_%ilX1D*b+Xs1B?G6?;?kKX{&Hq$w7r#sO
z$%Cy!2=H67>Ce*0mDGWKvir~nX6i^vcRq7$+?KoZPqx!N{VNcaRl(0<Iqw#KN}-1`
ztnT6O7xrM$WpCl_Ps8?IuQgoWIQ^Uc>}L3IGaI+0YkUI~FXJ0kyRobr6Z%tXo2HdN
zfk}z28%3v245(hItxF0#k#3+p;D<sVNNo&s4i{VFm;V89)p~ukD*yl=iwFb&00000
W{{{d;LjnLB00RI3000000000X>-j|h

literal 5727
zcmV-l7NF@LiwFb&00000{{{d;LjnLo7VTW=Z`;U{|LplIhzu5+U>%t^DSNmJT9>t>
zE284Oyg@)vWJ}^RMQTVYw&t(js^%$*mNlCbc={4c?CP%SI)Bw{<=0;qu16L@xVB>L
zzSaFKHn?AY{WT;D5)$7g8V+AyU0pd=Y+bec_tMqxy{mjUx>{QCgE`;&jz>(sAziJO
zT+{S}m9@Fr+L7HTzBd+D6kD+Pd*5BS#OV<~h<6(jwu1HA@|}U}6YaNtThliiZ#N~j
z<+=0Firv7!yeDA<e>I`Oz2)96S;Z3P3-z*X#YB5$STXazTc6=XhLf))^_o?fR>znz
zJj;o%E<2mRUTK`j@`|LWsv^nCWr^i@gB&95mE~LBE^;Fig{l#PasACUHc=tkAD@_L
z6+Cvgo_A}l2W}K=3(JcL6B5e-bu2(XX5{`wG#^Nq%_{JTxeC^>?gxXRI(s0yWH7dt
z+DK1k!s0Orae3E=51(Ppgor%Z-qs=g{ihj`8Ijlww;<U2QGYxtp07!`bp55K4+h1<
zjpaD_kY#c^!UBJ8slC6Np&38$+XPsyT_5!0+EY+?ur>&BQe8Q2WX(O&cH<QZQK01U
z+IpI9{n!O%z2!>Zz-{748}>)<XB~raf=55H)*BBdCKkK>Xd3!(GUyaz#Eb)T5n49d
zFHzb4m>qm3p|vC$$FeLlxBk?|-7fy^S)1{O_!9z75$`mSZLrM3eKOZ>1nLMJ;+TJK
z!8hWa8QDQdqEvXx^MXfn6T}2$hA+DpXKW)f4qdR4<-MchZa(CYg|+qKTPq^P5a<MA
zhDI#dWe_Zreciq@rpYJZ0=)kslwAZ?PW5cMzC*n9X~KZ4-PjCL?t3E}5D>ft7D<FP
zJrY4=tSn!XnWY!ZL0SVxCZQXIAQb2ae9LucsF+b4x_0ct9L?wtI-?o5+8T6iZo*(0
zT5B^TFxkntJrJ^aW+=54IH9#B=Gxk9(BL!(>29}R%|-e-+o-h#)8cNQOdG<AN(9bl
zdV&3bvmaN)H!<GuOBm$5{OuaHhy#(v(Qo!<n>{f>>m~+A`XHrngM13htOM}TAe_+$
zoyEenU9kPGwEP@4h{BtoNr<F0j%n;K!(h8H=bmLhm>y=7By|wSU513Fob0m^({L$0
zxzANU!%#40(|&7~I%A$m%twecbL&&hfIEULsbovyLm0%wSXmGsx2<ma{IeCh;7qA5
zX_U1;j`U%_Wg4B%@0Ib51?IzuM<4ZqC1xJ*o6p44`lId`AJKTsZMP!cCgYSQ_~aYS
zwAvr9r-5?|Yr{GWnqQ;OnJCzX2bt;t%nEM04zEIRi7U`hU@9ytORNZ)(KlIk@7dL6
z1yRLcwN|d@WDgtLWy4;D!J6hx;xs<ni#MOZI#2}?u@59F#=fdpVY%>_+hs-0$lo+5
z$r~mPAq}j61(ol6G0(U%0oW-BmyLNCU^qdZkGC;tz+Zpaz>?8$1NgR`M~(Z$Sy}=3
z<5T=0biq>$7=c1z+G__AIi@Ic1+f4w^reeN@F_9_-`johQgeIp#cRY-THDirXf`2X
zdV$k|RJB;7)qLqgLd`Hhx$Jl8qr7mgp_*vMXu_mRCi#g*_sF!}r4^=x={S@ifOqx^
z<8clqrpd6J?Nd4&i=c+bkiUbR-K+op#+XgHX1QxteGap!tgK2Bf5vQQ%y!0XXUulS
zZ2#_<?No3?Nq%A!M`QcP;kbnT8Oi-{B-bCcJA<DH%T>ouKwi};d${YXI5Q-^sN%q=
zs$ByHE}5^#Vi+K@Nr%*L!=r1Jc*q5h?+8%;fIfHl=?qc7F+}C1GeSKh)H6c;DG=&1
ztxlka2jaP_0BmGp1@^St@$8P-*r7?nFn}VAXB6Okowp)m2Wu#0#2Ad*Sv2~v4p5BR
zdgj-b8=9*%*bUuk4N*$9rVYluLB}x6mfpJSk9rK`tfXk~#jM6Dn2$^lgix=_Qw!t=
z3H5CN`j`wdhr(5osIvpeVJ1YaneGh@D3u7Zm8hE~P7jOn5<`a}9U;DP{Vf$YOG_M^
z5AS8Gbk?1<@u(gREJbQc+$QsFmTD9DV~i%c=Mk@*#ir=e#$=CfVEIsh7C~O@^5T&M
zDS%i=AP1N|L+bBGx8wJt_M5|rl#IGW`T*p1ot?&01yxKyEwKMgWGVK{pKfSU;#oR<
z6Fh=SF-W0i11Lcwpro_{HkyeM3XvSLiC1`Hk(J!EDnr1O<VK{MssZ3QA}9(FJH!6W
zyw|6FeKa%7@o4Z7oE#<xf9;k5OxmHNuv}sYvogeixzUdc4@nkn_l?&0zBB#EKvne}
zD|Ad0*$(`L>rip(H^*bqI&BITf;gXJrOxoOCD9ZXYR88tGqZVe3MvybK$-%+u>}HH
zdqI;X156mm^%<ff9!i;G13<$EbEHY3^AEi>GZCFnI`9va8}bCxsXuyIg>>f7ZyxT6
zKCW*jPXj-60h1GyG6KR*mL??J{<&d0Z4DHp1D<G5b4RE;2a{FF@K9d>Swy8LX^zJv
z{6`Lx^}UDRr(vr(FPE-)_L>);J1^&1UKAANne%sS0P_sq!r1znXZbvT&8yce|BP&k
zBCEV0i|4cabE9-V%RhHY=T7P8?38}g{fm|TQ{KN=*+209MJ@C5sr?T>wg0~NFTN>c
z6S!+mcowp$iolA(8M2)r+ZnQ*A=??U{kuc9Q?b!ccKzZ&(C>Tw;u-d5Aos(8-1okI
z@dA{q44;6zDl^V^FV1%_&iM3v_u>ps&+zo8fv0D9bJ`lk)f6b^_ST~w(D9du4knfx
z-VKN8yMMHBu#J9)dVm)(rX&iyElrB8@$hzRLRMuU`ZlXEy-ASR2OOmG)x_7ykPf0E
zou)B|gVYVsGwDd?ROZ@|(u);%;l*9-ZuMclVzL|RwJj)8owj+WPX^=J#|h^6@BX^n
z0UggmcdbLk_zY@VQV5mTrM4T{`n2-|>7A;8m+)z?r<WfJEX$-NR)?7OojxXBNGVd~
zxiOkVcta$6M_+!8-11zS>Y!*o7~A!XMDZol#(N<XFd<EfcdgoY&Bg}Sftl^}f>3nH
z%plp2^^|NN8K$H``00Et<Y&1S@{oq5_6NsaGQ9MO`eEP4>J^bJXrpcrwsv-LMG6Ya
zB*c8a%x)`i*f}GlX<PmhHGM}bQwk?w&FTK|%yRRR51Ry}5^c-iU`iIP+%Vc`JnAsL
z`&PI2?iNf3*i-^(L1LNxl=uybISV?TOO*v^6&4s$J#g%^_u1%3x)F_Mmi7knR&uY5
zI`@05`ox)*wgR9jy)6d8x?Ga5c%C(FlpHixghq~Sz`q;=EI^zdMKVY_#8xehxR~Je
zHE^KOcAd9NQQw~gn>5F29IrC;wNsj#qkGpHWLN%9z2JwiXWC0#;9uJ)&d@Tc%RXqm
z?8ba%A-sOr3PP}Xj80PE#^QDz{oKRMJ_{{BO8)PT%ChexVg!qrrX1oR+{}H#+Qi$C
z=ut%0b8puLcY&{%ze{j8@yWiAA|#xGUYh$S&_Ctg8gw^Ht48#)VKEX0#9ziMFdYtb
zN`zU)b^Y-i(P^>~AVkiy;1+fPFbFd$Ud~JsSrUb<@yO^*@AX;#zQbVjeg^-WgaL-_
zIIC|~N$a%_`~KmxBJI;5(661EFH=#8cxXKtc~>}x5Smy~WVp+<mEJ6Z2+1Jsevt6&
z54}MoUvtm?;+7W_Rb_=EH;54Z(bJ0d=n_I|4&lmO(vpwuEJW{Rn?!c#(gO?amFu_X
zP#RzDPn~5`*<@<hvzFC)_>uu*WKE*VY`wxbfzbQ&(;CsfJz(DLWz!qaPhV_(>Y}}S
zqcXewH)hvx?+s6wZ`pPN%PMmRC(K2phKCsYQ+2&{3SaQqK(}3|qVdV(gvrqI#RPIc
zf{X*vTkTVpu-hZOe*Wa{6bS-vH-UDr!g_LoyU@U#tgHxud1d#l-U$m}!KFQPuQQ!g
zb$`nI71Woc3|6?0PS|;CbiY6A?i8N-0vdbEhpITVHr<f?Sy`iV!Wx}z7yuG|ev#HG
zf<*4~E3|trCJV4zR!yt9;{F8fIhdh}u6uul=_vslQZzoaq6ZRIrq5289+x0}#rpb-
zo3{b`ff+rh>w)do?4WH_BcL`@L-QBoDby@5ps}+nZrC$U$vuEd;+3_gfv`V+I6Yzh
z+O@-AI?^jfEOM8AME)dr82ZBqK%4n*KA>fSLe3D&UNQP$mUbPNND{B;nG|7y&><mX
zMGT(R5U0YvTVdbz0xPa%j^2wWue&eic~Wtm0T}kMa{xjnSY!<qvcaB&IE51Sm#3$Y
zX*XJ7RI|B%zNP7rS%krw`QUnX@D%tA*!qEqwreJab{4toUd-G_4F=i-1_254XLin6
z@yFX2eVRrjGbBhAt=vszfs_1M-#|rrasfXjNBIo>A;t~#mh0p}1YN<7Sy`V6GbhAn
zvY%yCqvu6wKfo;L_jjm-)V$chl}@Lf|EVlWyy#`UchlY!*QL0c!BZp4#omg^P6n9#
zXkecOG^nR&$%umuvq@`iwt?)Fu>yUJmqngnI|lB8{TQ=(3u|n#kwZDVZnPrde@02q
zoOn$ffK_7X2ABw6hdRvMwJQ54eshFA__738zcHX72y=V5-~tG;834lT0Qo)!OhlF;
zSz34;0Zx!LhWpI>{yZcE4<=dxSu9@WUQB>QMAAd!7^57vOMTHv!U*Ncg0pn=7vcMW
zw~5m*c(hH=lPpjwV=XN_S)?=XtkBuL{Qz(a#94wyH=x6s)IS__PD`H9%(!ca<wBD3
zqReYp=nYM2$lWSs+K71#n_zr`h!xsrZ&zeE*foCETW`KJ9pXMt$in}hT`aPceZoe%
zU%rudFZ<ARXfVU1ab%KLMZ4-E?CI5^vJ&%85(Y^&sV@4r>p}PBdhj`|2Q&ev-SBgD
zVbtl-p5B~I_5R2-`m@f9sdu(=qwLKE^9V}cMi45HR_0sRbEr{#8sw1VA;To8Ga^*J
zIf#@_t@b-(VE*+bNz!bMBun(vdQsk0vv`d?(qk*8INiI-tUcc)=gbvQI?*d=tk!`>
zIk3Ptwh#jprZNfGWwyQxI4erO5~mazX=uGp&irtr`BDCxiX`t(Xb)zGd|%TOkLqIJ
zWY9i-3tY43$9cIttX9Ts`l6>~<^tVWc=aHR_cgn};{7}P{qf7283aovaG-$DTSwHt
z32NQ~Z`<Z?|LD89ZIq6s=d3;^EOb!nmxvmzZG<fuvv!~G$dKNash|OS>uo|eKE-@1
zQhFN#Rm}w#c}#QaV-SLt)0f~VQ<_S9wIW<5sXa-}c%DwzOTLz4FN1RU^0CrT?Z1lo
zCdsY<<SVGtjGY|CQcPTN7f5)&r<Z+R?v9W+-=&>g_{owCel2UElFzI#9xqyeB=Aku
z39LS>2wvkz)r_z3samt!1z-TcBy`_aHf|Tt)J_@gr?l>PI@D)LcT+t?d%mG_?F6?+
zdEJIo!XQ+%K;xx`2Zxf<TdKPyo?p^5Ft3w-fnvzvJB^nl;FWxx)4~RBz*r&49q&z*
zeXm;hYuRU~7<+ZUyvpo^c*_>4PM+k)hrMG-A73hX?q8VD_t1P9q@peI))KLfLl!i1
zq=vX}b(si&mkR*6fTi93VAh#ZSq8x)Zn2ue9#d?vtI|UjoX(O9ll0Oa-u+|a0RA-D
zO125NU=@7%<tA?ivJ0LJy}dw1PwxHEjDjrzEY1u9=bCxTHAICIncp<#iTRCr%Qggo
z7f+C-39@i3zH|MruO$ugmm6$T_*?A1qAn=&Tos9Cj>rNH`<EF7JN_Cm#D>c8awZ4c
z;CcCFF}}zib%LM@gr=y7jV7;3e}n&5(?h{Mw@IR!sxY%4xO0~)`uYT#rn+OQqUV8l
zrGyAZ)>!87Rv07w$1{ueg3;ist%)}&o<>1T3!qauN5~+S>(iIfPwAlwu8D10a=cZ$
zVn0d&R<$$7RcM{0z^hy+_m<GZ7}J((c`f`wVm)?S(E0|#!OpCF0z-eo7~l!h!MYK>
z?FmWKwO?4$%?Y6dEwki__0?S_&}$VAFFpEVRpVqy;SR^SY@8RHyu`|U4Gmw=2wbx%
zD8eC)kn?G(EDwb2A&pqjaGWIaihM{T<uoeW<Ro6MSx3%k6b=?*W#y1Y$!Vk}C#bSq
zBeR;*h_WDyAod|mv!LO4P5=ee=#b5c1W6Jg{AvW^xDpZGEaqx<&+!G3ph}$Dlxm0y
zU30KXQ&rR&qC(lMB8W{%tWhLaXq)9FffJfF8sKtuH(5av6{$v(T&`}F<C+rB)iKM}
zt*ES`@M4|Rxw^q>yxLT2c=5TmWdWp?_!?e(u51x9U$a(?e6DMW6IfA_YwYol!iywc
zVWm26{8)$#kx{!0XD&TRbl~=8NUuW|YL`W!PM_HXzA^#yFE(p9h}i_8%0X=^YV*hh
zu{xm%(NHr%$tFk#9KcXDBB|K~xiUdefyEJerAn`GN{#)R8NFJeZ?dd<jJ{dj0ESSF
zK<NaYWe;|&KqMTQP&Kn4bG6>iWfQ7q7Fmw14YO=Q)yx98MIDEHLe<PFC)5UMKH<R3
ziYTc^Syb&TLUFIlmTESkYG{#{>+(;o)&oQHa-C}z6Amm5S*6am`OFYp)zm!xZ#0J!
z-uu~vBCE@CZR{7MMQ&GRRjN^CL0V*ZS>%Npi3`#q&kKse9bL7^_6jeFH4a>mmP#zo
zH*12UAXQ4FoOo2$qTtFLS63tp(xUDvO|H(hbJC(7z^Z^fj*)Vud@eMP3Bd1JORn+b
zg0!%f3U%a|UU{*!ExB0};st4LEpRff4ssQwg|z@#9Mx4}EM67rl3qbt7>kEmT$3vb
z(!y9x09!{Wh1^!6AUBUGRVaX|0`Tn^C0}GsSvqEZAvcsLHN~Uc#bQmV4)PV0xtRn}
zJgQJZE)ru?I7+D$a;kOUwxG<7Bmgcy7JEW&BLJ5CG43K;mh+f3#R6qPayh1Xv9OUK
z@kiGzY$OOrZ6u1Nz4OQRE*3U|Bz<(v+(v>59m+9ErC3v}Ya9y7l8uz3o-7tNQYD(m
zT6fd&@M1DHF8b|@X{URk4`vtd-s^*lq+`6mliQ2r=<Z_g*OwRjKVrPt|J~#-SSBx2
zg=0njB8(txCCCN>chS0N=@+nbuEFsa6WuT_x;+j0%rNk8x@~RiKls6;f72UkqoD~l
z7gf2z-}DBWp%2>}chj2C9~7`Rtsee73=sNFYouWbxfFN>{$T}~zvLjLz)0JeV0Ut(
zw>17n?`n98Cfw*FO};7LJJG}jE2}r{+nYAf^Wsg{(71-c^EdRZwU+$NrVrm~GKc_w
zSq=Zq7J17y09xL>YeSN1u)@uBqHEkGe=~ir={Nf?gpgJQKE`Z%le|^IO)1Tr^i_t|
zHExOndF_4g@~(B!n_cVe>-OZj_wIT!xL~iiYsmIIcP&cS&E~ZVp~_yWYLkvY(d3%f
zf{Y^yI5K-FvjU8;*AOZ2kAMU0B?p-udR*m~{{wJOSzUv1001A02m}BC000301^_}s
R0stET0{{R300000008>47%2b%

diff --git a/public/biopet-tools/src/test/resources/unvep_online.vcf.gz.tbi b/public/biopet-tools/src/test/resources/unvep_online.vcf.gz.tbi
index bb43ff545f591dd276973e7158919e6d14c78f23..c52458d5601b25642ad72a90934955f280428dea 100644
GIT binary patch
delta 30
lcmYdIouDt46BQ82#wsU}dCY*TjjQ3PEJH}Ef#XD11ptas2;%?%

delta 30
mcmYdIouDrk{%PSRo^weGmwFOKctjbxtr@iEJh?T|RRI9If(y+6

diff --git a/public/biopet-tools/src/test/resources/unvepped.vcf b/public/biopet-tools/src/test/resources/unvepped.vcf
index 2ea7ff21e..95c2f1e17 100644
--- a/public/biopet-tools/src/test/resources/unvepped.vcf
+++ b/public/biopet-tools/src/test/resources/unvepped.vcf
@@ -1,11 +1,5 @@
 ##fileformat=VCFv4.1
 ##reference=file:///data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta
-##UnifiedGenotyperCommandLine=<ID=ApplyRecalibration,Version=3.1-1-g07a4bf8,Date="Wed Jul 09 06:04:23 CEST 2014",Epoch=1404878663625,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=1000 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false input=[(RodBinding name=input source=/data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/Child_7006504.ug.chrom_merged.vcf)] recal_file=(RodBinding name=recal_file source=/data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.snp.recal) tranches_file=/data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.snp.tranches out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub ts_filter_level=99.0 lodCutoff=null ignore_filter=null excludeFiltered=false mode=SNP filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
-##UnifiedGenotyperCommandLine=<ID=ApplyRecalibration,Version=3.1-1-g07a4bf8,Date="Wed Jul 09 06:07:15 CEST 2014",Epoch=1404878835287,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=1000 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false input=[(RodBinding name=input source=/data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.snp.recalibrated.vcf)] recal_file=(RodBinding name=recal_file source=/data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.indel.recal) tranches_file=/data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_ug/Child_7006504.indel.tranches out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub ts_filter_level=99.0 lodCutoff=null ignore_filter=null excludeFiltered=false mode=INDEL filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
-##UnifiedGenotyperCommandLine=<ID=UnifiedGenotyper,CommandLineOptions="analysis_type=UnifiedGenotyper input_file=[/data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Child_7006504.ready.bam, /data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Mother_7006508.ready.bam, /data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Father_7006506.ready.bam] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=[chrM] excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=250 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false genotype_likelihoods_model=BOTH pcr_error_rate=1.0E-4 computeSLOD=false annotateNDA=false pair_hmm_implementation=LOGLESS_CACHING min_base_quality_score=17 max_deletion_fraction=0.05 allSitePLs=false min_indel_count_for_genotyping=5 min_indel_fraction_per_sample=0.25 indelGapContinuationPenalty=10 indelGapOpenPenalty=45 indelHaplotypeSize=80 indelDebug=false ignoreSNPAlleles=false allReadsSP=false ignoreLaneInfo=false reference_sample_calls=(RodBinding name= source=UNBOUND) reference_sample_name=null sample_ploidy=2 min_quality_score=1 max_quality_score=40 site_quality_prior=20 min_power_threshold_for_calling=0.95 min_reference_depth=100 exclude_filtered_reference_sites=false output_mode=EMIT_VARIANTS_ONLY heterozygosity=0.001 indel_heterozygosity=1.25E-4 genotyping_mode=DISCOVERY standard_min_confidence_threshold_for_calling=20.0 standard_min_confidence_threshold_for_emitting=20.0 alleles=(RodBinding name= source=UNBOUND) max_alternate_alleles=6 input_prior=[] contamination_fraction_to_filter=0.0 contamination_fraction_per_sample_file=null p_nonref_model=EXACT_INDEPENDENT exactcallslog=null dbsnp=(RodBinding name=dbsnp source=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/dbsnp_137.hg19_nohap.vcf) comp=[] out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub onlyEmitSamples=[] debug_file=null metrics_file=null annotation=[] excludeAnnotation=[] filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Wed Jul 09 04:32:56 CEST 2014",Epoch=1404873176100,Version=3.1-1-g07a4bf8>
-##HaplotypeCallerCommandLine=<ID=ApplyRecalibration,Version=3.1-1-g07a4bf8,Date="Wed Jul 09 10:12:42 CEST 2014",Epoch=1404893562589,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=1000 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false input=[(RodBinding name=input source=/data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/Child_7006504.hc.chrom_merged.vcf)] recal_file=(RodBinding name=recal_file source=/data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.snp.recal) tranches_file=/data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.snp.tranches out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub ts_filter_level=99.0 lodCutoff=null ignore_filter=null excludeFiltered=false mode=SNP filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
-##HaplotypeCallerCommandLine=<ID=ApplyRecalibration,Version=3.1-1-g07a4bf8,Date="Wed Jul 09 10:15:50 CEST 2014",Epoch=1404893750488,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=1000 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false input=[(RodBinding name=input source=/data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.snp.recalibrated.vcf)] recal_file=(RodBinding name=recal_file source=/data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.indel.recal) tranches_file=/data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/recalibration_hc/Child_7006504.indel.tranches out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub ts_filter_level=99.0 lodCutoff=null ignore_filter=null excludeFiltered=false mode=INDEL filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
-##HaplotypeCallerCommandLine=<ID=HaplotypeCaller,CommandLineOptions="analysis_type=HaplotypeCaller input_file=[/data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Child_7006504.ready.bam, /data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Mother_7006508.ready.bam, /data/DIV5/KG/kg_wes_mr_devel/runs/trio_7006504_run_00/trio_7006504/phase2/bams/Father_7006506.ready.bam] showFullBamList=false read_buffer_size=null phone_home=AWS gatk_key=null tag=NA read_filter=[] intervals=[chrM] excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/ucsc.hg19_nohap.fasta nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=250 baq=OFF baqGapOpenPenalty=40.0 fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 disable_indel_quals=false emit_original_quals=false preserve_qscores_less_than=6 globalQScorePrior=-1.0 validation_strictness=SILENT remove_program_records=false keep_program_records=false sample_rename_mapping_file=null unsafe=null disable_auto_index_creation_and_locking_when_reading_rods=false num_threads=1 num_cpu_threads_per_data_thread=1 num_io_threads=0 monitorThreadEfficiency=false num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false generateShadowBCF=false variant_index_type=DYNAMIC_SEEK variant_index_parameter=-1 logging_level=INFO log_to_file=null help=false version=false out=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub no_cmdline_in_header=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub bcf=org.broadinstitute.sting.gatk.io.stubs.VariantContextWriterStub likelihoodCalculationEngine=PairHMM heterogeneousKmerSizeResolution=COMBO_MIN graphOutput=null bamOutput=null bam_compression=null disable_bam_indexing=null generate_md5=null simplifyBAM=null bamWriterType=CALLED_HAPLOTYPES dbsnp=(RodBinding name=dbsnp source=/data/DIV5/KG/references/gatk_bundle_2.5/hg19_nohap/dbsnp_137.hg19_nohap.vcf) dontTrimActiveRegions=false maxDiscARExtension=25 maxGGAARExtension=300 paddingAroundIndels=150 paddingAroundSNPs=20 comp=[] annotation=[ClippingRankSumTest, DepthPerSampleHC] excludeAnnotation=[SpanningDeletions, TandemRepeatAnnotator] heterozygosity=0.001 indel_heterozygosity=1.25E-4 genotyping_mode=DISCOVERY standard_min_confidence_threshold_for_calling=20.0 standard_min_confidence_threshold_for_emitting=20.0 alleles=(RodBinding name= source=UNBOUND) max_alternate_alleles=6 input_prior=[] contamination_fraction_to_filter=0.0 contamination_fraction_per_sample_file=null p_nonref_model=EXACT_INDEPENDENT exactcallslog=null kmerSize=[10, 25] dontIncreaseKmerSizesForCycles=false numPruningSamples=1 recoverDanglingHeads=false dontRecoverDanglingTails=false consensus=false emitRefConfidence=NONE GVCFGQBands=[5, 20, 60] indelSizeToEliminateInRefModel=10 min_base_quality_score=10 minPruning=2 gcpHMM=10 includeUmappedReads=false useAllelesTrigger=false useFilteredReadsForAnnotations=false phredScaledGlobalReadMismappingRate=45 maxNumHaplotypesInPopulation=128 mergeVariantsViaLD=false pair_hmm_implementation=LOGLESS_CACHING keepRG=null justDetermineActiveRegions=false dontGenotype=false errorCorrectKmers=false debug=false debugGraphTransformations=false dontUseSoftClippedBases=false captureAssemblyFailureBAM=false allowCyclesInKmerGraphToGeneratePaths=false noFpga=false errorCorrectReads=false kmerLengthForReadErrorCorrection=25 minObservationsForKmerToBeSolid=20 pcr_indel_model=CONSERVATIVE activityProfileOut=null activeRegionOut=null activeRegionIn=null activeRegionExtension=null forceActive=false activeRegionMaxSize=null bandPassSigma=null min_mapping_quality_score=20 filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Wed Jul 09 04:32:56 CEST 2014",Epoch=1404873176271,Version=3.1-1-g07a4bf8>
 ##INFO=<ID=AC,Number=A,Type=Integer,Description="Allele count in genotypes, for each ALT allele, in the same order as listed">
 ##INFO=<ID=AF,Number=A,Type=Float,Description="Allele Frequency, for each ALT allele, in the same order as listed">
 ##INFO=<ID=AN,Number=1,Type=Integer,Description="Total number of alleles in called genotypes">
@@ -130,5 +124,5 @@
 ##contig=<ID=chrUn_gl000247,length=36422>
 ##contig=<ID=chrUn_gl000248,length=39786>
 ##contig=<ID=chrUn_gl000249,length=38502>
-#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	Child_7006504	Father_7006506	Mother_7006508
+#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	Sample_101	Sample_102	Sample_103
 chr1	14599	.	T	A	35.13	VQSRTrancheSNP99.00to99.90	AC=4;AF=1.0;AN=4;BaseQRankSum=-0.736;DP=3;Dels=0.0;FS=0.0;HaplotypeScore=0.4161;MLEAC=4;MLEAF=1.0;MQ=20.77;MQ0=0;MQRankSum=0.736;QD=11.71;ReadPosRankSum=-0.736;VQSLOD=0.076;culprit=FS;GATKCaller=UG	GT:AD:DP:GQ:PL	./.:.:.:.:.	1/1:0,1:1:3:32,3,0	1/1:1,1:2:3:27,3,0
diff --git a/public/biopet-tools/src/test/resources/unvepped.vcf.gz b/public/biopet-tools/src/test/resources/unvepped.vcf.gz
index 836a5646314f87f60c541fd56d3fced47ab53294..f8d6fe3f8bfb1a1a4a4512eef76c53a8fc39c49a 100644
GIT binary patch
literal 2514
zcmV;@2`%;?iwFb&00000{{{d;LjnM`3AI^mbJ|D}{!D(w)aJt_wKiy8^oGM-vBo&m
zVjBZ-?)HmoKm({EBvul3e82vBB%}c`E0c4!{K1mw>F(*L`{^F({riO<h(#2ydD?ln
z>_4$4e)s-;EEXab;aqg&&+S&L<?)oadV>eE_4{?J{5feYdHU^fwhg^NJd&o_S}n2r
z7)C3;X>I4pys7S+3!bDLln-wEcb!jzUZ;C`cDr5AMBM3~O<y*mGYC_$6!BS4B=gwc
zq<$23e(45*2n3o(+b~6bh?XLZQjkv05crJ*pRZ7NI7K{@JCh>mN}z<V1&U%%#E2&-
z@RL+{zx?ySe7|D8A4ELW8t=#A&n;~D<;^y43!AuR%hM?30SYrg6fJgU5@}{GuRRrn
zL-oJ$MBI;g_-(RXR}5)3_s4lMk79uqakNHX{a_w_k3tl~a+GY>C>2R6t#^YIJtR#u
z77H{31L)5!5B$`x8t?s9m^2ELe0gZmn@w&<=tfAXWaV!v3Pz^{r@Pt2(RUvbR-jmL
z4|!sft}aj#t@8&7TJsmQ;!gt2ga{F&a0=;KL83Q#)!VPI3i0qo24M&9M#A7l&tZP;
z{i~u>1j*}h@5g*DucY2_=I~_{rrZzxaET@x9%dJyOGv3ol&@}kr$lLzMjNz=62CAc
zH{gS0{KSJpjs3~%E{|3aKu8xZlW%hlI<^9B6X|E)PgWxSM}ow6*jAo3X#}XG0&bHO
z^2tLppI7$rnQwwDg(orxU-tk)l*odDA}sOnAqAfObcJ||){(UP9HB%kV*#`zeOnbq
zD#i!lEEd88k;_Od7Oz`Z3xJ4Vg+KQN&=FKQY&XkAFE1nMS5ch60N+I5uc8Pb-iY`-
zPa5>Rf~=fJ>kWsXBEAIVrSGuvGmpKO^RL1W!JQ>UbOSNu(uG6+Z-!T=k>?GEN`YEd
zp6n|{A3)K^jMAiazQ?o)`T(0ip#YLkc<cu+is4|_^16HZ@vU?y_i+LX|3AA}6qJ3!
zM*45x$o=a%G`^Xj8@}1dB)ZQi^rfoAIK4YmG@(C49OZ0MRr<sIp#OS5c-r-YjDXiW
zSN@_e++JOGr-O&9$LY8`xP6=qrdO|0?{wuS=n3HBP-V{$^lbvA0<<#Q`hk~$5@tmX
zBo7ReQ)ePF|7PIQZfV7n(H$87`JUu7tMOzhGga@(eU-&Ecs%T^0`)C?PtIGd0`1ug
z^ap5rhf2qzZlSeVlR&rSU;)$x{VWBB7HhzABY;S^p$~Ibvbz=DUTN$~>w_Q8gKb8J
z<s|sA;>nlF_9kfz@|DmY+x;L1jHIdzl57+uf5%*NkjG`aJnmML>G;)5Dd2+Gag?m0
zI6aK)p~i>%$?)!kBp60Z6nW4<Vo7#cehV3S4%r4hxX8>W7^t?%?zaUGU`&IdT;!CP
z<)xe&`8G*F9IgHD=@ufijyzGB2Ikvf6Z`2Y@JX@_fn0MIeqQ;IBF{TceU4&?GBiDT
zwvayUsPtOv!dXu3IW+?oNv1Eq-LJi3lta<B63)xO&q(O+Bs(>bub@w(xiksBWGQ>-
z4G`Y%m$F*oN<Hf0D22;iX7hbql8n6Fk3wdqe>cAAPV@6b@05biPBH#GpOllbQ$q0~
zIh=9W<K7WN{**Vny-v4h^DSdK`jAiP+0=G(-yx<6w({Fq)=D7oh|A*LLH)66{<_*~
zwdSvj8En7P?rWKzkl}@r>f+I`SbSYav+`6O`2Up@dN4olkA)J8@-Qc-3r{SvTHOHw
zXL6LlsN=)PsRqIRU^u-RXTA)hXSu~HX*}<`@k>=6cxOC?G@q};;I?-)bls+5q!Ijc
zjY@sTTR<fdVb<Mb*GcnrC=gE~$Xw4<<TWgZa8Qnhk=I7&q-oohg+6uAck~IJ8%@J8
ztrL84!53_Rxbf@H@<!Sx4sKeu^J5<T%y8`NmjTJ+Q{xUCl-~x8JZWn@F`JkgcHxK7
z#1?j5b>mH_=p2zDppNUZCN@p;2kif}I4J8SCV6NwXXbDq_m_D?T*WtMfe4rBs)H@l
z#>c<o;&;MaVj33F$Pgt%vFlQs9?Q_uj_Vi%1dU@ErerWS8L_Qn8MBme3>TZk(&({D
z85@Hl!?KTM>{7;bF?B3U<FixBFpF9Y+&-3Zl?*1BLI5-|G)f_Anv~g0;}FJsLWW_P
zShIUfln`}H?6{^zNW~h17T2*IjgShqVN>RsOcO|~qHPe9V(MxlfXi^Z1~r*&Y9fit
za61^gCc!#mWw>p}ux-M0UYFqp)QICc8e61{wnf2flW1&_GT6+)#MOq8l(9CUc9>~t
z@c6gNB9qvLsmsRyD3_r!v}*9<GB6AS(*9AtQi|3p>=vWCcoqvtWdX#Wxf%<sSU{^5
z9M^W-ngy)7z=dkiEU=3O=8**eO5>7KEU+pIfW9!Rk+-YzHnuhRxrMw_k#`Mf*fsKQ
zbpx`5G!E?+5W_gyu?>|_v!DtywXinZaj~EZGc&NEEwf@l6=n))p|emfsKV@Esx8uT
z!4b?hGo4x!Rh${LJsn#*#eyo(j95DUDZ_dMG_iE4tri?14XmQeZBhV%t58Faax^}D
z35x}a)h$a~`$|-CyJI<~CX`B4F}%eH)p)E#70*-K#<f}%+uOusngmv&do7j_S5q8H
z)ZP=ttkzf6a4n2=O;U-fzH7T!m)fPM>IZTjH0sb&Z2|b4x^)ha1GFqn9xG9WmIHmH
zF0O>_ZOd{sC9XtEv>?f_Zsby;3N4r{YGb9qA`aD&o)T4H5$MGlu27;1EEk5Y8cAAW
z#i-@h1(mA7)Q0h`PLe3rw9GpAX$cfFT~=$IDV?Tn<WrI*Oq8+OfKp2lqf2WgZRMw<
z8*Y_k2@-|5ysq|Wi4hEzq|Q1s_Wi6wlc^~S$fYjwOkqS#Qma#8L}@KXjP326)a{)q
zjDVzTb(R=W2Tmw;lD5*xbO#3|xrdQmo5@UJ<d_+eU49<l-84pblg6Oe7+>`p-Ql!x
z|D`)@<VVhid~MgrU+**~*(bV3Y+!Xs)Fpi<dt(Db?z@d<W7_C8s5~|{{`%zJ=w5c%
zMYrF<P2-|_3)lN!={x5})27x%Z`7d|^7C;APPrHTN%rsX8!0?uz|I%h$3IY$eLTyw
z$`@&%#)b={BX0}kx$gTO9CDgAzBqi5UYO3`L`#FVb+LQ>+v!g()Jx3Hm+Qv$wB7Bs
zd!zRCeS0))G+WK~?yrGcxNV%_Hg3~45D7gqvPT#m5qJcA36I|W58^O~tZO0w03VA8
c1ONa4009360763o02=@U00000000000KE>*`v3p{

literal 5094
zcmV<C6B+CuiwFb&00000{{{d;LjnNF6YX4SbK5wR{;c^GsLY2=YHgYK;f!`m*5PPJ
zmnGSG^G>B$ln6;&QzVC^Y-@h~c7vxVT6Q+)2h~(fM6dxg`gr;Qh~)R*m$pNfUa&Ev
z#$BuXAl3Qbe*Zn7OBzsjK^yq`^5WvcG9&Y%-Mf=7{=B}(UWXSeGrBkCTi0@^A=KrI
z^@`UG*IS$Z#dZ-c>iKhhX@-#rgWtOL(x%omb-ifk)1c*THl}O!ZI?Fw=(QV7-*<Ks
zx-cDk9+;8sxo3AY2;t9$Sm)3A^A)F>(tN3%wath&-h8AM`Jb&rIGu3HB_~}9B58G|
zGa_)j^yaMNdy930mpDmNHAN8>K|U)?Ir4Fau<^!pO=lO{p@B+OC_!x3-$n)+r195h
z60W_+?$&Xx%uU}8qsG#7LP`Q^T40VP*vAO%e`v!65#p~sml|tt1LJ<25@xe|x{F^%
z=BhDl#yw&1hz2;k?ZS@-7&8{4PYY*j(cXS*DrAH-GJ-8Ac025ihI#vj1}ocLHJbf?
z-sqc_g^e7^{0Ia5y=DIXW{xJV=e8+Kxv^cakG+_H!M)9tQY*2QWryb6p=~=_(*PBU
zn>XguWa~yY80(xbJVDq*VPnu6-p)ExVtJ2lXl{H5I>sKm`eaO-gK@u;zoSOv8Oy+2
zp#Nf%?Y~=i4>T}Ww83*6N9N|=jZwFYe_osZ$fxd@LQq7zhQ!r5vb3L!jU9qHJd0Y!
z-&=@{XlH~AFQ8$fyy-aJqv3lI1)bqx_hOH2NJoJUUNW5<_TJ83S#)V`o#@I8Y5odq
z0y)Dn7W~rpR`I&7Zl;s?7t8{4|3leyA$U2lvthdyb&}S20=lsy!%Jr0d*g$G5G^oB
ztgPYC5He$Jx($V_9B&TVPVvn+u)P430{cL0*%nI`Ba8xj5xLN(G41s`!x@Cy25fEk
zfwv0GjS*1jY^Af^Q|dqSD)AKrp}C>P#`Jx3I7>pd+AVl<nfy){YHq=`xY{RU0cpiF
zg5Wb8Z*h;kAJ^11FyHVY@Upqw?FN>JFTw`To)-Q#Z85;=1|~<+kVx2G)`cS*4`S2{
zW~`yJT-pm8e7`FUKZgaP@&;HEGAYSpmiw!~+xo`bF&FoSgNPDC2es@~KpEtuzo!vP
zmqN>3X8IYX!gMz2wPr~$W=LW@LZ%s8mstkF5p+p3TTvI%AfnT?3Hfo=>L%?EW?(}w
zCAMT)*8VhX4tgzP+UfjRdheUyKFoNGQO8>$@<7}?P^Zxwc1PI6@-g$>nmT?uQv&eu
z6GPg@UvDR#bp>O?I0W6_VC^LIw!r~X-Gf^pOgG>oobe)nz_<WAd4|eyip)uXkFLRS
zdu<o~8gfgxXsvC>N*n9jReiAzybVK7YSkYW%Xgo_M^Fo5<q!1A-+c$=f?^RewyTo1
zq3hR~$U71R010cDh-vv_RY>P19(c|RR`q${VS)kjN85<j;ZJu}$0Aa<J$TyA!}?tu
zJ}nO-^eOrn*bqBYcmt)wSZpkSb_6YB4cP%f^|gbB@GCSt*V%pbQ1iv|tH+4Kq}nGb
zqU!*mj%T$1yq3$PzOP(>*pw9%i9nY%Wu<fj^`$W#j!80Pj94t+hsMHQnL#3$yoaI$
zgm1CN44s3E89Wy8Kas<CAyoAU5Z%kBd-Gr4pT(({y!<SS(?nU&)Dw$4vA7e9JF&PE
zi#xHnFCj`y<d^4lEYE*7vx}8Dal4<+?Rvv@r~eDFyYl;2aItciJtp=ILmB`vEQK)|
zY*&XDXXIN?H1Lr3B(JpZIaAjx%%tQq{Xi1?7wn4%pHDRQJJVP}KC##ni#@T}UxLN1
zl6nPpxTlW2_P|F5)?}w~HJW`O{vt4F5O`3Q@f-t!uX8Sug|~s?MyLHzJIzKH#sP_8
zyP3wdX$Qu71AfD>+CY{vt2O$g>wafCHCoNqhu-j-0A|Hyd@pAkyoz`P<;0;ZsxoK+
ze2~%JQ{W@xei~3XN}P2I4=~JttTmFgp#$ZDKspi&vpDEsP(dc_HNX)v9^2hAbyHa4
zx$N~GS|z<cn7+g8IK@(=CBkhw-=<I-6Qb#G+;kl36j5w~Aw3=M$@NVa3eeKaid|Mb
z;v@xH3jlJ!=u;klJG>g*4%_bzJ2EjA5=jH-ZCksBzzkG!0rS9qo5V3Tnzi~2lmw1-
z_q|6jDJCh@Y#=D;1k{vOz`7ngAy>(weze9@jkM$@RT&bdpf_^g#0)^tkyA0s*ctR@
z#$9vLYYt~qV>Il4f*^;^-aosQ2c32pC>$R<!YIAs!QB|gd4$A+wtLgo=&m#QM4+m=
zmKj(Es=V;rrEM{Fsuw3<(Ry7MECf;3$4r9ZWly3jOtg**S!SgDl!Ru+ZUC5qc_R}P
zF!z!s0RsdKob?%@JRb_k;RDfd!5s-GbpF+B%?#x0;|}}}$_;&j?kpaiv_e{Q*o_{p
zi8WVmHD5=5*Z>AEsznBb`&<S{xct)=?sa3JA{~gtDRXy-rn3;N3XX^R0?Z;Tv?Lgh
zApBQGlbd^uKV}P8ai%U^io&I=Ja?up@~Q$r^UN7NZV594Z((kI%kzGoyX1vSNq9y#
zU6d6;*7VbP|7oLiI`2R2lukRPUvsDQb8l^|7r*eWjrHQE-`c2D;dHY9(@*w)?5&OO
zOXB29nR}MRsStFUb|P^n5_ckTClYreaVHY@B|PbuytHwk?$5Zi@r)8DYWLHr-H*Ps
z@dCLky?zA~EA=?t*f`zTIC0q1jg1q9JyF<Sn8Kb2&+Ep>FSI~;Uu+!~0UdXRd|_<b
z!H2;hxiNqq_O{`lP!I6>#)O8Rvt^*z8V#;S2B0c|?DN+ncK;yu4+Kc*(ZG+;00&_X
zr%BG?OBMzgne0s(RK~`Vlj|6GO~zjCu9|~v#CSC<YFkjGI&I@abKD=zK8+FMfA~js
z3v4_K>`fCY#sipX#Yk0Fm)dr?Xihp$0PoBMf{d-#*Ue(1$Z;eou}#Qn*Xkkg0;I^L
z_UUjO;=Pje%YO0WbIY+A)Irg_pKdoZ8b)WNjkih0;6et9A6k`PsHZ-R13lZxwW9Eh
z%m8fYW<q@c!$i~zKA$d={3@48?vu1^{Ka!;M36tTIPAGty+WEMZP@jK)^3qstb&3v
z4iL{*=^Y0i_ti*g+NQfgOMhUMDWMZE=4Ai%%(SzT4~qn*5^js+U_zIz%rTAOXxJgw
zcdhRA%@w!~xTy@<g2oE>IgT4tbLMp%n;8q(Dhx1ScHp^ZZ^5xQ$wC?eSuK3Pt@t(?
z3+`L2`qY{fmI9<Hxnl;&x?0g7Z%?;zOdKqhr$)wYAik_Acz{~hjLD$r05@+*#zhn_
zy+HsCx0`GO753b*=O-B3;02Aa&r}IEhj+HwPcI6-_F5p)p2>9_0{^=Wqm(VvD)d3`
zr8n<W58<W7mKT7}BMg$9H|BTp*zapZ_E})MVf^oel!m<x>C{_B405Q2bTf+y(~q_R
zZH6J;%$;2q!UY}?zYBCXj>(>jD#SB+*9rEIVgHmxtJn2cW`*iS$6_Y*sk@5S;5vNK
zDF~*BYrCU4veS4WAc!;1yen7*&>%!qywVv1S)7Hf(Qw+C+%;#tyAHwZeSr8I2Og&F
zD6Ma1!RnQbJ@>FJhkf=6?7K+Zmzt;`JTRZ8*{*O7DKs|2aB8nMW^x}1G9;zA`$-~j
zzx2H$LFJF!J3{ssf1GNaK}yXbU7KgD<l}Z0viGb_!$n}TBaFtIII{`6<O&#&-QR^m
z2IeKyT&zj6KTB%C8H&IL5_bx~90Q<v#0&{CXllLLpT1j~zUz2qR5^J!i1)+0uby+5
z&dDEH@(du#K|BXR;w{rMhtpuC0B2-I*dLzl4HJ$~ny6xN|2}8YM3#ZKAs=mL;XQdS
z0bIBz;dVoS+eSEh(}~DkSSJA91jB^<o$fyOCvRV+Q*!B)=4yXNYmQ$>2Nb2kwe6P%
z7{A6iO&_*~Pc}q#e8583rP}vU@i%cp2w1%%Yx9AUIi)U1uh^C*o4%;+N9YACY==fj
z+>2$kaAsircWF={b}!q#n=l3(m!Ta&E^pdc@kq+oNJiN9rvTdU8Y@jB&nG_8wPYe(
zK)UzaAF(@qS>-VnNeCCLlEgoUF}7G<QO~X&u4(WeVa#OWsEIuAN(2}UtPbMIB6Hg;
zts}objxqSM24JEP6AF?r^9Sq|2tpq*BiL-?ypJ9UDQtm(?W1n%tX%I;-gD=$O-8!d
z^5}B;^6YsB03w=f<s-~;ST2i2D^4TSD^1SAhAhNWn)5hl2p%{vmT!VliMcXyzsGuh
z0RFo>|43~Y)LDT=eXt>`Yy0{S`kmK-CzKd_1G$`QQXG`Lg@N8NNW<C90@8-$Ei8h#
z6!uObkQ!%~BOI(6e>Yq2z7`H~9j`#)|IaFxDP+H5A>FTE$j!?!G#O0EATH-5KET_R
z2jN~{9ZDmSf5QeXMw9ZOUtABmFV};Iq#iH;PB_YQbz#`KW=zYNO`5%-G40JdFGBBZ
zZHMW35_trp1C_TBR_0sVv6xd_mgE5PfG{z1hLq_y2bGepB~8X7nEcyo62okTCM%58
zYE^E^C|==@WHUkXJ=aC*?fEX={qGMt88_2W>qi1;m=Oy+M~56>=K<hLJP;tT@dBO|
zCXduAq(+ijZ{wX9E;Kt}c~{cpwzRzXkg{aj!+GB_9$Fp@f^6)C--Fle#c@_H53AL5
zHhD2pQg?yvO#C1eMEjQAAMx%6|1R0fi|KnS;#p8Y*w<ke-xMuxLbPqOGg)?Gx($<e
z$?n%hz`_7!afz(a+=d~D!_MLfH*)N&W{C|fwvHdz(QDkdriC*qsA@KZ$YX-3k6r*)
zPCiYbP6?D=SE_IpLwgL(c*?=Xiyw=*mr^-A65dz4bU^l9z%D1`YpBy?5pR$gCoY8x
z0N$U;N7?;+tg<5tXRBPqA7)dK!N;N%*fibIWYC<&`-%2z5}xg1>_yC!!>m$5@)9Ix
zCakfr#lTl_u`5)%!kk}2)5+VIp)=MtxPkq<40FYgb5=^g@Wc#g=IVFl>95PBwsQLP
zY-F#}7Dt)ABEr)_%ALpkwDa{bn8m{t2R8Tsp?mzfNDI*y*=7#eg+-SvR})9vwYnq(
zs$&D?%!y#P*PnGJOqag*h)XPk#$!^AcV&6Noyjb|O2V!);H?NQ^574hE9G~v1+U;C
zpdFu(q!(Za`-;O<k8eb<97YljgOk2zU6OO5uBwVc{%DXV@&`HR>YOGkuh2yYU3d=9
zp8xY(X@m0%8edmb?aMUymgY$2&kPdV=ZHG6QT{sD;7+?j6S>ZdoSN#v)p>>2UR2|&
zyi%nInn2WaU8?i4EPuiL-^>pA?3qtu(=cW5h~Hkt1#uk=&K!0K+29pf<qzNU>3c!a
z1({QX3K=1n5qVt~Rq;?p%;wQGP5?pfP)5pScwUwSRXLQAGZ~H3d09{@#!)gEm4|^i
zMLm>JGZ|UuMf|sx4`sAWMp8sY0<{li^jwA)co7UxVM8tx5@lJGRH;HCo-YXDEgZgL
z^}LV^iJHu7x?CZY+nR?_bWPPNgmPnZswnBQRAERyw>BrpA}{I{7T_~;>zpV{s$5}7
zJ~Ow*^SUhXRg-1rRy9sl1*uBw%-rBLLDRL0xrEHxiU?ZELd9G{W^7601-;UYLS}0j
zaz~Pt3V-}7k0M!6Ik_qtKjkVzW>gNt8!Lz7I0$<qU{{B7vn!HVWzV#OQ0f5om-LDm
zq_l%ro<Y-9O+V5>DtFK!8!9@eX$SdW25?k`N?O`MDRlt&;-w?<YFS?8)e8UVsk~N_
z*EuNIN96VL0x*Rt6iPY>9Cxr{6*A#Shq9YRg|CcuKJ8F;v&8XSWtycO%5D}REUIS6
zI+We4@nU6?W*rXPtV*(Wbc(W{B`AAUXsM+g%8r%<r3!yCvmQ8FP^v;Z?{MI0K$WW4
z7E(v>WmiLu(kgVib<++xsw+xm?&qR8wrh$eSC}#v%@JOa1hGQnTr|gdQC0b)qvmL@
z3bIrYz`1Cl#0r95ksP_GS`g)>qq^n=SK;}pBAJWkbywB-s?g3v^Ll`MkUQqByet5o
zi~2DIggtL56>*%4=HAkvjvTYAAQhIS=oKlRi)P+}AmghjS1y`+3z)@GTjkCYG_eZw
za?#vb0@UIPtjI-kXX!w;j!24`uOv~?j~P|WiKz<o?U<yHqoyJsbHA85N|JTy=<HHH
zrd&n&a>>k1q9h$PsHo(?sEbD>)m%@lirnUsnUh4|<;QYQ%zOmIQaCod#1-{C=1nPQ
zSpY7_EHC9g5@q4&n7NNc@u-g^sjzn8*xIGsM*z}C$IN^rYOq5&CaLCQN>v+&T(aOJ
z^=KqZxsNoN0a@$AWHhLaN7GubU7K{ewPt@-ySZ)lYw?z|hR1fb_;{zb_qP^m`xpLe
z`+w8mH!PDi9@0Rs)oZg_vnJvWviAHYzSe9tr1#A({-a{=n?v|r{9xZW=jy7cyl;;i
z;(Pp(jF-*ccc<)+{m-S)Bmq#rXWs!~81@~Fjf!V#FpLU6xCTC^gU4gvv|+zfSNZq*
zXX@$H@u4(Ms4DN1li)^o`aVC;Y}{Vgu4k9c_GNo~d3|#^?$_!U^~>bH8h^oGa%cP{
ze<@x9rif=8YvQ3vfF^)WX!`B{00Hsncco|m03VA81ONa4009360763o02=@U00000
I0000005u)NzW@LL

diff --git a/public/biopet-tools/src/test/resources/unvepped.vcf.gz.tbi b/public/biopet-tools/src/test/resources/unvepped.vcf.gz.tbi
index 1895ffffc8cafd9cda37d3b9afb81a9a1bfe4034..2580d3da8a38a5516bf475e240870a07626279a2 100644
GIT binary patch
delta 24
fcma!#pP(<`))f%hmOF=sXGIc&@Bvo4iS`NrW&;Pl

delta 24
fcma!#pP(;b`7$7~E%*=*&x#}lHJ|X46YUiMZO;gp

diff --git a/public/biopet-tools/src/test/scala/nl/lumc/sasc/biopet/tools/BastyGenerateFastaTest.scala b/public/biopet-tools/src/test/scala/nl/lumc/sasc/biopet/tools/BastyGenerateFastaTest.scala
index d9d8c1bee..70e9aefd4 100644
--- a/public/biopet-tools/src/test/scala/nl/lumc/sasc/biopet/tools/BastyGenerateFastaTest.scala
+++ b/public/biopet-tools/src/test/scala/nl/lumc/sasc/biopet/tools/BastyGenerateFastaTest.scala
@@ -35,7 +35,7 @@ class BastyGenerateFastaTest extends TestNGSuite with MockitoSugar with Matchers
     val tmppath = tmp.getAbsolutePath
     tmp.deleteOnExit()
 
-    val arguments = Array("-V", chrQ_path, "--outputVariants", tmppath, "--sampleName", "Child_7006504", "--reference", chrQRef_path, "--outputName", "test")
+    val arguments = Array("-V", chrQ_path, "--outputVariants", tmppath, "--sampleName", "Sample_101", "--reference", chrQRef_path, "--outputName", "test")
     main(arguments)
   }
 
@@ -45,7 +45,7 @@ class BastyGenerateFastaTest extends TestNGSuite with MockitoSugar with Matchers
     val tmppath = tmp.getAbsolutePath
     tmp.deleteOnExit()
 
-    val arguments = Array("-V", chrQ_path, "--outputVariants", tmppath, "--bamFile", bam_path, "--sampleName", "Child_7006504", "--reference", chrQRef_path, "--outputName", "test")
+    val arguments = Array("-V", chrQ_path, "--outputVariants", tmppath, "--bamFile", bam_path, "--sampleName", "Sample_101", "--reference", chrQRef_path, "--outputName", "test")
     main(arguments)
   }
 
@@ -55,7 +55,7 @@ class BastyGenerateFastaTest extends TestNGSuite with MockitoSugar with Matchers
     val tmppath = tmp.getAbsolutePath
     tmp.deleteOnExit()
 
-    val arguments = Array("-V", chrQ_path, "--outputConsensus", tmppath, "--outputConsensusVariants", tmppath, "--bamFile", bam_path, "--sampleName", "Child_7006504", "--reference", chrQRef_path, "--outputName", "test")
+    val arguments = Array("-V", chrQ_path, "--outputConsensus", tmppath, "--outputConsensusVariants", tmppath, "--bamFile", bam_path, "--sampleName", "Sample_101", "--reference", chrQRef_path, "--outputName", "test")
     main(arguments)
   }
 
@@ -63,13 +63,13 @@ class BastyGenerateFastaTest extends TestNGSuite with MockitoSugar with Matchers
     val reader = new VCFFileReader(vepped, false)
     val record = reader.iterator().next()
 
-    val child = mock[Args]
-    when(child.sampleName) thenReturn "Child_7006504"
-    val father = mock[Args]
-    when(father.sampleName) thenReturn "Father_7006506"
+    val one = mock[Args]
+    when(one.sampleName) thenReturn "Sample_101"
+    val two = mock[Args]
+    when(two.sampleName) thenReturn "Sample_102"
 
-    getMaxAllele(record)(child) shouldBe "C-"
-    getMaxAllele(record)(father) shouldBe "CA"
+    getMaxAllele(record)(one) shouldBe "C-"
+    getMaxAllele(record)(two) shouldBe "CA"
 
   }
 
diff --git a/public/biopet-tools/src/test/scala/nl/lumc/sasc/biopet/tools/VcfFilterTest.scala b/public/biopet-tools/src/test/scala/nl/lumc/sasc/biopet/tools/VcfFilterTest.scala
index 80fe1980e..13d959122 100644
--- a/public/biopet-tools/src/test/scala/nl/lumc/sasc/biopet/tools/VcfFilterTest.scala
+++ b/public/biopet-tools/src/test/scala/nl/lumc/sasc/biopet/tools/VcfFilterTest.scala
@@ -65,20 +65,20 @@ class VcfFilterTest extends TestNGSuite with MockitoSugar with Matchers {
     val reader = new VCFFileReader(vepped, false)
     val record = reader.iterator().next()
 
-    hasGenotype(record, List(("Child_7006504", GenotypeType.HET))) shouldBe true
-    hasGenotype(record, List(("Child_7006504", GenotypeType.HOM_VAR))) shouldBe false
-    hasGenotype(record, List(("Child_7006504", GenotypeType.HOM_REF))) shouldBe false
-    hasGenotype(record, List(("Child_7006504", GenotypeType.NO_CALL))) shouldBe false
-    hasGenotype(record, List(("Child_7006504", GenotypeType.MIXED))) shouldBe false
+    hasGenotype(record, List(("Sample_101", GenotypeType.HET))) shouldBe true
+    hasGenotype(record, List(("Sample_101", GenotypeType.HOM_VAR))) shouldBe false
+    hasGenotype(record, List(("Sample_101", GenotypeType.HOM_REF))) shouldBe false
+    hasGenotype(record, List(("Sample_101", GenotypeType.NO_CALL))) shouldBe false
+    hasGenotype(record, List(("Sample_101", GenotypeType.MIXED))) shouldBe false
 
-    hasGenotype(record, List(("Mother_7006508", GenotypeType.HET))) shouldBe false
-    hasGenotype(record, List(("Mother_7006508", GenotypeType.HOM_VAR))) shouldBe false
-    hasGenotype(record, List(("Mother_7006508", GenotypeType.HOM_REF))) shouldBe true
-    hasGenotype(record, List(("Mother_7006508", GenotypeType.NO_CALL))) shouldBe false
-    hasGenotype(record, List(("Mother_7006508", GenotypeType.MIXED))) shouldBe false
+    hasGenotype(record, List(("Sample_103", GenotypeType.HET))) shouldBe false
+    hasGenotype(record, List(("Sample_103", GenotypeType.HOM_VAR))) shouldBe false
+    hasGenotype(record, List(("Sample_103", GenotypeType.HOM_REF))) shouldBe true
+    hasGenotype(record, List(("Sample_103", GenotypeType.NO_CALL))) shouldBe false
+    hasGenotype(record, List(("Sample_103", GenotypeType.MIXED))) shouldBe false
 
-    hasGenotype(record, List(("Mother_7006508", GenotypeType.HOM_REF), ("Child_7006504", GenotypeType.HET))) shouldBe true
-    hasGenotype(record, List(("Mother_7006508", GenotypeType.HET), ("Child_7006504", GenotypeType.HOM_REF))) shouldBe false
+    hasGenotype(record, List(("Sample_103", GenotypeType.HOM_REF), ("Sample_101", GenotypeType.HET))) shouldBe true
+    hasGenotype(record, List(("Sample_103", GenotypeType.HET), ("Sample_101", GenotypeType.HOM_REF))) shouldBe false
   }
 
   @Test def testMinQualScore() = {
@@ -142,42 +142,42 @@ class VcfFilterTest extends TestNGSuite with MockitoSugar with Matchers {
     val reader = new VCFFileReader(vepped, false)
     val record = reader.iterator().next()
 
-    mustHaveVariant(record, List("Child_7006504")) shouldBe true
-    mustHaveVariant(record, List("Child_7006504", "Father_7006506")) shouldBe true
-    mustHaveVariant(record, List("Child_7006504", "Father_7006506", "Mother_7006508")) shouldBe false
+    mustHaveVariant(record, List("Sample_101")) shouldBe true
+    mustHaveVariant(record, List("Sample_101", "Sample_102")) shouldBe true
+    mustHaveVariant(record, List("Sample_101", "Sample_102", "Sample_103")) shouldBe false
   }
 
   @Test def testSameGenotype() = {
     val reader = new VCFFileReader(vepped, false)
     val record = reader.iterator().next()
 
-    notSameGenotype(record, "Child_7006504", "Father_7006506") shouldBe false
-    notSameGenotype(record, "Child_7006504", "Mother_7006508") shouldBe true
-    notSameGenotype(record, "Father_7006506", "Mother_7006508") shouldBe true
+    notSameGenotype(record, "Sample_101", "Sample_102") shouldBe false
+    notSameGenotype(record, "Sample_101", "Sample_103") shouldBe true
+    notSameGenotype(record, "Sample_102", "Sample_103") shouldBe true
   }
 
   @Test def testfilterHetVarToHomVar() = {
     val reader = new VCFFileReader(vepped, false)
     val record = reader.iterator().next()
 
-    filterHetVarToHomVar(record, "Child_7006504", "Father_7006506") shouldBe true
-    filterHetVarToHomVar(record, "Child_7006504", "Mother_7006508") shouldBe true
-    filterHetVarToHomVar(record, "Father_7006506", "Mother_7006508") shouldBe true
+    filterHetVarToHomVar(record, "Sample_101", "Sample_102") shouldBe true
+    filterHetVarToHomVar(record, "Sample_101", "Sample_103") shouldBe true
+    filterHetVarToHomVar(record, "Sample_102", "Sample_103") shouldBe true
   }
 
   @Test def testDeNovo() = {
     val reader = new VCFFileReader(vepped, false)
     val record = reader.iterator().next()
 
-    denovoInSample(record, "Child_7006504") shouldBe false
-    denovoInSample(record, "Father_7006506") shouldBe false
-    denovoInSample(record, "Mother_7006508") shouldBe false
+    denovoInSample(record, "Sample_101") shouldBe false
+    denovoInSample(record, "Sample_102") shouldBe false
+    denovoInSample(record, "Sample_103") shouldBe false
   }
 
   @Test def testResToDom() = {
     val reader = new VCFFileReader(vepped, false)
     val record = reader.iterator().next()
-    val trio = new Trio("Child_7006504", "Father_7006506", "Mother_7006508")
+    val trio = new Trio("Sample_101", "Sample_102", "Sample_103")
 
     resToDom(record, List(trio)) shouldBe false
   }
@@ -185,7 +185,7 @@ class VcfFilterTest extends TestNGSuite with MockitoSugar with Matchers {
   @Test def testTrioCompound = {
     val reader = new VCFFileReader(vepped, false)
     val record = reader.iterator().next()
-    val trio = new Trio("Child_7006504", "Father_7006506", "Mother_7006508")
+    val trio = new Trio("Sample_101", "Sample_102", "Sample_103")
 
     trioCompound(record, List(trio))
   }
@@ -193,7 +193,7 @@ class VcfFilterTest extends TestNGSuite with MockitoSugar with Matchers {
   @Test def testDeNovoTrio = {
     val reader = new VCFFileReader(vepped, false)
     val record = reader.iterator().next()
-    val trio = new Trio("Child_7006504", "Father_7006506", "Mother_7006508")
+    val trio = new Trio("Sample_101", "Sample_102", "Sample_103")
 
     denovoTrio(record, List(trio))
   }
-- 
GitLab