SEQanswers

Go Back   SEQanswers > Applications Forums > Genomic Resequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
GATK error sridhar28 Bioinformatics 6 06-17-2013 12:53 AM
GATK IndelGenotyperV2 ERROR wanguan2000 Bioinformatics 2 07-26-2011 01:38 AM
GATK error poisson200 Bioinformatics 3 07-08-2011 10:45 PM
GATK error nvteja Bioinformatics 1 12-12-2010 08:47 AM
GATK GenomicAnnotator Error Michael Jones Bioinformatics 1 12-06-2010 05:41 PM

Reply
 
Thread Tools
Old 09-09-2013, 09:50 PM   #1
vishnuamaram
Member
 
Location: india

Join Date: Jun 2013
Posts: 42
Default GATK-usage-ERROR-Helpreg

Hi Seqanswers community,


I am doing human whole genome sequence analysis, of the data obtained Illumina Hiseq 2000.

This is the pipeline i followed & get errors in using GATK


a) BWA ------> SAM output
b) samtools ---> sortedindexedbamoutput
c) picard ------> for readgroup replacement & removing duplicates
d) GATK ----> for indelrealignment, BQSR, variant calling.


when i give the following command for GATK-indel realignment, i get the error.

U]cd used:-[/U]

java -Xmx30g -jar GenomeAnalysisTK-2.6-5-gba531bd/GenomeAnalysisTK.jar -T IndelRealigner \ -R INDEX_REF_GENOME_HG19/hg19.fa \ -I WG_PBMC_SAM9_CORDSORT_RMDUP.BAM -targetintervals intervalListFromRTC.intervals \ -o WG_PBMC_SAM9_coordsort_realigned.bam


the errors are :-



##### ERROR ------------------------------------------------------------------------------------------
##### ERROR A USER ERROR has occurred (version 2.6-5-gba531bd):
##### ERROR The invalid arguments or inputs must be corrected before the GATK can proceed
##### ERROR Please do not post this error to the GATK forum
##### ERROR
##### ERROR See the documentation (rerun with -h) for this tool to view allowable command-line arguments.
##### ERROR Visit our website and forum for extensive documentation and answers to
##### ERROR commonly asked questions http://www.broadinstitute.org/gatk
##### ERROR
##### ERROR MESSAGE: Argument with name '--targetIntervals' (-targetIntervals) is missing.
##### ERROR ------------------------------------------------------------------------------------------



kindly give me solutions, will be very thankful to you.

Please do let me know the proper usage of GATK.

Thanking you all,

Vishnu.
vishnuamaram is offline   Reply With Quote
Old 09-10-2013, 10:04 AM   #2
sheenams
Member
 
Location: Seattle, Washington

Join Date: Oct 2011
Posts: 15
Default

targetintervals intervalListFromRTC.intervals

needs to be:

targetIntervals intervalListFromRTC.intervals
sheenams is offline   Reply With Quote
Old 10-02-2013, 04:43 AM   #3
vishnuamaram
Member
 
Location: india

Join Date: Jun 2013
Posts: 42
Default

Hi sheenams,

Thank you very much.

I had tried using targetIntervals, but stil again it throws error as follows

##### ERROR MESSAGE: Couldn't read file Sample_WG-9/RAWFASTQ_CAT_FILES/intervalListFromRTC.intervals because The interval file does not exist.


How should i proceed exactly with GATK .

Guys and all members, i request all to help me using GATK with proper commands for indel realignment, BQSR and variant calling.

Thank you,
Vishnu.
vishnuamaram is offline   Reply With Quote
Old 10-08-2013, 08:15 PM   #4
whataBamBam
Member
 
Location: Italy

Join Date: May 2013
Posts: 27
Default

Hi Vishnu.. looks like the same error you had on the other thread.. sorry, I thought I'd replied to that but my post didn't register, internet here is quite flaky.

so you tried with relative pathname (i.e. just the file name) and it crashed with

##### ERROR MESSAGE: Couldn't read file Sample_WG-9/RAWFASTQ_CAT_FILES/intervalListFromRTC.intervals because The interval file does not exist.

It's looking in the current working directory for the file and finding it doesn't exist.

And you tried (from other thread) full path

and got

##### ERROR MESSAGE: Couldn't read file /data/odity/Project_Blood-GNPC-464/Sample_WG-9/RAWFASTQ_CAT_FILES/intervalListFromRTC.intervals because The interval file does not exist

