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?
vg giraffe -Z hprc-v1.1-mc-grch38.gbz -p -f ./small_sim.fq -o json--max-multimaps 10 >small_sim_aln_M10.json