This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Message seen when performing alignment

It's my first time doing alignment and I keep getting this message, is this normal?

enter image description here

bwa alignment
• 2,351 views • link
updated by GenoMax GenoMax • written by K • K

Please do not paste screenshots of plain text content, it is counterproductive. You can copy paste the content directly here (using the code formatting option shown below), or use a GitHub Gist if the content volume exceeds allowed length here.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

2 answers

Yes this is normal.

• 0 views • link
written by GenoMax GenoMax

Yes, normal and actually something that in 99.99% of cases you don't need. If you inspect the bwa mem C source code you'll see that you also cannot turn this off because the relevant lines miss an argument to respect the verbosity setting of bwa mem (which you should actually be able to set to make the tool less talky), so easiest is to simply ignore it.

For example, here is one line that you could turn off to be printed to screen when verbosity is 3 or higher:

https://github.com/lh3/bwa/blob/master/bwamem_pair.c#L91

And here is the stuff you see on screen without any verbosity switch:

https://github.com/lh3/bwa/blob/master/bwamem_pair.c#L97

• 0 views • link
written by ATpoint ATpoint

Okay thank you, i have another question if you don't mind. How long will it take if my sample size R1&R2 is 53 gb each

• 0 views • link
written by K • K

It will take several hours (4-5 h with 8-12 cores). The time would be dependent on kind of hardware you have access to and performance of the storage system.

• 0 views • link
written by GenoMax GenoMax

Oh if that's the case then maybe i am doing something wrong. It's been more than 15 hours since i started the process and still no sign of completing.

• 0 views • link
written by K • K

It depends on your storage, the memory available, the CPU speed, number of cores and to some extend on the experiment. I tend to remember that a 30x WGS sample took 12h on our high-performance cluster with 16 cores using a decent node. If the job runs (check by top) it is fine, get coffee and wait.

• 0 views • link
written by ATpoint ATpoint

My workstation processor is intel xeon with 128 gb ram and 10 TB storage, it's an old computer.

• 0 views • link
written by K • K

Long as you are using multiple cores (Xeon should have them) you are maximizing the resource use. If you have an older computer and spinning disks then things can take long(er).

• 0 views • link
written by GenoMax GenoMax

Log in to answer this question.

More posts like this