How to do damage calibration after vg map

Question

[puddingmeow516]

Hi all,

I am currently using vgtools on ancient DNA. My objective is to map reads to the variation graph and call variants using vg. In the past, I used bwa to do the alignments, and remove duplicates (Markduplicates) and calibrate PMD (MapDamage) on the BAM file. However, in the vg workflow, I put in two paired FASTQ file to do the mapping, and then I got a GAM file. I haven't figure out a way to put GAM into further modification.

Any suggestions on dealing with duplicates and PMD using vg?

Many thanks!

Monica

[4galaxy77]

It seems like you want to go from gam to bam using something like this: vg surject -x x.xg -b aln.gam > aln.bam

[puddingmeow516]

Yes exactly, thank you!