Comments (4)
Another accession with the same issue: AY395000.1
from vadr.
Thanks, Tomer. I'm working on a fix and will update you when I have more to report.
from vadr.
Thanks for reporting this, @taltman. This was due to a bug that is demonstrated if an entire internal CDS feature is deleted, as is the case with AY34999 and AY35000 (the hypothetical protein sars8a CDS is deleted).
I've fixed this in the develop
branch and the fix will go into the next release: v1.1.1 which should happen within the next few weeks. As a temporary workaround, using --skip_pv
with v-annotate.pl to avoid vadr exiting in error, but note that the sars8a CDS will be incorrectly annotated in the .ftr table with a length of '0' (start coordinate is 1 more than stop coordinate, e.g for AY394999: 27702 to 27701). With the fix, deleted features are not annotated in the .ftr file, but a new DELETION_OF_FEATURE error is reported.
The iss21-delcds
test in testfiles/github-issues/
exercises the fix for this bug.
from vadr.
Thanks for looking into this! I appreciate the work-around, and look forward to the upcoming version.
from vadr.
Related Issues (20)
- VADR doesn't annotate second segment of segmented CoV genome HOT 6
- VADR produces qualifier with invalid value HOT 4
- VADR predicted nested genes, prevents submission to ENA HOT 4
- Error in v-build.pl on example command HOT 2
- VADR doesn't detect annotated segments HOT 4
- Error: summed psi of split states in node 32146 not 1.0 but : 1.001000 HOT 8
- cdsstopp can be reported for stop codon that overlaps with replaced region (due to -r) HOT 1
- conda/mamba installation improvement HOT 5
- problem with git clone corona on bitbucket HOT 2
- Bug: parse_blast.pl incorrectly detects strand if alignment line length is 1 HOT 1
- Problem about fetching and parsing protein feature table file(s) HOT 9
- ERROR in utl_RunCommand() in vadr 1.5 installed via conda HOT 18
- Unable to build a model HOT 2
- ERROR in vdr_EutilsFetchToFile, problem fetching NC_063383.1 (undefined)
- ERROR in sqf_EslTranslateCdsToFastaFile, problem translating CDS feature, unable to find expected translated sequence in EU156171/EU156171.vadr.cds.esl-translate.2.fa: HOT 3
- Is there a database built for most viruses? HOT 2
- Bug in parse_blast.pl script causes max length deletion to be incorrectly reported in rare cases HOT 1
- Incorrect reporting of short frameshifts in 5' truncated CDS in rare cases HOT 1
- Bug: v-annotate.pl does not report cdsstopn for early in-frame stop codons in 3' truncated CDS if they are the final 3 nt of the input sequence HOT 2
- Bug: early stops misreported for multi-segment CDS 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 vadr.