BWA MEM anfter trimming
0
0
Entering edit mode
12 months ago
Edmond • 0

Hi there,

I have WES data. The read length is 35-101 bp. If I trim the sequence, I get 15-70 length sequence. I trimmed based and performed illumicaclip order on trimmomatic.

The question is as follows:

I should get a mean coverage of X100, (based on the analysis of the same fastq files). However I get 68X of coverage.

I think it could be related with the aligner: Mi idea is that when I perform the alignment on the trimmed data, I use BWA-MEM however is on the limit of the minimum read length. Should I modify the default parameters?

Because the Quality of the alignment according to qualimap is 40. I think this is a low BQ.

The per base sequence plot before and after are here:

enter image description here enter image description here

BWA Alignment MEM • 519 views
ADD COMMENT
1
Entering edit mode

Don't trim the initial base pairs as you are throwing away good data. See: https://sequencing.qcfail.com/articles/positional-sequence-bias-in-random-primed-libraries/ This is also true for tagmented libraries.

ADD REPLY
0
Entering edit mode

Thank you!

ADD REPLY

Login before adding your answer.

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