vg icon indicating copy to clipboard operation
vg copied to clipboard

the problem about 'vg surject+samtools'

Open xuxingyubio opened this issue 1 year ago • 7 comments

I'm very sorry to bother you.

I used 'vg mpmap' to get a '.gamp' file, and then I tried to project down to linear reference. The code like this:

vg surject -x "$GRAPH" -t 20 -m -b -S -A -F "$REF_PATHS" "$GAMP" > "$BAM" (vg paths -x $GRAPH" -L > REF_PATHS)

When I used 'samtools sort -@40 --output-fmt BAM "$BAM" -o ${BAM}_sort.bam‘, I met an error like this: [E::bam_read1] CIGAR and query sequence lengths differ for A00511:418:HWLJWDSXY:3:1170:3360:12320 samtools sort: truncated file. Aborting

What should I do to solve this problem?

xuxingyubio avatar Jul 22 '23 10:07 xuxingyubio