Advice on mapping reads to polyploid de-novo assembly
1
1
Entering edit mode
4.9 years ago
O.rka ▴ 740

My organism is a novel diatom so the assembly has been de-novo using PacBio reads. I'm getting a large number of reads mapping to more than one location in my genome. I suspect this organism is polyploid so the assembly may have "duplicate" regions. I ran minimap2 using minimap2 -PD -k19 -w19 -m200 -t8 scaffolds.fasta scaffolds.fasta > scaffolds-self.paf as suggested by the author and found large stretches of homology.

How should I proceed? Do I need to somehow dereplicate my polyploid assembly into a (pseudo?)haploid assembly? I've never dealt with this before so any help would be greatly appreciated.

My of my sample replicates do not have a consistent gene overlap because a lot of the reads are being mapped to multiple locations.

assembly RNA-Seq • 1.0k views
ADD COMMENT
0
Entering edit mode
2.5 years ago

Did you end up building a pangenome for this, or how did you proceed ?

ADD COMMENT

Login before adding your answer.

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