de novo variant SNP calling from barcoded single-cell data
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?
• 2,811 views
•
link
updated
by
niklas.lang •
•
written
by
dominicdhall •
•
written
by
dominicdhall •
1 answer
This Post was deleted.
• 0 views
•
link
written
by
niklas.lang •
Log in to answer this question.
More posts like this
-
bcftools not calling variants that appear high quality
written by bertb •
-
How does bcftools decide what sample name to assign when calling variants?
written by Paul •
-
Missing value in AD (allele depth) in VCF (eg: DP=8, AD=7,.)
written by Nicolas •
-
Undoing large genomic deletions using minimap and bcftools
written by Jamie •
-
How to force bcftools to call all variants
written by Hans
-
Window size in bcftools call accuracy
written by wenbinm •
-
Why Does Mpileup Skip some variants ?
written by sami •
-
Set bam read group from existing tag
written by jbalberge
-
why I used bcftools mpilup and call commond, no indel variance were genearte?
written by lemon •
-
Creating a cell by SNP matrix
written by dominicdhall •
Are you sure you have enough depth do to variant calling?
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.