Open
Description
Hi, my metagenome input files were metagenomes that were binned and used as input.
When running finding_overlapping_regions.sh, I get an error at the 'blastdbcmd' step.
It is saying that the start is not defined.
looking through the datasets,
ln 95:
"len=$length" '{if ($13 == "minus") print $2"\t"$10-len"
This results in '-1711' and causes an error in 'blastdbcmd' step.
Example of the blast output (tab input file):
NC_014248_1_genome_0:1000 3300047529_3_Ga0494058_00019 98.300 1000 13 1 1 1000 1284 289 0.0 1749 minus
NC_014248_1_genome_0:1000 3300047525_1_Ga0494053_00203 98.200 1000 14 1 1 1000 34610 35605 0.0 1744 plus
NC_014248_1_genome_0:1000 3300047525_1_Ga0494053_00209 97.000 1000 30 0 1 1000 30429 31428 0.0 1681 plus
NC_014248_1_genome_0:1000 3300047526_4_Ga0494054_00523 97.036 911 20 1 1 904 911 1 0.0 1526 minus
Metadata
Metadata
Assignees
Labels
No labels