Skip to content

Error at step 5 (blastdbcmd) of finding_overlapping_regions.sh #1

Open
@pyongheng

Description

@pyongheng

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

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions