Question: VCF combine only showing chromosome 1 in the output
1
gravatar for drsrideviunni
11 weeks ago by
drsrideviunni20 wrote:

Hi, I tried many times to combine 32 .vcf files in Galaxy using VCF Combine. In the output file I can see only a few lines (chromosome 1). I downladed the file and the downloaded file also contain only few lines in it. Using the edit attributes tool I edited the Number of comment lines to larger numbers, and also did unclick, so that it does not have any limit. But, everytime it showed up the same result. If you could help me with this it will be highly appreciated.

Thank you, Sri

snp galaxy • 179 views
ADD COMMENTlink modified 9 weeks ago by Guy Reeves1000 • written 11 weeks ago by drsrideviunni20

Please help me if any one has a solution for this. Hope this problem is common for everyone who uses galaxy for this application.

ADD REPLYlink written 10 weeks ago by drsrideviunni20
1
gravatar for drsrideviunni
10 weeks ago by
drsrideviunni20 wrote:

Hi, Thank you for the reply to my post. I sorted the .vcf files using vcf sort of Galaxy and ran the vcf combine. But the results are still the same. In my .vcf and sorted files I can see all the other chromosomes, but it is not showing up in the combined file. I really appreciate your help. Sri

ADD COMMENTlink written 10 weeks ago by drsrideviunni20

Where are you working? http://usegalaxy.org? If so a bug report can be sent in and we can review. This is the only report of a problem with the tool so far that didn't involve a memory problem (tool ran out of resources). https://galaxyproject.org/issues/#usage-problem-reporting

I am wondering if the tool is running out of memory in your case however some error or warning message should result from that situation now.

ADD REPLYlink modified 9 weeks ago • written 9 weeks ago by Jennifer Hillman Jackson22k

Yes, I am working on http://usegalaxy.org When I do the VCFcombine I do not get any error messages. So there is no bug report available. The program runs perfectly, and I get a single file as output which contains all my samples, but it shows only the top 130 lines (chromosome 1, not sure if it is complete data of chromosome 1). Thanks

ADD REPLYlink written 9 weeks ago by drsrideviunni20
1
gravatar for Guy Reeves
9 weeks ago by
Guy Reeves1000
Germany
Guy Reeves1000 wrote:

Hi Siri Can I sugest you try another tool. I am assuming that you are using usegalaxy.org. So I suggest you try "Combine Variants (Galaxy Version 0.0.4)" and see if you get the same result.

If you are using your own galaxy instance there are a number of other tools you could try Tell us what happens. Cheers Guy

ADD COMMENTlink modified 9 weeks ago • written 9 weeks ago by Guy Reeves1000

Agree - this can be an alternative depending on the desired output (merged vs combined).

ADD REPLYlink written 9 weeks ago by Jennifer Hillman Jackson22k

what is the difference between combine versus merge VCF?

ADD REPLYlink written 9 weeks ago by Guy Reeves1000

Hello, As suggested I did the Combine variants. I tried both .vcf files and also the sorted files. This created error messages.
Unlike VCFCombine (which produces a single file), Combine Variants program create separate files for each input file, which I do not understand why. Please suggest me if anything can be rectified from my end. Thank you Sri

ADD REPLYlink written 9 weeks ago by drsrideviunni20

The errors might be informative. Would you send one in as a bug report please? This will help get to the bottom of the problem (core tool issue vs input data problem vs usage clarification vs tool running out of memory). A shared history link can also be sent in from any history at http://usegalaxy.org, even if there are not red error datasets, directly to galaxy-bugs@lists.galaxyproject.org.

The simplest way to do this is: if the runs from the two tools are in different histories, use the one with the red error to send in a bug report and just include a shared history link to the other history in the comments. Note the datasets involved and make sure they are left in an undeleted state (inputs and outputs). Please also include a link to this Biostars post so we can associate the two. I didn't see any that look like yours so far, so correct me if one was sent already with this information - maybe by replying to the galaxy-bugs email address with the extra info, even if only a link to this post (if all runs are in that same history).

ADD REPLYlink modified 9 weeks ago • written 9 weeks ago by Jennifer Hillman Jackson22k

Hello, I have already sent the red error bug report yesterday. Also, I tried to share the history through galaxy-bugs@lists.galaxyproject.org. However, it shows invalid galaxy user ID. Pls suggest me if there is another email ID.

Thank you.

ADD REPLYlink written 9 weeks ago by drsrideviunni20

I just tried this tool and I get only one output file and one log file.

ADD REPLYlink written 9 weeks ago by Guy Reeves1000

Could you please guide me the steps so that I can also follow the same method.

Thank you.

ADD REPLYlink written 9 weeks ago by drsrideviunni20

Combining two VCFs I did this

https://www.dropbox.com/s/cm0z305ezc3yxyo/Screen%20Shot%202017-07-13%20at%2007.58.01.png?dl=0

and got

https://www.dropbox.com/s/idr09tsb5ieezx1/Screen%20Shot%202017-07-13%20at%2007.58.24.png?dl=0

!There us a difference between combine and merge !

ADD REPLYlink modified 9 weeks ago • written 9 weeks ago by Guy Reeves1000
0
gravatar for Jennifer Hillman Jackson
10 weeks ago by
United States
Jennifer Hillman Jackson22k wrote:

Hello,

Try using the tool VCFsort first then combine. It may be that your data is not ordered.

Thanks, Jen, Galaxy team

ADD COMMENTlink written 10 weeks ago by Jennifer Hillman Jackson22k

Hi, Thank you for the reply to my post. I sorted the .vcf files using vcf sort of Galaxy and ran the vcf combine. But the results are still the same. In my .vcf and sorted files I can see all the other chromosomes, but it is not showing up in the combined file. I really appreciate your help. Sri

ADD REPLYlink written 10 weeks ago by drsrideviunni20
0
gravatar for Guy Reeves
9 weeks ago by
Guy Reeves1000
Germany
Guy Reeves1000 wrote:

If you post the shared link to the history here that may also help me help you.

ADD COMMENTlink written 9 weeks ago by Guy Reeves1000
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: 114 users visited in the last hour