Comments (7)
Is this likely to change? As @nriddiford says, it seems a shame to have a 50% chance to incorporate a low quality base when the information is available to make the better call.
from tracy.
Just to add, I've also been wondering about this!
from tracy.
@tobiasrausch forgive me for pinging you, but are you intending to respond to this?
from tracy.
For tracy assemble
it's a simple majority vote. If you have, for instance, 3 traces and 2 support a gap -
and 1 a nucleotide then the gap -
is chosen. Ties are arbitrarily broken and tracy assemble
does not take into account qualities at the moment. For the pairwise case, tracy consensus
does use the qualities but gaps don't have any to begin with. Therefore, for tracy consensus
it depends on whether you use -i
or not.
from tracy.
@tobiasrausch Thanks for the clarification. How about in cases where you have 2 traces (like the image above). Is it just a 50:50 change to incorporate a low quality insertion?
from tracy.
Indeed, it's a 50:50 chance in theory but in order to make the algorithm deterministic the code currently favours nucleotides over gaps.
from tracy.
I think tracy consensus
in tracy v0.7.5 now properly handles the low-quality vs. high-quality base problem but low-quality insertions vs. gaps is still something I need to work on. Do you have some example traces that you can share with me where you think the insertion is incorrect? Thanks.
from tracy.
Related Issues (20)
- SCF version greater 2.9 required! HOT 4
- Only single-chromosome FASTA files are supported HOT 1
- Can tracy assemble trim option be disabled? HOT 1
- Decompose doesn't give bcf file give json instead HOT 3
- How can I provide a file with a list of input files to tracy? HOT 2
- Indigo HOT 4
- Cannot execute Tracy binary on mac HOT 5
- Bioconda tracy package runs very slowly where precompiled binary is fast HOT 2
- terminate called while attempting to run tracy on WSL HOT 1
- [feature request] Further options to control trimming with consensus HOT 5
- -p signal to noise ratio HOT 1
- Annotation error in the results generated sanger sequencer HOT 4
- `Couldn't anchor the Sanger trace in the selected reference genome.` error in tracy but not in Indigo HOT 3
- Issues Installing/Compiling on Mac HOT 2
- docker images differences between quay.io and dockerhub HOT 4
- Options for specifying overlap, 'fracmatch' in consensus subcommand
- Can tracy's consensus option be updated?
- Hemizygous calls HOT 2
- basecall problem HOT 2
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 tracy.