so maybe the path is wrong..

what happens when you type

ls /data/odity/Project_Blood-GNPC-464/Sample_WG-9/RAWFASTQ_CAT_FILES/intervalListFromRTC.intervals

at the command line

if it says..

ls: cannot access /data/odity/Project_Blood-GNPC-464/Sample_WG-9/RAWFASTQ_CAT_FILES/intervalListFromRTC.intervals: No such file or directory

then your path is wrong.

And you can find the right one by doing..

/data [HIT TAB KEY]

/data/o [HIT TAB KEY]

/data/odity/P [HIT TAB KEY]

etc.. unix will tab complete the paths for you and you get the path correct.
whataBamBam is offline   Reply With Quote
Old 10-08-2013, 10:13 PM   #5
vishnuamaram
Member
 
Location: india

Join Date: Jun 2013
Posts: 42
Default

Hey whataBamBam,

I had got your responses on other thread too.

Actually, i didnt generate interval file,

My question is how to generate an interval file.

Kindly help me in generating an interval file by giving a proper command usage.

Thank you in tons,
Vishnu.
vishnuamaram is offline   Reply With Quote
Old 10-09-2013, 12:57 AM   #6
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Quote:
Originally Posted by vishnuamaram View Post
Hey whataBamBam,

I had got your responses on other thread too.

Actually, i didnt generate interval file,

My question is how to generate an interval file.

Kindly help me in generating an interval file by giving a proper command usage.

Thank you in tons,
Vishnu.
Did you read my reply to the same question in the other thread? If the RealignerTargetCreator command is itself giving errors, then you should just post those.
dpryan is offline   Reply With Quote
Old 10-09-2013, 01:09 AM   #7
harryzs
Member
 
Location: Germany

Join Date: Dec 2010
Posts: 29
Default

Quote:
Originally Posted by vishnuamaram View Post
Hi Seqanswers community,


I am doing human whole genome sequence analysis, of the data obtained Illumina Hiseq 2000.

This is the pipeline i followed & get errors in using GATK


a) BWA ------> SAM output
b) samtools ---> sortedindexedbamoutput
c) picard ------> for readgroup replacement & removing duplicates
d) GATK ----> for indelrealignment, BQSR, variant calling.


when i give the following command for GATK-indel realignment, i get the error.

U]cd used:-[/U]

java -Xmx30g -jar GenomeAnalysisTK-2.6-5-gba531bd/GenomeAnalysisTK.jar -T IndelRealigner \ -R INDEX_REF_GENOME_HG19/hg19.fa \ -I WG_PBMC_SAM9_CORDSORT_RMDUP.BAM -targetintervals intervalListFromRTC.intervals \ -o WG_PBMC_SAM9_coordsort_realigned.bam


the errors are :-



##### ERROR ------------------------------------------------------------------------------------------
##### ERROR A USER ERROR has occurred (version 2.6-5-gba531bd):
##### ERROR The invalid arguments or inputs must be corrected before the GATK can proceed
##### ERROR Please do not post this error to the GATK forum
##### ERROR
##### ERROR See the documentation (rerun with -h) for this tool to view allowable command-line arguments.
##### ERROR Visit our website and forum for extensive documentation and answers to
##### ERROR commonly asked questions http://www.broadinstitute.org/gatk
##### ERROR
##### ERROR MESSAGE: Argument with name '--targetIntervals' (-targetIntervals) is missing.
##### ERROR ------------------------------------------------------------------------------------------



kindly give me solutions, will be very thankful to you.

Please do let me know the proper usage of GATK.

Thanking you all,

Vishnu.

You may follow this
http://www.broadinstitute.org/gatk/g...name=tutorials
harryzs is offline   Reply With Quote
Old 10-09-2013, 01:10 AM   #8
harryzs
Member
 
Location: Germany

Join Date: Dec 2010
Posts: 29
Default

Quote:
Originally Posted by harryzs View Post
-->(howto) Perform local realignment around indels
harryzs is offline   Reply With Quote
Old 10-09-2013, 01:13 AM   #9
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

The informative line of that error message is:
Quote:
##### ERROR MESSAGE: Argument with name '--targetIntervals' (-targetIntervals) is missing.
Your command used "--targetintervals" (note the capital "I" in the error message).
dpryan is offline   Reply With Quote
Reply

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 12:39 PM.


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