Question: Mark Duplicate Reads Error- "Value was put into PairInfoMap more than once"
2
gravatar for JP
2.6 years ago by
JP20
United States
JP20 wrote:

Hi everyone,

I am relatively new to galaxy and have been trying to set up a workflow for variant calling used paired end Illumina data (MiSeq).  I aligned the reads using BWA-MEM, but have hit a wall using the Mark Duplicate Reads tool.  I’ve repeatedly received the error “Value was put into PairInfoMap more than once”, and I don’t seem to be able to correct the problem.

I did use the –M option (Mark Shorter Split Hits as Secondary).  Here are the steps I took:

1) FASTQ Groomer on both Fastq Files

2) BWA-MEM with groomed FASTQ’s (-M)

3) Addorreplacereadgroups

4) Clean Sam (on .BAM file)

4) Paired read mate fixer

5) MarkDuplicate Reads (Failed)

I ran the analysis on the public galaxy server.  Any advice/suggestions would be greatly appreciated.

Thanks,


JP

bwa-mem gatk • 1.0k views
ADD COMMENTlink modified 2.6 years ago by Jennifer Hillman Jackson24k • written 2.6 years ago by JP20
0
gravatar for Jennifer Hillman Jackson
2.6 years ago by
United States
Jennifer Hillman Jackson24k wrote:

Hello,

If you would like to submit this as a bug report, we can examine the data in detail and provide feedback. Send on in from any of the red error datasets using the green bug icon. Please put a link to this biostars post in the comments so we can link the two.

Thanks, Jen, Galaxy team

ADD COMMENTlink written 2.6 years ago by Jennifer Hillman Jackson24k

The report is submitted, thank you for the help!

JP

ADD REPLYlink written 2.6 years ago by JP20
1

Hi JP, I do not see anything obviously wrong. But I did find a similar post from about 7 months ago. This was someone else, correct? https://www.biostars.org/p/130549/

The MAPQ error has been moved to a warning since that post was written. Perhaps this other fatal error needs to be migrated to be a warning as well. At this point I am opening a ticket for our team to review together. It can be followed here: https://trello.com/c/zhyLdScr

Thanks for reporting the problem and hopefully we can resolve it quickly. It is entirely possibly that I missed an input requirement, but we want to rule out a bug/usage issue if there is one. Jen, Galaxy team

ADD REPLYlink written 2.6 years ago by Jennifer Hillman Jackson24k

Hi Jen, 

You are correct, that was another user.....I came across that post as well.  Thank you very much for starting a ticket, I'm very interested to see the result of your review!

Best,


JP

ADD REPLYlink written 2.6 years ago by JP20
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: 88 users visited in the last hour