Skip to content

Error running novoBreak #10

@moldach

Description

@moldach

I'm getting an error trying to run novoBreak. It's in my path (placed export PATH=$PWD/nb_distribution/:$PATH in ~/.bashrc and sourced it. I've also made each of the scripts executable with chmod +x *

[moldach@cdr767 alignment]$ bash /scratch/moldach/bin/nb_distribution/run_novoBreak.sh /scratch/moldach/bin/nb_distribution/ ~/projects/def-mtarailo/common/indexes/WS265_wormbase/c_elegans.PRJNA13758.WS265.genomic.fa 470.sorted.dedupped.bam control.bam 1 novoBreak
[Tue May  5 17:04:43 2020]
Building reference kmer...
Calculating number of paired reads ...
[fillin_bitvec] processed   72635794 reads
There are 23315106 reference kmers and 59276238 non-reference kmers
Finished reference kmer building
[Tue May  5 17:06:58 2020]
Freed reference hash and begin building reads hash table...
[Tue May  5 17:06:58 2020]
/scratch/moldach/bin/nb_distribution/run_novoBreak.sh: line 26: 15470 Killed                  $novobreak -i $tumor_bam -c $normal_bam -r $ref -o kmer.stat
[M::bam2fq_mainloop] processed 0 reads
begin kmer2id ...
kmer2id takes 0 seconds
begin id2pair...
id2pair takes 0 seconds
begin sorting ids...
sorting ids takes 0 seconds
begin output results...
Outputting results takes 0 seconds
Finished
(standard_in) 1: parse error
*.txt
No such file or directory at /scratch/moldach/bin/nb_distribution/run_ssake.pl line 11.
awk: fatal: cannot open file `../group_reads/split/*.ssake.asm.out' for reading (No such file or directory)
[main] Version: 0.7.10-r806-dirty
[main] CMD: /scratch/moldach/bin/nb_distribution/bwa mem -t 1 -M /project/6013424/common/indexes/WS265_wormbase/c_elegans.PRJNA13758.WS265.genomic.fa ssake.fa
[main] Real time: 1.141 sec; CPU: 0.237 sec
x??
No such file or directory at /scratch/moldach/bin/nb_distribution/infer_bp_v4.pl line 13.

It give me this error but the file infer_bp_v4.pl is clearly there?

[moldach@cdr767 alignment]$ ll /scratch/moldach/bin/nb_distribution/
total 5008
-rwxr-x--- 1 moldach moldach  955882 May  5 15:17 bwa
-rwxr-x--- 1 moldach moldach    2774 May  5 15:17 fetch_discordant.pl
-rwxr-x--- 1 moldach moldach     980 May  5 15:17 filter_sv2.pl
-rwxr-x--- 1 moldach moldach    1436 May  5 15:17 filter_sv.bak.pl
-rwxr-x--- 1 moldach moldach    1171 May  5 15:17 filter_sv_icgc.pl
-rwxr-x--- 1 moldach moldach    1405 May  5 15:17 filter_sv.pl
-rwxr-x--- 1 moldach moldach    4004 May  5 15:17 group_bp_reads.pl
-rwxr-x--- 1 moldach moldach    4512 May  5 15:17 infer_bp.pl
-rwxr-x--- 1 moldach moldach    4897 May  5 15:17 infer_bp_v4.pl
-rwxr-x--- 1 moldach moldach   13282 May  5 15:17 infer_sv.pl
-rwxr-x--- 1 moldach moldach    1069 May  5 15:17 LICENSE
-rwxr-x--- 1 moldach moldach  702112 May  5 15:17 novoBreak
-rwxr-x--- 1 moldach moldach    5986 May  5 15:17 README.md
-rwxr-x--- 1 moldach moldach    2853 May  5 15:17 run_novoBreak.sh
-rwxr-x--- 1 moldach moldach    2276 May  5 15:17 run_ssake.pl
-rwxr-x--- 1 moldach moldach 3299767 May  5 15:17 samtools
-rwxr-x--- 1 moldach moldach   83364 May  5 15:17 SSAKE

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions