Hy guys,

I am new. I'm trying to annotate my bacterial genome. I ran Prodigal to find CDS, nucleotide seqs and amminoacids seqs. I blasted the nucleotide sequences against Uniprot by using blastx and the option number 6 as outfmt Then tryed to convert the tabular output from Blast vs Uniprot into gff by using this two scripts, blast_to_gff_wrapper.sh and blast_to_gff.py.

After that I sorted and indexed my gff file by using Tabix, as follow:

(grep ^"#" in.gff; grep -v ^"#" in.gff | sort -k1,1 -k4,4n) | bgzip > in-sorted.gff3.gz

tabix -p gff in-sorted.gff3.gz

When I loaded the gff.file into IGV I did not see anything. Do you have any idea?

Thank you in advance.

Best,

Salvatore