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 :
but I know why.
I get errors in the alignement step about some file or directory not found, as you can see :
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 :
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 !
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 :
Code:
Warning: Encountered reference sequence with only gaps
I get errors in the alignement step about some file or directory not found, as you can see :
Code:
-bash-4.1$ ./bs_seeker2-align.py -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 "./bs_seeker2-align.py", line 367, in <module> options.cut_format File "/XXXXXXX/BS-seeker2/bs_align/bs_rrbs.py", line 682, in bs_rrbs 'output_file' : CG2A} ]) File "/XXXXXXX/BS-seeker2/bs_utils/utils.py", 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/subprocess.py", line 394, in __init__ errread, errwrite) File "/local/miniconda3/envs/python-2.7.15/lib/python2.7/subprocess.py", line 1047, in _execute_child raise child_exception [B]OSError: [Errno 2] No such file or directory[/B]
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 :
Code:
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 would be very grateful if someone could help me !
Have a nice day !