Question: MPileup output errror ALT <X> QUAL 0
0
gravatar for edouglasrobertson
20 months ago by
edouglasrobertson • 0 wrote:

When using MPileup I received an error in that all the output was for chr10 (even when my file only contained info for chr20), the ALT position was always X, and the quality score was 0. Therefore for some reason the variant calling failed but I don't know why. Can anyone help?

mpileup • 598 views
ADD COMMENT • link • modified 20 months ago by Jennifer Hillman Jackson ♦ 25k • written 20 months ago by edouglasrobertson • 0

There might be a reference genome mismatch problem. Double check that the same genome was used for the mapping step as was assigned for the option "Using reference genome" on the Mpileup tool form. If using a Custom reference genome, ensure that it was first run through the tool NormalizeFasta with the option to wrap at 80 bases and clip the title line at the first whitespace (or is in that format already).

If you cannot determine the problem and are either working at http://usegalaxy.org or can reproduce the problem there, an email with a shared history link can be sent to galaxy-bugs@lists.galaxyproject.org. Please also include a link to this post.

Thanks, Jen, Galaxy team

ADD REPLY • link written 20 months ago by Jennifer Hillman Jackson ♦ 25k

I've sent a link to my history. Its not a custom reference so I don't see what the problem is.

ADD REPLY • link written 20 months ago by edouglasrobertson • 0
0
gravatar for Jennifer Hillman Jackson
20 months ago by
United States
Jennifer Hillman Jackson ♦ 25k wrote:

Hello,

The input fastq sequence quality is very, very low. It may not be usable. Run FastQC to review.

The quality impacts mapping (very few sequences map at all) plus all downstream steps. GATK tools failing is a known - these are deprecated and not recommended. Alternative variant detection methods are covered in the tutorials.

Thanks, Jen, Galaxy team

ADD COMMENT • link written 20 months ago by Jennifer Hillman Jackson ♦ 25k
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: 175 users visited in the last hour