SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
allele information from samtools mpileup contents fatakias Bioinformatics 2 11-07-2012 08:49 AM
References/information for Extended BAQ algorithm in SAMtools??? Genomics101 Bioinformatics 0 07-31-2012 08:02 AM
samtools mpileup to get vcf files including GL information xiangfeiloulan Bioinformatics 0 07-01-2012 06:57 PM
samtools question - redefine mate information after applying filters caswater Bioinformatics 0 03-07-2012 07:37 PM
about problems run into for samtools henry Bioinformatics 1 09-09-2009 01:08 AM

Reply
 
Thread Tools
Old 09-26-2012, 07:35 AM   #1
Rcucullatus
Junior Member
 
Location: UK

Join Date: Aug 2012
Posts: 7
Default Problems with outputting GT information using samtools

Hi everyone,

I'm trying to calculate heterozigosity for one individual genome using vcftools. For that, I need to have GT information in the FORMAT field for all sites. However, using the command "samtools mpileup -C50 -uDgf ref.fa output.sorted.bam | bcftools view -cg - > output.vcf" in samtools, only adds the GT information for variable sites relative to the reference, which is apparently a problem for vcftools. Anyone has any idea how to tell samtools to include the GT field for all sites? I thought the -g option in both mpilup and bcfview would take care of it, but apparently not.

Thanks!
Rcucullatus 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 06:49 AM.


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