Question: 3 peaks in per base coverage plot
gravatar for m.eitel
4 days ago by
m.eitel0 wrote:

I have problems interpreting the results from our assemblies. We have sequenced an animal with nanopore and Illumina reads. I first performed Canu v1.7.1 assemblies of subsets of the nanopore reads to see effects on the outcome. Here I generated individual assemblies with all reads >1kb, >2kb, etc. up to 20kb (with a jump from 15kb to 20kb straight). Afterwards, I corrected each Canu assembly with the Illumina reads each with 10 iterations of Pilon error correction (using Bowtie2 mapping). Then I used Bedtools to get histograms for each mapping from the bam files. The results of the per-base coverage plot can be seen on the attached graph (cov on the x- and millions of bases on the y-axis; the number in brackets behind the sets is the rough nanopore coverage based on the 1kb assembly length).

In the plot one can see coverage ratios of 0.25, 0.5 and 1. First I was interpreting the graph as a triploid organism, but this does not fit to the ratio. One would expect 0.33-0.66-1 for a triploid.

I also performed a mapping of the corrected & trimmed >1kb nanopore reads (generated by Canu) to the raw >1kb read assembly with bwa and also plotted a per base coverage.

Any Idea what is going on here? Any help in interpreting the results would be highly appreciated. This drives me nuts already...

Thanks for your help.


Illumina reads mapped to the corrected nanopore assemblies corrected & trimmed nanopore reads >1kb mapped to the raw nanopore assembly

next-gen assembly genome • 159 views
ADD COMMENTlink modified 3 days ago • written 4 days ago by m.eitel0

You may want to use k-mer analysis with your raw data to estimate ploidy.

ADD REPLYlink written 3 days ago by genomax51k

Hi there - unfortuntely your graph doesn't seem to be attached. Can you share more insights into the type of organism your assembling. I know some unicellular production organisms that became tetraploid but lost parts of their genome again - however in that case you might expect a peak or shoulder at .75, too.

ADD REPLYlink written 4 days ago by Carambakaracho360

Thanks for your comment. I added links to the graphs. Please le me know what you think.


ADD REPLYlink written 3 days ago by m.eitel0

How to add images to a Biostars post

This way they will be parsed in-line with your text.

ADD REPLYlink modified 3 days ago • written 3 days ago by genomax51k

thanks. however the link should work

ADD REPLYlink written 3 days ago by m.eitel0

While the links work, people are sometimes leery of clicking random links. Having the graphic inserted inline with the post makes things simple to follow.

ADD REPLYlink modified 3 days ago • written 3 days ago by genomax51k

updated the post to include images

ADD REPLYlink written 3 days ago by m.eitel0
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: 962 users visited in the last hour