How to split alignment into smaller chunks with node range
Hi,
I have a GFA file(graph) and GAM file(alignment). I want to split them into smaller chunks with node range.
I have splited GFA file into smaller chunks with command vg chunk -r 1527845:2392738 -x test.gfa -O gfa > chunk303.gfa, but I got an error when I splited GAM file with command vg chunk -a gfa_mapped_sort.gam -r 1527845:2392738 > chunk303_mapped.gam.
I don't know whether my command is right or not.
Thank you in advance!
• 989 views
•
link
updated
by
Jordan M Eizenga •
•
written
by
Wenhai •
•
written
by
Wenhai •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Extracting individual chromosomes from a .gfa graph
written by samuel.michel •
-
Trouble building combined pangenome graph with circularized genomes + embedded ORF1 paths
written by Floren •
-
Using vg gamsort with naive sorting algorithm
written by AshleeThomson •
-
vg rna index
written by Yumin •
-
Speed up vg call ?
written by colindaven
-
how to get path IDs from vg call with snarls
written by cqm19860612 •
-
Using VG to genotype a GAM using a graph constructed from GFA
written by Michael •
-
Issues calling variants with vg when using other tools upstream
written by brian.nadon •
-
Vg call in a specific region
written by arunkumarramrsh •
-
VG view command hangs with no result (VG Toolkit)
written by Diomfeas •
Can you be more specific about the error that you encountered? Was there any output to stderr? I think the error might be that your GAM was saved to a file in your working directory rather than emitted through stdout. Did you check the working directory?