Skip to content

Instantly share code, notes, and snippets.

@muhammadsohailraza
Last active August 6, 2017 07:08
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save muhammadsohailraza/3b371464c95c2ffde353f59879df6742 to your computer and use it in GitHub Desktop.
Save muhammadsohailraza/3b371464c95c2ffde353f59879df6742 to your computer and use it in GitHub Desktop.
Question: Switch the position of REF/ALT alleles based on AA (ancestral allele) annotations in VCF
##fileformat=VCFv4.2
##ALT=<ID=NON_REF,Description="Represents any possible alternative allele at this location">
##FILTER=<ID=LowQual,Description="Low quality">
##FILTER=<ID=VQSRTrancheINDEL99.00to99.90,Description="Truth sensitivity tranche level for INDEL model at VQS Lod: -4.0415 <= x < -0.3987">
##FILTER=<ID=VQSRTrancheINDEL99.90to100.00+,Description="Truth sensitivity tranche level for INDEL model at VQS Lod < -39528.47">
##FILTER=<ID=VQSRTrancheINDEL99.90to100.00,Description="Truth sensitivity tranche level for INDEL model at VQS Lod: -39528.47 <= x < -4.0415">
##FILTER=<ID=VQSRTrancheSNP99.90to100.00+,Description="Truth sensitivity tranche level for SNP model at VQS Lod < -240139.7406">
##FILTER=<ID=VQSRTrancheSNP99.90to100.00,Description="Truth sensitivity tranche level for SNP model at VQS Lod: -240139.7406 <= x < -8.9721">
##FORMAT=<ID=AD,Number=R,Type=Integer,Description="Allelic depths for the ref and alt alleles in the order listed">
##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Approximate read depth (reads with MQ=255 or with bad mates are filtered)">
##FORMAT=<ID=GQ,Number=1,Type=Integer,Description="Genotype Quality">
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
##FORMAT=<ID=MIN_DP,Number=1,Type=Integer,Description="Minimum DP observed within the GVCF block">
##FORMAT=<ID=PGT,Number=1,Type=String,Description="Physical phasing haplotype information, describing how the alternate alleles are phased in relation to one another">
##FORMAT=<ID=PID,Number=1,Type=String,Description="Physical phasing ID information, where each unique ID within a given sample (but not across samples) connects records within a phasing group">
##FORMAT=<ID=PL,Number=G,Type=Integer,Description="Normalized, Phred-scaled likelihoods for genotypes as defined in the VCF specification">
##FORMAT=<ID=RGQ,Number=1,Type=Integer,Description="Unconditional reference genotype confidence, encoded as a phred quality -10*log10 p(genotype call is wrong)">
##FORMAT=<ID=SB,Number=4,Type=Integer,Description="Per-sample component statistics which comprise the Fisher's Exact Test to detect strand bias.">
##GATKCommandLine.ApplyRecalibration.2=<ID=ApplyRecalibration,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null read_filter=[] disable_read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/share_bio/unisvx3/zengchq_group/sohail/Ref/ref-1000G/hs37d5.fa 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 refactor_NDN_cigar_string=false fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 static_quantized_quals=null round_down_quantized=false 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 no_cmdline_in_header=false sites_only=false never_trim_vcf_format_field=false bcf=false bam_compression=null simplifyBAM=false disable_bam_indexing=false generate_md5=false num_threads=15 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 reference_window_stop=0 phone_home= gatk_key=null tag=NA logging_level=INFO log_to_file=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/GATK/FINAL-JOINT/_LOG/log//PHASE4-VQSR-INDEL-ApplyRECALIB.log help=false version=false input=[(RodBinding name=input source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/GATK/FINAL-JOINT//VQSR_PHASE4_snp99.5_raw_indel_2017JUN30.vcf)] recal_file=(RodBinding name=recal_file source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/GATK/FINAL-JOINT//Recal_VQSR_INDEL_PHASE4_JOINT.recal) tranches_file=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/GATK/FINAL-JOINT/tranches_VQSR_INDEL_PHASE4_JOINT.tranches out=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/GATK/FINAL-JOINT/VQSR_PHASE4_snp99.5_indel99.0_2017JUN30.vcf ts_filter_level=99.0 useAlleleSpecificAnnotations=false lodCutoff=null ignore_filter=null ignore_all_filters=false excludeFiltered=false mode=INDEL filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Sun Jul 02 19:09:37 HKT 2017",Epoch=1498993777828,Version=3.6-0-g89b7209>
##GATKCommandLine.ApplyRecalibration=<ID=ApplyRecalibration,CommandLineOptions="analysis_type=ApplyRecalibration input_file=[] showFullBamList=false read_buffer_size=null read_filter=[] disable_read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/share_bio/unisvx3/zengchq_group/sohail/Ref/ref-1000G/hs37d5.fa 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 refactor_NDN_cigar_string=false fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 static_quantized_quals=null round_down_quantized=false 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 no_cmdline_in_header=false sites_only=false never_trim_vcf_format_field=false bcf=false bam_compression=null simplifyBAM=false disable_bam_indexing=false generate_md5=false num_threads=15 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 reference_window_stop=0 phone_home= gatk_key=null tag=NA logging_level=INFO log_to_file=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/GATK/FINAL-JOINT/_LOG/log//PHASE4-VQSR-SNP-ApplyRECALIB.log help=false version=false input=[(RodBinding name=input source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/GATK/FINAL-JOINT//PHASE4_JOINT_snp_indels_ALL_2017JUNE28-new.vcf)] recal_file=(RodBinding name=recal_file source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/GATK/FINAL-JOINT//Recal_VQSR_SNP_PHASE4_JOINT.recal) tranches_file=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/GATK/FINAL-JOINT/tranches_VQSR_SNP_PHASE4_JOINT.tranches out=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/GATK/FINAL-JOINT/VQSR_PHASE4_snp99.5_raw_indel_2017JUN30.vcf ts_filter_level=99.5 useAlleleSpecificAnnotations=false lodCutoff=null ignore_filter=null ignore_all_filters=false excludeFiltered=false mode=SNP filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Sun Jul 02 12:07:34 HKT 2017",Epoch=1498968454049,Version=3.6-0-g89b7209>
##GATKCommandLine.CombineVariants=<ID=CombineVariants,CommandLineOptions="analysis_type=CombineVariants input_file=[] showFullBamList=false read_buffer_size=null read_filter=[] disable_read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/share_bio/unisvx3/zengchq_group/sohail/Ref/ref-1000G/hs37d5.fa 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 refactor_NDN_cigar_string=false fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 static_quantized_quals=null round_down_quantized=false 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 no_cmdline_in_header=false sites_only=false never_trim_vcf_format_field=false bcf=false bam_compression=null simplifyBAM=false disable_bam_indexing=false generate_md5=false num_threads=5 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 reference_window_stop=0 phone_home= gatk_key=null tag=NA logging_level=INFO log_to_file=null help=false version=false variant=[(RodBindingCollection [(RodBinding name=variant source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/NO-CALL/CONCORD//PHASE4_snp99.5_CHR_HWE_BA_NOCALL_HiCOV-LoCoV-Concord.vcf)]), (RodBindingCollection [(RodBinding name=variant2 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/NO-CALL//PHASE4_snp99.5_CHR_HWE_BA_FILTER-HiCoVSAM-NOCALL_10.vcf)])] out=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/NO-CALL/COMBINE/PHASE4_snp99.5-CHR_HWE_BA_NOCALL_Combine.vcf genotypemergeoption=UNIQUIFY filteredrecordsmergetype=KEEP_IF_ANY_UNFILTERED multipleallelesmergetype=BY_TYPE rod_priority_list=null printComplexMerges=false filteredAreUncalled=false minimalVCF=false excludeNonVariants=false setKey=set assumeIdenticalSamples=false minimumN=1 suppressCommandLineHeader=false mergeInfoWithMaxAC=false filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Fri Jul 14 15:56:27 HKT 2017",Epoch=1500018987671,Version=3.6-0-g89b7209>
##GATKCommandLine.GenotypeGVCFs=<ID=GenotypeGVCFs,CommandLineOptions="analysis_type=GenotypeGVCFs input_file=[] showFullBamList=false read_buffer_size=null read_filter=[BadCigar] disable_read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/share_bio/unisvx3/zengchq_group/sohail/Ref/ref-1000G/hs37d5.fa 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 refactor_NDN_cigar_string=false fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 static_quantized_quals=null round_down_quantized=false 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 no_cmdline_in_header=false sites_only=false never_trim_vcf_format_field=false bcf=false bam_compression=null simplifyBAM=false disable_bam_indexing=false generate_md5=false num_threads=20 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 reference_window_stop=0 phone_home= gatk_key=null tag=NA logging_level=INFO log_to_file=null help=false version=false variant=[(RodBindingCollection [(RodBinding name=variant source=/share_bio/unisvx3/zengchq_group/sohail/_TIBETAN/GATK/VAR_CALL//RAW_ALI_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant2 source=/share_bio/unisvx3/zengchq_group/sohail/_TIBETAN/GATK/VAR_CALL//RAW_DINGRI_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant3 source=/share_bio/unisvx3/zengchq_group/sohail/_TIBETAN/GATK/VAR_CALL//RAW_QJ16_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant4 source=/share_bio/unisvx3/zengchq_group/sohail/_TIBETAN/GATK/VAR_CALL//RAW_QJ23_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant5 source=/share_bio/unisvx3/zengchq_group/sohail/_TIBETAN/GATK/VAR_CALL//RAW_QJ28_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant6 source=/share_bio/unisvx3/zengchq_group/sohail/_TIBETAN/GATK/VAR_CALL//RAW_QJ36_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant7 source=/share_bio/unisvx3/zengchq_group/sohail/_TIBETAN/GATK/VAR_CALL//RAW_QJ37_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant8 source=/share_bio/unisvx3/zengchq_group/sohail/_TIBETAN/GATK/VAR_CALL//RAW_QJ46_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant9 source=/share_bio/unisvx3/zengchq_group/sohail/_TIBETAN/GATK/VAR_CALL//RAW_QJ48_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant10 source=/share_bio/unisvx3/zengchq_group/sohail/_TIBETAN/GATK/VAR_CALL//RAW_QJ50_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant11 source=/share_bio/unisvx3/zengchq_group/sohail/_TIBETAN/GATK/VAR_CALL//RAW_QJ6_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant12 source=/share_bio/unisvx3/zengchq_group/sohail/_TIBETAN/GATK/VAR_CALL//RAW_QJ7_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant13 source=/share_bio/unisvx3/zengchq_group/sohail/_SGDP-SAMPLES/GATK/VAR_CALL//RAW_LP6005441-DNA_A09_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant14 source=/share_bio/unisvx3/zengchq_group/sohail/_SGDP-SAMPLES/GATK/VAR_CALL//RAW_LP6005441-DNA_B09_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant15 source=/share_bio/unisvx3/zengchq_group/sohail/_SGDP-SAMPLES/GATK/VAR_CALL//RAW_LP6005441-DNA_D12_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant16 source=/share_bio/unisvx3/zengchq_group/sohail/_SGDP-SAMPLES/GATK/VAR_CALL//RAW_LP6005442-DNA_G01_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant17 source=/share_bio/unisvx3/zengchq_group/sohail/_SGDP-SAMPLES/GATK/VAR_CALL//RAW_LP6005442-DNA_H01_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant18 source=/share_bio/unisvx3/zengchq_group/sohail/_SGDP-SAMPLES/GATK/VAR_CALL//RAW_LP6005443-DNA_E09_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant19 source=/share_bio/unisvx3/zengchq_group/sohail/_SGDP-SAMPLES/GATK/VAR_CALL//RAW_LP6005443-DNA_G01_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant20 source=/share_bio/unisvx3/zengchq_group/sohail/_SGDP-SAMPLES/GATK/VAR_CALL//RAW_LP6005443-DNA_H01_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant21 source=/share_bio/unisvx3/zengchq_group/sohail/_SGDP-SAMPLES/GATK/VAR_CALL//RAW_LP6005592-DNA_E05_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant22 source=/share_bio/unisvx3/zengchq_group/sohail/_SGDP-SAMPLES/GATK/VAR_CALL//RAW_LP6005592-DNA_F05_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant23 source=/share_bio/unisvx3/zengchq_group/sohail/_SGDP-SAMPLES/GATK/VAR_CALL//RAW_LP6005592-DNA_G05_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant24 source=/share_bio/unisvx3/zengchq_group/sohail/_SGDP-SAMPLES/GATK/VAR_CALL//RAW_LP6005619-DNA_A01_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant25 source=/p200/zengchq_group/sohail/_SHERPA/GATK/VAR_CALL//RAW_SRR1101035_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant26 source=/p200/zengchq_group/sohail/_SHERPA/GATK/VAR_CALL//RAW_SRR1101036_snp_indels.g.vcf)]), (RodBindingCollection [(RodBinding name=variant27 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022051.g.vcf)]), (RodBindingCollection [(RodBinding name=variant28 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022052.g.vcf)]), (RodBindingCollection [(RodBinding name=variant29 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022053.g.vcf)]), (RodBindingCollection [(RodBinding name=variant30 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022054.g.vcf)]), (RodBindingCollection [(RodBinding name=variant31 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022055.g.vcf)]), (RodBindingCollection [(RodBinding name=variant32 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022056.g.vcf)]), (RodBindingCollection [(RodBinding name=variant33 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022057.g.vcf)]), (RodBindingCollection [(RodBinding name=variant34 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022058.g.vcf)]), (RodBindingCollection [(RodBinding name=variant35 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022059.g.vcf)]), (RodBindingCollection [(RodBinding name=variant36 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022060.g.vcf)]), (RodBindingCollection [(RodBinding name=variant37 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022061.g.vcf)]), (RodBindingCollection [(RodBinding name=variant38 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022062.g.vcf)]), (RodBindingCollection [(RodBinding name=variant39 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022063.g.vcf)]), (RodBindingCollection [(RodBinding name=variant40 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022064.g.vcf)]), (RodBindingCollection [(RodBinding name=variant41 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022065.g.vcf)]), (RodBindingCollection [(RodBinding name=variant42 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022066.g.vcf)]), (RodBindingCollection [(RodBinding name=variant43 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022067.g.vcf)]), (RodBindingCollection [(RodBinding name=variant44 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022068.g.vcf)]), (RodBindingCollection [(RodBinding name=variant45 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022069.g.vcf)]), (RodBindingCollection [(RodBinding name=variant46 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022070.g.vcf)]), (RodBindingCollection [(RodBinding name=variant47 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022071.g.vcf)]), (RodBindingCollection [(RodBinding name=variant48 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022072.g.vcf)]), (RodBindingCollection [(RodBinding name=variant49 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022073.g.vcf)]), (RodBindingCollection [(RodBinding name=variant50 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022074.g.vcf)]), (RodBindingCollection [(RodBinding name=variant51 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022075.g.vcf)]), (RodBindingCollection [(RodBinding name=variant52 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022076.g.vcf)]), (RodBindingCollection [(RodBinding name=variant53 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022077.g.vcf)]), (RodBindingCollection [(RodBinding name=variant54 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022078.g.vcf)]), (RodBindingCollection [(RodBinding name=variant55 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022079.g.vcf)]), (RodBindingCollection [(RodBinding name=variant56 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022080.g.vcf)]), (RodBindingCollection [(RodBinding name=variant57 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022107.g.vcf)]), (RodBindingCollection [(RodBinding name=variant58 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022108.g.vcf)]), (RodBindingCollection [(RodBinding name=variant59 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022109.g.vcf)]), (RodBindingCollection [(RodBinding name=variant60 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022110.g.vcf)]), (RodBindingCollection [(RodBinding name=variant61 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022111.g.vcf)]), (RodBindingCollection [(RodBinding name=variant62 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022112.g.vcf)]), (RodBindingCollection [(RodBinding name=variant63 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022113.g.vcf)]), (RodBindingCollection [(RodBinding name=variant64 source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-3/GATK/VAR_CALL//RAW_AAGC022114.g.vcf)]), (RodBindingCollection [(RodBinding name=variant65 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_103RQ.g.vcf)]), (RodBindingCollection [(RodBinding name=variant66 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_109EP.g.vcf)]), (RodBindingCollection [(RodBinding name=variant67 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_16EJ.g.vcf)]), (RodBindingCollection [(RodBinding name=variant68 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_31CE.g.vcf)]), (RodBindingCollection [(RodBinding name=variant69 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_40JI.g.vcf)]), (RodBindingCollection [(RodBinding name=variant70 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_43AM.g.vcf)]), (RodBindingCollection [(RodBinding name=variant71 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_50JB.g.vcf)]), (RodBindingCollection [(RodBinding name=variant72 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_52AA.g.vcf)]), (RodBindingCollection [(RodBinding name=variant73 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_53JR.g.vcf)]), (RodBindingCollection [(RodBinding name=variant74 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_55EC.g.vcf)]), (RodBindingCollection [(RodBinding name=variant75 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_56AM.g.vcf)]), (RodBindingCollection [(RodBinding name=variant76 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_63PF.g.vcf)]), (RodBindingCollection [(RodBinding name=variant77 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_72DA.g.vcf)]), (RodBindingCollection [(RodBinding name=variant78 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_76AS.g.vcf)]), (RodBindingCollection [(RodBinding name=variant79 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_77JN.g.vcf)]), (RodBindingCollection [(RodBinding name=variant80 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_7PCE.g.vcf)]), (RodBindingCollection [(RodBinding name=variant81 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_82JC.g.vcf)]), (RodBindingCollection [(RodBinding name=variant82 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_93BS.g.vcf)]), (RodBindingCollection [(RodBinding name=variant83 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_96JP.g.vcf)]), (RodBindingCollection [(RodBinding name=variant84 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL//RAW_98AC.g.vcf)]), (RodBindingCollection [(RodBinding name=variant85 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_QUECHUA/_GATK/VAR-CALL//RAW_LP6005677-DNA_E01.g.vcf)]), (RodBindingCollection [(RodBinding name=variant86 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_QUECHUA/_GATK/VAR-CALL//RAW_LP6005677-DNA_F01.g.vcf)]), (RodBindingCollection [(RodBinding name=variant87 source=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_QUECHUA/_GATK/VAR-CALL//RAW_LP6005519-DNA_G02.g.vcf)]), (RodBindingCollection [(RodBinding name=variant88 source=/share_bio/unisvx3/zengchq_group/sohail/_ETHNIC/_GATK/8_VAR-CALL//RAW_LP6005441-DNA_B07.g.vcf)]), (RodBindingCollection [(RodBinding name=variant89 source=/share_bio/unisvx3/zengchq_group/sohail/_ETHNIC/_GATK/8_VAR-CALL//RAW_LP6005441-DNA_C05.g.vcf)]), (RodBindingCollection [(RodBinding name=variant90 source=/share_bio/unisvx3/zengchq_group/sohail/_ETHNIC/_GATK/8_VAR-CALL//RAW_LP6005441-DNA_C08.g.vcf)]), (RodBindingCollection [(RodBinding name=variant91 source=/share_bio/unisvx3/zengchq_group/sohail/_ETHNIC/_GATK/8_VAR-CALL//RAW_LP6005441-DNA_D04.g.vcf)]), (RodBindingCollection [(RodBinding name=variant92 source=/share_bio/unisvx3/zengchq_group/sohail/_ETHNIC/_GATK/8_VAR-CALL//RAW_LP6005441-DNA_D05.g.vcf)]), (RodBindingCollection [(RodBinding name=variant93 source=/share_bio/unisvx3/zengchq_group/sohail/_ETHNIC/_GATK/8_VAR-CALL//RAW_LP6005441-DNA_D08.g.vcf)]), (RodBindingCollection [(RodBinding name=variant94 source=/share_bio/unisvx3/zengchq_group/sohail/_ETHNIC/_GATK/8_VAR-CALL//RAW_LP6005441-DNA_F12.g.vcf)]), (RodBindingCollection [(RodBinding name=variant95 source=/share_bio/unisvx3/zengchq_group/sohail/_ETHNIC/_GATK/8_VAR-CALL//RAW_LP6005442-DNA_B01.g.vcf)]), (RodBindingCollection [(RodBinding name=variant96 source=/share_bio/unisvx3/zengchq_group/sohail/_ETHNIC/_GATK/8_VAR-CALL//RAW_LP6005443-DNA_A02.g.vcf)]), (RodBindingCollection [(RodBinding name=variant97 source=/share_bio/unisvx3/zengchq_group/sohail/_ETHNIC/_GATK/8_VAR-CALL//RAW_LP6005443-DNA_B01.g.vcf)]), (RodBindingCollection [(RodBinding name=variant98 source=/share_bio/unisvx3/zengchq_group/sohail/_ETHNIC/_GATK/8_VAR-CALL//RAW_LP6005443-DNA_B02.g.vcf)]), (RodBindingCollection [(RodBinding name=variant99 source=/share_bio/unisvx3/zengchq_group/sohail/_ETHNIC/_GATK/8_VAR-CALL//RAW_LP6005443-DNA_E01.g.vcf)]), (RodBindingCollection [(RodBinding name=variant100 source=/share_bio/unisvx3/zengchq_group/sohail/_ETHNIC/_GATK/8_VAR-CALL//RAW_LP6005443-DNA_F01.g.vcf)]), (RodBindingCollection [(RodBinding name=variant101 source=/share_bio/unisvx3/zengchq_group/sohail/_ETHNIC/_GATK/8_VAR-CALL//RAW_LP6005592-DNA_D03.g.vcf)]), (RodBindingCollection [(RodBinding name=variant102 source=/share_bio/unisvx3/zengchq_group/sohail/_ETHNIC/_GATK/8_VAR-CALL//RAW_SS6004469.g.vcf)]), (RodBindingCollection [(RodBinding name=variant103 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP31_4.g.vcf)]), (RodBindingCollection [(RodBinding name=variant104 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP31_5.g.vcf)]), (RodBindingCollection [(RodBinding name=variant105 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP32_2.g.vcf)]), (RodBindingCollection [(RodBinding name=variant106 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP38_2.g.vcf)]), (RodBindingCollection [(RodBinding name=variant107 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP38_3.g.vcf)]), (RodBindingCollection [(RodBinding name=variant108 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP42_2.g.vcf)]), (RodBindingCollection [(RodBinding name=variant109 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP45_2.g.vcf)]), (RodBindingCollection [(RodBinding name=variant110 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP45_3.g.vcf)]), (RodBindingCollection [(RodBinding name=variant111 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP47_2.g.vcf)]), (RodBindingCollection [(RodBinding name=variant112 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP47_3.g.vcf)]), (RodBindingCollection [(RodBinding name=variant113 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP49_2.g.vcf)]), (RodBindingCollection [(RodBinding name=variant114 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP49_3.g.vcf)]), (RodBindingCollection [(RodBinding name=variant115 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP50_2.g.vcf)]), (RodBindingCollection [(RodBinding name=variant116 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP50_3.g.vcf)]), (RodBindingCollection [(RodBinding name=variant117 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP55_2.g.vcf)]), (RodBindingCollection [(RodBinding name=variant118 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP55_3.g.vcf)]), (RodBindingCollection [(RodBinding name=variant119 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP58_2.g.vcf)]), (RodBindingCollection [(RodBinding name=variant120 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP58_3.g.vcf)]), (RodBindingCollection [(RodBinding name=variant121 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP59_2.g.vcf)]), (RodBindingCollection [(RodBinding name=variant122 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP64_2.g.vcf)]), (RodBindingCollection [(RodBinding name=variant123 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP64_3.g.vcf)]), (RodBindingCollection [(RodBinding name=variant124 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP67_5.g.vcf)]), (RodBindingCollection [(RodBinding name=variant125 source=/p200/zengchq_group/sohail/_HM-SAMPLES//RAW_HMP72_3.g.vcf)]), (RodBindingCollection [(RodBinding name=variant126 source=/share_bio/unisvx3/zengchq_group/sohail/_HM-SAMPLES/GATK/VAR-CALL//RAW_HMP67_4.g.vcf)])] out=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/GATK/FINAL-JOINT/PHASE4_JOINT_snp_indels_ALL_2017JUNE28-new.vcf includeNonVariantSites=false uniquifySamples=false annotateNDA=false heterozygosity=0.001 indel_heterozygosity=1.25E-4 standard_min_confidence_threshold_for_calling=30.0 standard_min_confidence_threshold_for_emitting=30.0 max_alternate_alleles=6 max_num_PL_values=100 input_prior=[] sample_ploidy=2 annotation=[] group=[StandardAnnotation] dbsnp=(RodBinding name=dbsnp source=/share_bio/unisvx3/zengchq_group/sohail/Ref//dbsnp_138.b37.vcf) filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Wed Jun 28 12:27:02 CST 2017",Epoch=1498624022719,Version=3.6-0-g89b7209>
##GATKCommandLine.HaplotypeCaller=<ID=HaplotypeCaller,CommandLineOptions="analysis_type=HaplotypeCaller input_file=[/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/BQSR//50JB_BQSR.bam] showFullBamList=false read_buffer_size=null read_filter=[BadCigar] disable_read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/share_bio/unisvx3/zengchq_group/sohail/Ref/ref-1000G/hs37d5.fa nonDeterministicRandomSeed=false disableDithering=false maxRuntime=-1 maxRuntimeUnits=MINUTES downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=500 baq=OFF baqGapOpenPenalty=40.0 refactor_NDN_cigar_string=false fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 static_quantized_quals=null round_down_quantized=false 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 no_cmdline_in_header=false sites_only=false never_trim_vcf_format_field=false bcf=false bam_compression=null simplifyBAM=false disable_bam_indexing=false generate_md5=false num_threads=1 num_cpu_threads_per_data_thread=50 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=LINEAR variant_index_parameter=128000 reference_window_stop=0 phone_home= gatk_key=null tag=NA logging_level=INFO log_to_file=/share_bio/unisvx3/zengchq_group/sohail/_ANDEANS/_GATK/VAR_CALL/_LOG/log//HC-4.log help=false version=false likelihoodCalculationEngine=PairHMM heterogeneousKmerSizeResolution=COMBO_MIN dbsnp=(RodBinding name=dbsnp source=/share_bio/unisvx3/zengchq_group/sohail/Ref//dbsnp_138.b37.vcf) dontTrimActiveRegions=false maxDiscARExtension=25 maxGGAARExtension=300 paddingAroundIndels=150 paddingAroundSNPs=20 comp=[] annotation=[StrandBiasBySample] excludeAnnotation=[ChromosomeCounts, FisherStrand, StrandOddsRatio, QualByDepth] group=[StandardAnnotation, StandardHCAnnotation] debug=false useFilteredReadsForAnnotations=false emitRefConfidence=GVCF bamOutput=null bamWriterType=CALLED_HAPLOTYPES emitDroppedReads=false disableOptimizations=false annotateNDA=false heterozygosity=0.001 indel_heterozygosity=1.25E-4 standard_min_confidence_threshold_for_calling=-0.0 standard_min_confidence_threshold_for_emitting=-0.0 max_alternate_alleles=6 max_num_PL_values=100 input_prior=[] sample_ploidy=2 genotyping_mode=DISCOVERY alleles=(RodBinding name= source=UNBOUND) contamination_fraction_to_filter=0.0 contamination_fraction_per_sample_file=null p_nonref_model=null exactcallslog=null output_mode=EMIT_VARIANTS_ONLY allSitePLs=true gcpHMM=10 pair_hmm_implementation=VECTOR_LOGLESS_CACHING pair_hmm_sub_implementation=ENABLE_ALL always_load_vector_logless_PairHMM_lib=false phredScaledGlobalReadMismappingRate=45 noFpga=false sample_name=null kmerSize=[10, 25] dontIncreaseKmerSizesForCycles=false allowNonUniqueKmersInRef=false numPruningSamples=1 recoverDanglingHeads=false doNotRecoverDanglingBranches=false minDanglingBranchLength=4 consensus=false maxNumHaplotypesInPopulation=128 errorCorrectKmers=false minPruning=2 debugGraphTransformations=false allowCyclesInKmerGraphToGeneratePaths=false graphOutput=null kmerLengthForReadErrorCorrection=25 minObservationsForKmerToBeSolid=20 GVCFGQBands=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 70, 80, 90, 99] indelSizeToEliminateInRefModel=10 min_base_quality_score=10 includeUmappedReads=false useAllelesTrigger=false doNotRunPhysicalPhasing=false keepRG=null justDetermineActiveRegions=false dontGenotype=false dontUseSoftClippedBases=false captureAssemblyFailureBAM=false errorCorrectReads=false pcr_indel_model=CONSERVATIVE maxReadsInRegionPerSample=10000 minReadsPerAlignmentStart=10 mergeVariantsViaLD=false activityProfileOut=null activeRegionOut=null activeRegionIn=null activeRegionExtension=null forceActive=false activeRegionMaxSize=null bandPassSigma=null maxProbPropagationDistance=50 activeProbabilityThreshold=0.002 min_mapping_quality_score=20 filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Wed Mar 29 14:14:43 HKT 2017",Epoch=1490768083716,Version=3.6-0-g89b7209>
##GATKCommandLine.SelectVariants.2=<ID=SelectVariants,CommandLineOptions="analysis_type=SelectVariants input_file=[] showFullBamList=false read_buffer_size=null read_filter=[] disable_read_filter=[] intervals=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, X, Y] excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/share_bio/unisvx3/zengchq_group/sohail/Ref/ref-1000G/hs37d5.fa 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 refactor_NDN_cigar_string=false fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 static_quantized_quals=null round_down_quantized=false 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 no_cmdline_in_header=false sites_only=false never_trim_vcf_format_field=false bcf=false bam_compression=null simplifyBAM=false disable_bam_indexing=false generate_md5=false num_threads=5 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 reference_window_stop=0 phone_home= gatk_key=null tag=NA logging_level=INFO log_to_file=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/SELECT-VAR/_LOG/log//PHASE4-JOINT-SelectSNPs-CHR.log help=false version=false variant=(RodBinding name=variant source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/SELECT-VAR//PHASE4_snp99.5_only_2017JUN30.vcf) discordance=(RodBinding name= source=UNBOUND) concordance=(RodBinding name= source=UNBOUND) out=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/SELECT-VAR/PHASE4_snp99.5_only_CHR_2017JUN30.vcf sample_name=[] sample_expressions=null sample_file=null exclude_sample_name=[] exclude_sample_file=[] exclude_sample_expressions=[] selectexpressions=[] invertselect=false excludeNonVariants=false excludeFiltered=false preserveAlleles=false removeUnusedAlternates=false restrictAllelesTo=ALL keepOriginalAC=false keepOriginalDP=false mendelianViolation=false invertMendelianViolation=false mendelianViolationQualThreshold=0.0 select_random_fraction=0.0 remove_fraction_genotypes=0.0 selectTypeToInclude=[SNP] selectTypeToExclude=[] keepIDs=null excludeIDs=null fullyDecode=false justRead=false maxIndelSize=2147483647 minIndelSize=0 maxFilteredGenotypes=2147483647 minFilteredGenotypes=0 maxFractionFilteredGenotypes=1.0 minFractionFilteredGenotypes=0.0 maxNOCALLnumber=2147483647 maxNOCALLfraction=1.0 setFilteredGtToNocall=false ALLOW_NONOVERLAPPING_COMMAND_LINE_SAMPLES=false forceValidOutput=false filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Wed Jul 05 11:42:16 HKT 2017",Epoch=1499226136876,Version=3.6-0-g89b7209>
##GATKCommandLine.SelectVariants.3=<ID=SelectVariants,CommandLineOptions="analysis_type=SelectVariants input_file=[] showFullBamList=false read_buffer_size=null read_filter=[] disable_read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/share_bio/unisvx3/zengchq_group/sohail/Ref/ref-1000G/hs37d5.fa 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 refactor_NDN_cigar_string=false fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 static_quantized_quals=null round_down_quantized=false 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 no_cmdline_in_header=false sites_only=false never_trim_vcf_format_field=false bcf=false bam_compression=null simplifyBAM=false disable_bam_indexing=false generate_md5=false num_threads=5 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 reference_window_stop=0 phone_home= gatk_key=null tag=NA logging_level=INFO log_to_file=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/SELECT-BA/_LOG/log//Phase4-Select-BA.log help=false version=false variant=(RodBinding name=variant source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/PLINK/SELECT-SNP//PHASE4_snp99.5_only_CHR_HWE_P20new.recode.vcf) discordance=(RodBinding name= source=UNBOUND) concordance=(RodBinding name= source=UNBOUND) out=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/SELECT-BA/PHASE4_snp99.5_CHR_HWE_BA_FILTER.vcf sample_name=[] sample_expressions=null sample_file=null exclude_sample_name=[] exclude_sample_file=[] exclude_sample_expressions=[] selectexpressions=[] invertselect=false excludeNonVariants=false excludeFiltered=false preserveAlleles=false removeUnusedAlternates=false restrictAllelesTo=BIALLELIC keepOriginalAC=false keepOriginalDP=false mendelianViolation=false invertMendelianViolation=false mendelianViolationQualThreshold=0.0 select_random_fraction=0.0 remove_fraction_genotypes=0.0 selectTypeToInclude=[SNP] selectTypeToExclude=[] keepIDs=null excludeIDs=null fullyDecode=false justRead=false maxIndelSize=2147483647 minIndelSize=0 maxFilteredGenotypes=2147483647 minFilteredGenotypes=0 maxFractionFilteredGenotypes=1.0 minFractionFilteredGenotypes=0.0 maxNOCALLnumber=2147483647 maxNOCALLfraction=1.0 setFilteredGtToNocall=false ALLOW_NONOVERLAPPING_COMMAND_LINE_SAMPLES=false forceValidOutput=false filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Sun Jul 09 23:57:54 HKT 2017",Epoch=1499615874196,Version=3.6-0-g89b7209>
##GATKCommandLine.SelectVariants.4=<ID=SelectVariants,CommandLineOptions="analysis_type=SelectVariants input_file=[] showFullBamList=false read_buffer_size=null read_filter=[] disable_read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/share_bio/unisvx3/zengchq_group/sohail/Ref/ref-1000G/hs37d5.fa 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 refactor_NDN_cigar_string=false fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 static_quantized_quals=null round_down_quantized=false 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 no_cmdline_in_header=false sites_only=false never_trim_vcf_format_field=false bcf=false bam_compression=null simplifyBAM=false disable_bam_indexing=false generate_md5=false num_threads=5 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 reference_window_stop=0 phone_home= gatk_key=null tag=NA logging_level=INFO log_to_file=null help=false version=false variant=(RodBinding name=variant source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/SELECT-BA//PHASE4_snp99.5_CHR_HWE_BA_FILTER.vcf) discordance=(RodBinding name= source=UNBOUND) concordance=(RodBinding name= source=UNBOUND) out=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/SELECT-BA/SEP-SAMPLE/PHASE4_snp99.5_CHR_HWE_BA_FILTER-LoCoVSamples.vcf sample_name=[] sample_expressions=null sample_file=[/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/SELECT-BA/SEP-SAMPLE/LoCOV-SAMPLELIST.txt] exclude_sample_name=[] exclude_sample_file=[] exclude_sample_expressions=[] selectexpressions=[] invertselect=false excludeNonVariants=false excludeFiltered=false preserveAlleles=false removeUnusedAlternates=false restrictAllelesTo=BIALLELIC keepOriginalAC=false keepOriginalDP=false mendelianViolation=false invertMendelianViolation=false mendelianViolationQualThreshold=0.0 select_random_fraction=0.0 remove_fraction_genotypes=0.0 selectTypeToInclude=[SNP] selectTypeToExclude=[] keepIDs=null excludeIDs=null fullyDecode=false justRead=false maxIndelSize=2147483647 minIndelSize=0 maxFilteredGenotypes=2147483647 minFilteredGenotypes=0 maxFractionFilteredGenotypes=1.0 minFractionFilteredGenotypes=0.0 maxNOCALLnumber=2147483647 maxNOCALLfraction=1.0 setFilteredGtToNocall=false ALLOW_NONOVERLAPPING_COMMAND_LINE_SAMPLES=false forceValidOutput=false filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Mon Jul 10 11:32:01 HKT 2017",Epoch=1499657521554,Version=3.6-0-g89b7209>
##GATKCommandLine.SelectVariants.5=<ID=SelectVariants,CommandLineOptions="analysis_type=SelectVariants input_file=[] showFullBamList=false read_buffer_size=null read_filter=[] disable_read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/share_bio/unisvx3/zengchq_group/sohail/Ref/ref-1000G/hs37d5.fa 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 refactor_NDN_cigar_string=false fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 static_quantized_quals=null round_down_quantized=false 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 no_cmdline_in_header=false sites_only=false never_trim_vcf_format_field=false bcf=false bam_compression=null simplifyBAM=false disable_bam_indexing=false generate_md5=false num_threads=5 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 reference_window_stop=0 phone_home= gatk_key=null tag=NA logging_level=INFO log_to_file=null help=false version=false variant=(RodBinding name=variant source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/SELECT-BA/SEP-SAMPLE//PHASE4_snp99.5_CHR_HWE_BA_FILTER-LoCoVSamples.vcf) discordance=(RodBinding name= source=UNBOUND) concordance=(RodBinding name=concordance source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/NO-CALL//PHASE4_snp99.5_CHR_HWE_BA_FILTER-HiCoVSAM-NOCALL_10.vcf) out=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/NO-CALL/CONCORD/PHASE4_snp99.5_CHR_HWE_BA_NOCALL_HiCOV-LoCoV-Concord.vcf sample_name=[] sample_expressions=null sample_file=null exclude_sample_name=[] exclude_sample_file=[] exclude_sample_expressions=[] selectexpressions=[] invertselect=false excludeNonVariants=false excludeFiltered=false preserveAlleles=false removeUnusedAlternates=false restrictAllelesTo=ALL keepOriginalAC=false keepOriginalDP=false mendelianViolation=false invertMendelianViolation=false mendelianViolationQualThreshold=0.0 select_random_fraction=0.0 remove_fraction_genotypes=0.0 selectTypeToInclude=[] selectTypeToExclude=[] keepIDs=null excludeIDs=null fullyDecode=false justRead=false maxIndelSize=2147483647 minIndelSize=0 maxFilteredGenotypes=2147483647 minFilteredGenotypes=0 maxFractionFilteredGenotypes=1.0 minFractionFilteredGenotypes=0.0 maxNOCALLnumber=2147483647 maxNOCALLfraction=1.0 setFilteredGtToNocall=false ALLOW_NONOVERLAPPING_COMMAND_LINE_SAMPLES=false forceValidOutput=false filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Fri Jul 14 09:51:06 HKT 2017",Epoch=1499997066797,Version=3.6-0-g89b7209>
##GATKCommandLine.SelectVariants.6=<ID=SelectVariants,CommandLineOptions="analysis_type=SelectVariants input_file=[] showFullBamList=false read_buffer_size=null read_filter=[] disable_read_filter=[] intervals=[1] excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/share_bio/unisvx3/zengchq_group/sohail/Ref/ref-1000G/hs37d5.fa 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 refactor_NDN_cigar_string=false fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 static_quantized_quals=null round_down_quantized=false 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 no_cmdline_in_header=false sites_only=false never_trim_vcf_format_field=false bcf=false bam_compression=null simplifyBAM=false disable_bam_indexing=false generate_md5=false num_threads=5 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 reference_window_stop=0 phone_home= gatk_key=null tag=NA logging_level=INFO log_to_file=null help=false version=false variant=(RodBinding name=variant source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/NO-CALL/COMBINE//PHASE4_snp99.5-CHR_HWE_BA_NOCALL_Combine_RN1.vcf) discordance=(RodBinding name= source=UNBOUND) concordance=(RodBinding name= source=UNBOUND) out=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/EPO-AA/SELECT-CHR/PHASE4-CHR1.vcf sample_name=[] sample_expressions=null sample_file=null exclude_sample_name=[] exclude_sample_file=[] exclude_sample_expressions=[] selectexpressions=[] invertselect=false excludeNonVariants=false excludeFiltered=false preserveAlleles=false removeUnusedAlternates=false restrictAllelesTo=ALL keepOriginalAC=false keepOriginalDP=false mendelianViolation=false invertMendelianViolation=false mendelianViolationQualThreshold=0.0 select_random_fraction=0.0 remove_fraction_genotypes=0.0 selectTypeToInclude=[SNP] selectTypeToExclude=[] keepIDs=null excludeIDs=null fullyDecode=false justRead=false maxIndelSize=2147483647 minIndelSize=0 maxFilteredGenotypes=2147483647 minFilteredGenotypes=0 maxFractionFilteredGenotypes=1.0 minFractionFilteredGenotypes=0.0 maxNOCALLnumber=2147483647 maxNOCALLfraction=1.0 setFilteredGtToNocall=false ALLOW_NONOVERLAPPING_COMMAND_LINE_SAMPLES=false forceValidOutput=false filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Thu Jul 20 19:22:19 HKT 2017",Epoch=1500549739085,Version=3.6-0-g89b7209>
##GATKCommandLine.SelectVariants.7=<ID=SelectVariants,Version=3.6-0-g89b7209,Date="Tue Aug 01 14:48:54 HKT 2017",Epoch=1501570134220,CommandLineOptions="analysis_type=SelectVariants input_file=[] showFullBamList=false read_buffer_size=null read_filter=[] disable_read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/share_bio/unisvx3/zengchq_group/sohail/Ref/ref-1000G/hs37d5.fa 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 refactor_NDN_cigar_string=false fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 static_quantized_quals=null round_down_quantized=false 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 no_cmdline_in_header=false sites_only=false never_trim_vcf_format_field=false bcf=false bam_compression=null simplifyBAM=false disable_bam_indexing=false generate_md5=false num_threads=5 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 reference_window_stop=0 phone_home= gatk_key=null tag=NA logging_level=INFO log_to_file=null help=false version=false variant=(RodBinding name=variant source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/EPO-AA/FILL-AA//PHASE4_snp99.5-CHR_HWE_BA_NOCALL_Combine_RN1-AA.vcf) discordance=(RodBinding name= source=UNBOUND) concordance=(RodBinding name= source=UNBOUND) out=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/SELECTION/FREQ/SHP/PHASE4_snp99.5-CHR_FINAL_SHERPA-AA.vcf sample_name=[] sample_expressions=null sample_file=[/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/_POP-SPECIFIC-SNP/SHP/SHP.txt] exclude_sample_name=[] exclude_sample_file=[] exclude_sample_expressions=[] selectexpressions=[] invertselect=false excludeNonVariants=false excludeFiltered=false preserveAlleles=false removeUnusedAlternates=false restrictAllelesTo=BIALLELIC keepOriginalAC=false keepOriginalDP=false mendelianViolation=false invertMendelianViolation=false mendelianViolationQualThreshold=0.0 select_random_fraction=0.0 remove_fraction_genotypes=0.0 selectTypeToInclude=[SNP] selectTypeToExclude=[] keepIDs=null excludeIDs=null fullyDecode=false justRead=false maxIndelSize=2147483647 minIndelSize=0 maxFilteredGenotypes=2147483647 minFilteredGenotypes=0 maxFractionFilteredGenotypes=1.0 minFractionFilteredGenotypes=0.0 maxNOCALLnumber=2147483647 maxNOCALLfraction=1.0 setFilteredGtToNocall=false ALLOW_NONOVERLAPPING_COMMAND_LINE_SAMPLES=false forceValidOutput=false filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false">
##GATKCommandLine.SelectVariants=<ID=SelectVariants,CommandLineOptions="analysis_type=SelectVariants input_file=[] showFullBamList=false read_buffer_size=null read_filter=[] disable_read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL interval_padding=0 reference_sequence=/share_bio/unisvx3/zengchq_group/sohail/Ref/ref-1000G/hs37d5.fa 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 refactor_NDN_cigar_string=false fix_misencoded_quality_scores=false allow_potentially_misencoded_quality_scores=false useOriginalQualities=false defaultBaseQualities=-1 performanceLog=null BQSR=null quantize_quals=0 static_quantized_quals=null round_down_quantized=false 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 no_cmdline_in_header=false sites_only=false never_trim_vcf_format_field=false bcf=false bam_compression=null simplifyBAM=false disable_bam_indexing=false generate_md5=false num_threads=5 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 reference_window_stop=0 phone_home= gatk_key=null tag=NA logging_level=INFO log_to_file=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/SELECT-VAR/_LOG/log//PHASE4-JOINT-SelectSNPs.log help=false version=false variant=(RodBinding name=variant source=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/GATK/FINAL-JOINT//VQSR_PHASE4_snp99.5_indel99.0_2017JUN30.vcf) discordance=(RodBinding name= source=UNBOUND) concordance=(RodBinding name= source=UNBOUND) out=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/SELECT-VAR/PHASE4_snp99.5_only_2017JUN30.vcf sample_name=[] sample_expressions=null sample_file=null exclude_sample_name=[] exclude_sample_file=[] exclude_sample_expressions=[] selectexpressions=[] invertselect=false excludeNonVariants=false excludeFiltered=false preserveAlleles=false removeUnusedAlternates=false restrictAllelesTo=ALL keepOriginalAC=false keepOriginalDP=false mendelianViolation=false invertMendelianViolation=false mendelianViolationQualThreshold=0.0 select_random_fraction=0.0 remove_fraction_genotypes=0.0 selectTypeToInclude=[SNP] selectTypeToExclude=[] keepIDs=null excludeIDs=null fullyDecode=false justRead=false maxIndelSize=2147483647 minIndelSize=0 maxFilteredGenotypes=2147483647 minFilteredGenotypes=0 maxFractionFilteredGenotypes=1.0 minFractionFilteredGenotypes=0.0 maxNOCALLnumber=2147483647 maxNOCALLfraction=1.0 setFilteredGtToNocall=false ALLOW_NONOVERLAPPING_COMMAND_LINE_SAMPLES=false forceValidOutput=false filter_reads_with_N_cigar=false filter_mismatching_base_and_quals=false filter_bases_not_stored=false",Date="Wed Jul 05 11:14:17 HKT 2017",Epoch=1499224457660,Version=3.6-0-g89b7209>
##GVCFBlock0-1=minGQ=0(inclusive),maxGQ=1(exclusive)
##GVCFBlock1-2=minGQ=1(inclusive),maxGQ=2(exclusive)
##GVCFBlock10-11=minGQ=10(inclusive),maxGQ=11(exclusive)
##GVCFBlock11-12=minGQ=11(inclusive),maxGQ=12(exclusive)
##GVCFBlock12-13=minGQ=12(inclusive),maxGQ=13(exclusive)
##GVCFBlock13-14=minGQ=13(inclusive),maxGQ=14(exclusive)
##GVCFBlock14-15=minGQ=14(inclusive),maxGQ=15(exclusive)
##GVCFBlock15-16=minGQ=15(inclusive),maxGQ=16(exclusive)
##GVCFBlock16-17=minGQ=16(inclusive),maxGQ=17(exclusive)
##GVCFBlock17-18=minGQ=17(inclusive),maxGQ=18(exclusive)
##GVCFBlock18-19=minGQ=18(inclusive),maxGQ=19(exclusive)
##GVCFBlock19-20=minGQ=19(inclusive),maxGQ=20(exclusive)
##GVCFBlock2-3=minGQ=2(inclusive),maxGQ=3(exclusive)
##GVCFBlock20-21=minGQ=20(inclusive),maxGQ=21(exclusive)
##GVCFBlock21-22=minGQ=21(inclusive),maxGQ=22(exclusive)
##GVCFBlock22-23=minGQ=22(inclusive),maxGQ=23(exclusive)
##GVCFBlock23-24=minGQ=23(inclusive),maxGQ=24(exclusive)
##GVCFBlock24-25=minGQ=24(inclusive),maxGQ=25(exclusive)
##GVCFBlock25-26=minGQ=25(inclusive),maxGQ=26(exclusive)
##GVCFBlock26-27=minGQ=26(inclusive),maxGQ=27(exclusive)
##GVCFBlock27-28=minGQ=27(inclusive),maxGQ=28(exclusive)
##GVCFBlock28-29=minGQ=28(inclusive),maxGQ=29(exclusive)
##GVCFBlock29-30=minGQ=29(inclusive),maxGQ=30(exclusive)
##GVCFBlock3-4=minGQ=3(inclusive),maxGQ=4(exclusive)
##GVCFBlock30-31=minGQ=30(inclusive),maxGQ=31(exclusive)
##GVCFBlock31-32=minGQ=31(inclusive),maxGQ=32(exclusive)
##GVCFBlock32-33=minGQ=32(inclusive),maxGQ=33(exclusive)
##GVCFBlock33-34=minGQ=33(inclusive),maxGQ=34(exclusive)
##GVCFBlock34-35=minGQ=34(inclusive),maxGQ=35(exclusive)
##GVCFBlock35-36=minGQ=35(inclusive),maxGQ=36(exclusive)
##GVCFBlock36-37=minGQ=36(inclusive),maxGQ=37(exclusive)
##GVCFBlock37-38=minGQ=37(inclusive),maxGQ=38(exclusive)
##GVCFBlock38-39=minGQ=38(inclusive),maxGQ=39(exclusive)
##GVCFBlock39-40=minGQ=39(inclusive),maxGQ=40(exclusive)
##GVCFBlock4-5=minGQ=4(inclusive),maxGQ=5(exclusive)
##GVCFBlock40-41=minGQ=40(inclusive),maxGQ=41(exclusive)
##GVCFBlock41-42=minGQ=41(inclusive),maxGQ=42(exclusive)
##GVCFBlock42-43=minGQ=42(inclusive),maxGQ=43(exclusive)
##GVCFBlock43-44=minGQ=43(inclusive),maxGQ=44(exclusive)
##GVCFBlock44-45=minGQ=44(inclusive),maxGQ=45(exclusive)
##GVCFBlock45-46=minGQ=45(inclusive),maxGQ=46(exclusive)
##GVCFBlock46-47=minGQ=46(inclusive),maxGQ=47(exclusive)
##GVCFBlock47-48=minGQ=47(inclusive),maxGQ=48(exclusive)
##GVCFBlock48-49=minGQ=48(inclusive),maxGQ=49(exclusive)
##GVCFBlock49-50=minGQ=49(inclusive),maxGQ=50(exclusive)
##GVCFBlock5-6=minGQ=5(inclusive),maxGQ=6(exclusive)
##GVCFBlock50-51=minGQ=50(inclusive),maxGQ=51(exclusive)
##GVCFBlock51-52=minGQ=51(inclusive),maxGQ=52(exclusive)
##GVCFBlock52-53=minGQ=52(inclusive),maxGQ=53(exclusive)
##GVCFBlock53-54=minGQ=53(inclusive),maxGQ=54(exclusive)
##GVCFBlock54-55=minGQ=54(inclusive),maxGQ=55(exclusive)
##GVCFBlock55-56=minGQ=55(inclusive),maxGQ=56(exclusive)
##GVCFBlock56-57=minGQ=56(inclusive),maxGQ=57(exclusive)
##GVCFBlock57-58=minGQ=57(inclusive),maxGQ=58(exclusive)
##GVCFBlock58-59=minGQ=58(inclusive),maxGQ=59(exclusive)
##GVCFBlock59-60=minGQ=59(inclusive),maxGQ=60(exclusive)
##GVCFBlock6-7=minGQ=6(inclusive),maxGQ=7(exclusive)
##GVCFBlock60-70=minGQ=60(inclusive),maxGQ=70(exclusive)
##GVCFBlock7-8=minGQ=7(inclusive),maxGQ=8(exclusive)
##GVCFBlock70-80=minGQ=70(inclusive),maxGQ=80(exclusive)
##GVCFBlock8-9=minGQ=8(inclusive),maxGQ=9(exclusive)
##GVCFBlock80-90=minGQ=80(inclusive),maxGQ=90(exclusive)
##GVCFBlock9-10=minGQ=9(inclusive),maxGQ=10(exclusive)
##GVCFBlock90-99=minGQ=90(inclusive),maxGQ=99(exclusive)
##GVCFBlock99-2147483647=minGQ=99(inclusive),maxGQ=2147483647(exclusive)
##INFO=<ID=AA,Number=1,Type=String,Description="Ancestral Allele, ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/pilot_data/technical/reference/ancestral_alignments/README">
##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">
##INFO=<ID=BaseQRankSum,Number=1,Type=Float,Description="Z-score from Wilcoxon rank sum test of Alt Vs. Ref base qualities">
##INFO=<ID=ClippingRankSum,Number=1,Type=Float,Description="Z-score From Wilcoxon rank sum test of Alt vs. Ref number of hard clipped bases">
##INFO=<ID=DB,Number=0,Type=Flag,Description="dbSNP Membership">
##INFO=<ID=DP,Number=1,Type=Integer,Description="Approximate read depth; some reads may have been filtered">
##INFO=<ID=DS,Number=0,Type=Flag,Description="Were any of the samples downsampled?">
##INFO=<ID=END,Number=1,Type=Integer,Description="Stop position of the interval">
##INFO=<ID=ExcessHet,Number=1,Type=Float,Description="Phred-scaled p-value for exact test of excess heterozygosity">
##INFO=<ID=FS,Number=1,Type=Float,Description="Phred-scaled p-value using Fisher's exact test to detect strand bias">
##INFO=<ID=HaplotypeScore,Number=1,Type=Float,Description="Consistency of the site with at most two segregating haplotypes">
##INFO=<ID=InbreedingCoeff,Number=1,Type=Float,Description="Inbreeding coefficient as estimated from the genotype likelihoods per-sample when compared against the Hardy-Weinberg expectation">
##INFO=<ID=MLEAC,Number=A,Type=Integer,Description="Maximum likelihood expectation (MLE) for the allele counts (not necessarily the same as the AC), for each ALT allele, in the same order as listed">
##INFO=<ID=MLEAF,Number=A,Type=Float,Description="Maximum likelihood expectation (MLE) for the allele frequency (not necessarily the same as the AF), for each ALT allele, in the same order as listed">
##INFO=<ID=MQ,Number=1,Type=Float,Description="RMS Mapping Quality">
##INFO=<ID=MQRankSum,Number=1,Type=Float,Description="Z-score From Wilcoxon rank sum test of Alt vs. Ref read mapping qualities">
##INFO=<ID=NEGATIVE_TRAIN_SITE,Number=0,Type=Flag,Description="This variant was used to build the negative training set of bad variants">
##INFO=<ID=POSITIVE_TRAIN_SITE,Number=0,Type=Flag,Description="This variant was used to build the positive training set of good variants">
##INFO=<ID=QD,Number=1,Type=Float,Description="Variant Confidence/Quality by Depth">
##INFO=<ID=RAW_MQ,Number=1,Type=Float,Description="Raw data for RMS Mapping Quality">
##INFO=<ID=ReadPosRankSum,Number=1,Type=Float,Description="Z-score from Wilcoxon rank sum test of Alt vs. Ref read position bias">
##INFO=<ID=SOR,Number=1,Type=Float,Description="Symmetric Odds Ratio of 2x2 contingency table to detect strand bias">
##INFO=<ID=VQSLOD,Number=1,Type=Float,Description="Log odds of being a true variant versus being false under the trained gaussian mixture model">
##INFO=<ID=culprit,Number=1,Type=String,Description="The annotation which was the worst performing in the Gaussian mixture model, likely the reason why the variant was filtered out">
##INFO=<ID=set,Number=1,Type=String,Description="Source VCF for the merged record in CombineVariants">
##VcfRenameSamplesCmdLine=-o /share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/NO-CALL/COMBINE//PHASE4_snp99.5-CHR_HWE_BA_NOCALL_Combine_RN1.vcf -f /share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/NO-CALL/COMBINE//renamesample.txt -E /share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/VAR-FILTER/NO-CALL/COMBINE//PHASE4_snp99.5-CHR_HWE_BA_NOCALL_Combine.vcf
##VcfRenameSamplesHtsJdkHome=lib/com/github/samtools/htsjdk/2.9.1/htsjdk-2.9.1.jar
##VcfRenameSamplesHtsJdkVersion=2.9.1
##VcfRenameSamplesVersion=02fa35350f66938cbbbf14b590aec49b8f67e514
##contig=<ID=1,length=249250621>
##contig=<ID=2,length=243199373>
##contig=<ID=3,length=198022430>
##contig=<ID=4,length=191154276>
##contig=<ID=5,length=180915260>
##contig=<ID=6,length=171115067>
##contig=<ID=7,length=159138663>
##contig=<ID=8,length=146364022>
##contig=<ID=9,length=141213431>
##contig=<ID=10,length=135534747>
##contig=<ID=11,length=135006516>
##contig=<ID=12,length=133851895>
##contig=<ID=13,length=115169878>
##contig=<ID=14,length=107349540>
##contig=<ID=15,length=102531392>
##contig=<ID=16,length=90354753>
##contig=<ID=17,length=81195210>
##contig=<ID=18,length=78077248>
##contig=<ID=19,length=59128983>
##contig=<ID=20,length=63025520>
##contig=<ID=21,length=48129895>
##contig=<ID=22,length=51304566>
##contig=<ID=X,length=155270560>
##contig=<ID=Y,length=59373566>
##contig=<ID=MT,length=16569>
##contig=<ID=GL000207.1,length=4262>
##contig=<ID=GL000226.1,length=15008>
##contig=<ID=GL000229.1,length=19913>
##contig=<ID=GL000231.1,length=27386>
##contig=<ID=GL000210.1,length=27682>
##contig=<ID=GL000239.1,length=33824>
##contig=<ID=GL000235.1,length=34474>
##contig=<ID=GL000201.1,length=36148>
##contig=<ID=GL000247.1,length=36422>
##contig=<ID=GL000245.1,length=36651>
##contig=<ID=GL000197.1,length=37175>
##contig=<ID=GL000203.1,length=37498>
##contig=<ID=GL000246.1,length=38154>
##contig=<ID=GL000249.1,length=38502>
##contig=<ID=GL000196.1,length=38914>
##contig=<ID=GL000248.1,length=39786>
##contig=<ID=GL000244.1,length=39929>
##contig=<ID=GL000238.1,length=39939>
##contig=<ID=GL000202.1,length=40103>
##contig=<ID=GL000234.1,length=40531>
##contig=<ID=GL000232.1,length=40652>
##contig=<ID=GL000206.1,length=41001>
##contig=<ID=GL000240.1,length=41933>
##contig=<ID=GL000236.1,length=41934>
##contig=<ID=GL000241.1,length=42152>
##contig=<ID=GL000243.1,length=43341>
##contig=<ID=GL000242.1,length=43523>
##contig=<ID=GL000230.1,length=43691>
##contig=<ID=GL000237.1,length=45867>
##contig=<ID=GL000233.1,length=45941>
##contig=<ID=GL000204.1,length=81310>
##contig=<ID=GL000198.1,length=90085>
##contig=<ID=GL000208.1,length=92689>
##contig=<ID=GL000191.1,length=106433>
##contig=<ID=GL000227.1,length=128374>
##contig=<ID=GL000228.1,length=129120>
##contig=<ID=GL000214.1,length=137718>
##contig=<ID=GL000221.1,length=155397>
##contig=<ID=GL000209.1,length=159169>
##contig=<ID=GL000218.1,length=161147>
##contig=<ID=GL000220.1,length=161802>
##contig=<ID=GL000213.1,length=164239>
##contig=<ID=GL000211.1,length=166566>
##contig=<ID=GL000199.1,length=169874>
##contig=<ID=GL000217.1,length=172149>
##contig=<ID=GL000216.1,length=172294>
##contig=<ID=GL000215.1,length=172545>
##contig=<ID=GL000205.1,length=174588>
##contig=<ID=GL000219.1,length=179198>
##contig=<ID=GL000224.1,length=179693>
##contig=<ID=GL000223.1,length=180455>
##contig=<ID=GL000195.1,length=182896>
##contig=<ID=GL000212.1,length=186858>
##contig=<ID=GL000222.1,length=186861>
##contig=<ID=GL000200.1,length=187035>
##contig=<ID=GL000193.1,length=189789>
##contig=<ID=GL000194.1,length=191469>
##contig=<ID=GL000225.1,length=211173>
##contig=<ID=GL000192.1,length=547496>
##contig=<ID=NC_007605,length=171823>
##contig=<ID=hs37d5,length=35477943>
##reference=file:///share_bio/unisvx3/zengchq_group/sohail/Ref/ref-1000G/hs37d5.fa
##source=SelectVariants
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT AAGC022072 AAGC022073 AAGC022074 AAGC022075 AAGC022076 LP6005592-DNA_E05 LP6005592-DNA_F05 SHERPA-001 SHERPA-002
1 10575 . C G 67.27 PASS AA=.;AC=0;AF=0.00;AN=18;DP=85;set=Intersection GT:AD:DP:GQ:PL 0/0:1,0:1:3:0,3,28 0/0:11,0:11:30:0,30,450 0/0:7,0:7:21:0,21,226 0/0:9,0:9:24:0,24,360 0/0:6,0:6:18:0,18,168 0/0:26,0:26:60:0,60,900 0/0:7,0:7:21:0,21,218 0/0:10,0:10:27:0,27,405 0/0:8,0:8:24:0,24,246
1 12882 . C G 660.23 PASS AA=c;AC=0;AF=0.00;AN=18;DP=143;set=Intersection GT:AD:DP:GQ:PL 0/0:7,0:7:21:0,21,253 0/0:11,0:11:33:0,33,377 0/0:12,0:12:36:0,36,417 0/0:23,0:23:66:0,66,990 0/0:16,0:16:13:0,13,470 0/0:29,0:29:81:0,81,839 0/0:24,0:24:63:0,63,945 0/0:15,0:15:45:0,45,481 0/0:6,0:6:18:0,18,193
1 13079 . C G 1594.47 PASS AA=c;AC=3;AF=0.167;AN=18;DP=340;set=Intersection GT:AD:DP:GQ:PL 0/1:10,4:14:78:78,0,226 0/1:32,4:36:9:9,0,744 0/0:29,3:32:4:0,4,710 0/0:37,0:37:99:0,109,1206 0/0:56,4:60:61:0,61,1501 0/0:33,0:33:99:0,99,1155 0/1:51,8:59:41:41,0,1276 0/0:55,0:55:0:0,0,1269 0/0:14,0:14:42:0,42,425
1 17730 . C A 9050.45 PASS AA=-;AC=2;AF=0.111;AN=18;DP=230;set=Intersection GT:AD:DP:GQ:PGT:PID:PL 0/0:28,0:28:48:.:.:0,48,696 0/0:31,0:31:93:.:.:0,93,959 0/0:28,2:30:1:0|1:17722_A_G:0,1,1171 0/0:46,0:46:99:.:.:0,113,1434 0/0:35,0:35:69:.:.:0,69,1054 0/0:27,0:27:46:.:.:0,46,884 0/1:9,4:13:99:0|1:17722_A_G:141,0,617 0/1:7,2:9:63:0|1:17722_A_G:63,0,320 0/0:11,0:11:14:.:.:0,14,274
1 49272 rs370116346 G A 397.31 PASS AA=N;AC=4;AF=0.286;AN=14;DP=180;set=Intersection GT:AD:DP:GQ:PL 1/1:0,3:3:9:95,9,0 0/0:25,2:27:28:0,28,742 0/0:31,0:31:90:0,90,1350 1/1:0,4:4:12:125,12,0 0/0:36,4:40:1:0,1,939 0/0:36,0:36:81:0,81,1215 ./.:0,0:0:.:0,0,0 0/0:39,0:39:83:0,83,1138 ./.:0,0:0:.:0,0,0
1 936210 rs3121569 C A 124552 PASS AA=A;AC=17;AF=0.944;AN=18;DP=224;set=Intersection GT:AD:DP:GQ:PGT:PID:PL 1/1:0,37:37:99:.:.:1098,108,0 1/1:0,23:23:69:1|1:936194_A_G:987,69,0 1/1:0,29:29:85:.:.:860,85,0 1/1:0,28:28:83:.:.:911,83,0 1/1:0,20:20:60:.:.:673,60,0 1/1:0,24:24:75:1|1:936194_A_G:1071,75,0 1/1:0,29:29:87:.:.:1011,87,0 0/1:7,7:14:99:.:.:206,0,144 1/1:1,19:20:34:.:.:514,34,0
##fileformat=VCFv4.1
##fileDate=06082017_14h14m53s
##source=SHAPEIT2.v837
##log_file=/share_bio/unisvx3/zengchq_group/sohail/_PHASE-4/VAR-PROCESS/_PHASING/DENOVO/VCF/_LOG/log//PHASE4-CHR-1.log
##FORMAT=<ID=GT,Number=1,Type=String,Description="Phased Genotype">
##INFO=<ID=AA,Number=1,Type=String,Description="Ancestral Allele, ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/pilot_data/technical/reference/ancestral_alignments/README">
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT 103-RQ 109-EP 16-EJ 31-CE 40-JI 43-AM 50-JB 52-AA
1 10575 . C G . PASS AA=. GT 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0
1 12882 . C G . PASS AA=c GT 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0
1 13079 . C G . PASS AA=c GT 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0
1 13110 . G A . PASS AA=g GT 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0
1 13116 rs201725126 T G . PASS AA=t GT 1|1 0|0 0|0 0|0 1|1 0|0 0|0
1 13118 rs200579949 A G . PASS AA=a GT 1|1 0|0 0|0 0|0 1|1 0|0 0|0
1 13143 . G C . PASS AA=g GT 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0
1 13178 . G A . PASS AA=g GT 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0
1 13273 . G C . PASS AA=g GT 0|0 0|0 1|1 1|1 1|1 0|1 0|1 0|1
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment