vg
vg copied to clipboard
the problem about 'vg surject+samtools'
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?