About "xx_h1_sorted.bam" file about ervcaller HOT 4 CLOSED

Can you show the commands that you are using?

"_h1_sorted.bam" should be the whole-genome alignment file if you are using the fastq reads as the input. After you successfully got the VCF file, you can consider to either delete or keep it for other usages.

Xun

from ervcaller.

Hi Xun.

Thanks for quick reply.

I didn't use fastq reads as the input.

perl $ERVcaller/ERVcaller_v1.4.pl
-i ${Sample_name}
-f .bam
-H $Human_genome
-T ${ERVdatabase}/Human_TE_library.fa
-I ${Input_bam_dir}
-O ${Output_dir}
-r 150
-t 16
-S 40
-BWA_MEM
-G 2>> ${Output_dir}/${Sample_name}.$Tool.log

from ervcaller.

Hi,

I may answer this question in another comment. Firstly you can check if your input BAM file is from BAM MEM; Secondly you can check if you sucessfully get the VCF file (with called TE loci) or the job may fail; Thirdly you also can run a 1KGP low-coverage sample HG00100 as a technical control, which I just ran successfully in few days.

Best,
Xun

from ervcaller.

Thank you.

I'll check and try.

Best,
Oh

from ervcaller.