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
0
Entering edit mode

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

0
Entering edit mode

0
Entering edit mode

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

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?

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.

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.

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.

0
Entering edit mode