Getting reads contributing to a velvet assembly - contig
1
0
Entering edit mode
5.1 years ago

Dear all,

I have a fasta file of contigs generated after assembling the reads using Velvet. Can anyone help me with how to get information on the original read ids , contributing to each of the contig using the outputs from Velvet?

Thank you

Assembly Mapping • 1.5k views
ADD COMMENT
0
Entering edit mode

I am not sure how to do this using Velvet output but you can always align reads back to the contigs and extract this information from the sam or bam file.

ADD REPLY
1
Entering edit mode
5.1 years ago
lexnederbragt ★ 1.3k

if you use the -amos_file yes flag for velvetg, you'll get an afg file which contains this information. The viewer Tablet can be used to inspect the results. There may be suitable parsers/converters for the afg file that help you get what you need.

ADD COMMENT

Login before adding your answer.

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