Miniasm Assembler producing smaller size of reads
0
0
Entering edit mode
7.4 years ago

Hi, all.

I am trying to assemble CCS reads from PacBio with Miniasm. We used a bait library to enrich sequences of interest and expect that assembled reads would be about a few Mb. Canu was able to give what I wanted. But the outputs from Miniasm are all smaller than 0.5Mb (when converted to fasta). Has anyone tried Miniasm to assemble reads of a small size? How would I need to improve this?

Here is what I did:

minimap -Sw5 -L100 -m0  CCS_reads.fasta CCS_reads.fasta > CCS_reads.paf

miniasm -f CCS_reads.fasta CCS_reads.paf > reads.gfa

awk '/^S/{print ">"$2"\n"$3}' reads.gfa > reads.utg.fa

Thank you in advance.

Assembly • 1.7k views
ADD COMMENT

Login before adding your answer.

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