user reference vs. built-in index
Hi. I created a workflow to map IGA reads using bowtie and generate a pileup at the end. The workflow is FASTQ Groomer -> Map with Bowtie for Illumina -> SAM-to-BAM -> Generate pileup. I ran the workflow successfully using the built-in mm9 index as reference for both bowtie and the pileup generation. Then I changed the workflow to use a reference file I uploaded (fasta format) and now I get an error on the pileup generation step: An error occurred running this job: /The output file is empty. Your input file may have had no matches, or there may be an error with your input file or settings./ The prior SAM-to-BAM step shows a sizeable BAM file being generated. So I am wondering if the Generate pileup tool requires the reference in a special format. V. Patel Dept. of Human Genetics Emory University School of Medicine
----- Original Message -----
Hi. I created a workflow to map IGA reads using bowtie and generate a pileup at the end. The workflow is FASTQ Groomer -> Map with Bowtie for Illumina -> SAM-to-BAM -> Generate pileup. I ran the workflow successfully using the built-in mm9 index as reference for both bowtie and the pileup generation. Then I changed the workflow to use a reference file I uploaded (fasta format) and now I get an error on the pileup generation step:
An error occurred running this job: The output file is empty. Your input file may have had no matches, or there may be an error with your input file or settings.
The prior SAM-to-BAM step shows a sizeable BAM file being generated. So I am wondering if the Generate pileup tool requires the reference in a special format.
V. Patel Dept. of Human Genetics Emory University School of Medicine
I am a new Galaxy user but it seems to me that the problem lies with using a built-in index vs. your own fasta file that needs an index built for it. I suspect that the parameters for the workflow are set to use built-in indexes. My suggestion is to manually run the steps instead of relying on the workflow. -- Rick Westerman westerman@purdue.edu Bioinformatics specialist at the Genomics Facility. Phone: (765) 494-0505 FAX: (765) 496-7255 Department of Horticulture and Landscape Architecture 625 Agriculture Mall Drive West Lafayette, IN 47907-2010 Physically located in room S049, WSLR building
On 5/13/10 1:14 PM, Rick Westerman wrote:
----- Original Message -----
Hi. I created a workflow to map IGA reads using bowtie and generate a pileup at the end. The workflow is FASTQ Groomer -> Map with Bowtie for Illumina -> SAM-to-BAM -> Generate pileup. I ran the workflow successfully using the built-in mm9 index as reference for both bowtie and the pileup generation. Then I changed the workflow to use a reference file I uploaded (fasta format) and now I get an error on the pileup generation step:
An error occurred running this job: The output file is empty. Your input file may have had no matches, or there may be an error with your input file or settings.
The prior SAM-to-BAM step shows a sizeable BAM file being generated. So I am wondering if the Generate pileup tool requires the reference in a special format.
V. Patel Dept. of Human Genetics Emory University School of Medicine
I am a new Galaxy user but it seems to me that the problem lies with using a built-in index vs. your own fasta file that needs an index built for it. I suspect that the parameters for the workflow are set to use built-in indexes. My suggestion is to manually run the steps instead of relying on the workflow.
I changed the workflow to use my own uploaded reference file, not the built-in index, but judging from the nomenclature I suspect I may have to process (index?) the reference file first into a format usable by the pileup generation tool. However I do not see how I can do that.
I changed the workflow to use my own uploaded reference file, not the built-in index, but judging from the nomenclature I suspect I may have to process (index?) the reference file first into a format usable by the pileup generation tool. However I do not see how I can do that.
Galaxy should do this automatically. At least it did it for me yesterday (admittedly using SOLiD data instead of Illumina. I still recommend manually running each step instead of relying on a workflow. -- Rick Westerman westerman@purdue.edu Bioinformatics specialist at the Genomics Facility. Phone: (765) 494-0505 FAX: (765) 496-7255 Department of Horticulture and Landscape Architecture 625 Agriculture Mall Drive West Lafayette, IN 47907-2010 Physically located in room S049, WSLR building
The indexing should be automatic. We'll need to look at it in more detail. -- jt (composed on my phone) On May 13, 2010, at 1:22 PM, Viren Patel <vcpatel@emory.edu> wrote:
On 5/13/10 1:14 PM, Rick Westerman wrote:
----- Original Message -----
Hi. I created a workflow to map IGA reads using bowtie and generate a pileup at the end. The workflow is FASTQ Groomer -> Map with Bowtie for Illumina -> SAM-to-BAM -> Generate pileup. I ran the workflow successfully using the built-in mm9 index as reference for both bowtie and the pileup generation. Then I changed the workflow to use a reference file I uploaded (fasta format) and now I get an error on the pileup generation step:
An error occurred running this job: The output file is empty. Your input file may have had no matches, or there may be an error with your input file or settings.
The prior SAM-to-BAM step shows a sizeable BAM file being generated. So I am wondering if the Generate pileup tool requires the reference in a special format.
V. Patel Dept. of Human Genetics Emory University School of Medicine
I am a new Galaxy user but it seems to me that the problem lies with using a built-in index vs. your own fasta file that needs an index built for it. I suspect that the parameters for the workflow are set to use built-in indexes. My suggestion is to manually run the steps instead of relying on the workflow.
I changed the workflow to use my own uploaded reference file, not the built-in index, but judging from the nomenclature I suspect I may have to process (index?) the reference file first into a format usable by the pileup generation tool. However I do not see how I can do that. <vcpatel.vcf> _______________________________________________ galaxy-user mailing list galaxy-user@lists.bx.psu.edu http://lists.bx.psu.edu/listinfo/galaxy-user
Hi, You will not have to generate the index separately (you actually can't do this--if you're not using one of the built-in indexes, you have to supply a fasta file to be indexed on each mapping run). Changing to an uploaded reference should work within your workflow, just as it should work if you ran each step outside of the workflow. I can look into what's going wrong here if you will share both your workflow and your history with me (kpvincent@bx.psu.edu). To share your workflow: From the workflow menu, click on the arrow next to your workflow and select Share or Publish. Then click the Share with a user button near the bottom of the page and enter my email address in the following box. To share your history: From the history pane, click on the Options button in the upper right and select Share or Publish. Then click the Share with a user button and enter my email address in the following box. Regards, Kelly On Thu 13 May 2010, at 1:22 PM, Viren Patel wrote:
On 5/13/10 1:14 PM, Rick Westerman wrote:
----- Original Message -----
Hi. I created a workflow to map IGA reads using bowtie and generate a pileup at the end. The workflow is FASTQ Groomer -> Map with Bowtie for Illumina -> SAM-to-BAM -> Generate pileup. I ran the workflow successfully using the built-in mm9 index as reference for both bowtie and the pileup generation. Then I changed the workflow to use a reference file I uploaded (fasta format) and now I get an error on the pileup generation step:
An error occurred running this job: The output file is empty. Your input file may have had no matches, or there may be an error with your input file or settings.
The prior SAM-to-BAM step shows a sizeable BAM file being generated. So I am wondering if the Generate pileup tool requires the reference in a special format.
V. Patel Dept. of Human Genetics Emory University School of Medicine
I am a new Galaxy user but it seems to me that the problem lies with using a built-in index vs. your own fasta file that needs an index built for it. I suspect that the parameters for the workflow are set to use built-in indexes. My suggestion is to manually run the steps instead of relying on the workflow.
I changed the workflow to use my own uploaded reference file, not the built-in index, but judging from the nomenclature I suspect I may have to process (index?) the reference file first into a format usable by the pileup generation tool. However I do not see how I can do that. <vcpatel.vcf>_______________________________________________ galaxy-user mailing list galaxy-user@lists.bx.psu.edu http://lists.bx.psu.edu/listinfo/galaxy-user
participants (4)
-
James Taylor -
Kelly Vincent -
Rick Westerman -
Viren Patel