swan_scan, Killed!
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)
-
repo owner -
repo owner - edited description
- changed status to resolved
- Log in to comment
Dear Dr. Varma,
I think it might be memory problem. Can you try with '-c 22' see if it works.
THanks.
Charlie