Question: VariantAnnotator GATK Error
0
gravatar for cvu
4.8 years ago by
cvu130
India
cvu130 wrote:

I' m trying to get rs ID# in my .vcf file which i have generated through samtools. i am using GATK VariantAnnotator to get rs ID# in my vcf file.

java -jar GenomeAnalysisTK-3.1-1/GenomeAnalysisTK.jar -R ambda_virus.fa -T VariantAnnotator -I local_sorted.bam -o output.vcf --variant sim_variants.vcf -L sim_variants.vcf --dbsnp 00-All.vcf

but getting this error:

ERROR MESSAGE: Fasta index file /home/bio1/Desktop/Shivangi/BOWTIE/bowtie2-2.2.3/example/reference/lambda_virus.fa.fai.fai for reference /home/bio1/Desktop/Shivangi/BOWTIE/bowtie2-2.2.3/example/reference/lambda_virus.fa.fai does not exist. Please see http://gatkforums.broadinstitute.org/discussion/1601/how-can-i-prepare-a-fasta-file-to-use-as-reference for help creating it.

please let me know, these command line arguments are okay?

 

 

 

 

snp alignment assembly sequence • 2.2k views
ADD COMMENTlink modified 4.8 years ago • written 4.8 years ago by cvu130
1

Firstly, the command you posted wouldn't create that error message (you have a typo). Secondly, I suspect that you used -R lambda_virus.fa.fai rather than -R lambda_virus.fa, since that would indeed cause that error. Try running the command again like that and see if the error still occurs (if so, just copy and paste the exact command you used into a new comment).
 

ADD REPLYlink written 4.8 years ago by Devon Ryan88k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1138 users visited in the last hour