reads of different Sequencing platform to call variants
0
0
Entering edit mode
4.9 years ago
liu9827885 • 0

Hi, I'm new for WGS. And there were some problems about the reads of different sequencing platform. One sample were sequenced by the illumina with PE150. Meanwhile, it was used to sequence by the BGI with PE100. I'm not sure could I combine the bam of the two different reads after bwa-mem and sort ? Because they have the same base quality standard. Thanks a lot !

SNP alignment sequencing • 1.0k views
ADD COMMENT
1
Entering edit mode

Illumina is a company that offers multiple sequencing platforms so by the Illumina is not informative. BGI is a service provides. You should find out 1) on which platform (HiSeq, NExtseq, Novaseq...) the first and second sample have been sequenced. If thes platforms are similar probably merging is ok. Alternatively, call variants on both datasets and only take reproducable variants for downstream analysis.

ADD REPLY
0
Entering edit mode

Thanks for your reply. I check the platform of illumina and bgi is HiSeq-2500 and BEGISEQ-500. But the different length of PE can merge directly?

ADD REPLY
0
Entering edit mode

You could merge it but I wouldn't unless I had a very good reason to do so. Both are radically different seq technologies. There have been a number of comparisons which show they are quite similar, but you'll still find some issues around different protocols eg amplification biases, different TAs, software eg base callers and so on. They will each deal with certain motifs differently.

If you can, I would call SNVs on each using the same methods and use each of them as technical variates, both of which are not able to detect the true biological variation 100% accurately.

ADD REPLY
0
Entering edit mode

Oh, and @ATpoint - BGI and MGI sequencers do exist these days (BGI is not just a service provider). Apparently costs are about 50% of Illumina, which makes for some welcome competition.

ADD REPLY
0
Entering edit mode

Ok, did not know that, my bad.

ADD REPLY

Login before adding your answer.

Traffic: 2498 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6