Biostar Test Site

This is site is used for testing only. Visit: https://www.biostars.org to ask a question.

How to do damage calibration after vg map
0
0
Entering edit mode
3 months ago

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

vgtool vg PMD aDNA • 103 views
0
Entering edit mode

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

0
Entering edit mode

Yes exactly, thank you!