Collapsing my paired-end reads before or after mapping ?
0
0
Entering edit mode
11 months ago
yolek64754 ▴ 30

Hi everyone,

I am about to map some Paired-End (PE) reads (aDNA) and I don't know what is the pro and cons about those two methods:

  1. Run fastp (to remove the left-over adapters) -> map my PE dataset (using BWA) -> single BAM -> remove duplicates using Picard -> Mapping Quality filtering
  2. Run Clip&Merge from EAGER (that collapse reads and trim adapters, specialised for aDNA) -> map my new SE dataset (using BWA) -> remove duplicates using Picard -> Mapping Quality filtering

My question is; is there any advantage for any of those two methods ?

Thanks a lot.

paired-end mapping single-end collapse • 522 views
ADD COMMENT
0
Entering edit mode

aDNA = amplicon? If EAGER is a special pipeline for such data then use that.

ADD REPLY
0
Entering edit mode

My bad, should not use abbreviations. aDNA stand for ancient DNA, 30-120bp damaged reads.

ADD REPLY
0
Entering edit mode

I would still recommend using EAGER since it seems to be designed for the special type of data you have.

ADD REPLY
0
Entering edit mode

seconding Genomax on this -- why would you not want to use EAGER? Are there any specific reasons why you think alternative ways might be better suited for the task at hand?

ADD REPLY
0
Entering edit mode

I meant, correct me if I am wrong but if we collapse reads before mapping, we can loose some information from sequencing errors ? I aDNA most of the time the reads are short so the reads are fully sequenced in Left and Right.

ADD REPLY
1
Entering edit mode

EAGER must be taking that into account if you are following a standard protocol. That said every dataset has its own characteristics. You may need to start with #2 and if that does not work then try #1. I would still suggest this route.

ADD REPLY
0
Entering edit mode

This would be the point where you could dig into the documentation of EAGER to find out how the collapsing is done. Presumably they are leveraging the fact that both reads tend to cover the same piece of a fragment.

ADD REPLY
0
Entering edit mode

what is your insert size?

ADD REPLY

Login before adding your answer.

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