Question: How to Convert Illumina reports to Vcf files?
0
gravatar for vikrantvicky.1987
2.6 years ago by
vikrantvicky.19870 wrote:

Hi All,

I am new to CNV analysis and Vcf file conversion.

I have an Illumina text report and I am looking for a solution how can I convert them into a vcf file(with BAF and LRR) for copy number variation analysis by using bcftools/cnv

my input file looks like:

[Header]
GSGT Version 1.9.4
Processing Date 3/6/2014 10:39 AM
Content humancore-12v1-0_a.bpm
Num SNPs 298930
Total SNPs 298930
Num Samples 24
Total Samples 24
File 7 of 24
[Data]

Sample Index Sample ID Sample Name SNP Index SNP Name Chr Position GT Score GC Score Allele1 - AB Allele2 - AB Allele1 - Top Allele2 - Top Allele1 - Forward Allele2 - Forward Allele1 - Design Allele2 - Design Theta R X Raw Y Raw X Y B Allele Freq Log R Ratio SNP Aux SNP ILMN Strand Top Genomic Sequence Customer Strand

Thanks in Advance

Vikrant

snp • 3.2k views
ADD COMMENTlink modified 9 months ago by freeseek40 • written 2.6 years ago by vikrantvicky.19870

Hi Igor,

Thank you for the answer. I have already tried Plink to get vcf files by using above threads and generated vcf files.

Output given below:

##fileformat=VCFv4.2
##fileDate=20160806
##source=PLINKv1.90
##INFO=<ID=PR,Number=0,Type=Flag,Description="Provisional reference allele, may not be based on real reference genome">
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
#CHROM  POS ID  REF ALT QUAL    FILTER  INFO    FORMAT  7_BG01V
1   752566  rs3094315   T   .   .   .   PR  GT  0/0
1   752721  rs3131972   C   .   .   .   PR  GT  0/0
1   798959  rs11240777  G   .   .   .   PR  GT  0/0
1   838555  rs4970383   T   G   .   .   PR  GT  0/1
1   846808  rs4475691   T   C   .   .   PR  GT  0/1

But I need vcf file with BAF and LRR values in it for further analysis.
Do you have any idea how can I do that?

Thanks in advance.

Regards, Vikrant

ADD REPLYlink modified 9 months ago by finswimmer11k • written 2.6 years ago by vikrantvicky.19870

Hello Vikrant:

I have the same situation. I received the Illumina text report. I was asked to convert into vcf. Did you use plink to convert the Illumina text report into PED then VCF? If so, how did it go?

Thanks Jenny

ADD REPLYlink written 2.2 years ago by jtftyl20
1
gravatar for igor
2.6 years ago by
igor7.6k
United States
igor7.6k wrote:

You may have to convert to plink ped and then vcf.

See previous discussion on converting to ped here:

And then ped to vcf:

ADD COMMENTlink modified 2.6 years ago • written 2.6 years ago by igor7.6k
0
gravatar for freeseek
9 months ago by
freeseek40
freeseek40 wrote:

I wrote gtc2vcf, a high performing bcftools plugin to convert Illumina GTC files to VCF on Linux. It can also be made to work with IDAT files. Some features are currently missing (as of 2018-07-11) but I plan to include them in the future if there is demand.

ADD COMMENTlink written 9 months ago by freeseek40
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 879 users visited in the last hour