Bowtie2 hg19 reference for gatk MuTect
3
0
Entering edit mode
5.0 years ago
varsha619 ▴ 90

Hello, I understand that the suggested aligner to use with GATK is bwa. If I want to use Bowtie2 as the aligner, which reference file should I be using? The reference in GATK bundle (Homo_sapiens_assembly19.fasta) does not seem to work with Bowtie2 and using hg19 Bowtie2 index from http://bowtie-bio.sourceforge.net/bowtie2/index.shtml gives a compatibility error with dbsnp_138.hg19.vcf file while running GATK MuTect.

ERROR MESSAGE: Input files dbsnp_138.hg19.vcf and reference have incompatible contigs: No overlapping contigs found.

Is there any way to fix this other than switching to the bwa aligner? Is there a different reference or dbsnp vcf file I could use? Thank you for your help.

gatk bowtie2 • 3.0k views
ADD COMMENT
0
Entering edit mode

The reference in GATK bundle (Homo_sapiens_assembly19.fasta) does not seem to work with Bowtie2

What make you think this reference doesn't work with Bowtie2? Do you encounter any errors? What were the steps you run?

ADD REPLY
0
Entering edit mode

This is the error I get with the fasta file I downlaoded from - ftp://gsapubftp-anonymous@ftp.broadinstitute.org/bundle/hg19/

bowtie2 -p 4 -x ucsc.hg19.fasta -1 R1_001.fastq.gz -2 R2_001.fastq.gz | samtools view -bS - > 001_out.bam

(ERR): "ucsc.hg19.fasta" does not exist or is not a Bowtie 2 index

Exiting now ...

[samopen] no @SQ lines in the header.

[sam_read1] missing header? Abort!

Do I need to create new indices to be able to use this reference file from the GATK bundle with bowtie2?

ADD REPLY
0
Entering edit mode

If you download the sequence from one source and the index from a different source, it is not surprising that they are not be compatible.

ADD REPLY
0
Entering edit mode

I did keep the reference files consistent through the steps. I meant to say that I tried using both files separately. I can successfully align the files using hg19 Bowtie2 indices but then I get a compatibility error in MuTect with dbsnp_138.hg19.vcf file while running GATK MuTect. If I try to use the reference in GATK bundle to align instead, I get the error: file does not exist or is not a Bowtie 2 index.

ADD REPLY
0
Entering edit mode

I have another question.

  • I want to compare GATK/BWA-MEM + Variant Caller vs GATK/BOWTIE2 + Variant Caller. For that purpose, can I use the GATK interleaved fastq file generated according to their Best Practices or it is necessary to put the original 1 2 paired fastq files?
  • This is confusing for me because in Bowtie2 manual there is the option --interleaved that accepts a single interleaved fastq as input but I'm not sure if this is the same file that generated by GATK or not.
  • Using the former option I can achieve a more similar result but I'm afraid of using the wrong input.
ADD REPLY
1
Entering edit mode

can I use the GATK interleaved fastq file generated according to their Best Practices

GATK Best Practices does not generate FASTQs. FASTQs are the input and come from the user.

ADD REPLY
0
Entering edit mode

Update: Yes, the Interleaved fastq file generated from unmapped BAMs (uBAMs) can be used as input for Bowtie2

ADD REPLY
2
Entering edit mode
5.0 years ago
igor 13k

Since GATK does not provide a Bowtie2 index, the index does not exist. Thus, you need to generate it yourself.

ADD COMMENT
0
Entering edit mode

Yes this fixed the issue, my bad. Thank you!

ADD REPLY
1
Entering edit mode

If an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted. You can accept more than one if they work.
Upvote|Bookmark|Accept

ADD REPLY
1
Entering edit mode
5.0 years ago
h.mon 35k

Do I need to create new indices to be able to use this reference file from the GATK bundle with bowtie2?

Yes, creating the bowtie2 indices would help.

ADD COMMENT
0
Entering edit mode

Worked after I created the bowtie2 indices for the GATK bundle hg19 genome file, thank you.

ADD REPLY
0
Entering edit mode
5.0 years ago

Input files dbsnp_138.hg19.vcf and reference have incompatible contigs

Okay, so did you do the obvious thing and check to see that chromosome names are the same between the two files?

ADD COMMENT
0
Entering edit mode

Yes both the files have the format chr1, chr2, ...

ADD REPLY

Login before adding your answer.

Traffic: 2323 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6