I am trying to do something very simple here, but due to memory limitations I am stuck, and can no longer use Velvet's Columbus feature.
I had the reads from a shotgun sequencing of a bacteria.
I performed genome assembly using velvet (kmer size 31) and received the contigs
Now I want to perform assembly once again on those reads, but this time, reference-based assembly, using those Velvet contigs as reference.
How can I do that in an easy way?
What I tried to do was this: using BWA to index the contigs, align the contigs with my reads and then using bwa sampe to get the aligned sam file. Sorted the sam file. used velveth with -reference option on the contigs file and -sam on the sorted sam file, and I received a Segmentation Fault (core dumped) error... Which as far as I have searched for, is related to memory.. Since this is an urgent task and I need to finish it soon. I wonder whether there are any other ways to do reference-based assembly or is there an online version of Velvet? I tried some on galaxy sites, but none offered the -reference option.
I had the reads from a shotgun sequencing of a bacteria.
I performed genome assembly using velvet (kmer size 31) and received the contigs
Now I want to perform assembly once again on those reads, but this time, reference-based assembly, using those Velvet contigs as reference.
How can I do that in an easy way?
What I tried to do was this: using BWA to index the contigs, align the contigs with my reads and then using bwa sampe to get the aligned sam file. Sorted the sam file. used velveth with -reference option on the contigs file and -sam on the sorted sam file, and I received a Segmentation Fault (core dumped) error... Which as far as I have searched for, is related to memory.. Since this is an urgent task and I need to finish it soon. I wonder whether there are any other ways to do reference-based assembly or is there an online version of Velvet? I tried some on galaxy sites, but none offered the -reference option.
Comment