##fileformat=VCFv4.1 ##FORMAT= ##FORMAT= ##FORMAT= ##FORMAT= ##FORMAT= ##INFO= ##INFO= ##INFO= ##INFO= ##INFO= ##INFO= ##INFO= ##INFO= ##INFO= ##INFO= ##INFO= ##INFO= ##INFO= ##INFO= ##INFO= ##INFO= ##UnifiedGenotyper="analysis_type=UnifiedGenotyper input_file=[bam/M10478.conc.on.pos.realigned.bam, bam/M10475.conc.on.pos.realigned.bam, bam/M10500.conc.on.pos.realigned.bam, bam/M128215.conc.on.pos.realigned.bam] read_buffer_size=null phone_home=STANDARD read_filter=[] intervals=null excludeIntervals=null interval_set_rule=UNION interval_merging=ALL reference_sequence=/m/cphg-quinlan/cphg-quinlan/shared/genomes/hg19/bwa/gatk/hg19_gatk.fa rodBind=[] nonDeterministicRandomSeed=false downsampling_type=BY_SAMPLE downsample_to_fraction=null downsample_to_coverage=250 baq=OFF baqGapOpenPenalty=40.0 performanceLog=null useOriginalQualities=false defaultBaseQualities=-1 validation_strictness=SILENT unsafe=null num_threads=10 num_cpu_threads=null num_io_threads=null num_bam_file_handles=null read_group_black_list=null pedigree=[] pedigreeString=[] pedigreeValidationType=STRICT allow_intervals_with_unindexed_bam=false logging_level=INFO log_to_file=null help=false genotype_likelihoods_model=BOTH p_nonref_model=EXACT heterozygosity=0.0010 pcr_error_rate=1.0E-4 genotyping_mode=DISCOVERY output_mode=EMIT_VARIANTS_ONLY standard_min_confidence_threshold_for_calling=30.0 standard_min_confidence_threshold_for_emitting=30.0 computeSLOD=false alleles=(RodBinding name= source=UNBOUND) min_base_quality_score=17 max_deletion_fraction=0.05 multiallelic=false max_alternate_alleles=5 min_indel_count_for_genotyping=5 indel_heterozygosity=1.25E-4 indelGapContinuationPenalty=10.0 indelGapOpenPenalty=45.0 indelHaplotypeSize=80 bandedIndel=false indelDebug=false ignoreSNPAlleles=false dbsnp=(RodBinding name= source=UNBOUND) out=org.broadinstitute.sting.gatk.io.stubs.VCFWriterStub NO_HEADER=org.broadinstitute.sting.gatk.io.stubs.VCFWriterStub sites_only=org.broadinstitute.sting.gatk.io.stubs.VCFWriterStub debug_file=null metrics_file=null annotation=[] excludeAnnotation=[] filter_mismatching_base_and_quals=false" ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##contig= ##reference=file:///m/cphg-quinlan/cphg-quinlan/shared/genomes/hg19/bwa/gatk/hg19_gatk.fa ##SnpEffVersion="SnpEff 3.0g (build 2012-08-31), by Pablo Cingolani" ##SnpEffCmd="SnpEff -i vcf -o vcf GRCh37.66 test4.vep.vcf " ##INFO= ##INFO= #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT M10475 M10478 M10500 M128215 chr10 1142208 . T C 3404.3 . AC=8;AF=1.00;AN=8;DP=122;Dels=0.00;FS=0.000;HRun=0;HaplotypeScore=2.6747;MQ=36.00;MQ0=0;QD=27.90;EFF=DOWNSTREAM(MODIFIER||||208|WDR37|protein_coding|CODING|ENST00000436154|),INTRON(MODIFIER||||494|WDR37|protein_coding|CODING|ENST00000263150|),INTRON(MODIFIER||||494|WDR37|protein_coding|CODING|ENST00000358220|),STOP_LOST(HIGH|MISSENSE|Tga/Cga|*250R|249|WDR37|protein_coding|CODING|ENST00000381329|exon_10_1142110_1142566);CSQ=intron_variant|||ENSG00000047056|WDR37|ENST00000263150||||-/494|protein_coding|1,downstream_gene_variant|||ENSG00000047056|WDR37|ENST00000436154||||-/209|protein_coding|1,intron_variant|||ENSG00000047056|WDR37|ENST00000358220||||-/494|protein_coding|1,stop_lost|Tga/Cga|*/R|ENSG00000047056|WDR37|ENST00000381329|9/9|||250/249|protein_coding|1 GT:AD:DP:GQ:PL 1/1:1,37:39:87.16:940,87,0 1/1:0,29:29:78.20:899,78,0 1/1:0,24:24:66.14:729,66,0 1/1:0,30:30:75.18:836,75,0 chr10 48003992 . C T 1047.87 . AC=4;AF=0.50;AN=8;BaseQRankSum=-0.053;DP=165;Dels=0.00;FS=6.377;HRun=0;HaplotypeScore=4.3830;MQ=20.94;MQ0=0;MQRankSum=-0.368;QD=9.53;ReadPosRankSum=1.346;EFF=NON_SYNONYMOUS_CODING(MODERATE|MISSENSE|tGt/tAt|C540Y|610|ASAH2C|protein_coding|CODING|ENST00000420079|exon_10_48003968_48004056),NON_SYNONYMOUS_CODING(MODERATE|MISSENSE|tGt/tAt|C552Y|622|ASAH2C|protein_coding|CODING|ENST00000426610|exon_10_48003968_48004056);CSQ=missense_variant|tGt/tAt|C/Y|ENSG00000072444|ASAH2C|ENST00000420079|16/17|benign(0)|tolerated(1)|542/612|protein_coding|1,missense_variant|tGt/tAt|C/Y|ENSG00000072444|ASAH2C|ENST00000426610|17/18|benign(0)|tolerated(1)|552/622|protein_coding|1 GT:AD:DP:GQ:PL 1/1:0,28:28:56.88:461,57,0 0/1:20,18:38:99:317,0,370 0/1:24,20:44:99:319,0,361 0/0:55,0:55:99:0,105,865 chr10 52004315 . T C 40.11 . AC=4;AF=1.00;AN=4;DP=2;Dels=0.00;FS=0.000;HRun=0;HaplotypeScore=0.0000;MQ=37.00;MQ0=0;QD=20.06;EFF=INTRON(MODIFIER||||240|ASAH2|protein_coding|CODING|ENST00000374028|),INTRON(MODIFIER||||726|ASAH2|protein_coding|CODING|ENST00000329428|),INTRON(MODIFIER||||745|ASAH2|protein_coding|CODING|ENST00000447815|),INTRON(MODIFIER||||780|ASAH2|protein_coding|CODING|ENST00000395526|),UPSTREAM(MODIFIER||||622|ASAH2|protein_coding|CODING|ENST00000443575|);CSQ=intron_variant|||ENSG00000188611|ASAH2|ENST00000329428||||-/726|protein_coding|1,upstream_gene_variant|||ENSG00000188611|ASAH2|ENST00000443575||||-/622|protein_coding|1,intron_variant|||ENSG00000188611|ASAH2|ENST00000447815||||-/745|protein_coding|1,intron_variant|||ENSG00000188611|ASAH2|ENST00000395526||||-/780|protein_coding|1 GT:AD:DP:GQ:PL ./. ./. 1/1:0,1:1:3.01:36,3,0 1/1:0,1:1:3.01:37,3,0 chr10 52497529 . G C 33.61 . AC=4;AF=1.00;AN=4;DP=2;Dels=0.00;FS=0.000;HRun=0;HaplotypeScore=0.0000;MQ=31.58;MQ0=0;QD=16.81;EFF=UPSTREAM(MODIFIER||||160|ASAH2B|protein_coding|CODING|ENST00000185907|),UPSTREAM(MODIFIER||||160|ASAH2B|protein_coding|CODING|ENST00000374007|),UPSTREAM(MODIFIER||||165|ASAH2B|protein_coding|CODING|ENST00000374006|),UPSTREAM(MODIFIER|||||ASAH2B|processed_transcript|CODING|ENST00000483649|);CSQ=upstream_gene_variant|||ENSG00000204147|ASAH2B|ENST00000374007||||-/160|protein_coding|1,upstream_gene_variant|||ENSG00000204147|ASAH2B|ENST00000483649|||||processed_transcript|1,upstream_gene_variant|||ENSG00000204147|ASAH2B|ENST00000374006||||-/165|protein_coding|1,upstream_gene_variant|||ENSG00000204147|ASAH2B|ENST00000185907||||-/160|protein_coding|1 GT:AD:DP:GQ:PL ./. 1/1:0,1:1:3.01:37,3,0 1/1:0,1:1:3.01:30,3,0 ./. chr10 126678092 . G A 89.08 . AC=1;AF=0.13;AN=8;BaseQRankSum=-3.120;DP=185;Dels=0.00;FS=3.490;HRun=0;HaplotypeScore=3.3843;MQ=25.32;MQ0=0;MQRankSum=6.568;QD=2.02;ReadPosRankSum=-5.871;EFF=DOWNSTREAM(MODIFIER||||708|ZRANB1|protein_coding|CODING|ENST00000359653|),DOWNSTREAM(MODIFIER|||||CTBP2|processed_transcript|CODING|ENST00000395705|),STOP_GAINED(HIGH|NONSENSE|Caa/Taa|Q445*|445|CTBP2|protein_coding|CODING|ENST00000337195|exon_10_126676421_126678267),STOP_GAINED(HIGH|NONSENSE|Caa/Taa|Q445*|445|CTBP2|protein_coding|CODING|ENST00000411419|exon_10_126677794_126678267),STOP_GAINED(HIGH|NONSENSE|Caa/Taa|Q445*|445|CTBP2|protein_coding|CODING|ENST00000494626|exon_10_126677794_126678267),STOP_GAINED(HIGH|NONSENSE|Caa/Taa|Q445*|445|CTBP2|protein_coding|CODING|ENST00000531469|exon_10_126677794_126678267),STOP_GAINED(HIGH|NONSENSE|Caa/Taa|Q513*|513|CTBP2|protein_coding|CODING|ENST00000334808|exon_10_126677704_126678267),STOP_GAINED(HIGH|NONSENSE|Caa/Taa|Q985*|985|CTBP2|protein_coding|CODING|ENST00000309035|exon_10_126677507_126678267);CSQ=stop_gained|Caa/Taa|Q/*|ENSG00000175029|CTBP2|ENST00000531469|11/11|||445/445|protein_coding|1,stop_gained|Caa/Taa|Q/*|ENSG00000175029|CTBP2|ENST00000309035|9/9|||985/985|protein_coding|1,downstream_gene_variant|||ENSG00000019995|ZRANB1|ENST00000359653||||-/708|protein_coding|1,stop_gained|Caa/Taa|Q/*|ENSG00000175029|CTBP2|ENST00000494626|11/11|||445/445|protein_coding|1,stop_gained|Caa/Taa|Q/*|ENSG00000175029|CTBP2|ENST00000337195|11/11|||445/445|protein_coding|1,stop_gained|Caa/Taa|Q/*|ENSG00000175029|CTBP2|ENST00000334808|9/9|||513/513|protein_coding|1,stop_gained|Caa/Taa|Q/*|ENSG00000175029|CTBP2|ENST00000411419|11/11|||445/445|protein_coding|1,downstream_gene_variant|||ENSG00000175029|CTBP2|ENST00000395705|||||processed_transcript|1 GT:AD:DP:GQ:PL 0/0:64,3:67:99:0,165,1505 0/0:11,1:12:7.31:0,7,240 0/0:52,10:62:54.97:0,55,1263 0/1:35,9:44:99:125,0,693 chr10 135210791 . T C 65.41 . AC=4;AF=0.50;AN=8;BaseQRankSum=2.054;DP=11;Dels=0.00;FS=0.000;HRun=0;HaplotypeScore=0.2489;MQ=35.12;MQ0=0;MQRankSum=0.248;QD=16.35;ReadPosRankSum=-1.001;EFF=INTRON(MODIFIER||||226|MTG1.1|protein_coding|CODING|ENST00000537620|),INTRON(MODIFIER||||283|MTG1.1|protein_coding|CODING|ENST00000432508|),INTRON(MODIFIER||||334|MTG1.1|protein_coding|CODING|ENST00000317502|),INTRON(MODIFIER||||339|MTG1|protein_coding|CODING|ENST00000468317|),INTRON(MODIFIER|||||MTG1.1|processed_transcript|CODING|ENST00000460848|),INTRON(MODIFIER|||||MTG1.1|processed_transcript|CODING|ENST00000473735|),INTRON(MODIFIER|||||MTG1.1|processed_transcript|CODING|ENST00000477902|),INTRON(MODIFIER|||||MTG1.1|processed_transcript|CODING|ENST00000495014|),INTRON(MODIFIER|||||MTG1.1|processed_transcript|CODING|ENST00000498334|),INTRON(MODIFIER|||||MTG1.1|processed_transcript|CODING|ENST00000498790|),UPSTREAM(MODIFIER|||||MTG1.1|processed_transcript|CODING|ENST00000492266|);CSQ=intron_variant|||ENSG00000148824|MTG1|ENST00000317502||||-/334|protein_coding|1,upstream_gene_variant|||ENSG00000148824|MTG1|ENST00000492266|||||retained_intron|1,intron_variant&nc_transcript_variant|||ENSG00000148824|MTG1|ENST00000460848|||||retained_intron|1,intron_variant|||ENSG00000254536|RP11-108K14.8|ENST00000468317||||-/339|protein_coding|1,intron_variant|||ENSG00000148824|MTG1|ENST00000477902||||-/293|protein_coding|1,intron_variant&nc_transcript_variant|||ENSG00000148824|MTG1|ENST00000473735|||||retained_intron|1,intron_variant&nc_transcript_variant|||ENSG00000148824|MTG1|ENST00000498790|||||retained_intron|1,intron_variant|||ENSG00000148824|MTG1|ENST00000432508||||-/283|protein_coding|1,intron_variant&nc_transcript_variant|||ENSG00000148824|MTG1|ENST00000495014|||||retained_intron|1 GT:AD:DP:GQ:PL 0/0:4,0:4:9:0,9,84 1/1:0,3:3:6.02:74,6,0 1/1:0,1:1:3.01:37,3,0 0/0:3,0:3:9.02:0,9,100 chr10 135336656 . G A 38.34 . AC=4;AF=1.00;AN=4;DP=2;Dels=0.00;FS=0.000;HRun=4;HaplotypeScore=0.0000;MQ=37.00;MQ0=0;QD=19.17;EFF=INTRON(MODIFIER||||151|SPRN|protein_coding|CODING|ENST00000541506|),INTRON(MODIFIER||||493|CYP2E1|protein_coding|CODING|ENST00000463117|),INTRON(MODIFIER||||693|RP11-108K14.4.1|protein_coding|CODING|ENST00000356567|),INTRON(MODIFIER|||||RP11-108K14.4.1|retained_intron|CODING|ENST00000488261|),UPSTREAM(MODIFIER||||305|CYP2E1|protein_coding|CODING|ENST00000418356|),UPSTREAM(MODIFIER||||355|CYP2E1|protein_coding|CODING|ENST00000421586|),UPSTREAM(MODIFIER||||493|CYP2E1|protein_coding|CODING|ENST00000252945|),UPSTREAM(MODIFIER||||85|CYP2E1|protein_coding|CODING|ENST00000541261|),UPSTREAM(MODIFIER|||||CYP2E1|processed_transcript|CODING|ENST00000477500|),UPSTREAM(MODIFIER|||||CYP2E1|processed_transcript|CODING|ENST00000480558|);CSQ=upstream_gene_variant|||ENSG00000130649|CYP2E1|ENST00000418356||||-/306|protein_coding|1,intron_variant&nc_transcript_variant|||ENSG00000214279|RP11-108K14.4|ENST00000488261|||||retained_intron|1,intron_variant|||ENSG00000203772|SPRN|ENST00000541506||||-/151|protein_coding|1,upstream_gene_variant|||ENSG00000130649|CYP2E1|ENST00000541261||||-/86|protein_coding|1,downstream_gene_variant|||ENSG00000268925|AL161645.2|ENST00000599428||||-/20|protein_coding|1,upstream_gene_variant|||ENSG00000130649|CYP2E1|ENST00000421586||||-/356|protein_coding|1,intron_variant|||ENSG00000130649|CYP2E1|ENST00000463117||||-/493|protein_coding|1,upstream_gene_variant|||ENSG00000130649|CYP2E1|ENST00000252945||||-/493|protein_coding|1,upstream_gene_variant|||ENSG00000130649|CYP2E1|ENST00000541080||||-/43|nonsense_mediated_decay|1,upstream_gene_variant|||ENSG00000130649|CYP2E1|ENST00000477500|||||processed_transcript|1,upstream_gene_variant|||ENSG00000130649|CYP2E1|ENST00000480558|||||processed_transcript|1,intron_variant&nc_transcript_variant|||ENSG00000214279|RP11-108K14.4|ENST00000356567|||||unitary_pseudogene|1 GT:AD:DP:GQ:PL ./. 1/1:0,1:1:3.01:37,3,0 ./. 1/1:0,1:1:3.01:34,3,0 chr10 135369532 . T C 122.62 . AC=2;AF=0.25;AN=8;BaseQRankSum=2.118;DP=239;Dels=0.00;FS=5.194;HRun=2;HaplotypeScore=5.7141;MQ=36.02;MQ0=0;MQRankSum=0.082;QD=2.31;ReadPosRankSum=-0.695;EFF=DOWNSTREAM(MODIFIER|||||SYCE1|processed_transcript|CODING|ENST00000460441|),DOWNSTREAM(MODIFIER|||||SYCE1|processed_transcript|CODING|ENST00000482127|),EXON(MODIFIER|||||CYP2E1|retained_intron|CODING|ENST00000368520|),EXON(MODIFIER|||||SYCE1|processed_transcript|CODING|ENST00000479535|),INTRON(MODIFIER||||151|SPRN|protein_coding|CODING|ENST00000541506|),NON_SYNONYMOUS_CODING(MODERATE|MISSENSE|aAg/aGg|K147R|282|SYCE1|protein_coding|CODING|ENST00000368517|exon_10_135369485_135369551),NON_SYNONYMOUS_CODING(MODERATE|MISSENSE|aAg/aGg|K147R|282|SYCE1|protein_coding|CODING|ENST00000432597|exon_10_135369485_135369551),NON_SYNONYMOUS_CODING(MODERATE|MISSENSE|aAg/aGg|K183R|318|SYCE1|protein_coding|CODING|ENST00000303903|exon_10_135369485_135369551),NON_SYNONYMOUS_CODING(MODERATE|MISSENSE|aAg/aGg|K183R|351|SYCE1|protein_coding|CODING|ENST00000343131|exon_10_135369485_135369551);CSQ=missense_variant|aAg/aGg|K/R|ENSG00000171772|SYCE1|ENST00000303903|9/13|benign(0.001)|tolerated(1)|183/318|protein_coding|1,missense_variant|aAg/aGg|K/R|ENSG00000171772|SYCE1|ENST00000432597|10/14|benign(0.001)|tolerated(1)|147/282|protein_coding|1,missense_variant|aAg/aGg|K/R|ENSG00000171772|SYCE1|ENST00000343131|9/13|benign(0.002)|tolerated(1)|183/351|protein_coding|1,intron_variant|||ENSG00000203772|SPRN|ENST00000541506||||-/151|protein_coding|1,non_coding_exon_variant&nc_transcript_variant|||ENSG00000130649|CYP2E1|ENST00000368520|6/6||||retained_intron|1,missense_variant|aAg/aGg|K/R|ENSG00000171772|SYCE1|ENST00000368517|9/13|benign(0.001)|tolerated(1)|147/282|protein_coding|1,non_coding_exon_variant&nc_transcript_variant|||ENSG00000171772|SYCE1|ENST00000479535|7/11||||retained_intron|1 GT:AD:DP:GQ:PL 0/0:40,0:40:63.06:0,63,628 0/1:27,5:32:49.03:49,0,658 0/1:15,5:21:99:115,0,261 0/0:145,0:146:99:0,319,3399 chr16 72057435 . C T 572.98 . AC=1;AF=0.13;AN=8;BaseQRankSum=-2.270;DP=260;Dels=0.00;FS=0.000;HRun=0;HaplotypeScore=4.5319;MQ=36.53;MQ0=0;MQRankSum=0.596;QD=8.07;ReadPosRankSum=0.927;EFF=NON_SYNONYMOUS_CODING(MODERATE|MISSENSE|Cgg/Tgg|R346W|395|DHODH|protein_coding|CODING|ENST00000219240|exon_16_72057373_72057532);CSQ=missense_variant|Cgg/Tgg|R/W|ENSG00000102967|DHODH|ENST00000219240|8/9|probably_damaging(0.956)|deleterious(0)|346/395|protein_coding|1,non_coding_exon_variant&nc_transcript_variant|||ENSG00000102967|DHODH|ENST00000571392|3/4||||retained_intron|1,downstream_gene_variant|||ENSG00000102967|DHODH|ENST00000572003|||||retained_intron|1,downstream_gene_variant|||ENSG00000102967|DHODH|ENST00000573843|||||retained_intron|1,downstream_gene_variant|||ENSG00000102967|DHODH|ENST00000573922|||||processed_transcript|1,intron_variant|||ENSG00000102967|DHODH|ENST00000574309||||-/193|protein_coding|1,missense_variant|Cgg/Tgg|R/W|ENSG00000102967|DHODH|ENST00000572887|8/9|probably_damaging(0.946)|deleterious(0)|344/393|protein_coding|1 GT:AD:DP:GQ:PL 0/1:37,33:71:99:609,0,666 0/0:56,0:56:99:0,135,1543 0/0:67,1:68:99:0,120,1460 0/0:65,0:65:99:0,180,1982