ruanjue / wtdbg2

Redbean: A fuzzy Bruijn graph approach to long noisy reads assembly
GNU General Public License v3.0
497 stars 91 forks source link

contig size doubled when doing local assembly #270

Closed volcano1998 closed 1 month ago

volcano1998 commented 1 month ago

I'm using wtdbg2 to perform by phase block local assembly for the human genome using Nanopore reads(48x). I used the phase block length as the estimated genome size. However, the final fasta file size is usually doubled compared to the input genome size.

Here is my example code for a 4.6M phase block.

wtdbg2.pl -t 10 -x ont -g 4613k -o PS128766436_128766069_133378765_hp1 PS128766436_128766069_133378765_hp1.fq

And the resulting fasta file is 14M.

Any thoughts why this might happen and how to fix it?

Thank you!

volcano1998 commented 1 month ago

Sorry I think found my problem. I have some bug in my input.