Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Samtools-mpileup aslihan Bioinformatics 3 01-18-2017 10:53 PM
samtools mpileup: -B or not to -B? yenhuahuang1 Bioinformatics 4 03-19-2013 03:21 AM
samtools mpileup Ling1 Bioinformatics 1 05-24-2012 10:09 AM
Many Ns in samtools mpileup mikheyev Bioinformatics 5 03-22-2012 09:56 AM
samtools mpileup hv4 Bioinformatics 2 10-29-2011 01:21 PM

Thread Tools
Old 03-02-2011, 04:39 AM   #1
Location: London

Join Date: Jan 2010
Posts: 65
Default samtools mpileup

Hello all,

I'm using samtools-0.1.12a for snp calling from mutiple samples, after remove the INDELs ( varFilter applied), I find some snp have "N" as ref allele

1 145833120 . N T 11.10 . AF1=0.9402;CI95=0.0625,1;DP=4;DP4=0,0,0,2;GENE=NBPF1;MQ=49;TX=NBPF1:uc001emp.3:in12;TYPE=INTRONIC GT:GQ:PL 1/1:9:0,0,0 1/1:9:0,0,0 1/1:9:0,0,0 1/1:15:42,6,0 1/1:9:0,0,0 1/1:9:0,0,0 1/1:9:0,0,0 1/1:9:0,0,0

Why these loci with "N" as ref are treated as SNPs?

Many thanks
bair is offline   Reply With Quote
Old 03-02-2011, 11:37 PM   #2
Nils Homer
nilshomer's Avatar
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285

Looking in hg19, there is an N base in the reference.
nilshomer is offline   Reply With Quote
Old 03-03-2011, 03:38 AM   #3
Senior Member
Location: Germany

Join Date: Oct 2008
Posts: 415

You can exclude these positions using the -N switch in bcftools view

-N skip sites where REF is not A/C/G/T
colindaven is offline   Reply With Quote

n ref in snps

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 03:43 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO