Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Error during bowtie2-align run , (ERR): bowtie2-align exited with value 139 harijay Bioinformatics 2 02-08-2017 01:24 AM
Bowtie2 alignment error MichielMunckhof Bioinformatics 7 02-06-2015 03:29 AM
Help- Bowtie2-alignment error Votinhkiem90 Bioinformatics 1 03-20-2014 01:17 AM
bowtie2 gapped alignment Bulak Bioinformatics 0 03-11-2013 07:39 AM
Bowtie2 'Could not open alignment output file' error Lcontami Bioinformatics 1 02-28-2013 03:14 AM

Thread Tools
Old 10-10-2018, 02:16 AM   #1
Junior Member
Location: France

Join Date: Oct 2018
Posts: 1
Default BS-seeker2 alignment error (bowtie2)

Hi everyone ,

I'm trying to align some of ours RRBS data on human genome hg38 with BS-seeker2.
Indexing went good without any error except the warning :
Warning: Encountered reference sequence with only gaps
but I know why.

I get errors in the alignement step about some file or directory not found, as you can see :

-bash-4.1$ ./ -i /YYYYYYY/NNRD204.fastq.gz --aligner=bowtie2 -o PostCT_6m_bs-seeker2.bam -g ../../sequence/hg38.fa -r -L 20 -U 250 -t Y

     BS-Seeker2 v2.1.5 - Dec. 21, 2017

[2018-10-10 11:24:42] Mode: Bowtie2, local alignment
[2018-10-10 11:24:42] Filter for tag XS: #(mCH)/#(all CH)>50.00% and #(mCH)>5
[2018-10-10 11:24:42] Temporary directory: /tmp/bs_seeker2_PostCT_6m_bs-seeker2.bam_-bowtie2-local-TMP-FAvs7M
[2018-10-10 11:24:42] Reduced Representation Bisulfite Sequencing: True
[2018-10-10 11:24:42] Single end
[2018-10-10 11:24:42] Aligner command: None/bowtie2 --local --quiet -p 2 -D 50 --norc --sam-nohead  -k 2 -x %(reference_genome)s -f -U %(input_file)s -S %(output_file)s
[2018-10-10 11:24:42] ----------------------------------------------
[2018-10-10 11:24:42] Read filename: /YYYYYYYY/NNRD204.fastq.gz
[2018-10-10 11:24:42] The first base (for mapping): 1
[2018-10-10 11:24:42] The  last base (for mapping): 200
[2018-10-10 11:24:42] Path for short reads aligner: None/bowtie2 --local --quiet -p 2 -D 50 --norc --sam-nohead  -k 2 -x %(reference_genome)s -f -U %(input_file)s -S %(output_file)s

