swan_scan, Killed!

Issue #6 resolved
G Varma created an issue

Hi,

I was trying to use SWAN for my Whole Genome dataset, I used following command line. I was able to run "swan_stat" successfully, however while running "swan_scan" it got killed, any suggestions that I can avoid this error?

#-----------command line used----------------------------------

$SWAN_BIN/swan_scan -c a -n /home/user/Desktop/swan-s/data/human_g1k_v37.fasta.gap.bed /media/softwares/human_g1k_v37.fasta /media/HCM_WGS_BAM/B00H7EW_1-2-3_sort_mark.bam -Info: swan_scan vesion: 0.2.0 -Info: invoking command: /usr/lib/R/bin/exec/R --slave --no-restore --file=/home/user/R/x86_64-pc-linux-gnu-library/3.2/swan/bin/swan_scan --args -c a -n /home/user/Desktop/swan-s/data/human_g1k_v37.fasta.gap.bed /media/softwares/human_g1k_v37.fasta /media/HCM_WGS_BAM/B00H7EW_1-2-3_sort_mark.bam all chrs: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 X Y -Info: doing scan, gtk=123.578s,gmk=(rss,vsz) 3276264 3545308kb -Info: apply_chr: 1 , rg_files: 1 [1] "/media/HCM_WGS_BAM/B00H7EW_1-2-3_sort_mark.bam"

-Info: apply_chr: 1 , sp_prefix: [1] "/media/HCM_WGS_BAM/B00H7EW_1-2-3_sort_mark.1" -Info: apply_chr: 1 , rg_prefix: 1 [1] "/media/HCM_WGS_BAM/B00H7EW_1-2-3_sort_mark.1.rg1"

=Info apply_scan: group 1 of 1 groups within sample 1 of 1 samples =Info apply_scan: rg_file /media/HCM_WGS_BAM/B00H7EW_1-2-3_sort_mark.bam =Info apply_scan: rg_prefix /media/HCM_WGS_BAM/B00H7EW_1-2-3_sort_mark.1.rg1 =Info apply_scan: rg_out /media/HCM_WGS_BAM/B00H7EW_1-2-3_sort_mark.1.rg1 =Info: find rg isize modes at 19 332 with densities 4.827486e-05 0.005674746 ==Info: scan_bam: t_start= 1 t_end= 1000000 ; 1 -th of 250 Killed

#------------------------------------------------

Regards, Varma

Comments (2)

  1. Charlie Xia repo owner

    Dear Dr. Varma,

    I think it might be memory problem. Can you try with '-c 22' see if it works.

    THanks.

    Charlie

  2. Log in to comment