Question: Suggestions for variant callers/ problem with BCF view
gravatar for murphykatharyne
16 months ago by
murphykatharyne0 wrote:

Samtools mpileup only outputs genotypes fields for samples when a BCF output is selected. The BCF view tool in galaxy is supposed to convert mpileup's BCF variants file to VCF variants file, but unfortunately does not change the binary encoding. I am attempting to call variants for two drosophila samples against the latest dm6 reference, but other variant callers do not seem to offer the same version reference genome as an index file or do not have the option to output genotypes. I am specifically interested in determining homozygous vs heterozygous variant sites. Any suggestions on variant callers that are fit for this specific task would be greatly appreciated.

bcf galaxy tools vfc • 395 views
ADD COMMENTlink modified 16 months ago by Jennifer Hillman Jackson25k • written 16 months ago by murphykatharyne0
gravatar for Jennifer Hillman Jackson
16 months ago by
United States
Jennifer Hillman Jackson25k wrote:


The reference genome dm6 should be indexed for most variant callers at Galaxy Main

GATK is one exception and those tools are deprecated. GEMINI is also a noted exception; at Galaxy Main, only hg19 is indexed/supported.

If the tool you want to use does (not including GATK or GEMINI) do not have dm6 indexed - please comment back about which tool so we can address that if possible. Meanwhile, a custom reference genome can be used with nearly all other tools if needed, this is how:

As an alternative, pileup format (instead of VCF) can be directly given to Varscan to compute the consensus genotype. After, tools in the group VCF Manipulation can be used to further analyze those results. Or you can use Freebayes in one step (input the mapped BAM data to produce a VCF with genotype information).

Galaxy tutorials:

Thanks, Jen, Galaxy team

ADD COMMENTlink modified 16 months ago • written 16 months ago by Jennifer Hillman Jackson25k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 175 users visited in the last hour