mismatch sites from GAM/GAF
I would like to know if mismatch sites can be extracted directly from GAM/GAF files.
vg
• 188 views
•
link
updated
by
colindaven
844
• |
written
by
zhengluo
0
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Is gaf.gz supported by vg pack and gamsort?
written by gulin 1I am inspired by the article: https://jmonlong.github.io/manu-vggafannot/. I am not working well with gam format. So I am thinking about build a pipline based on …
-
GATK SNP calling from BAM
written by zhengluo 0I used vg giraffe to align long reads to a graph, generating a GAM file. Then I converted the GAM to a BAM file. I …
-
Does vg stats support gaf file?
written by Wang Cong 2Hi, I am trying to get alignment stats from a gaf file using vg stats -a <gaf> and got a long error. May I ask …
-
Can a read be aligned more than once?
written by guntul 4In linear reference(sam,bam) or in pangenomic alignment(gam,gaf), is it possible to a read that is aligned twice or more? Can we observe that a read …
-
Difference between gam and sam/bam files?
written by Uveyik 8Hi everyone, I am trying to understand pangenomic data/file formats. In the vg's descripton of file format, it says gam is similar to sam/bam and …
-
vg convert does not preserve the order of the reads
written by Yihang 0When I use `vg convert` (v.1.39) to convert a `gam` file to a `gaf` file (or vice versa), I find that the order of the …
-
How to calculate amino acid coding changes and predict coding outcomes
written by yoser4 1Hello everyone. I have a snp vcf file of 99 samples from which I screened for my gene of interest and extracted 14 loci in …
-
Converting gam to gaf
written by hmg 2I noticed it is now possible to output alignments in gaf files instead of gam, which is great! How difficult would it be to convert …
-
Merging .gam files
written by elizabethmareeross 3Can anyone point me towards a method to merge gam files that are created by using ' vg map ' to map illumina reads to …
-
How to annotate splicing sites, ESEs and ESSs in a gene
written by Aurelie MLB 36Hello, Given a big list of genes, I would like to annotate the positions of splicing sites, ESEs (Exonic splicing enhancers) and ESSs (Exoninc Splicing …
I don't think this is possible yet. Not many SNP callers are available for pangenomes. I think most will work in the VCF reference space after vg surject.
See https://github.com/colindaven/awesome-pangenomes