Question: VariantAnnotator GATK Error
gravatar for cvu
5.2 years ago by
cvu140 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 for help creating it.

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





snp alignment assembly sequence • 2.3k views
ADD COMMENTlink modified 5.2 years ago • written 5.2 years ago by cvu140

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 5.2 years ago by Devon Ryan91k
Please log in to add an answer.


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