sequenza-utils bam2seqz command is not terminating.

67 views
Skip to first unread message

TANAY BISWAS

unread,
May 17, 2020, 1:52:01 PM5/17/20
to Sequenza User Group
Hi
I'm using sequenza to identify CNVs. I'm now in the step to generate seqz file. First of all, although i'm giving tumor and normal bam files as an input, it is showing [mpileup] after giving the command. The output file is generating but the command is not terminating. The command is running since 9hrs. And if there is not any problem as the output file is generating, can I speed up this function?
Here's my command:
@tanay:~$ sequenza-utils bam2seqz -n ~/normal.recal.bam -t ~/tumor.recal.bam --fasta /home/tanay/hg19.fa -gc ~/copynumber/hg19.gc50Base.wig.gz -o ~/copynumber/out6.seqz.gz
[mpileup] 1 samples in 1 input files
[mpileup] 1 samples in 1 input files

See, it is showing [mpileup] and running till now. I've also attached the output file, just showing few lines of the main seqz.gz file. Please help me out of this.

Regards,
Tanay
out6.seqz.ods
Reply all
Reply to author
Forward
0 new messages