Fq Column In Mpileup Vcf Output.

Hello the line that I have is :

22 22716000 . T A,X 76 . DP=23;AF1=1;AC1=2;DP4=0,0,8,3;MQ=57;FQ=-60 GT:PL:GQ 1/1:109,33,0,109,33,109:63

I'm wondering what FQ means as well as AF1 and DP4

thanks

mpileup vcf

• 3,778 views • link

written 15.0 years ago by Tiffani •

1 answer

Consensus quality. If positive, FQ equals the phred-scaled probability of there being two or more different alleles. If negative, FQ equals the minus phred-scaled probability of all chromosomes being identical. Notably, given one sample, FQ is positive at hets and negative at homs.

• 1 views • link

written 15.0 years ago by Tiffani •

More posts like this

VCF file REF column only contains N
written 10.2 years ago by genie66 •

I am using samtools mpileup to call variant. Command used: samtools mpileup -u -f "/gpfs/home/g/o/Homo_sapiens.GRCh37.75.dna.primary_assembly.fa" "/gopisiva/Mason/Single_cell/SC_Stressed/*.bam" | bcftools view -bvcg - > SC_Stressed.raw.bcf" bcftools view …
Hom vs Het Variant calls
written 11.7 years ago by hellbio

Hi all, I have a variant detected in four samples using samtools-0.1.19 as shown below In one of the sample it is identified as "het" …
Why Is The Variant File Produced By Samtools So Large? It Is About ~34G.
written 13.0 years ago by Jordan

Hi, I have a <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file of size 62G. When I do the variant analysis by <a href='http://samtools.sourceforge.net/'>samtools</a>, it gives a huge bcf file …
Meaning Of The Initial Column Of Filtered Out By Vcfutils.Pl Varfilter -P
written 13.7 years ago by bioinformatica2005 •

Hello all, Which are the possible values and meanings of the first column in the output (STDERR) of the variants that are obtained when the …
Extracting Specific Reads From Vcf File
written 14.1 years ago by Bnfoguy •

Hello, I am trying to extract indels with coverage ~1 and they are in a VCF file. How do I write a python script or …
Strand Information On Vcf
written 14.1 years ago by Jirapong •

My mpileup output looks like this. <pre><code>#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT xyzzy.W1g7rI9gs2.bam X 533 . C G 25 . DP=42;VDB=0.0033;AF1=0.5;AC1=1;DP4=20,0,18,0;MQ=20;FQ=26.8;PV4=1,7.1e-22,1,1 GT:PL:GQ …
I Did A Samtools Pileup And Ended In A Vcf Format Like This
written 14.3 years ago by satish.nair •

Hi Everyone I was doing the <a href='http://samtools.sourceforge.net/'>samtools</a> mpileup piping with bcftools using this command ("samtools mpileup -uf hg19.fa SingleCell5.chrM.bam | bcftools view -cg -v …
Mpileup And Snp Calling
written 14.5 years ago by User 7317 •

I am reading the mathematical on Samtools Algorithms to understand how mpileup works. I ran mpileup per each sample and multiples samples. /share/bin/samtools-0.1.16/samtools mpileup -Euf …
How To Filter The Snps/Indel Results From Mpileup?
written 15.0 years ago by Haiping •

Hi. I trid to use mpileup to identified the SNP/Indels. The commands I used were: <pre><code>samtools fillmd –bAr sample.sorted.bam ref.fa > sample.sorted.baq.bam samtools mpileup -uf …
How To Filter The Snps/Indel Results From Mpileup?
written 15.0 years ago by User 7551 •

Hi. I trid to use mpileup to identified the SNP/Indels. The commands I used were: samtools fillmd –bAr sample.sorted.bam ref.fa > sample.sorted.baq.bam samtools mpileup -ugf …