Comments (4)
I have added some information.
from nanocaller.
I have updated the usage docs.
from nanocaller.
Hi, thanks for making this tool available. I'd like to try it, but I'm having trouble understanding some of the options and getting it to run.
First, regarding indel mode. In your preprint you say:
"We use WhatsHap to phase variants, allowing the algorithm to re-genotype variants. In the phasing step, all variant calls made by NanoCaller, both in high-confidence regions and outside confidence regions are used. Any candidate site whose genotype changes to homozygous reference is considered as non-variantcall. Afterwards, we use WhatsHap again to phase alignment reads in order to create pileup images for indels."
But, in the usage information you say the bam file must be phased if indel mode is chosen.
The description in the preprint makes me think NanoCaller does the phasing, but the usage information suggests I have to do it. Can you clarify?
Second, it appears I can pass a bed file to NanoCaller. Can you describe what this is doing? Does it override -chrom?
Finally, can you provide an example dataset for running this program so I can make sure I'm using it correctly?
Thanks!
from nanocaller.
The default NanoCaller workflow, which is the mode 'both', goes like this:
bam file > SNP calling by NanoCaller > SNP phasing by whatshap > Reads phasing by whatshap > Indel calling by NanoCaller on reads phased in previous step > Combine indel and SNP calls
If you choose 'indel' mode, then we do not have any SNPs that we can use to phase the reads in order to call indels. This option is allowed so that if a user prefers to use a different tool than whatshap to carry out SNPs and reads phasing, or if they already have phased reads that they believe are phased using some other reliable method, then they can bypass the whatshap subroutine by NanoCaller and supply their own phased bam file.
Another reason for this option is that if a user runs NanoCaller in 'both' mode, but the program crashes during indel calling stage for some likely scenarios such as memory exceeding system limits, then they can resume NanoCaller at the indel calling stage, instead of having to do the entire SNP calling and phasing all over again.
We have removed the option for bed files for now. We found for whole genome variant calling that using a bed file would get rid of potential variants near the boundaries of the interval in a bed file. So we found it better to leave it up to a user to filter variants with a bed file using a more sophisticated tool like rtg vcffilter
. But I can see having an option like bed file can be useful if you you dont want to spend extra time in calling variants that are outside the intervals in a bed file, and submitting individual calls for each line in a bed file can be cumbersome. We will revisit this issue in our team. Do you have any preferences or thoughts about this?
We provide an example here: https://github.com/WGLab/NanoCaller/tree/master/sample
You can see the command we used in the file: https://github.com/WGLab/NanoCaller/blob/master/sample/sample_call
from nanocaller.
Related Issues (20)
- High coverage genotyping outputs only G HOT 10
- Using NanoCaller with BED files and Indel Calling Gives Empty VCF files HOT 3
- bcf tools error and error in whatshap HOT 4
- docker build fails due to missing NanoCaller_WGS file HOT 2
- Models compatible with Guppy6.x.x HOT 1
- Failed installation HOT 2
- bam file issue!
- what do you want to fill in the β--chromβ option? I read the instructions and there is no document. HOT 22
- header error in variant_calls.snps.phrased.vcf.gz HOT 6
- Unnormal results ! HOT 8
- AD (allelic depth) field HOT 5
- TypeError HOT 2
- Supporting reads HOT 3
- conda install incompatibilities HOT 2
- Models for R10.4.1 dorado
- Module not found tqtm after install with conda HOT 1
- Error phasing large chromosomes (550mb chromosome) due to inability to store as tbi index HOT 3
- Applications HOT 12
- indels HOT 2
- keras import issue by conda install
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 nanocaller.