Question: Aligning Mb length contigs to a reference genome
0
gravatar for sanjitsbatra
4.6 years ago by
United States
sanjitsbatra0 wrote:

I am using Bowtie2 for aligning contigs with N50 ~5Mb to a reference. I am getting the following error:

bowtie2-align died with signal 6 (ABRT) (core dumped)

 

I am using a machine with 100 GB RAM, and 24 cores, so I think memory should not be an issue. Is there any fix for this issue?

 

Alternatively, I am considering BWA-mem, BLAST, BLAT, SOAP2 and Mummer for the above task. Which one of these is most appropriate, given that I prefer an end-to-end alignment for each contig, instead of an accurate local chimeric alignment?

 

Thanks!

alignment • 3.2k views
ADD COMMENTlink modified 4.6 years ago by thiagomafra70 • written 4.6 years ago by sanjitsbatra0
1

I probably saw this too late, probably helps someone looking for an alternative now.
bwasw/bwamem were designed for this particularly, so it is pretty fast.. Speed of blat would decrease when the query size is more than 200kb. Blast would be highly accurate, but slow. Mummer should be able to do it (I have no experience with it). Lastz would be accurate too but speed would be low, and it is not for multiple cores.

ADD REPLYlink modified 4.1 years ago • written 4.1 years ago by Rohit1.4k

this is helpful for me, however I was wondering if I use blast, if gives me alignments as blocks, what if I need to obtain a general similarity % coverage and score? also the same for bwa mem, is does not give me this kind of statistics if I am not mistaken? Thanks in advance

ADD REPLYlink written 3.9 years ago by User000270

The similarity %, coverage and score must be easy to do based on alignment hits. You could look for calculate these by yourself or use extended version of blast with extended columns. Something like
blastn -task megablast -query query.fa -db database -out hits -perc_identity 90 -outfmt "6 qseqid sseq
id pident length mismatch evalue bitscore qcovs"

 

ADD REPLYlink written 3.9 years ago by Rohit1.4k
1
gravatar for thiagomafra
4.6 years ago by
thiagomafra70
Brazil
thiagomafra70 wrote:

Try using the software CONTIGuator. Take a look in More Info (in green). I used the web version as well as Command-line version with success.

http://contiguator.sourceforge.net/

[ ]'s

ADD COMMENTlink written 4.6 years ago by thiagomafra70

I tried the contiguator. However, since it uses blast+'s blastn within it, it turns out to be very slow(>24 hours) for contigs of one human chromosome. Is there a faster and/or better alternative?

ADD REPLYlink written 4.6 years ago by sanjitsbatra0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 662 users visited in the last hour