3 peaks in per base coverage plot
0
0
Entering edit mode
6.4 years ago
m.eitel • 0

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.

Michael

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

genome Assembly next-gen • 2.3k views
ADD COMMENT
1
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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 REPLY
0
Entering edit mode

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

-Michael

ADD REPLY
0
Entering edit mode

How to add images to a Biostars post

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

ADD REPLY
0
Entering edit mode

thanks. however the link should work

ADD REPLY
1
Entering edit mode

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 REPLY
0
Entering edit mode

updated the post to include images

ADD REPLY

Login before adding your answer.

Traffic: 1895 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