First, I calculated some basic information about my bam file which is derived from alignment after duplicates Marking. They are: Mean Insert Size: 478.5; Minimum Insert Size is 100; Maximum Insert Size is 121827452. Standard Deviation is 21.6.

Then, I use the "Mean Insert Size- 4Standard Deviation" as the minimum parameter and "Mean Insert Size+4Standard Deviation" as the maximum one to convert the bam file into DIVET format which is input format for VariationHunter.

At last, I run the ./VariationHunter_SC:

Please enter the minimum paired-end insert size: 100
Please enter the maximum paired-end insert size: 121827452
Please enter the pre-processing mapping prune probability: 1
Please enter the minimum support for a cluster: 2

I am learning on NGS. I don't know where the problem is. I also tried other parameters in the command line, like setting a smaller maximum insert size. But it still failed. I think I need to change some settings in order to use the tool with non-human genome. Right? Any help will be appreciated!

Thank you!