Manuel Holtgrewe
Manuel Holtgrewe
Also requires updates to baserow ETL scripts, thus would not be as helpful as needed right now. Postponing.
Seqvars mehari has separation of UTR exonic/intronic, needs to be ported to SVs.
VEP can handle it so should we...
Actually, `MaybeRefOrGetter` should be used and `toValue(value)` can be used to resolve this.
@mrlubos As always, I'm grateful for the work you put into this wonderful library. Is there any timeline regarding a fix to this?
As a follow-up, I now get incorrect ploidy calls and subsequent ploidy in the CNV call files for ~10 out of ~300 exome data sets enriched with Agilent V8 chemistry....
Curious, all transcripts of `HGNC:29279` are affected by this but UCSC shows RefSeq transcript alignments.
OK, this appears to be legit by the data ... ``` # grep 'NM_001377448|NM_001291324' ~/mehari-data-tx/pass-1/txs.bin.zst.report Skipping transcript NM_001291324.3 because of missing stop codon in translated CDS Skipping transcript NM_001377448.1 because...
There are 139 such transcripts, see below. Most probably, we cannot easily resolve this for GRCh37. ``` # grep 'because of missing stop codon' ~/mehari-data-tx/pass-1/txs.bin.zst.report | wc -l 139 #...
I agree, we need to find a solution with either pre-compute for all exon or retrieve online from BAM.