-
Notifications
You must be signed in to change notification settings - Fork 18
Description
Dear @durrantmm ,
Thank you for a great resource!
I am interested in running the MGEfinder on metagenome contigs/short sequences and using the mgefinder wholegenome command for the same.
I am trying to find an IS element (IS1169) that ISfinder finds in my sequence of interest, but I am not able to find the IS element using MGEfinder
Code-
mgefinder wholegenome X71444.1.fna GCA_021405735.1_ASM2140573v1_genomic.fna --output_prefix test_GCA_021405735.1_ASM2140573v1_X71444.1
The files-
GCA_021405735.1 https://www.ncbi.nlm.nih.gov/datasets/genome/GCA_021405735.1/
X71444.1 https://www.ncbi.nlm.nih.gov/nuccore/X71444.1/
Output-
Empty files
Any suggestions on how to find a specific IS element of interest?
Another question-
ISFinder method is very different from MGEfinder. MGEfinder is more specific and looks for not just direct repeats but also inverted repeats. Whereas ISFinder blast search may suggest the presence of an IS element, but the inverted repeats may or maynot be present. Is there any benchmarking with the ISFinder? Sorry, if it seems naive, still trying to understand my data.
Looking forward to your reply!