Question: mpileup calling every ALT <X>
0
gravatar for rquigl01
22 months ago by
rquigl0110
rquigl0110 wrote:

Im using the mpileup application on galaxy and every base is being called with the ALT <x>. It has 110,000,000 lines as well which seems to be too much. I'm just trying to call SNPs and INDELS. Here is a screenshot. https://gyazo.com/f70ca31ebd7d47421b4e702ae3cb91f3 Does anybody know what I am doing wrong?

galaxy mpileup • 900 views
ADD COMMENTlink modified 22 months ago by Jennifer Hillman Jackson24k • written 22 months ago by rquigl0110
1
gravatar for Jennifer Hillman Jackson
22 months ago by
United States
Jennifer Hillman Jackson24k wrote:

Hello,

This prior post at Biostars.org has help for the same issue: https://www.biostars.org/p/161981/. See Devon's reply (currently the second answer).

In Galaxy: When using MPilieup, output BCF (instead of VCF) then convert to VCF excluding non-ATGC bases with the tool bcftools view.

Hopefully this helps! Jen, Galaxy team

ADD COMMENTlink written 22 months ago by Jennifer Hillman Jackson24k

on bcftools view, there is an option to skip sites where the REF field is not A/C/G/T, but my problem is in the ALT column. What option on bcftools are you referring to?

ADD REPLYlink written 22 months ago by rquigl0110
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 16.09
Traffic: 113 users visited in the last hour