Question: VCF file REF column only contains N
1
gravatar for genie66
3.3 years ago by
genie6620
United States
genie6620 wrote:

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 SC_Stressed.raw.bcf -v vcfutils.pl varFilter -D100 > SC_Stressed.flt100.vcf

OUTPUT:

INFO= <ID=PR,Number=1,Type=Integer,Description="# permutations yielding a smaller PCHI2.">
CHROM   POS ID  REF ALT QUAL    FILTER  INFO    FORMAT SRR1481944_Aligned.sorted.bam

    chr1    16961   .   N   T   39.8    .   DP=2;VDB=5.960000e-02;AF1=1;AC1=2;DP4=0,0,2,0;MQ=60;FQ=-33;AF1=1;AC1=2;FQ=-33   GT:PL:GQ    1/1:71,6,0:10
    chr1    16962   .   N   C   39.8    .   DP=2;VDB=5.960000e-02;AF1=1;AC1=2;DP4=0,0,2,0;MQ=60;FQ=-33;AF1=1;AC1=2;FQ=-33   GT:PL:GQ    1/1:71,6,0:10
    chr1    16963   .   N   G   39.8    .   DP=2;VDB=5.960000e-02;AF1=1;AC1=2;DP4=0,0,2,0;MQ=60;FQ=-33;AF1=1;AC1=2;FQ=-33   GT:PL:GQ    1/1:71,6,0:10
    chr1    16964   .   N   T   39.8    .   DP=2;VDB=5.960000e-02;AF1=1;AC1=2;DP4=0,0,2,0;MQ=60;FQ=-33;AF1=1;AC1=2;FQ=-33   GT:PL:GQ    1/1:71,6,0:10
mpileup samtools ref vcf • 1.3k views
ADD COMMENTlink modified 3.3 years ago by jackfrost219970 • written 3.3 years ago by genie6620
0
gravatar for jackfrost2199
3.3 years ago by
Washington D.C.
jackfrost219970 wrote:

If you run your bcftools line like:

bcftools view SC_Stressed.raw.bcf -v | vcfutils.pl varFilter -D100 > SC_Stressed.flt100.vcf

Do you get a better result? It looks like the pipe to vcfutils.pl might have been missed. Otherwise everything looks fine to me so its strange it is reporting things with a read depth of 2 and not a minimum of 100.

ADD COMMENTlink written 3.3 years ago by jackfrost219970

When I add pipe, I am getting nothing after the header description. Do you find any mistake in the command?

ADD REPLYlink written 3.3 years ago by genie6620

I think everything looks good. What is the size of "SC_Stressed.raw.bcf" after you do the mpileup stage but before the bcftools?

ADD REPLYlink written 3.3 years ago by jackfrost219970
1

Thanks @jackfrost this following command worked for me. Hope it help someone

samtools mpileup -uf hg19.fa SRR14819$i.sorted.bam | bcftools call -mv > SRR14819$i.Control.vcf
ADD REPLYlink written 3.3 years ago by genie6620
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: 771 users visited in the last hour