Failed to read BGZF block data at offset 4088217654 expected 4336 bytes; hread returned 1914

[abcyulongwang]

Hello! I've successfully gotten the .gamp file of the transcriptome file with vg mpmap and he seems to have no problem. But when I run rpvg, he has the following error:

Random number generator seed: 1679924110 Fragment length distribution parameters found in alignment (mean: 225.014, standard deviation: 224.567) Loaded graph and GBWT (7.58161 seconds, 10.123 GB) [E::bgzf_read_block] Failed to read BGZF block data at offset 5048343298 expected 8709 bytes; hread returned 6944 terminate called after throwing an instance of 'std::runtime_error' what(): [vg::io::MessageIterator] obsolete, invalid, or corrupt input at message 330847274859937 group 330837660586273

In fact, I successfully ran rpvg once, but when I used the same script but different .gamp files, he got the error above.

This is my code: docker run --rm \ -v /home/yulong/Pan_genome/Pan_transcript/2023.3.21submit:/input \ quay.io/jsibbesen/rpvg:latest \ rpvg \ -t 10 -g /input/vg_rpvg.spliced.xg -p /input/vg_rpvg.haplotx.gbwt -f /input/vg_rpvg.txorigin.tsv -a /input/sample1.vg.gamp -o /input/sample --inference-model haplotype-transcripts Looking forward to your useful help！ Sincerely yulong

[jeizenga]

Could you also provide the vg mpmap command you used?

[jeizenga]

You can send material to joeizeng@gmail.com. Just to be clear: is this a separate issue from the one posted above by @abcyulongwang ?

[CarlosAmadeo7]

I have a similar issue. Did you solve it?