Question: Wgs Using Solid And Illumina Platform
gravatar for Ashutosh Pandey
8.6 years ago by
Ashutosh Pandey12k wrote:

Hi Everyone,

I have a genomic data for a specific strain of mouse. The data consists of short reads generated from both the SOLiD (~20X) and Illumina (~20X) platforms. I have already aligned the reads on the reference genome using BWA (Illumina) and Lifescope (SOLiD) respectively. Following are my questions:

1) What should be my approach for calling variants. Should I call the variants separately for BAM files generated from two different platforms or I can merge these two BAM files and then call for the variants ?

There are many tools which are platform dependent such as Dindel. But tools such as samtools SNP, Indel caller OR GATK are platform independent and in this case the two bam files (SOLiD and Illumina) could be merged and used to call the variants. The colour space reads have already been transformed into base space in the BAM file so I don't think there should be a problem.

2) If Dindel uses BAM files to call for short indels then why it can't be used on SOLiD bam files with transformed nucleotide sequence ? OR why it is not recommended?

3) What should be my approach ? Doing platform specific analysis including calling variants and then choose the common variants ? OR to create a high coverage BAM files by merging two independent BAM files from different platforms and use GATK/samtools ?



calling variant • 2.1k views
ADD COMMENTlink written 8.6 years ago by Ashutosh Pandey12k
gravatar for JC
8.6 years ago by
JC12k wrote:

Hi Ashutosh,

1) and 3) I might consider analyzing both samples separately and then recover the union (not the intersection) of high quality SNPs, indels or structural variants. This is because each technology have their own bias and specific advantages.

2) Maybe you need to recompute or adjust the mapping quality score and the quality scores per base to avoid errors.

Also check this paper: it's a 2009 paper and tech have improve but it can give you an idea.

ADD COMMENTlink modified 8.6 years ago • written 8.6 years ago by JC12k

Thanks JC. Your suggestion makes sense. Thanks for the paper too.

ADD REPLYlink written 8.6 years ago by Ashutosh Pandey12k

your welcome, good luck

ADD REPLYlink written 8.6 years ago by JC12k
Please log in to add an answer.


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