How to obtain full alignment results?
Hi, When I prepare to align my short reads to the human pan-genome graph genome, the result is only a path, 1. what I want is to provide me all the sequences that have a mismatch and full alignment with this reads?
- And how can I know the annotation information of the sequences that are aligned to the reads, such as to
exon regions, cds regions or intronregions, from the results?
Shell:
vg giraffe -Z hprc-v1.1-mc-grch38.gbz -p -f ./small_sim.fq -o json--max-multimaps 10 >small_sim_aln_M10.json
Results:
• 628 views
•
link
written
by
nrqstudent •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
VG Crash with HPRC V2.0 Graph and Indexes
written by Alexandria A. •
-
Removing PCR duplicates in gam file
written by gulin201909 •
-
GBWT Construction Error: "Sources X and Y both have node Z" When Building GBZ from GFA
written by Sauers •
-
GRCh38 reference paths in human pangenome
written by Rugare •
-
Unexpected behavior of `vg surject`
written by Ilya •
-
What path list should be used to surject GAM to GRCh38/CHM13 reference genome
written by Hanc •
-
vg autoindex seems stuck after a week of running
written by StevieP •
-
Merging alignment graphs vs merging BAMs
written by Daniel •
-
display vcf after giraffe alignment on IGV
written by Michal •
-
Indexing the human pangenome draft
written by lisle.mose •
