Question: bwa -mem log messages
4
gravatar for aleka
3.7 years ago by
aleka100
United Kingdom
aleka100 wrote:

Hello,

I am using bwa mem to aligh paired ends and I get the following errors:

[M::process] read 1057202 sequences (80000122 bp)...
[M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (2, 440646, 0, 1)
[M::mem_pestat] skip orientation FF as there are not enough pairs
[M::mem_pestat] analyzing insert size distribution for orientation FR...
[M::mem_pestat] (25, 50, 75) percentile: (100, 120, 148)
[M::mem_pestat] low and high boundaries for computing mean and std.dev: (4, 244)
[M::mem_pestat] mean and std.dev: (125.80, 35.21)
[M::mem_pestat] low and high boundaries for proper pairs: (1, 292)
[M::mem_pestat] skip orientation RF as there are not enough pairs
[M::mem_pestat] skip orientation RR as there are not enough pairs
[M::mem_process_seqs] Processed 1056888 reads in 132.524 CPU sec, 132.576 real sec

Please let me know if you have any idea of what this might be.

Best.

sequencing alignment next-gen • 2.5k views
ADD COMMENTlink modified 3.7 years ago by Istvan Albert ♦♦ 81k • written 3.7 years ago by aleka100
3
gravatar for Istvan Albert
3.7 years ago by
Istvan Albert ♦♦ 81k
University Park, USA
Istvan Albert ♦♦ 81k wrote:

None of these are error messages.

Different library preparations and sequencing instruments may read DNA fragments in different orientations.

bwa tries all these orientations to detect your fragment sizes.

So these messages just say that all your reads seem to be in FR orientation - which is the most common library prep.

ADD COMMENTlink modified 3.7 years ago • written 3.7 years ago by Istvan Albert ♦♦ 81k
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: 1139 users visited in the last hour