I have assembled e. coli genomes in which I want to search for a list of genes.
I have access to a script which takes a multifasta file of genes as the database and makes a blast comparison. The output is either crunch files for each genome that can be viewed or a list of hits. However the list of hits does not contain details on which gene from the multifasta file is hit.
So I want a quick way to search for multiple genes at the same time in one or more genomes and I want a list of what genes are found, length of query and the e-value, does anyone know how to do this?
Cheers!!
I have access to a script which takes a multifasta file of genes as the database and makes a blast comparison. The output is either crunch files for each genome that can be viewed or a list of hits. However the list of hits does not contain details on which gene from the multifasta file is hit.
So I want a quick way to search for multiple genes at the same time in one or more genomes and I want a list of what genes are found, length of query and the e-value, does anyone know how to do this?
Cheers!!