Question: de novo variant SNP calling from barcoded single-cell data
1
gravatar for dominicdhall
2.8 years ago by
dominicdhall40
dominicdhall40 wrote:

I have a BAM file containing sorted, aligned reads from a number of single cells. Each read has an associated cell barcode (CB tag). I would like to consider all aligned reads and produce a VCF file of SNP variants from my single cell data. I have looked over the documentation for bcftools mpileup and found useful options such as --ignore-RG to ignore the readgroup when assigning a genotype probability. I therefore wanted to run something like:

bcftools mpileup -Ou --ignore-RG --no-reference -b <input.bam> | bcftools call -mv -Oz -o calls.vcf

My question is: will mpileup group reads by barcode when calculating genotype probability? Will this method successfully output SNP variants present within my aligned reads?

rna-seq single-cell vcf • 1.3k views
ADD COMMENTlink modified 21 months ago by niklas.lang20 • written 2.8 years ago by dominicdhall40
1

Are you sure you have enough depth do to variant calling?

ADD REPLYlink written 2.8 years ago by swbarnes29.6k

I didn't produce the data set but I have been assured that there is the necessary depth to work out a sample of variant SNPs present within each cell.

ADD REPLYlink written 2.8 years ago by dominicdhall40
0
gravatar for niklas.lang
21 months ago by
niklas.lang20
niklas.lang20 wrote:

This Post was deleted.

ADD COMMENTlink modified 21 months ago • written 21 months ago by niklas.lang20
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2241 users visited in the last hour
_