How to do damage calibration after vg map
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
• 1,291 views
•
link
written
by
puddingmeow516 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
GATK SNP calling from BAM
written by zhengluo •
-
Removing PCR duplicates in gam file
written by gulin201909 •
-
Does vg giraffe output secondary and supplimentary alignments?
written by lushjia •
-
VG : No reference-sense paths available in the graph; falling back to generic paths.
written by Hang •
-
vg giraffe
written by puddingmeow516 •
-
Tools to Construct variation graph from contigs
written by Microbiomolecule •
-
vg call on augmented graph blank
written by kingcohn •
-
Segmentation fault (core dumped) during vg surject
written by LilithElina •
-
Extracting specific regions(ex, chr1:100-200 of hg38) from whole genome .vg and .gam file
written by hd00ljy3 •
-
VG view command hangs with no result (VG Toolkit)
written by Diomfeas •
It seems like you want to go from gam to bam using something like this:
vg surject -x x.xg -b aln.gam > aln.bamYes exactly, thank you!