SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
some error message from BFAST match lei Bioinformatics 6 11-23-2015 04:42 AM
GATK: -I parameter in RealignerTargetCreator jorge Bioinformatics 2 10-18-2011 04:10 AM
GATK RealignerTargetCreator -B and -D options efoss Bioinformatics 1 10-03-2011 09:13 AM
Error during Bfast match step James W. MacDonald Bioinformatics 2 01-05-2011 07:15 AM
454IsotigsLayout.txt - contigs don't match to any isotig? dschika Bioinformatics 9 10-11-2010 05:58 AM

Reply
 
Thread Tools
Old 06-27-2011, 07:20 AM   #1
Hkins552
Member
 
Location: Baltimore MD

Join Date: Jun 2011
Posts: 18
Question GATK RealignerTargetCreator Error = Contigs do not match

I'm trying to do a local realignment using GATK and seem to be stuck at the first step. I've created my ref.dict file and then tried the following script:

$ java -jar gatk -T RealignerTargetCreator -I input.bam -R genome.fasta -o output.intervals

An error was generated which said the contigs in my reference fasta and input bam do not match. Ex: Chrom 1 vs. 1

I attempted to use Picardtools ReorderSam command to at least order them correctly, but it will not recognize that "1" in the input.bam represents "chrom 1" in the reference. Is there a way to rename the contigs on one of these files such that ReorderSam can be used?

Thanks for your input!
Hkins552 is offline   Reply With Quote
Old 07-01-2011, 05:59 AM   #2
Hkins552
Member
 
Location: Baltimore MD

Join Date: Jun 2011
Posts: 18
Default Contigs Do Not Match

Still running into the same error. Any input would be greatly appreciated In the mean time, I will likely write a script to rename the contigs in the file and see if that pans out.
Hkins552 is offline   Reply With Quote
Reply

Tags
contigs, gatk, picard tools, reordersam

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 07:04 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO