Comments (5)
Hi @mbhall88! Thank you very much for the feedback, I hadn't heard of nanoq. I'll see about get these swapped out and getting an update pushed in the next day or two.
Please don't hesitate to provide any other suggestions for improvements. Thanks again!
from dragonflye.
Replaced filtlong with nanoq 5c5d2ac It'll be in the next version (v1.0.1)
Thanks again for the feedback, please feel free to reopen.
[dragonflye] Filter reads based on length and quality
[dragonflye] Running: nanoq --min_length 10000 --fastx READS.sub.fq.gz --detail 2>&1 1> READS.fq | sed 's/^/[nanoq] /' | tee -a dragonflye.log
[nanoq]
[nanoq] Nanoq Read Summary
[nanoq] ====================
[nanoq]
[nanoq] Number of reads: 9,545
[nanoq] Number of bases: 186,524,370
[nanoq] N50 read length: 20,856
[nanoq] Longest read: 99,830
[nanoq] Shortest read: 10,001
[nanoq] Mean read length: 19,541
[nanoq] Median read length: 17,507
[nanoq] Mean read quality: 13.54
[nanoq] Median read quality: 13.77
[nanoq]
from dragonflye.
Hi @mbhall88
I flipped the two methods, so as of v1.0.4 (https://github.com/rpetit3/dragonflye/releases/tag/v1.0.4) read length filtering is done before read depth reduction.
Thanks again for the feedback!
from dragonflye.
Sorry, and one other question/comment. Would it not make more sense to filter the reads before subsampling? Otherwise, if you subsample to say 60x and half the reads are less than the minimum read length, you end up with a filtered fastq that is much less than 60x?
from dragonflye.
Yeah, totally agree. I think it might have been to reduce the filtlong run time, but with nanoq that won't be an issue.
I'll get the switched, and haven't in next version update
Thanks again for the great feedback!
from dragonflye.
Related Issues (20)
- Failed to run medaka consensus. - ModelStoreTF exception <class 'NotImplementedError'> HOT 14
- Memory default HOT 5
- Memory calculation HOT 3
- Medaka v1.7.3 HOT 8
- Feature request: dragonflye 1.1.N default to flye 2.9.2? HOT 3
- mismatch between model names valid for dragonflye 1.1.0 and medaka 1.8.0 HOT 14
- medaka fails to open model file for r1041_e82_400bps_sup_g615 HOT 9
- Niche QoL inquiry: polypolish acceptance of nonstandard file affixes? HOT 4
- Feature request: replicon rotation HOT 12
- flye log HOT 4
- Coverage of "0"? HOT 3
- Polypolish error HOT 6
- Citation HOT 2
- Typo HOT 1
- Citation for Nanoq in readme - suggested update HOT 1
- Polypolish error HOT 7
- Missed plasmid in reoriented.fa
- Execution error with Rasusa v1.0.0 in conda environment HOT 2
- how to check software in dragonflye?
- Can't locate FindBin.pm in @INC HOT 1
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 dragonflye.