Open xuxingyubio opened 1 year ago
Hi, thanks for writing. Sorry for the delayed response. Could you share the command-line you used? The error happens because the chromosome "chr1" (found in one of the annotations) is not present as a path in the graph or it has another name in the graph.
I use 'vg paths -Lx
Recently, I encountered a new problem, how do you handle transcripts that span multiple chromosomes when processing gff files?
The error is like this:
[vg rna] Parsing graph file ...
[vg rna] Converting graph format ...
[vg rna] Graph and GBWT index parsed in 460.391 seconds, 8.89319 GB
[vg rna] Adding transcript splice-junctions and exon boundaries to graph ...
vg: src/transcriptome.cpp:718: int32_t vg::Transcriptome::parse_transcripts(std::vector
ERROR: Signal 6 occurred. VG has crashed. Visit https://github.com/vgteam/vg/issues/new/choose to report a bug. Please include this entire error log in your bug report!
When pangenome is used for transcriptome analysis, I want to know which of the graph you choose, full graph, clip graph, or filter graph?
Recently, I encountered a new problem, how do you handle transcripts that span multiple chromosomes when processing gff files?
vg rna
does not work for transcripts that span multiple chromosomes. This also the reason for why you see this error:
vg: src/transcriptome.cpp:718: int32_t vg::Transcriptome::parse_transcripts(std::vectorvg::Transcript, uint32_t, std::istream*, const bdsg::PositionOverlay&, const gbwt::GBWT&, bool) const: Assertion `transcript->chrom == chrom' failed.
When pangenome is used for transcriptome analysis, I want to know which of the graph you choose, full graph, clip graph, or filter graph?
For the HPRC RNA-seq analysis we used a filtered graph.
I'm very sorry to bother you
Recently, I choosed alt contigs on a chromosome to construct a pangenome.
When I used a gbz file constructed by minigraph cactus, gencode. v38. primary.gff, and alt.annotation. gff3 as input,the following error occurred: [vg rna] Adding transcript splice-junctions and exon boundaries to graph ... ERROR: Chromomsome path "chr1" not found in graph or haplotypes index (line 8).
How can I deal this problem?