Do I need to go back and filter my long-reads?
0
1
Entering edit mode
14 days ago
eesiribloom ▴ 80

A while back I asked about filtering basecalls from dorado. I was taken up with other projects so am only returning to this particular dataset now.

I took people's advice and did not filter my data and aligned to human reference genome GRCh38. My most recent samples have higher basecall quality and better alignment but the older samples have more of a tail (see violin plot below).

Is this good enough or should I filter my data e.g. for basecall quality > Q9 ?

Plot of percent identity of long-read alignments

alignment nanopore filtering QC ONT • 187 views
ADD COMMENT
0
Entering edit mode

Is this good enough

You need to make the call. We don't know what you are planning to do with this data downstream.

older samples have more of a tail

That is a characteristic of the sample (assuming same experimental procedure and flowcell/pore, dorado model was used). You could track it as a batch in case that has any effect on the analysis later.

ADD REPLY

Login before adding your answer.

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