Getting assembly from sam/bam file
1
0
Entering edit mode
3.0 years ago

My assembly (asm.fa) is fragmented, so I align it against the reference (ref.fa). How can I use the sam output to get a new fasta file for my assembly, where the contigs are put together in scaffolds?

mininmap -ax asm5 ref.fa asm.fa > aln.sam

samtools sort aln.sam > aln.sorted.sam



Thank you and Best Regards,

Priyanka

Assembly alignment • 2.2k views
ADD COMMENT
1
Entering edit mode
3.0 years ago

I think it's probably best if you use RaGOO for that - https://github.com/malonge/RaGOO

It runs minimap2 using your contigs and a reference assembly and then orders the contigs based on the alignment. You may have to rerun your alignments, I don't know whether you can give it the sam file directly.

ADD COMMENT

Login before adding your answer.

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