Question: Bacterial genome annotation using annovar
0
gravatar for Sondes Kalboussi
21 months ago by
Belgium
Sondes Kalboussi0 wrote:

Hello, am writing my own pipeline in python in order to annotate bacterial genome MTB, am new in this field and am a bit lost , I converted my VCF to appropriate annovar input format, then i got stack, i have to use dbSNP to annotate SNP and hrv37 as the reference genome for annotation but dont really know the correct command format or what i really need to provide more. I read the manual but it is not really helping me. Anyone with experience in using Annovar to annotate bacterial genomes? Thanks in advance

ADD COMMENTlink modified 21 months ago • written 21 months ago by Sondes Kalboussi0
1

I am assuming that you read this manual , to be able to get help here; you need first to try and write down what did you do and what error/s or where exactly you stumble.

ADD REPLYlink modified 21 months ago • written 21 months ago by Medhat8.2k

Thanks for the clarification, actually i still preparing the files that i would need to run table_annovar.pl, my output must have gene location and gene function annotation so am using the followin g command

Perl {}/table_annovar.pl --vcfinput {} {} --outfile {} --buildver H37RV --nastring . --protocol refGene,dbsnp --operation g,r,f --vcfinput""".format(self.annovar,input_vcf,MTBdb)

Am a bit lost cause am not sure if the databases am using are correct ones for gene location and gene function annotation

ADD REPLYlink written 20 months ago by Sondes Kalboussi0
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: 1253 users visited in the last hour