mismatch sites from GAM/GAF
I would like to know if mismatch sites can be extracted directly from GAM/GAF files.
vg
• 175 views
•
link
updated
by
colindaven
•
written
by
zhengluo •
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 •I 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 •I 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 •Hi, 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 •In 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 •Hi 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 •When 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 •Hello 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 •I 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 •Can 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 MLBHello, 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