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