SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
GATK variant calling for multiple exome sehrrot Bioinformatics 2 04-24-2012 05:09 AM
GATK variant calling on uniquely mapped reads? pravee1216 Bioinformatics 5 03-26-2012 09:19 AM
GATK excludes some samples for cohort variant calling liu_xt005 Bioinformatics 2 02-01-2012 12:58 PM
minimum depth variant calling samtools/gatk m_elena_bioinfo Bioinformatics 1 12-06-2011 09:31 AM
variant calling kjaja Bioinformatics 1 11-04-2011 08:16 AM

Reply
 
Thread Tools
Old 06-28-2012, 01:36 PM   #1
pooja_singh
Junior Member
 
Location: Austria

Join Date: Jun 2012
Posts: 2
Default Problems using GATK for variant calling

Hi,

I am busy working with Eucalyptus transcriptome data generated using the Illumina platform. I mapped the paired end reads using TopHat and cufflinks. Then I defined read groups using Picard tools and sorted and indexed the BAM files using SAM tools. Now I am trying to use GATK to call variants using the BAM files as input. However, i can't carry out the step where you 'realign around known indels' since there aren't any for this species. So when I skip this step and go straight to calling SNP and indels i get an unknown error:

old_5354, scaffold_5355, scaffold_5356, scaffold_5357, scaffold_5358, scaffold_5359, scaffold_5360, scaffold_5361, scaffold_5362, scaffold_5363, scaffold_5364, scaffold_5365, scaffold_5366, scaffold_5367, scaffold_5368, scaffold_5369, scaffold_5370, scaffold_5371, scaffold_5372, scaffold_5373, scaffold_5374, scaffold_5375, scaffold_5376, scaffold_5377, scaffold_5378, scaffold_5379]
##### ERROR --------------------------------------------------------------------------

I really am not sure how to fix this since this is the first time that I am working with RNAseq data. Does anyone have any solutions?
pooja_singh is offline   Reply With Quote
Old 07-07-2012, 11:23 AM   #2
jflowers
Member
 
Location: New York, NY

Join Date: Oct 2011
Posts: 41
Default

You don't need to have known indels. The --known argument in the RealignerTargetCreator and IndelRealigner is optional. You should however run both the target creator and realigner for best alignment results.
jflowers is offline   Reply With Quote
Reply

Tags
gatk, rnaseq, tophat, transcript mapping

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 09:21 PM.


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