priors for calculting genotype likelihood using RNAseq data · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

priors for calculting genotype likelihood using RNAseq data

Can anyone point me to a literature or explain the need of population priors during genotype likelihood estimation using RNAseq data? I did get results connected to DNAseq. I am using bcftools mpileup (v1.4) for likelihood estimation and then genotype calling and variant calling using bcftools call (v1.4).

rna-seq

• 1,269 views • link

written 9.2 years ago by prasundutta87

Why would variant calling on RNAseq be any different than that on DNAseq data in this respect?

• 0 views • link

written 9.2 years ago by Devon Ryan

0 answers

No answers yet.

Log in to answer this question.

More posts like this

Why does bcftools call from different max-depth (-d) value in mpileup give inconstant GT values?
written 3.5 years ago by magnolia •

Hi, I want to call variants from a BAM file with bcftools. To use `call` command, I have to use `mpileup` command to create a …
Does bcftools call ignore duplicate markings?
written 4.8 years ago by Egelbets •

I am using Picard [MarkDuplicates][1] to mark duplicates in my position sorted BAM. Then, I variant call using bcftools [mpileup][2] and bcftools [call][3] like so: …
bcftools mpileup before bcftools call
written 4.8 years ago by Egelbets •

I want to variant call using [bcftools call][1]. However, when researching, I see a lot of people running [bcftools mpileup][2] before doing the actual variant …
some good practices doubts
written 8.7 years ago by prasundutta87

Hi, I have a paired end DNA sequence alignment BAM file (got from using BWA -M) and would like to know if the points mentioned …
RNAseq variant filtering criteria
written 9.0 years ago by prasundutta87

Hello, Can anyone guide me or share what parameters should be considered for filtering RNAseq variants and why? I haven't got any published material on …
mpileup BAQ score
written 9.4 years ago by prasundutta87

Hi, I was calling variants from RNAseq data using samtools mpileup|bcftools call. I am aware that the BAQ (base alignment quality) strategy is invoked by …
variant annotation issue
written 9.4 years ago by prasundutta87

Hi, I have recently performed a variant annotation using Snpeff in water buffalo. The variants were called through bcftools call program (multiallelic caller) from Samtools …
difference in genotypes in a multisample vcf file
written 9.5 years ago by prasundutta87

Hi, I have two individuals of the same animal (one male and one female water buffalo). I have RNAseq data from both the animals and …
Bcftools magic string
written 9.7 years ago by blur

Hi, I ran mpileup and then wanted to transform the output by bcftools using this line: bcftools call -mv - > Sequence.vcf I got these …
Output AO using bcftools
written 10.3 years ago by L. A. Liggett

I calling variants from sorted bam files using `samtools mpileup` and then `bcftools call -cv input > output.vcf`. But in the resulting variant file there …