Comments (2)
I think this is a similar issue I outlined in more detail here artic-network/fieldbioinformatics#125 with the artic minion step.
You could examine your vcf files and see if it seems similar.
Essentially it seems there are overlapping variants in the fail.vcf and pass.vcf which causes preconsensus fasta file to be masked and thus bcftools complains when applying the "pass" variant as it expects the reference but instead finds the masked "N".
Typically I have seen this happen in indels when an amplicon drops out - when medaka variant is run on the pool containing the dropped amplicon it prodcues variants that should be filtered out - but dont seem to be.
I was wondering if it was something odd/bug happening with the longshot annotation and dropped amplicons (as it doesnt happen when medaka annotate is used via the --no-longshot flag in minion) causing the artifact variants to be kept from the dropped amplicon. Unforunately I havent had to time to look further.
One way around the issue is to run minion manually using the --no-longshot flag or use bcftools norm with the --check-ref x
flag to exclude those sites.
from viralrecon.
@dejenie21 I think that @Sam-Sims might be onto something.
BTW, the second script you run seems to have issue accessing file on github given the errors
from viralrecon.
Related Issues (20)
- [bug] Fails to find FASTQ files when stored on S3
- epi2me-labs/wf-cnv failed during analysis HOT 7
- Allow viralrecon to take gtf as input annotation file
- primer_set is not taking param
- Artic v5 mismatched primer names in artic-ncov2019 repo cause certain amplicons to be erroneously filtered/removed
- Add QIAseq DIRECT SARS-CoV-2 Kit amplicons
- Argument input-fasta is missing in NFCORE_VIRALRECON:ILLUMINA:CONSENSUS_BCFTOOLS:CONSENSUS_QC:NEXTCLADE_RUN" HOT 4
- Problem installing viralrecon
- Adding "aggregate" and "plot" methods in the freyja subworkflow HOT 4
- Error running version 2.6.0 with Nanopore data in process NFCORE_VIRALRECON:NANOPORE:ARTIC_MINION
- Unable to download the python script HOT 3
- Split authors in generated RO crate
- Adding "--grouplineages" parameter in the nf-core/viralrecon HOT 3
- Make cutadapt primer's position ext.arg
- Properly deal with multiqc in the config files before the next release HOT 3
- Temp file problem in VARIANTS_IVAR:BCFTOOLS_SORT HOT 1
- Allow skipping `freyja boot` HOT 1
- `MOSDEPTH_AMPLICON` doesn't run in `dev` branch HOT 1
- Non-SCV2 amplicon run returns consensus genomes with no low-coverage masking HOT 4
- nf-core/viralrecon run halted due to R version clash 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 viralrecon.