Comments (10)
I have also just encountered the issue with low read quality using V14 chemistry (sequencing_MIN114_DNA_e8_2_400T.toml) using readfish 0.0.10dev2, however the unblock read lengths looked fine. Did the enrichment results look as expected using this fix?
from readfish.
I'm about to test the enrichment over the next few days. Will post an update after that.
from readfish.
OK, so with this new run even setting break_reads_after = 1.0 didn't fix the issue. We had poor quality and unblocked read lengths were well over 1kb... Maybe we need a new GPU...
from readfish.
OK - whilst I would still advise one second for running at present, I am not convinced that the quality issue here is related to the break reads setting. The quality that @rainwala is showing above is lower than I would expect. As an example, here is a run we currently have on promethION which is running adaptive with break reads at 1 second using kit14:
(This is the modal quality score).
You can see the unblock peak is fine -
And the unblock read length N50 is around 900 bases.
We are unblocking a lot of the data here which is why you can't see any on target reads at all.
If I zoom in on the approrpriate lengths we see this:
Having said all the above, I do note that I see some strange differences in the sequenced reads when measuring read quality. I've not noticed specifically that this impacts enrichment and I would actually look at the aligned data and not the quality scores as reported by minknow.
My conclusion with the observation of @rainwala is that this might be a more fundamental library quality issue. @samhorsfield96 can you provide any more information on what you have seen?
Also can you report minKNOW versions here too?
from readfish.
@mattloose I think it isn't necessarily our library because we get much higher Q scores with readfish / adaptive sampling (we are trying both to try to get anything to work!) turned off. It may be something related to this:
https://community.nanoporetech.com/posts/low-q-score-with-new-kit-1
from readfish.
That is an interesting post. Reading my comment above I realise it isn't quite clear - in my comment about a fundemental library quality issue I really meant it's something other than the change in break_reads length.
We also see some strange changes in reported read quality - I can't comment if recalling the data using the super model will resolve the problem but I can imagine that to be plausible.
What I can say is that I didn't think I had seen the quality drop off in real world data espeically. But I will look more.
from readfish.
Could you both check the accuracy of on target and off target reads after mapping to the reference and look at the quality that way?
from readfish.
Hi @mattloose sorry for not getting back sooner. Alignment looks fine using pass and fail reads - I haven't got formal measures of quality, however when looking at enrichment by abundance, values are as expected even using failed reads. I've noticed that the high volume of failed reads doesn't happen with every run either, and when running with the exact same setup and library I get different outcomes (sometimes many reads fail, sometimes the run looks normal). Not sure if this info helps your end, but it looks to me that using 'failed' reads doesn't actually impact alignment or calculation of enrichment.
from readfish.
This issue is stale because it has been open 30 days with no activity. Remove stale label or comment or this will be closed in 5 days.
from readfish.
Closing as stale.
from readfish.
Related Issues (20)
- Readfish and boss-runs in 2024 HOT 3
- readfish with just depletion? HOT 4
- Not being able to finish the tutorial/test HOT 6
- Provide additional TOML files which provide examples of host depletion and other use cases for readfish.
- location of readfish output HOT 4
- toml_live handling issue
- Depletion setting and unexpected behavior HOT 2
- Question about Readfish basecalling performance for messy samples HOT 6
- Playback issue HOT 14
- Problem communicating with basecaller after MinKnow update HOT 10
- Detect MinKNOW version and warn. HOT 1
- Circular Binary Segmentation HOT 3
- Having only one region in toml caused issue HOT 12
- chore: update github workflow actions HOT 1
- Handshake failed with fatal error SSL_ERROR_SSL: error:1000007d:SSL routines:OPENSSL_internal:CERTIFICATE_VERIFY_FAILED HOT 4
- comparison between MinKNOW and Readfish HOT 4
- Error with flongle flow cell splitting HOT 4
- chore: update FAQ
- next release candidate checklist
- High count of no-barcode that map to my targets HOT 15
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from readfish.