...
| Code Block |
|---|
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT FBS1-LNCAP-RNA chr1 14542 rs1045833 A G 71.84 . AC=2;AF=1.00;AN=2;DB;DP=3;ExcessHet=3.0103;FS=0.000;MLEAC=1;MLEAF=0.500;MQ=23.70;QD=23.95;SOR=2.833;ANN=G|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*133A>G|||||133|,G|downstream_gene_variant|MODIFIER|MIR68MIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*2827T>C|||||2827|,G|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|10/10|n.1254-41T>C|||||| GT:AD:DP:GQ:PL 1/1:0,3:3:9:85,9,0 chr1 14574 rs28503599 A G 121.05 . AC=2;AF=1.00;AN=2;DB;DP=10;ExcessHet=3.0103;FS=0.000;MLEAC=2;MLEAF=1.00;MQ=22.52;QD=12.10;SOR=4.804;ANN=G|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*165A>G|||||165|,G|downstream_gene_variant|MODIFIER|MIR68MIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*2795T>C|||||2795|,G|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|10/10|n.1254-73T>C|||||| GT:AD:DP:GQ:PL 1/1:0,10:10:27:135,27,0 chr1 14599 rs531646671 T A 301.02 . AC=2;AF=1.00;AN=2;DB;DP=7;ExcessHet=3.0103;FS=0.000;MLEAC=2;MLEAF=1.00;MQ=22.00;QD=25.36;SOR=4.174;ANN=A|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*190T>A|||||190|,A|downstream_gene_variant|MODIFIER|MIR685 chr1 14604 rs541940975 AMIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*2770A>T|||||2770|,A|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|10/10|n.1254-98A>T|||||| GGT:AD:DP:GQ:PGT:PID:PL:PS 1|1:0,7:7:21:1|1:14599_T_A:315,21,0:14599 chr1 14604 rs541940975 A G 356.05 . AC=2;AF=1.00;AN=2;DB;DP=9;ExcessHet=3.0103;FS=0.000;MLEAC=2;MLEAF=1.00;MQ=22.23;QD=28.73;SOR=3.056;ANN=G|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*195A>G|||||195|,G|downstream_gene_variant|MODIFIER|MIR685 chr1 14610 . T C 356.05 . MIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*2765T>C|||||2765|,G|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|10/10|n.1254-103T>C|||||| GT:AD:DP:GQ:PGT:PID:PL:PS 1|1:0,9:9:27:1|1:14599_T_A:370,27,0:14599 chr1 14610 . T C 356.05 . AC=2;AF=1.00;AN=2;DP=9;ExcessHet=3.0103;FS=0.000;MLEAC=2;MLEAF=1.00;MQ=22.23;QD=30.97;SOR=3.056;ANN=C|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*201T>C|||||201|,C|downstream_gene_variant|MODIFIER|MIR6859-1|ENSG000 chr11|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*2759A>G|||||2759|,C|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|10/10|n.1254-109A>G|||||| 14653 rs62635297GT:AD:DP:GQ:PGT:PID:PL:PS 1|1:0,9:9:27:1|1:14599_T_A:370,27,0:14599 chr1 14653 rs62635297 C T 1816.06 . AC=2;AF=1.00;AN=2;DB;DP=67;ExcessHet=3.0103;FS=0.000;MLEAC=2;MLEAF=1.00;MQ=34.03;QD=27.94;SOR=2.412;ANN=T|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*244C>T|||||244|,T|downstream_gene_variant|MODIFIER|MIR68 chr1 14677 rs201327123 G MIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*2716G>A|||||2716|,T|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|10/10|n.1254-152G>A|||||| GT:AD:DP:GQ:PL 1/1:0,65:65:99:1830,195,0 chr1 14677 rs201327123 G A 926.64 . AC=1;AF=0.500;AN=2;BaseQRankSum=2.83;DB;DP=113;ExcessHet=3.0103;FS=1.545;MLEAC=1;MLEAF=0.500;MQ=37.82;MQRankSum=-1.398e+00;QD=8.35;ReadPosRankSum=0.962;SOR=0.852;ANN=A|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcrip chr1transcript||n.*268G>A|||||268|,A|downstream_gene_variant|MODIFIER|MIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*2692C>T|||||2692|,A|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|10/10|n.1254-176C>T|||||| 16257 GT:AD:DP:GQ:PL 0/1:63,48:111:99:934,0,1266 chr1 16257 rs11489794 G C 49.64 . AC=1;AF=0.500;AN=2;BaseQRankSum=-1.136e+00;DB;DP=23;ExcessHet=3.0103;FS=0.000;MLEAC=1;MLEAF=0.500;MQ=34.10;MQRankSum=0.381;QD=2.16;ReadPosRankSum=0.00;SOR=1.022;ANN=C|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript chr1 16288 rs113141985 C G 48.64 . AC=1;AF=0.500;AN=2;BaseQRankSum=-1.150e+00;||n.*1848G>C|||||1848|,C|downstream_gene_variant|MODIFIER|MIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*1112C>G|||||1112|,C|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|8/10|n.1068-310C>G|||||| GT:AD:DP:GQ:PL 0/1:18,5:23:57:57,0,384 chr1 16288 rs113141985 C G 48.64 . AC=1;AF=0.500;AN=2;BaseQRankSum=-1.150e+00;DB;DP=9;ExcessHet=3.0103;FS=0.000;MLEAC=1;MLEAF=0.500;MQ=32.24;MQRankSum=-1.350e-01;QD=6.08;ReadPosRankSum=-1.029e+00;SOR=0.307;ANN=G|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_ chr1 16298 rs62636498 C T 112.14 . AC=2;AF=1.00;AN=2;DB;DP=4;ExcessHet=3.0103;FS=0.000;MLEAC=2;MLEAF=1.00;MQ=22.25;QD=28.04;SOR=3.258;ANN=Ttranscript||n.*1879C>G|||||1879|,G|downstream_gene_variant|MODIFIER|DDX11L1MIR6859-1|ENSG00000223972ENSG00000278267|transcript|ENST00000456328ENST00000619216.21|processed_transcriptmiRNA||n.*1889C>T1081G>C|||||18891081|,TG|downstream_geneintron_variant|MODIFIER|MIR6 chr1 136048 rs371677125|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|8/10|n.1067+319G>C|||||| GT:AD:DP:GQ:PL 0/1:5,3:8:56:56,0,107 chr1 16298 rs62636498 C T 125112.14 . AC=2;AF=1.00;AN=2;DB;DP=4;ExcessHet=3.0103;FS=0.000;MLEAC=2;MLEAF=1.00;MQ=4022.0025;QD=3128.2904;SOR=03.693258;ANN=T|upstreamdownstream_gene_variant|MODIFIER|RP11-34P13.15|ENSG00000268903DDX11L1|ENSG00000223972|transcript|ENST00000494149ENST00000456328.2|processed_pseudogenetranscript||n.-153G>A*1889C>T|||||1531889|,T|downstream_gene_variant|MODIFIER|MIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*1071G>A|||||1071|,T|intron_variant|MODIFIER|RP chr1 136573 .|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|8/10|n.1067+309G>A|||||| GT:AD:DP:GQ:PL T 1/1:0,4:4:12:126,12,0 chr1 136048 rs371677125 C T 113125.9714 . AC=2;AF=1.00;AN=2;DB;DP=74;ExcessHet=3.0103;FS=0.000;MLEAC=2;MLEAF=1.00;MQ=4940.6200;QD=1931.0029;SOR=10.329693;ANN=CT|TFupstream_binding_sitegene_variant|LOWMODIFIER|RP11-34P13.15|ENSG00000268903|Egr1transcript|MA0162ENST00000494149.2|processed_pseudogene||n.136573A>G-153G>A|||||153|,CT|upstreamdownstream_gene_variant|MODIFIER|RP11-34P13.1514|ENSG00000268903ENSG00000239906|transcript|ENST00000494149.2ENST00000493797.1|antisense||n.*3742G>A|||||3742|,T|downstream_gene_variant|MODIFIER|CICP27|ENSG00000233750|transcript|ENST00000442987.3|processed_pseudogene||n chr1 136962 rs373582709 C T 118.68 ..*1212C>T|||||1212|,T|downstream_gene_variant|MODIFIER|RP11-34P13.16|ENSG00000269981|transcript|ENST00000595919.1|processed_pseudogene||n.*1634G>A|||||1634|,T|intergenic_region|MODIFIER|RP11-34P13.15-RP11-34P13.16|ENSG00000268903-ENSG00000269981|intergenic_region|ENSG00000268903-ENSG00000269981|||n.136048C>T|||||| AC=1;AF=0.500;AN=2;BaseQRankSum=-8.870e-01;DB;DP=8GT:AD:DP:GQ:PGT:PID:PL:PS 1|1:0,4:4:12:1|1:136048_C_T:139,12,0:136048 chr1 136573 . T C 113.97 . AC=2;AF=1.00;AN=2;DP=7;ExcessHet=3.0103;FS=0.000;MLEAC=12;MLEAF=01.50000;MQ=5049.36;MQRankSum=0.489;62;QD=1619.95;ReadPosRankSum=0.489;00;SOR=01.330329;ANN=TC|TF_binding_site_variant|MODIFIERLOW|||Egr1|MA0162.2|||n.136962G>A136573A>G||||||,TC|upstream_gene_variant|MODIFIER|RP1 chr1 137622 rs376555721 G A 54.66 . AC=1;AF=0.500;AN=2;BaseQRankSum=0.674;DB;DP=4;ExcessHet=3.0103;FS=6.021;MLEAC=1;MLEAF=0.500;MQ=20.00;MQRankSum=0.00;QD=13.67;ReadPosRankSum=-3.190e-01;SOR=2.788;ANN=A|upstreamRP11-34P13.15|ENSG00000268903|transcript|ENST00000494149.2|processed_pseudogene||n.-678A>G|||||678|,C|downstream_gene_variant|MODIFIER|RP11-34P13.1514|ENSG00000268903ENSG00000239906|transcript|ENST00000494149.2|processed_pseudo chr1 137825 rs147252685 G A 622.06 . AC=2;AF=1.00;AN=2;DB;DP=26;ExcessHet=3.0103;FS=0.000;MLEAC=2;MLEAF=1.00;MQ=36.46;QD=24.88;SOR=0.770;ANN=A|upstreamENST00000493797.1|antisense||n.*3217A>G|||||3217|,C|downstream_gene_variant|MODIFIER|RP11-34P13.13|ENSG00000241860|transcript|ENST00000484859.1|antisense||n.*4901A>G|||||4901|,C|downstream_gene_variant|MODIFIER|RP11-34P13.15|ENSG00000268903CICP27|ENSG00000233750|transcript|ENST00000494149ENST00000442987.23|processed_pseudogene||n.-1930C>T*1737T>C|||||19301737|,AC|downstream_gene_variant|MODIFIER chr1 138156 rs370691115 G T |RP11-34P13.16|ENSG00000269981|transcript|ENST00000595919.1|processed_pseudogene||n.*1109A>G|||||1109|,C|intergenic_region|MODIFIER|RP11-34P13.15-RP11-34P13.16|ENSG00000268903-ENSG00000269981|intergenic_region|ENSG00000268903-ENSG00000269981|||n.136573T>C|||||| GT:AD:DP:GQ:PL 1/1:0,6:6:18:128,18,0 chr1 136962 rs373582709 C 197.04 T 118.68 . AC=21;AF=10.00500;AN=2;BaseQRankSum=-8.870e-01;DB;DP=8;ExcessHet=3.0103;FS=0.000;MLEAC=21;MLEAF=10.00500;MQ=22.00;QD=24.63=50.36;MQRankSum=0.489;QD=16.95;ReadPosRankSum=0.489;SOR=0.693330;ANN=T|TF_binding_site_variant|MODIFIER|||Egr1|MA0162.2|||n.136962G>A||||||,T|upstream_gene_variant|MODIFIER|RP11-34P13.15|ENSG00000268903|transcript|ENST00000494149.2|processed_pseudogene||n.-2261C>A1067G>A|||||22611067|,T|upstreamdownstream_gene_variant|MODIFIER|RP chr1 138817 rs556938922 T C 945.06 . AC=2;AF=1.00;AN=2;DB;DP=45;ExcessHet=3.0103;FS=0.000;MLEAC=2;MLEAF=1.00;MQ=35.16;QD=21.98;SOR=0.739;ANN=C|upstreamRP11-34P13.14|ENSG00000239906|transcript|ENST00000493797.1|antisense||n.*2828G>A|||||2828|,T|downstream_gene_variant|MODIFIER|RP11-34P13.13|ENSG00000241860|transcript|ENST00000484859.1|antisense||n.*4512G>A|||||4512|,T|downstream_gene_variant|MODIFIER|RP11-34P13.15|ENSG00000268903CICP27|ENSG00000233750|transcript|ENST00000494149ENST00000442987.23|processed_pseudogene||n.-2922A>G*2126C>T|||||29222126|,CT|upstreamdownstream_gene_variant|MODIFIER|R chr1 184246 . TRP11-34P13.16|ENSG00000269981|transcript|ENST00000595919.1|processed_pseudogene||n.*720G>A|||||720|,T|intergenic_region|MODIFIER|RP11-34P13.15-RP11-34P13.16|ENSG00000268903-ENSG00000269981|intergenic_region|ENSG00000268903-ENSG00000269981|||n.136962C>T|||||| GT:AD:DP:GQ:PL 0/1:1,6:8:10:126,0,10 chr1 137622 rs376555721 G CA 5254.8466 . AC=21;AF=10.00500;AN=2;BaseQRankSum=0.674;DB;DP=4;ExcessHet=3.0103;FS=06.000021;MLEAC=1;MLEAF=0.500;MQ=41.80;QD=17.61=20.00;MQRankSum=0.00;QD=13.67;ReadPosRankSum=-3.190e-01;SOR=12.179788;ANN=CA|downstreamupstream_gene_variant|MODIFIER|FO538757RP11-34P13.215|ENSG00000279928ENSG00000268903|transcript|ENST00000624431ENST00000494149.12|proteinprocessed_codingpseudogene||cn.*88T>C-1727C>T|||||881727|,CA|downstream_gene_variant|MODIFIER|FO538757RP11-34P13.1|ENSG000002 chr1 185194 . G C 91.64 . AC=1;AF=0.500;AN=2;BaseQRankSum=-1.345e+00;DP=10;ExcessHet=3.0103;FS=4.260;MLEAC=1;MLEAF=0.500;MQ=24.41;MQRankSum=-1.150e+00;QD=11.45;ReadPosRankSum=0.00;SOR=2.833;ANN=C|3_prime_UTR_variant|MODIFIER|FO538757.1|ENSG00000279457|transcript|ENST00000623083.3|protein_coding|11/11|c.*23 |
Adding MAF information using bcftools
The following command line adds MAF information to either annotated or non-annotated VCF files.
| Code Block |
|---|
#non-annotated file:
bcftools +fill-tags HaplotypeCaller_sampleID.vcf > HaplotypeCaller_sampleID_tags.vcf
#VEP annotated file:
bcftools +fill-tags HaplotypeCaller_sampleID_VEP.vcf > HaplotypeCaller_sampleID_VEP_tags.vcf
#snpEff annotated file:
bcftools +fill-tags HaplotypeCaller_sampleID_snpEff.vcf > HaplotypeCaller_sampleID_snpEff_tags.vcf |
For the example shown above, find the added MAF information to the VCF annotated using snpEff.
| Code Block |
|---|
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT FBS1-LNCAP-RNA chr1 14542 rs1045833 A G 71.84 . AC=2;AF=1;AN=2;DB;DP=314|ENSG00000239906|transcript|ENST00000493797.1|antisense||n.*2168C>T|||||2168|,A|downstream_gene_variant|MODIFIER|RP11-34P13.13|ENSG00000241860|transcript|ENST00000484859.1|antisense||n.*3852C>T|||||3852|,A|downstream_gene_variant|MODIFIER|CICP27|ENSG00000233750|transcript|ENST00000442987.3|processed_pseudogene||n.*2786G>A|||||2786|,A|downstream_gene_variant|MODIFIER|RP11-34P13.16|ENSG00000269981|transcript|ENST00000595919.1|processed_pseudogene||n.*60C>T|||||60|,A|intergenic_region|MODIFIER|RP11-34P13.15-RP11-34P13.16|ENSG00000268903-ENSG00000269981|intergenic_region|ENSG00000268903-ENSG00000269981|||n.137622G>A|||||| GT:AD:DP:GQ:PL 0/1:1,3:4:13:62,0,13 chr1 137825 rs147252685 G A 622.06 . AC=2;AF=1.00;AN=2;DB;DP=26;ExcessHet=3.0103;FS=0.000;MLEAC=12;MLEAF=01.500;MQ=2336.746;QD=2324.9588;SOR=20.833770;ANN=GA|downstreamupstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972RP11-34P13.15|ENSG00000268903|transcript|ENST00000456328ENST00000494149.2|processed_transcriptpseudogene||n.*133A>G-1930C>T|||||1331930|,GA|downstream_gene_variant|MODIFIER|MIR6859-1|ENSG0 chr1 14574 rs28503599 A G 121.05 . AC=2;AF=1;AN=2;DB;DP=10;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=22.52;QD=12.1;SOR=4.804;ANN=GRP11-34P13.14|ENSG00000239906|transcript|ENST00000493797.1|antisense||n.*1965C>T|||||1965|,A|downstream_gene_variant|MODIFIER|RP11-34P13.13|ENSG00000241860|transcript|ENST00000484859.1|antisense||n.*3649C>T|||||3649|,A|downstream_gene_variant|MODIFIER|DDX11L1CICP27|ENSG00000223972ENSG00000233750|transcript|ENST00000456328ENST00000442987.23|processed_transcriptpseudogene||n.*165A>G2989G>A|||||1652989|,G|downstream_geneA|non_coding_transcript_exon_variant|MODIFIER|MIR6859-1|ENSG00RP11-34P13.16|ENSG00000269981|transcript|ENST00000595919.1|processed_pseudogene|1/1|n.141C>T|||||| GT:AD:DP:GQ:PL 1/1:0,25:25:75:636,75,0 chr1 14599138156 rs370691115 rs531646671 G T AT 301197.0204 . AC=2;AF=1.00;AN=2;DB;DP=78;ExcessHet=3.0103;FS=0.000;MLEAC=2;MLEAF=1.00;MQ=22.00;QD=2524.3663;SOR=40.174693;ANN=AT|downstreamupstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972RP11-34P13.15|ENSG00000268903|transcript|ENST00000456328ENST00000494149.2|processed_transcriptpseudogene||n.*190T>A-2261C>A|||||1902261|,AT|downstreamupstream_gene_variant|MODIFIER|MIR6859-1|ENSG00000 chr1 14604 rs541940975 A G 356.05 . AC=2;AF=1;AN=2;DB;DP=9;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=22.23;QD=28.73;SOR=3.056;ANN=GRP11-34P13.16|ENSG00000269981|transcript|ENST00000595919.1|processed_pseudogene||n.-191C>A|||||191|,T|downstream_gene_variant|MODIFIER|RP11-34P13.14|ENSG00000239906|transcript|ENST00000493797.1|antisense||n.*1634C>A|||||1634|,T|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972RP11-34P13.13|ENSG00000241860|transcript|ENST00000456328ENST00000484859.21|processed_transcriptantisense||n.*195A>G3318C>A|||||1953318|,GT|downstream_gene_variant|MODIFIER|MIR6859-1|ENSG00 chr1 14610 . CICP27|ENSG00000233750|transcript|ENST00000442987.3|processed_pseudogene||n.*3320G>T|||||3320|,T|intergenic_region|MODIFIER|RP11-34P13.16-RP11-34P13.14|ENSG00000269981-ENSG00000239906|intergenic_region|ENSG00000269981-ENSG00000239906|||n.138156G>T|||||| GT:AD:DP:GQ:PL 1/1:0,8:8:24:211,24,0 chr1 138817 rs556938922 T C 356945.0506 . AC=2;AF=1.00;AN=2;DB;DP=945;ExcessHet=3.0103;FS=0.000;MLEAC=2;MLEAF=1.00;MQ=2235.2316;QD=3021.9798;SOR=30.056739;ANN=C|downstreamupstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972RP11-34P13.15|ENSG00000268903|transcript|ENST00000456328ENST00000494149.2|processed_pseudogene||n.-2922A>G|||||2922|,C|upstream_gene_variant|MODIFIER|RP11-34P13.16|ENSG00000269981|transcript|ENST00000595919.1|processed_pseudogene||n.*201T>C-852A>G|||||201852|,C|downstream_gene_variant|MODIFIER|MIR6859RP11-134P13.14|ENSG00000278267|t chr1 14653 rs62635297 C T 1816.06 . AC=2;AF=1;AN=2;DB;DP=67;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=34.03;QD=27.94;SOR=2.412;ANN=TENSG00000239906|transcript|ENST00000493797.1|antisense||n.*973A>G|||||973|,C|downstream_gene_variant|MODIFIER|RP11-34P13.13|ENSG00000241860|transcript|ENST00000484859.1|antisense||n.*2657A>G|||||2657|,C|downstream_gene_variant|MODIFIER|DDX11L1CICP27|ENSG00000223972ENSG00000233750|transcript|ENST00000456328ENST00000442987.23|processed_transcriptpseudogene||n.*244C>T3981T>C|||||2443981|,TC|downstreamintergenic_gene_variantregion|MODIFIER|MIR6859-1|ENSG0 chr1 14677 rs201327123|RP11-34P13.16-RP11-34P13.14|ENSG00000269981-ENSG00000239906|intergenic_region|ENSG00000269981-ENSG00000239906|||n.138817T>C|||||| GT:AD:DP:GQ:PL 1/1:0,43:43:99:959,126,0 chr1 184246 . GT AC 92652.6484 . AC=12;AF=01.500;AN=2;BaseQRankSum=2.83;DB;DP=1134;ExcessHet=3.0103;FS=10.545000;MLEAC=1;MLEAF=0.5500;MQ=3741.8280;MQRankSumQD=-117.398;QD=8.35;ReadPosRankSum=0.962;61;SOR=01.852179;ANN=AC|downstream_gene_variant|MODIFIER|DDX11L1FO538757.2|ENSG00000223972ENSG00000279928|transcript|ENST00000456328ENST00000624431.21|processedprotein_transcriptcoding||n.*26 chr1 16257 rs11489794 G C 49.64 . AC=1;AF=0.5;AN=2;BaseQRankSum=-1.136;DB;DP=23;ExcessHet=3.0103;FS=0;MLEAC=1;MLEAF=0.5;MQ=34.1;MQRankSum=0.381;QD=2.16;ReadPosRankSum=0;SOR=1.022;ANN=C|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*1848G>C|||| chr1 16288 rs113141985 Cc.*88T>C|||||88|,C|downstream_gene_variant|MODIFIER|FO538757.1|ENSG00000279457|transcript|ENST00000623083.3|protein_coding||c.*971A>G|||||679|,C|downstream_gene_variant|MODIFIER|MIR6859-2|ENSG00000273874|transcript|ENST00000612080.1|miRNA||n.*3645A>G|||||3645|,C|intergenic_region|MODIFIER|FO538757.2-FO538757.1|ENSG00000279928-ENSG00000279457|intergenic_region|ENSG00000279928-ENSG00000279457|||n.184246T>C|||||| GT:AD:DP:GQ:PL 1/1:0,3:3:9:66,9,0 chr1 185194 . G C 48 91.64 . AC=1;AF=0.5500;AN=2;BaseQRankSum=-1.15;DB345e+00;DP=910;ExcessHet=3.0103;FS=04.260;MLEAC=1;MLEAF=0.5500;MQ=3224.2441;MQRankSum=-01.135150e+00;QD=611.0845;ReadPosRankSum=-10.02900;SOR=02.307833;ANN=GC|downstream3_prime_geneUTR_variant|MODIFIER|DDX11L1FO538757.1|ENSG00000223972ENSG00000279457|transcript|ENST00000456328ENST00000623083.23|processedprotein_transcriptcoding|11/11|n.*1879C> chr1c.*23C>G|||||23|,C|downstream_gene_variant|MODIFIER|FO538757.2|ENSG00000279928|transcript|ENST00000624431.1|protein_coding||c.*1036G>C|||||1036|,C|downstream_gene_variant|MODIFIER|MIR6859-2|ENSG00000273874|transcript|ENST00000612080.1|miRNA||n.*2697C>G|||||2697| 16298 rs62636498 C T 112.14 . AC=2;AF=1;AN=2;DB;DP=4;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=22.25;QD=28.04;SOR=3.258;ANN=T|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*1889C>T|||||1889|,T|downstream_gene_variant|MODIFIER|MIR6859-1|ENSG chr1 136048 rs371677125 C TGT:AD:DP:GQ:PL 0/1:3,5:8:65:99,0,65 |
Adding MAF information using bcftools
The following command line adds MAF information to either annotated or non-annotated VCF files.
| Code Block |
|---|
#non-annotated file:
bcftools +fill-tags HaplotypeCaller_sampleID.vcf > HaplotypeCaller_sampleID_tags.vcf
#VEP annotated file:
bcftools +fill-tags HaplotypeCaller_sampleID_VEP.vcf > HaplotypeCaller_sampleID_VEP_tags.vcf
#snpEff annotated file:
bcftools +fill-tags HaplotypeCaller_sampleID_snpEff.vcf > HaplotypeCaller_sampleID_snpEff_tags.vcf |
For the example shown above, find the added MAF information to the VCF annotated using snpEff.
| Code Block |
|---|
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT FBS1-LNCAP-RNA chr1 14542 rs1045833 A G 12571.1484 . AC=2;AF=1;AN=2;DB;DP=43;ExcessHet=3.0103;FS=0;MLEAC=21;MLEAF=10.5;MQ=4023.7;QD=3123.2995;SOR=02.693833;ANN=TG|upstreamdownstream_gene_variant|MODIFIER|RP11-34P13.15|ENSG00000268903DDX11L1|ENSG00000223972|transcript|ENST00000494149ENST00000456328.2|processed_pseudogenetranscript||n.-153G>A*133A>G|||||153133|,TG|downstream_gene_variant|MODIFIER|RP11-34P13.14|E chr1 136573 . T C 113.97 . AC=2;AF=1;AN=2;DP=7;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=49.62;QD=19;SOR=1.329;ANN=C|TF_binding_site_variant|LOW|||Egr1|MA0162.2|||n.136573A>G||||||,C|upstream_gene_variant|MODIFIER|RP11-34P13.15|ENSG00000268903|transcript|ENST00000494149.2|processed_pseudogene||n.-678A>G||||| chr1 136962 rs373582709 CMIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*2827T>C|||||2827|,G|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|10/10|n.1254-41T>C||||||;NS=1;MAF=0;AC_Het=0;AC_Hom=2;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PL 1/1:0,3:3:9:85,9,0 chr1 14574 rs28503599 A G 121.05 . AC=2;AF=1;AN=2;DB;DP=10;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=22.52;QD=12.1;SOR=4.804;ANN=G|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*165A>G|||||165|,G|downstream_gene_variant|MODIFIER|MIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*2795T>C|||||2795|,G|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|10/10|n.1254-73T>C||||||;NS=1;MAF=0;AC_Het=0;AC_Hom=2;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PL 1/1:0,10:10:27:135,27,0 chr1 14599 rs531646671 T A 301.02 . AC=2;AF=1;AN=2;DB;DP=7;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=22;QD=25.36;SOR=4.174;ANN=A|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*190T>A|||||190|,A|downstream_gene_variant|MODIFIER|MIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*2770A>T|||||2770|,A|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|10/10|n.1254-98A>T||||||;NS=1;MAF=0;AC_Het=0;AC_Hom=2;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PGT:PID:PL:PS 1|1:0,7:7:21:1|1:14599_T_A:315,21,0:14599 chr1 14604 rs541940975 A G 356.05 . AC=2;AF=1;AN=2;DB;DP=9;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=22.23;QD=28.73;SOR=3.056;ANN=G|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*195A>G|||||195|,G|downstream_gene_variant|MODIFIER|MIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*2765T>C|||||2765|,G|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|10/10|n.1254-103T>C||||||;NS=1;MAF=0;AC_Het=0;AC_Hom=2;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PGT:PID:PL:PS 1|1:0,9:9:27:1|1:14599_T_A:370,27,0:14599 chr1 14610 . T C 356.05 . AC=2;AF=1;AN=2;DP=9;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=22.23;QD=30.97;SOR=3.056;ANN=C|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*201T>C|||||201|,C|downstream_gene_variant|MODIFIER|MIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*2759A>G|||||2759|,C|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|10/10|n.1254-109A>G||||||;NS=1;MAF=0;AC_Het=0;AC_Hom=2;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PGT:PID:PL:PS 1|1:0,9:9:27:1|1:14599_T_A:370,27,0:14599 chr1 14653 rs62635297 C T 1816.06 . AC=2;AF=1;AN=2;DB;DP=67;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=34.03;QD=27.94;SOR=2.412;ANN=T|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*244C>T|||||244|,T|downstream_gene_variant|MODIFIER|MIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*2716G>A|||||2716|,T|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|10/10|n.1254-152G>A||||||;NS=1;MAF=0;AC_Het=0;AC_Hom=2;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PL 1/1:0,65:65:99:1830,195,0 chr1 14677 rs201327123 G A 926.64 . AC=1;AF=0.5;AN=2;BaseQRankSum=2.83;DB;DP=113;ExcessHet=3.0103;FS=1.545;MLEAC=1;MLEAF=0.5;MQ=37.82;MQRankSum=-1.398;QD=8.35;ReadPosRankSum=0.962;SOR=0.852;ANN=A|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*268G>A|||||268|,A|downstream_gene_variant|MODIFIER|MIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*2692C>T|||||2692|,A|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|10/10|n.1254-176C>T||||||;NS=1;MAF=0.5;AC_Het=1;AC_Hom=0;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PL 0/1:63,48:111:99:934,0,1266 chr1 16257 rs11489794 G C 49.64 . AC=1;AF=0.5;AN=2;BaseQRankSum=-1.136;DB;DP=23;ExcessHet=3.0103;FS=0;MLEAC=1;MLEAF=0.5;MQ=34.1;MQRankSum=0.381;QD=2.16;ReadPosRankSum=0;SOR=1.022;ANN=C|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*1848G>C|||||1848|,C|downstream_gene_variant|MODIFIER|MIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*1112C>G|||||1112|,C|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|8/10|n.1068-310C>G||||||;NS=1;MAF=0.5;AC_Het=1;AC_Hom=0;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PL 0/1:18,5:23:57:57,0,384 chr1 16288 rs113141985 C G 48.64 . AC=1;AF=0.5;AN=2;BaseQRankSum=-1.15;DB;DP=9;ExcessHet=3.0103;FS=0;MLEAC=1;MLEAF=0.5;MQ=32.24;MQRankSum=-0.135;QD=6.08;ReadPosRankSum=-1.029;SOR=0.307;ANN=G|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*1879C>G|||||1879|,G|downstream_gene_variant|MODIFIER|MIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*1081G>C|||||1081|,G|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|8/10|n.1067+319G>C||||||;NS=1;MAF=0.5;AC_Het=1;AC_Hom=0;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PL 0/1:5,3:8:56:56,0,107 chr1 16298 rs62636498 C T 112.14 . AC=2;AF=1;AN=2;DB;DP=4;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=22.25;QD=28.04;SOR=3.258;ANN=T|downstream_gene_variant|MODIFIER|DDX11L1|ENSG00000223972|transcript|ENST00000456328.2|processed_transcript||n.*1889C>T|||||1889|,T|downstream_gene_variant|MODIFIER|MIR6859-1|ENSG00000278267|transcript|ENST00000619216.1|miRNA||n.*1071G>A|||||1071|,T|intron_variant|MODIFIER|WASH7P|ENSG00000227232|transcript|ENST00000488147.1|unprocessed_pseudogene|8/10|n.1067+309G>A||||||;NS=1;MAF=0;AC_Het=0;AC_Hom=2;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PL 1/1:0,4:4:12:126,12,0 chr1 136048 rs371677125 C T 125.14 . AC=2;AF=1;AN=2;DB;DP=4;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=40;QD=31.29;SOR=0.693;ANN=T|upstream_gene_variant|MODIFIER|RP11-34P13.15|ENSG00000268903|transcript|ENST00000494149.2|processed_pseudogene||n.-153G>A|||||153|,T|downstream_gene_variant|MODIFIER|RP11-34P13.14|ENSG00000239906|transcript|ENST00000493797.1|antisense||n.*3742G>A|||||3742|,T|downstream_gene_variant|MODIFIER|CICP27|ENSG00000233750|transcript|ENST00000442987.3|processed_pseudogene||n.*1212C>T|||||1212|,T|downstream_gene_variant|MODIFIER|RP11-34P13.16|ENSG00000269981|transcript|ENST00000595919.1|processed_pseudogene||n.*1634G>A|||||1634|,T|intergenic_region|MODIFIER|RP11-34P13.15-RP11-34P13.16|ENSG00000268903-ENSG00000269981|intergenic_region|ENSG00000268903-ENSG00000269981|||n.136048C>T||||||;NS=1;MAF=0;AC_Het=0;AC_Hom=2;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PGT:PID:PL:PS 1|1:0,4:4:12:1|1:136048_C_T:139,12,0:136048 chr1 136573 . T C 113.97 . AC=2;AF=1;AN=2;DP=7;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=49.62;QD=19;SOR=1.329;ANN=C|TF_binding_site_variant|LOW|||Egr1|MA0162.2|||n.136573A>G||||||,C|upstream_gene_variant|MODIFIER|RP11-34P13.15|ENSG00000268903|transcript|ENST00000494149.2|processed_pseudogene||n.-678A>G|||||678|,C|downstream_gene_variant|MODIFIER|RP11-34P13.14|ENSG00000239906|transcript|ENST00000493797.1|antisense||n.*3217A>G|||||3217|,C|downstream_gene_variant|MODIFIER|RP11-34P13.13|ENSG00000241860|transcript|ENST00000484859.1|antisense||n.*4901A>G|||||4901|,C|downstream_gene_variant|MODIFIER|CICP27|ENSG00000233750|transcript|ENST00000442987.3|processed_pseudogene||n.*1737T>C|||||1737|,C|downstream_gene_variant|MODIFIER|RP11-34P13.16|ENSG00000269981|transcript|ENST00000595919.1|processed_pseudogene||n.*1109A>G|||||1109|,C|intergenic_region|MODIFIER|RP11-34P13.15-RP11-34P13.16|ENSG00000268903-ENSG00000269981|intergenic_region|ENSG00000268903-ENSG00000269981|||n.136573T>C||||||;NS=1;MAF=0;AC_Het=0;AC_Hom=2;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PL 1/1:0,6:6:18:128,18,0 chr1 136962 rs373582709 C T 118.68 . AC=1;AF=0.5;AN=2;BaseQRankSum=-0.887;DB;DP=8;ExcessHet=3.0103;FS=0;MLEAC=1;MLEAF=0.5;MQ=50.36;MQRankSum=0.489;QD=16.95;ReadPosRankSum=0.489;SOR=0.33;ANN=T|TF_binding_site_variant|MODIFIER|||Egr1|MA0162.2|||n.136962G>A||||||,T|upstream_gene_variant|MODIFIER|RP11-34P13.15|ENSG00000268903|transcript|ENST00000494149.2|processed_pseudogene||n.-1067G>A|||||1067|,T|downstream_gene_variant|MODIFIER|RP11-34P13.14|ENSG00000239906|transcript|ENST00000493797.1|antisense||n.*2828G>A|||||2828|,T|downstream_gene_variant|MODIFIER|RP11-34P13.13|ENSG00000241860|transcript|ENST00000484859.1|antisense||n.*4512G>A|||||4512|,T|downstream_gene_variant|MODIFIER|CICP27|ENSG00000233750|transcript|ENST00000442987.3|processed_pseudogene||n.*2126C>T|||||2126|,T|downstream_gene_variant|MODIFIER|RP11-34P13.16|ENSG00000269981|transcript|ENST00000595919.1|processed_pseudogene||n.*720G>A|||||720|,T|intergenic_region|MODIFIER|RP11-34P13.15-RP11-34P13.16|ENSG00000268903-ENSG00000269981|intergenic_region|ENSG00000268903-ENSG00000269981|||n.136962C>T||||||;NS=1;MAF=0.5;AC_Het=1;AC_Hom=0;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PL 0/1:1,6:8:10:126,0,10 chr1 137622 rs376555721 G A 54.66 . AC=1;AF=0.5;AN=2;BaseQRankSum=0.674;DB;DP=4;ExcessHet=3.0103;FS=6.021;MLEAC=1;MLEAF=0.5;MQ=20;MQRankSum=0;QD=13.67;ReadPosRankSum=-0.319;SOR=2.788;ANN=A|upstream_gene_variant|MODIFIER|RP11-34P13.15|ENSG00000268903|transcript|ENST00000494149.2|processed_pseudogene||n.-1727C>T|||||1727|,A|downstream_gene_variant|MODIFIER|RP11-34P13.14|ENSG00000239906|transcript|ENST00000493797.1|antisense||n.*2168C>T|||||2168|,A|downstream_gene_variant|MODIFIER|RP11-34P13.13|ENSG00000241860|transcript|ENST00000484859.1|antisense||n.*3852C>T|||||3852|,A|downstream_gene_variant|MODIFIER|CICP27|ENSG00000233750|transcript|ENST00000442987.3|processed_pseudogene||n.*2786G>A|||||2786|,A|downstream_gene_variant|MODIFIER|RP11-34P13.16|ENSG00000269981|transcript|ENST00000595919.1|processed_pseudogene||n.*60C>T|||||60|,A|intergenic_region|MODIFIER|RP11-34P13.15-RP11-34P13.16|ENSG00000268903-ENSG00000269981|intergenic_region|ENSG00000268903-ENSG00000269981|||n.137622G>A||||||;NS=1;MAF=0.5;AC_Het=1;AC_Hom=0;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PL 0/1:1,3:4:13:62,0,13 chr1 137825 rs147252685 G A 622.06 . AC=2;AF=1;AN=2;DB;DP=26;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=36.46;QD=24.88;SOR=0.77;ANN=A|upstream_gene_variant|MODIFIER|RP11-34P13.15|ENSG00000268903|transcript|ENST00000494149.2|processed_pseudogene||n.-1930C>T|||||1930|,A|downstream_gene_variant|MODIFIER|RP11-34P13.14|ENSG00000239906|transcript|ENST00000493797.1|antisense||n.*1965C>T|||||1965|,A|downstream_gene_variant|MODIFIER|RP11-34P13.13|ENSG00000241860|transcript|ENST00000484859.1|antisense||n.*3649C>T|||||3649|,A|downstream_gene_variant|MODIFIER|CICP27|ENSG00000233750|transcript|ENST00000442987.3|processed_pseudogene||n.*2989G>A|||||2989|,A|non_coding_transcript_exon_variant|MODIFIER|RP11-34P13.16|ENSG00000269981|transcript|ENST00000595919.1|processed_pseudogene|1/1|n.141C>T||||||;NS=1;MAF=0;AC_Het=0;AC_Hom=2;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PL 1/1:0,25:25:75:636,75,0 chr1 138156 rs370691115 G T 118197.6804 . AC=12;AF=0.51;AN=2;BaseQRankSum=-0.887;DB;DP=8;ExcessHet=3.0103;FS=0;MLEAC=12;MLEAF=0.51;MQ=50.36;MQRankSum=0.48922;QD=16.95;ReadPosRankSum=0.48924.63;SOR=0.33693;ANN=T|TF_bindingupstream_sitegene_variant|MODIFIER|RP11-34P13.15|ENSG00000268903|Egr1transcript|MA0162ENST00000494149.2|processed_pseudogene||n.136962G>A-2261C>A|||||2261|,T|upstream_gene_variant|MODIFIER|RP11-34P13.15|EN chr1 137622 rs376555721 G A 54.66 . AC=1;AF=0.5;AN=2;BaseQRankSum=0.674;DB;DP=4;ExcessHet=3.0103;FS=6.021;MLEAC=1;MLEAF=0.5;MQ=20;MQRankSum=0;QD=13.67;ReadPosRankSum=-0.319;SOR=2.788;ANN=A|upstreamRP11-34P13.16|ENSG00000269981|transcript|ENST00000595919.1|processed_pseudogene||n.-191C>A|||||191|,T|downstream_gene_variant|MODIFIER|RP11-34P13.15|ENSG00000268903|transcript|ENST00000494149.2|processed_pseudogene||n.-1727C chr1 137825 rs147252685 G A 622.06 . AC=2;AF=1;AN=2;DB;DP=26;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=36.46;QD=24.88;SOR=0.77;ANN=A|upstream.14|ENSG00000239906|transcript|ENST00000493797.1|antisense||n.*1634C>A|||||1634|,T|downstream_gene_variant|MODIFIER|RP11-34P13.13|ENSG00000241860|transcript|ENST00000484859.1|antisense||n.*3318C>A|||||3318|,T|downstream_gene_variant|MODIFIER|CICP27|ENSG00000233750|transcript|ENST00000442987.3|processed_pseudogene||n.*3320G>T|||||3320|,T|intergenic_region|MODIFIER|RP11-34P13.16-RP11-34P13.1514|ENSG00000268903|transcript|ENST00000494149.2|processed_pseudogeneENSG00000269981-ENSG00000239906|intergenic_region|ENSG00000269981-ENSG00000239906|||n.-1930C>T138156G>T|||||1930|,A|downstream_gene_variant|MODIFIER|RP11-34P13|;NS=1;MAF=0;AC_Het=0;AC_Hom=2;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PL 1/1:0,8:8:24:211,24,0 chr1 138156138817 rs370691115rs556938922 GT TC 197945.0406 . AC=2;AF=1;AN=2;DB;DP=845;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=2235.16;QD=2421.6398;SOR=0.693;ANN=T|upstream=0.739;ANN=C|upstream_gene_variant|MODIFIER|RP11-34P13.15|ENSG00000268903|transcript|ENST00000494149.2|processed_pseudogene||n.-2922A>G|||||2922|,C|upstream_gene_variant|MODIFIER|RP11-34P13.16|ENSG00000269981|transcript|ENST00000595919.1|processed_pseudogene||n.-852A>G|||||852|,C|downstream_gene_variant|MODIFIER|RP11-34P13.1514|ENSG00000268903ENSG00000239906|transcript|ENST00000494149ENST00000493797.21|processed_pseudogeneantisense||n.-2261C>A*973A>G|||||2261973|,TC|upstreamdownstream_gene_variant|MODIFIER|RP11-34P13.16|E chr1 138817 rs556938922 T C 945.06 . AC=2;AF=1;AN=2;DB;DP=45;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=35.16;QD=21.98;SOR=0.739;ANN=C|upstream_gene_variant34P13.13|ENSG00000241860|transcript|ENST00000484859.1|antisense||n.*2657A>G|||||2657|,C|downstream_gene_variant|MODIFIER|CICP27|ENSG00000233750|transcript|ENST00000442987.3|processed_pseudogene||n.*3981T>C|||||3981|,C|intergenic_region|MODIFIER|RP11-34P13.15|ENSG00000268903|transcript|ENST00000494149.2|processed_pseudogene16-RP11-34P13.14|ENSG00000269981-ENSG00000239906|intergenic_region|ENSG00000269981-ENSG00000239906|||n.-2922A>G138817T>C|||||2922|,C|upstream_gene_variant|MODIFIER|RP11-34P13.|;NS=1;MAF=0;AC_Het=0;AC_Hom=2;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PL 1/1:0,43:43:99:959,126,0 chr1 184246 . T C 52.84 . AC=2;AF=1;AN=2;DP=4;ExcessHet=3.0103;FS=0;MLEAC=1;MLEAF=0.5;MQ=41.8;QD=17.61;SOR=2;DP=4;ExcessHet=3.0103;FS=0;MLEAC=1;MLEAF=0.5;MQ=41.8;QD=17.61;SOR=1.179;ANN=C|downstream_gene_variant1.179;ANN=C|downstream_gene_variant|MODIFIER|FO538757.2|ENSG00000279928|transcript|ENST00000624431.1|protein_coding||c.*88T>C|||||88|,C|downstream_gene_variant|MODIFIER|FO538757.1|ENSG00000279457|transcript|ENST00000623083.3|protein_coding||c.*971A>G|||||679|,C|downstream_gene_variant|MODIFIER|MIR6859-2|ENSG00000273874|transcript|ENST00000612080.1|miRNA||n.*3645A>G|||||3645|,C|intergenic_region|MODIFIER|FO538757.2-FO538757.1|ENSG00000279928-ENSG00000279457|transcript|ENST00000624431.1|protein_coding||c.*88T>Cintergenic_region|ENSG00000279928-ENSG00000279457|||n.184246T>C|||||88|,C|downstream_gene_variant|MODIFIER|FO538757.1|ENSG00000279457|tran|;NS=1;MAF=0;AC_Het=0;AC_Hom=2;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PL 1/1:0,3:3:9:66,9,0 chr1 185194 . G C 91.64 . AC=1;AF=0.5;AN=2;BaseQRankSum=-1.345;DP=10;ExcessHet=3.0103;FS=4.26;MLEAC=1;MLEAF=0.5;MQ=24.41;MQRankSum=-1.15;QD=11.45;ReadPosRankSum=0;SOR=2.833;ANN=C|3_prime_UTR_variant|MODIFIER|FO538757.1|ENSG00000279457|transcript|ENST00000623083.3|protein_coding|11/11|c.*23C>G|||||23|,C|downstream_gene_variant|MODIFIER|FO538757.1|ENSG00000279457|transcript|ENST00000623083.3|protein_coding|11/11|c.*23C>G|FO538757.2|ENSG00000279928|transcript|ENST00000624431.1|protein_coding||c.*1036G>C|||||1036|,C|downstream_gene_variant|MODIFIER|MIR6859-2|ENSG00000273874|transcript|ENST00000612080.1|miRNA||n.*2697C>G|||||23|,C|dow2697|;NS=1;MAF=0.5;AC_Het=1;AC_Hom=0;AC_Hemi=0;HWE=1;ExcHet=1 GT:AD:DP:GQ:PL 0/1:3,5:8:65:99,0,65 |