Problem in Step 2

[ceromova]

Hello,

I runned ERVcaller and I found this problem:

Step 2: Detecting TE insertions...

~~~~~ the input bam file was indexed
sh: extractSoftclipped: command not found

What shoul I do?
Thank you.

[xunchen85]

Hi, you may need to first compile the SE-MEI tool under the Scripts folder (tar vxzf SE-MEI-master.tar.gz & cd SE-MEI & make) and then use "export" command to add it to the $PATH.

Thanks,
Xun

[ceromova]

Ok. Solved. Thank you. Another problem showed up after that and it is the following one:
[E::bwa_idx_load_from_disk] fail to locate the index files

I aligned the sequenes with BWA and the files of the indexation are: .fna, .amb, .ann, .bwt, .fai, .pac, .sa

What could be the problem in here?

[xunchen85]

it could due to the empty output. ERVcaller did a reciprocal alignment, if the output is empty the next alignment will fail and show this error.

You could show the file sizes of your input, ref (especially the chr IDs), and generated outputs. Maybe I could help to double check.

Thanks,
Xun

[ceromova]

Hi,

About the error [E::bwa_idx_load_from_disk] fail to locate the index files

I downloaded the ref from here: GCA_000001405.15_GRCh38_no_alt_analysis_set.fna.gz

The file sizes of my input are 40Gb.

The generated outputs are an empty VCF file and a temp folder with the following information:

[xunchen85]

It looks like the error happened at an earlier time than the reciprocal alignment. The human reference file is okay.

Maybe you could show the full log information till the error for the debugging? It also depends on the ERV references that you are using, for example, no TE insertions of your ERV reference were detected.

Have you successfully run the test data yet which could confirm your successful installation of ERVcaller?

Xun