Identify parent of each read in a GAF
I have a .gfa created by running the fasta files of two genomes through the Cactus-Minigraph pipeline. I am aligning PacBio Hifi reads to that reference using GraphAligner. Column 6 of the resulting GAF file lists the segments for the alignment path of a given read.
I expect to see reads that come exclusively from either genome, and also 'recombined' reads. What is an efficient way to classify each read as coming Genome A or Genome B, or recombined/mixed? Relatedly, how do I get the list of segment IDs belonging to each genome so that I can use it to 'decode' the parent from the segment IDs in column 6 of the GAF file?
• 1,116 views
•
link
updated
by
Jordan M Eizenga •
•
written
by
cfourps •
•
written
by
cfourps •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Structural Variant identification on pangenome graphs
written by Panos
-
Type of Pangenome Type/program to use for fungal isolates of the same species, but have major trans…
written by Drew •
-
Rank in rGFA files
written by tsnorri •
-
How to generate a GAF that shows alignments to all component genomes in graph
written by peterdfields •
-
Right method to annotate genes on pangenome graph
written by Quanyu •
-
Extract alignment between two samples from a GFA?
written by cmdcolin
-
Hifiasm coverage rd:i:0
written by ccruz •
-
What is the node_id in vg giraffe output (gaf or gam)?
written by Hendricks27 •
-
vg rna index
written by Yumin •
-
Issues calling variants with vg when using other tools upstream
written by brian.nadon •
It's not exactly clean, but one option is
vg paths -Awith-pset to get the two input genomes individually as GAF, thenvg pack -don each genome path GAF to get a table of nodes that are included.