Comments (2)
It is present in the snp_stats file because it was picked up as a candidate site, but not included in the vcf file because it was determined to be false positive. NanoCaller calculated probability of presence of A base =0.1847 which is too low for a variant call. Are you using Nanopore reads or PacBio? It might help to zoom out on IGV to see the surrounding 1-2000 bp for a better understanding of why this was regarded as false positive.
We are planning to release an update which allows you to get a different snapshot of the bam file than IGV, similar to the one in Fig1 of our biorxiv paper. It will show you only the high alternative allele frequency sites and skips other bases, and this would allow you to see if an allele might be false or not.
from nanocaller.
If I filter for only primary read mappings, then it calls the variant fine:
snps.vcf:
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT SAMPLE
Chr_22 42131531 . G A 38.810 PASS . GT:DP:FQ 0/1:45:0.4444
snp_stats:
pos,ref,prob_GT,prob_A,prob_G,prob_T,prob_C,DP,freq
42131531,G,0.9283,0.5908,0.9420,0.0000,0.0000,45,0.4444
I'm trying to understand why it calculates the probability of presence of A base to be so low when the supplementary reads are included.
I'm guessing including the supplementary reads introduces too much 'noise' in the surrounding area?
I'm using ONT data that has been corrected with Netcat. I just wanted to compare the variant callings with and without error correction.
from nanocaller.
Related Issues (20)
- Runtime HOT 7
- failed log report HOT 2
- Error: Resource temporarily unavailable HOT 4
- "-phase_bam" doesn't work HOT 1
- provide the pre-trained model in other websites HOT 2
- ONT Indel calling on polyploid plant
- multisample SNP calling? HOT 7
- Link to HiFi trained models HOT 2
- AttributeError: tuple object has no attribute ndims HOT 1
- The output data is now only an args file (previously several vcf files were generated).
- index file .fai error HOT 2
- NanoCaller for microbial genomes
- 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
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.