Comments (3)
Hi @mmlian05,
Thank you for your kind words :o)
Very large SV can sometime overlap too many features locations. It appears that depending on the visualisation program used (spreadsheet programs mostly) this annotation can be truncated. In order to avoid such embarrassing glitch and maybe also because overlapping so many features is already a problem, AnnotSV restricts the number of overlapping reported features.
Here are the features in question:
- RE_gene
- po_B_gain_someG_source
- po_B_gain_someG_coord
- po_B_loss_someG_source
- po_B_loss_someG_coord
If you want to change this behavior, you have to hack the code in 2 files:
(replace 19 or 49 by the desired maximum)
$ANNOTSV/share/tcl/AnnotSV/AnnotSV-benignsv.tcl:
lappend L_poBenignText($SVtoAnn) "[join [lrange $L_benign_source($SVtoAnn,$svtype,allG) 0 19] ";"]..."
lappend L_poBenignText($SVtoAnn) "[join [lrange $L_benign_coord($SVtoAnn,$svtype,allG) 0 19] ";"]..."
lappend L_poBenignText($SVtoAnn) "[join [lrange $L_benign_source($SVtoAnn,$svtype,someG) 0 19] ";"]..."
lappend L_poBenignText($SVtoAnn) "[join [lrange $L_benign_coord($SVtoAnn,$svtype,someG) 0 19] ";"]..."
$ANNOTSV/share/tcl/AnnotSV/AnnotSV-regulatoryelements.tcl:
set g_re($SV) "[join [lrange $g_re($SV) 0 49] ";"]..."
Best,
Véronique
from annotsv.
Hi Véronique,
Thank you for the advice. Will try them out.
Wondering for user's own annotation should i look into the following *.tcl and line to modify the code?
$ANNOTSV/share/tcl/AnnotSV/AnnotSV-userBED.tcl:
set L_userBEDAnn_${i}($SVtoAnn) [lrange [set L_userBEDAnn_${i}($SVtoAnn)] 0 50]
Thank you! :)
Best regards,
Michelle
from annotsv.
Wondering for user's own annotation should i look into the following *.tcl and line to modify the code?
$ANNOTSV/share/tcl/AnnotSV/AnnotSV-userBED.tcl
Absolutely, yes
from annotsv.
Related Issues (20)
- Unusual large memory needed for very large sized coupled BNDs HOT 4
- Why is there never a class 2 in ACMG_class? HOT 1
- CytoBand appear... HOT 3
- MAF in gnomAD HOT 1
- Is there a way to include ensembl_ids and/or hgnc_ids along with gene_name? HOT 3
- New long-read population reference HOT 1
- Question about reverse strand genes HOT 7
- Rogue double quote in PRDM10 HOT 6
- custom bed annotation. overlap within n bases HOT 2
- Multiple column and tab errors with CNVkit vcf output HOT 4
- left gene HOT 15
- Mapping of structural variant to intron or coding regions HOT 2
- exomiser error: connect failed connection refused HOT 2
- vcf errors HOT 1
- AnnotSV passing wrong arguments to variant convert HOT 5
- Keep AnnotSV_ID consistent with SV_chrom HOT 10
- frameshift between full and split mode HOT 5
- UCSC link failed to open the correct position HOT 5
- AnnotSV" line 283 HOT 3
- VCF for false positive deletion and for compund etherozygosity
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 annotsv.