Comments (1)
Hello!
If I an not wrong, this is mostly due to ShoRAH not (yet) being designed for multiplex PCR protocols, like the typical ARTIC v3/v4 that we use for wastewater.
ShoRAH assumes randomly distributed reads (e.g.: like in some fragmented protocols, that sequence 50~70bp, used for clinical swabs), so it positions its windows at regular intervals.
That's not the case with the typical ~250bp pair ended sequencing that we use with SARS-CoV-2 wastewater, there all reads start at the extremities of the multiplex PCR amplicons, and end also all roughly at similar places.
This has two consequences:
- ShoRAH can place windows at position along the genome where all the reads are only R1 from one amplicon. This will trigger the message that you report, but should be problematic.
- ShoRAH can also place windows at the border between two amplicons, and no single read will occupy >85% of the windows, causing ShoRAH to complain about coverage within that window, even if there is enough coverage according to the alignment.
We have PhD and Master students working on improving ShoRAH (making it multiplex PCR amplicon aware, and correctly matching the windows to the boundaries of the amplicons - this will completely negate the problems mentionned above) which will be released in the comming weeks.
from v-pipe.
Related Issues (20)
- Haplotype percentage or frequency estimation
- The order of jobs does not reflect the order of execution. Any help? HOT 1
- wrong url for release download in quick install HOT 2
- Running with tutorial dataset giving all empty file HOT 6
- Getting snake error HOT 23
- please document haploclique.extra_parameters HOT 3
- SARS-CoV-2 workflow comparison - kindly check if your work is represented correctly HOT 1
- Conflict in quality_assurance.smk with custom datadir-path
- frameshift_deletion_checks, simBench & Haplotype callers HOT 6
- [UX issue] unknown keywords in config file HOT 1
- [UX issue] multiple GFF HOT 1
- predicthaplo does not work on sars-cov-2 tutorial HOT 1
- hiv-tutorial works with haploclique, but fails with predicthaplo
- hiv-tutorial fails with haplotype_reconstruction: savage, and SAVAGE command output is dumped into terminal instead of log
- Error in rule hmm_align: HIV in Vpipe HOT 2
- Missing input files for rule generate_web_visualization
- vpipe broken? ModuleNotFoundError: No module named 'pandas' HOT 2
- Activating snv in sars tutorial leads to error HOT 2
- Snakemake dry run 404 HOT 9
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 v-pipe.