Building an indexed genome file for GATK tools

Hi,

Does anyone have any pointers for indexing a reference sequence for use with GATK tools? I have indexed reference sequences that I use for BWA and minimap2, but I can’t get the reference sequences to appear with GATK. I receive an error when I use the “GATK-sorted picard indexes builder”.

Thanks!

1 Like

Hi @jsalem

Nearly all GATK tools, including the Data Manager, as wrapped for Galaxy, are considered deprecated. These were all based on earlier releases of GATK.

Currently, the only GATK4 tool that has been wrapped is:

For native genome fasta files to be accessible to this tool (option: Choose the source for the reference list) – adding the genome with the fetch DM is probably enough, although to really make new genomes useful there is a short-list of core recommended indexes. All have Data Mangers. You can certainly run more DMs after those (GATK4 Mutect2 won’t need more… but other tools can). See the topic below for help – the same process applies to all genomes you plan to index, not just the one referenced in that particular Q&A.

Hi @jennaj,

Thanks for getting back to me. I was hoping to use GATK indel realigner, which requires a reference sequence. Do you know of any other good realigner tools?

1 Like

Hi – Review the LoFreq tools