Extract nucleotide from read at SNP position
I have read alignments in SAM format, I have read fastqs, I have genome fasta, and I have a VCF file with locations of SNPs. Is there an easy way to pull out the sequenced nucleotide for each read at the SNPs?
I'm not interested in calling SNPs (this is already done, I have my VCF) but I would like to know what nucleotide a read had at particular SNPs
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Tools to call structural variants in haploid yeast from WGS data
written by Tasneem Fetian •
-
Extract allelic depths in a BAM
written by mikothebichon •
-
Extract VCF of tumor variants from a list of SNPs
written by vctrm67 •
-
How to determine what allele a sample has from FASTQ data?
written by psyc_biostars •
-
Finding some specific SNPs in aligned sam file
written by hasani.iut6
-
Do I need to reconstruct haplotypes from SNP data to calculate nucleotide diversity?
written by suzyhocking •
-
Identify point mutations from each read in sam files
written by DVA
-
How can I filter SNP for 30 samples in vcf files?
written by sun.nation
-
How to discover new mutations in a derived strain relative to a reference genome?
written by vladimir.gritsenko
-
Is There A Simple Sam Viewer, That Can View Copy/Pasted Sam Formatted Alignment Info?
written by jnfass •
The following posts are relevant. There is also a python code in the first post that works with the pileup format.
SAM has read sequences for each alignment as well as the alignment position -- you can calculate offset into the read that corresponds to your SNP's location and look at the base
I am trying to do the same. Did you find a solution?