Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: Run Quiver for hybrid assembly?
1
gravatar for marcela.uliano
4.8 years ago by
marcela.uliano • 60
European Union
marcela.uliano • 60 wrote:

Hey guys,

I have a hybrid assembly (illumina + pacbio subreads - DBG2OCL assembler) for an 1.5Gb eukaryotic genome, and now I'm in doubt if I should quiver this assembly or not.

Once, as I read, quiver calls consensus based on the pacbio quality of the aligned subreads to the draft genome, what it would do to the regions that are covered only by my Illumina contigs in the hybrid assembly?

Does it make sense to run quiver in a hybrid assembly at all?

Thanks a lot for the help! =)

alignment assembly genome • 2.3k views
ADD COMMENT • link •
modified 4.8 years ago by rhall • 160 • written 4.8 years ago by marcela.uliano • 60

Tagging Dr. rhall

ADD REPLY • link written 4.8 years ago by GenoMax ♦ 95k
2
gravatar for rhall
4.8 years ago by
rhall • 160
United States
rhall • 160 wrote:

As it's a hybrid assembly, my guess is that you don't have enough coverage of pacbio data to quiver correct (<30x), otherwise it would probably make more sense to do a pacbio only assembly? If you do have ~50x pacbio data then the pacbio coverage is likely more even than the illumina, as pacbio has no GC bias, so quiver correcting will probably work. If you have contigs with no pacbio coverage, even with only a small amount of pacbio data then I wouldn't have a lot of confidence in them, given no GC bias and that it is easier to map longer reads, all contigs should have coverage. If you don't have enough pacbio data to quiver correct then I would simply use the illumina data for correction, using pilon pilon.
ADD COMMENT • link written 4.8 years ago by rhall • 160
Please log in to add an answer.

Similar posts • Search »

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