Should I wait for finishing the BWA MEM alignment?
0
0
Entering edit mode
5 weeks ago
oğuzhan • 0

I am trying to align two paired-reads each are 20 gb in fastq. gz format, with 2 cpu and 16 gb ram. It has taken 16 hours for now and I think ıt will not finish today. Do you have any suggestion to speed up the bwa mem alingment?

Bwa mem. • 483 views
ADD COMMENT
2
Entering edit mode

are they 20gb in zipped size? If so that are gigantic files and yes that will take some time indeed, especially given the resources you assigned to the process.

ADD REPLY
1
Entering edit mode

The author of BWA-MEM recently published chromap [REF] which is much faster than BWA-MEM with similar accuracy. You may want to consider trying it depending on your data type.

ADD REPLY
2
Entering edit mode

context: " an ultrafast method for aligning and preprocessing high throughput chromatin profiles." (copy/pasted) and "is over 10 times faster than traditional workflows on bulk ChIP-seq/Hi-C profiles and than 10x Genomics’ CellRanger v2.0.0 pipeline on single-cell ATAC-seq profiles." (copy/pasted). Chromap is good for chromatin studies, not a general purpose faster aligner (compared to BWA-mem and Bowtie2) , as I understand from the manuscript.

ADD REPLY
0
Entering edit mode

Since the output bam file keeps expanding (it is now more than 200gb), the process is possibly continuing without any crash. Thanks for the chromap suggestion. But I deduced from the manuel of choromap that it is not suitable for wgs (my sequences are wgs derived). Therefore it seems the only way is to wait for me. I hope it will not take whole week.

ADD REPLY
0
Entering edit mode

If the BAM file is growing, that's a sure sign your process is humming along and simply taking time. Next time, if you're curious, you might try timing the alignment on the first million lines of your fastq files to get an estimate on how long it might take to complete given your hardware etc.

ADD REPLY

Login before adding your answer.

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