Question: Where to get SNPs and Indels?
0
gravatar for PaulW
8 months ago by
PaulW60
PaulW60 wrote:

I downloaded some files from Broad Institute but Realigner Target Creator won't work with these files because the sort order is different to my BAM file.

/hg19/Mills_and_1000G_gold_standard.indels.hg19.sites.vcf /hg19/dbsnp_138.hg19.excluding_sites_after_129.vcf /hg19/ucsc.hg19.fasta

I can't sort the BAM file to match these files because in Galaxy the REFERENCE parameter is hidden for the SortSAM tool.

Where can I get these files in the correct sort order?

realigner target creator • 249 views
ADD COMMENTlink modified 8 months ago • written 8 months ago by PaulW60
0
gravatar for Jennifer Hillman Jackson
8 months ago by
United States
Jennifer Hillman Jackson23k wrote:

Hello,

Assign the database "hg19" to the BAM datasets. This is how certain tools interpret the reference genome. I am assuming that you have installed this reference genome on your server already or it exists on the server you are working on.

Click on the pencil icon for the dataset to make the change. Next time, you can assign the database in the Upload tool if that is used.

Sorting help: https://github.com/jennaj/support-prior-qa/wiki/Sort-your-inputs

Note: GATK tools at http://usegalaxy.org are not natively indexed for hg19 and jobs run with a hg19 custom reference genome will often fail for memory. Instead, hg_g1k_v37 is the human build indexed (1000 Genomes version). Also - the tools here are deprecated and may not work - we will be hiding them soon so that they are not used. Most other tools (Samtools, Picard, Mapping and Variant analysis tools, etc) are indexed for hg19.

Thanks, Jen, Galaxy team

ADD COMMENTlink written 8 months ago by Jennifer Hillman Jackson23k
0
gravatar for PaulW
8 months ago by
PaulW60
PaulW60 wrote:

Thanks for your help. It turns out that Picard ReorderSAM allows sorting the BAM to match the reference file from Broad which then makes everything work fine :-)

ADD COMMENTlink written 8 months ago by PaulW60
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: 103 users visited in the last hour