[2018-10-10 11:24:42] Reference genome library path: /XXXXXX/BS-seeker2/bs_utils/reference_genomes/hg38.fa_rrbs_20_250_bowtie2
[2018-10-10 11:24:42] Un-directional library
[2018-10-10 11:24:42] Number of mismatches allowed: 4
[2018-10-10 11:24:42] --------------------------------
[2018-10-10 11:24:42] Start reading and trimming the input sequences
[2018-10-10 11:24:55] Processing read file: /tmp/bs_seeker2_PostCT_6m_bs-seeker2.bam_-bowtie2-local-TMP-FAvs7M/NNRD204.fastq.gz-s-1
[2018-10-10 11:25:09] Processing input is done
[2018-10-10 11:25:09] Start mapping
[2018-10-10 11:25:09] Starting commands:
[2018-10-10 11:25:09] Launched: None/bowtie2 --local --quiet -p 2 -D 50 --norc --sam-nohead  -k 2 -x /XXXXXX/BS-seeker2/bs_utils/reference_genomes/hg38.fa_rrbs_20_250_bowtie2/W_C2T -f -U /tmp/bs_seeker2_PostCT_6m_bs-seeker2.bam_-bowtie2-local-TMP-FAvs7M/Trimmed_C2T.fa.tmp-9757942 -S /tmp/bs_seeker2_PostCT_6m_bs-seeker2.bam_-bowtie2-local-TMP-FAvs7M/W_C2T_m4.mapping.tmp-9757942
[2018-10-10 11:25:09] Launched: None/bowtie2 --local --quiet -p 2 -D 50 --norc --sam-nohead  -k 2 -x /XXXXXX/BS-seeker2/bs_utils/reference_genomes/hg38.fa_rrbs_20_250_bowtie2/C_C2T -f -U /tmp/bs_seeker2_PostCT_6m_bs-seeker2.bam_-bowtie2-local-TMP-FAvs7M/Trimmed_C2T.fa.tmp-9757942 -S /tmp/bs_seeker2_PostCT_6m_bs-seeker2.bam_-bowtie2-local-TMP-FAvs7M/C_C2T_m4.mapping.tmp-9757942
[2018-10-10 11:25:09] Launched: None/bowtie2 --local --quiet -p 2 -D 50 --norc --sam-nohead  -k 2 -x /XXXXXXX/BS-seeker2/bs_utils/reference_genomes/hg38.fa_rrbs_20_250_bowtie2/W_G2A -f -U /tmp/bs_seeker2_PostCT_6m_bs-seeker2.bam_-bowtie2-local-TMP-FAvs7M/Trimmed_G2A.fa.tmp-9757942 -S /tmp/bs_seeker2_PostCT_6m_bs-seeker2.bam_-bowtie2-local-TMP-FAvs7M/W_G2A_m4.mapping.tmp-9757942
[2018-10-10 11:25:09] Launched: None/bowtie2 --local --quiet -p 2 -D 50 --norc --sam-nohead  -k 2 -x /XXXXXXX/BS-seeker2/bs_utils/reference_genomes/hg38.fa_rrbs_20_250_bowtie2/C_G2A -f -U /tmp/bs_seeker2_PostCT_6m_bs-seeker2.bam_-bowtie2-local-TMP-FAvs7M/Trimmed_G2A.fa.tmp-9757942 -S /tmp/bs_seeker2_PostCT_6m_bs-seeker2.bam_-bowtie2-local-TMP-FAvs7M/C_G2A_m4.mapping.tmp-9757942

Traceback (most recent call last):
  File "./", line 367, in <module>
  File "/XXXXXXX/BS-seeker2/bs_align/", line 682, in bs_rrbs
    'output_file' : CG2A} ])
  File "/XXXXXXX/BS-seeker2/bs_utils/", line 332, in run_in_parallel
    for i, proc in enumerate([subprocess.Popen(args = shlex.split(cmd), stdout = stdout) for cmd, stdout in commands]):
  File "/local/miniconda3/envs/python-2.7.15/lib/python2.7/", line 394, in __init__
    errread, errwrite)
  File "/local/miniconda3/envs/python-2.7.15/lib/python2.7/", line 1047, in _execute_child
    raise child_exception
OSError: [Errno 2] No such file or directory

So I tried to launch bowtie2 alignements with the commands just preciding the error (without the --quiet flag) and it went good, no error to declare :
bowtie2 --local -p 2 -D 50 --norc --sam-nohead  -k 2 -x /XXXXXX/BS-seeker2/bs_utils/reference_genomes/hg38.fa_rrbs_20_250_bowtie2/W_C2T -f -U /tmp/bs_seeker2_PostCT_6m_bs-seeker2.bam_-bowtie2-local-TMP-FAvs7M/Trimmed_C2T.fa.tmp-9757942 -S /tmp/bs_seeker2_PostCT_6m_bs-seeker2.bam_-bowtie2-local-TMP-FAvs7M/W_C2T_m4.mapping.tmp-9757942
I don't get why it raises me this error. I tried on local and on a distant cluster .
I would be very grateful if someone could help me !

Have a nice day !
Hyzteria is offline   Reply With Quote

bowtie2, bs-seeker2, error

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 05:35 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO