Comments (14)
I think a useful analysis to inform this discussion would be:
a) map human reads to the ncov reference to see what maps
b) map coronavirus reads to human to see what gets lost
(NB: in real amplicon sequencing sets very few reads (if any) will be human so a) will be a massive overestimate of what happens in a real experiment)
from covid-19-signal.
This was resolved in commit 7b580ba right?
from covid-19-signal.
I believe we wanted confirmation to be through HISAT2 alignment against the human genome? This looks like the removal of non-SARS-CoV-2 reads, but not necessarily the confirmation step.
from covid-19-signal.
Considering we are throwing away any reads that don't map to the reference with the new "core" (see current branch) we can probably dispense with this entirely.
from covid-19-signal.
There was discussion in one of the groups that step 5 (removal of non-SARS-CoV-2) be via alignment to the SARS-CoV-2 genome but step 7 verification was alignment against human genome as this combination would absolutely pass ethics/privacy requirements and provide full confidence to our healthcare colleagues.
from covid-19-signal.
Just chatting about that on the cancogen call. Sounds like: map raw sorted reads to human reference and remove any reads that map. Then go ahead with the remaining reads for trimming.
The human removed but otherwise raw reads can then be uploaded to SRA?
I can add that and remake the PR.
from covid-19-signal.
Are we sure removing based on alignment to human won't exclude some legitimate SARS-CoV-2 data? Don't want to create a region of false low coverage.
from covid-19-signal.
Thinking about it a bit more human WGS won't be representative of the type of off-target sequences we might see so I'm going off the idea of a) a bit. b) is worth doing though
from covid-19-signal.
Yeah, not sure about benefit of a. B seems useful, as that's the issue we're worried about. Might be good to try b with a diverse sample of strains from GISAID, to get a sense of this
from covid-19-signal.
Yeah the thought with a) was to test whether mapping to human is necessary to remove any potential host reads. It would be far simpler if we can just map to coronavirus and discard everything else, I wanted to address that with a) but realized it won't really answer the question
from covid-19-signal.
Yes, agree, a) doesn't address that question.
from covid-19-signal.
Just throwing this in here so its all together.
- Need to consider mapping scores: galaxyproject/SARS-CoV-2#49
- Could do with a set of SRA archives across the nextstrain tree for doing this analysis (among other QC): e.g. create a composite reference and see if there is a good threshold for distinguishing host contamination.
from covid-19-signal.
I took the one wuhan scheme illumina sample I had to hand and ran BWA-MEM versus a composite human + viral reference.
Of those reads which mapped to viral and human contigs (~200 or 0.02%) the distribution the respective mapping qualities looking like this:
So most of the small number of problematic reads are a clear viral hit and lower quality human hit (as Torsten suggests in the linked thread).
If we just take those multihit reads with a MAPQ>=30 to the human reference, we are left with 13 (0.002%) reads:
We could save a whole 4 of them by comparing the map-qualities between human and viral. The remaining 11 reads with equally good hits to viral and human aren't likely to majorly affect the viral consensus or variant calling.
I could grab a bunch more SRAs and do this across a lot more samples but honestly, we are probably fine just using BWA-MEM MAPQ>=30 in the host removal stage and calling it a good'un.
from covid-19-signal.
This excellent, please make sure it is clear in the documentation, including the supporting data.
from covid-19-signal.
Related Issues (20)
- Lineage Assignment Failure preventing postprocess HOT 1
- Switch remaining default output to freebayes vs iVar
- Primer Scheme Correction HOT 1
- Revisit output layout
- Add Continuous Integration tests and linting
- Fix containerization
- Convert rules for wrapped tools to official wrapper
- support for pangolin v3.0.1 HOT 2
- ValidationError: 'sample' is a required property HOT 2
- Incorporate scorpio calls into report
- Galaxy workflow?
- Support IonTorrent .bam input
- Cleanup
- RuleException in line 549 of: filepath/Snakefile: 'Params' object has no attribute 'unlabelled_output_dir' HOT 2
- Update npm nextclade to the proper nextclade_cli HOT 1
- SARS-CoV-2 workflow comparison - kindly check if your work is represented correctly HOT 2
- gunzip silently fails when running get_data_dependencies.sh HOT 1
- Pipeline fails at run_lineage_assignment when trying to use Nextclade version 3.0.0 HOT 4
- Dependency issues on new install with ivar conda environment HOT 5
- Artic-ncov2019 primer scheme V5.3.2 HOT 3
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 covid-19-signal.