SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Varscan2 multi-variant samples VCF output with empty sample call FORMAT column? BhariD Bioinformatics 0 04-30-2015 01:56 PM
Does anybody know how to convert dbsnp to .vcf (variant calling format)? alexbmp Bioinformatics 13 04-05-2015 10:24 PM
what is DP4 in VCF (Variant Call Format) and AF1? gcrdb Bioinformatics 8 01-20-2015 09:37 AM
Asking for vcf file when calling SNP using GATK likebiology Bioinformatics 3 07-07-2014 11:59 PM
Help figuring out variant call format (vcf) bmtb Bioinformatics 1 12-11-2012 12:09 AM

Reply
 
Thread Tools
Old 06-14-2017, 10:32 AM   #1
MGCBrown
Member
 
Location: Canada

Join Date: May 2017
Posts: 15
Question samtools SNP calling pipeline outputs VCard file instead of variant call format (VCF)

Hello,

I used this workflow to call SNPs on bam files that were aligned to a reference de novo transcriptome assembly:

samtools mpileup -uf normalized_Trinity.fasta aln2.bam aln2.bam | bcftools view -bvcg \
- > var.raw.bcf

bcftools view var.raw.bcf | vcfutils.pl varFilter -D100 > var.flt.vcf



However, the final var.flt.vcf output is a VCard file in my Bitvise ssh client.

Any insight as to why this might be would be greatly appreciated!

Thanks,

Mike
MGCBrown is offline   Reply With Quote
Old 06-14-2017, 11:44 AM   #2
HESmith
Senior Member
 
Location: Bethesda MD

Join Date: Oct 2009
Posts: 498
Default

It's not actually VCard format. It's just that VCard uses the same file extension (.vcf) as VCF, and your computer is configured to read the former by default. Open the VCF in a flat text editor and you'll see.
HESmith is offline   Reply With Quote
Old 06-14-2017, 04:38 PM   #3
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

For a Windows text editor, I recommend Notepad++. Though normally I find Excel slightly easier for VCFs, if they are small enough.
Brian Bushnell is offline   Reply With Quote
Old 06-15-2017, 05:55 AM   #4
MGCBrown
Member
 
Location: Canada

Join Date: May 2017
Posts: 15
Default

Thanks for the help! I see that when I open my ssh client on a different computer, the files are no longer using the VCard extension.

Thanks again.
MGCBrown is offline   Reply With Quote
Reply

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 11:22 PM.


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