Comments (4)
Hi Tunc,
The ICGC indel calls were made from a much earlier version (~18 mo ago) that did not fill in the genotype fields, but instead put the tumor / normal read counts in the INFO field. This has been fixed for some time, but the ICGC data is frozen in the old format. For visualization in IGV, it does seem that removing the headers that mark empty genotype fields is the best way to go.
from svaba.
Dear Jeremiah,
Thank you for the rapid response.
that did not fill in the genotype fields, but instead put the tumor / normal read counts in the INFO field
I see this info kept by putting N/T prefix in to the measure. Since I am also interested in the homogeneity of the mutation, by which measure should I calculate if that mutation is homozygous or heterozygous. or Is it any possible ?
Best regards,
Tunc.
from svaba.
Hi Tunc,
Apologies for the slow response. Genotype calculations in SvABA (was snowman) are calculated using eqn 2 of https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3198575/
However the version you have did not have genotype calculations, so indeed you would have to implement eqn 2 manually, using the REF and ALT counts.
from svaba.
Closing since this is outdated relative to much newer versions.
from svaba.
Related Issues (20)
- ERROR: Unable to open index file: ucsc.hg19.fa, yet it's in the current directory. HOT 2
- SV and indel size cutoff
- Recommendation of QUAL? HOT 2
- SvABA not reading regions file HOT 4
- Extraction of variant types
- --case-bam is unrecognized
- Is it possible to recognize complex structural variants with SvABA? HOT 2
- "Duplicate" kgID when doing annotation. HOT 2
- Is it possible to infer duplications/deletions/translocations/insertions/inversions from SvABA output? HOT 2
- How can I find contig info of "~~01D"? HOT 2
- How to input multiple case bams? HOT 1
- No SV output with metagenomics data HOT 1
- Request to use Python instead of C++ and Java
- option --max-reads-mate-region (-M) doesn't work in v1.2
- "Caught stoi/stod/stof error" in breakpoints file
- Installation issue HOT 10
- [E::bwa_idx_load_from_disk] fail to locate the index files HOT 1
- Number of sequences in BAM header mismatches reference HOT 5
- The target SV in *.svaba.sv.vcf is not found in *.alignments.txt.gz HOT 4
- recommended parameters when one chr has higher coverage than the others HOT 1
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 svaba.