Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
What if you disable the -I and -X option, just run a basic bowtie2 run and see what happens?
-
That's odd. It must be the shell script that you are using that is causing the error.
Did you get the pre-compiled binary or compiled from source yourself?
Leave a comment:
-
So if I run the code below, it works. I had to take out my options (I want to use "--minins 200" and "--maxins 600").
/(path)/bowtie2 -x /(path)/genomeIndex -1 /(path)/(reads1).fastq -2 /(path)/(reads2).fastq -S (alignment).sam
If I run it with my desired options (code below):
/(path)/bowtie2 --minins 200 --maxins 600 -x /(path)/genomeIndex -1 /(path)/(reads1).fastq -2 /(path)/(reads2).fastq -S (alignment).sam
I get the error:
Extra parameter(s) specified: "—-maxins", "600"
Note that if <mates> files are specified using -1/-2, a <singles> file cannot
also be specified. Please run bowtie separately for mates and singles.
Error: Encountered internal Bowtie 2 exception (#1)
Command: /(path)/bowtie2-align-s --wrapper basic-0 --minins 200 -x /(path)/genomeIndex -S (alignment).sam -1 (reads1).fastq -2 (reads2).fastq --maxins 600
(ERR): bowtie2-align exited with value 1
@GenoMax, I am using OS X version 10.9.5.
I also tried updating my bowtie2 from 2.2.3 to 2.2.5 and that has not helped.
Leave a comment:
-
Tested on a Mac with latest bowtie (v.2.2.5) and the lambda test data that comes with bowtie. No errors.
@millrick10: Are you using the latest Mac binary?
Leave a comment:
-
@GenoMax and @millrick10
It does look like a Mac system which may cause problems.
BTW, my bowtie --version (on Linux) looks like:
/group/bioinfo/apps/apps/bowtie2-2.2.5/bowtie2-align-s version 2.2.5
64-bit
Built on localhost.localdomain
Thu Mar 5 18:34:52 EST 2015
Compiler: gcc version 4.1.2 20080704 (Red Hat 4.1.2-54)
Options: -O3 -m64 -msse2 -funroll-loops -g3 -DPOPCNT_CAPABILITY
Sizeof {int, long, long long, void*, size_t, off_t}: {4, 8, 8, 8, 8, 8}
Leave a comment:
-
This is strange. Not the 'Note that ...' part. That message (as per the source code) will be output if extra parameters are found. The strange part is why bowtie2 thinks you have extra parameters. What happens if you take out the minins and maxins parameters?
BTW, the source code says "... TODO: Argument parsing is very, very flawed. ..." which, having looked at it, is a statement I agree with.
Leave a comment:
-
@westerman
Results of "bowtie2 --version":
"64-bit
Built on pc-201305311418.ebalto.jhmi.edu
Thu May 29 12:48:36 EDT 2014
Compiler: gcc version 4.2.1 (Based on Apple Inc. build 5658) (LLVM build 2336.11.00)
Options: -O3 -m64 -msse2 -funroll-loops -g3 -DPOPCNT_CAPABILITY
Sizeof {int, long, long long, void*, size_t, off_t}: {4, 8, 8, 8, 8, 8}"
Leave a comment:
-
@westerman
Sure thing, it is below.
My code:
/(path)/bowtie2 --minins 200 --maxins 600 -x /(path)/genomeIndex -1 (reads1).fastq -2 (reads2).fastq -S (alignment).sam
The error:
Extra parameter(s) specified: "--minins", "200", "--maxins", "600"
Note that if <mates> files are specified using -1/-2, a <singles> file cannot
also be specified. Please run bowtie separately for mates and singles.
Error: Encountered internal Bowtie 2 exception (#1)
Command: /(path)/bowtie2-align-s --wrapper basic-0 -x /(path)/bowtie2-2.2.3/indexes/(genomeIndex) -S (alignment).sam -1 (reads1).fastq -2 (reads2).fastq --minins 200 --maxins 600
(ERR): bowtie2-align exited with value 1
Also, I have to apologize, in the code that I originally posted, I did have double dashes with I and X. I fixed it though, and it did not get rid of the error.
Leave a comment:
-
@westerman
Nope, they are single dashes in "-I" and "-X." For good measure, I also swapped "--minins" for "-I" and "--maxins" for "-X" and it still gave me the same error.
Thanks for the idea though.
Leave a comment:
-
I think you have double dashes for '-I' and '-X'. They should be single dashes or the long name with double dash; e.g., '--minins' and '--maxins'. Take care of that problem first and then let us know if there are problems with '-1' and '-2'.
Leave a comment:
-
Using bowtie2 with min and max fragment fragment
I am pretty new to bioinformatics, and I was hoping that I could get some help with a bowtie2 script that I am putting together.
I am aligning ChIP-seq paired-end reads (2 samples: Input control and ChIP sample) to a reference and I want to include a minimum fragment size of 200bp and maximum fragment size of 600bp. A pseudocode version of my script is below.
/(path)/bowtie2 —I 200 —X 600 -x /(path)/genomeIndex -1 (chip/IN).fastq -2 (chip/IN).fastq -S alignment.sam
When I run this code, I get the error:
"Extra parameter(s) specified: "—I", "200", "—X", "600"
Note that if <mates> files are specified using -1/-2, a <singles> file cannot
also be specified. Please run bowtie separately for mates and singles.
Error: Encountered internal Bowtie 2 exception (#1)"
and:
"(ERR): bowtie2-align exited with value 1"
I am not specifying any singles files, only mate pair files, so I am thoroughly confused. Any suggestions?
Latest Articles
Collapse
-
by seqadmin
Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...-
Channel: Articles
10-18-2024, 07:11 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 11-01-2024, 06:09 AM
|
0 responses
22 views
0 likes
|
Last Post
by seqadmin
11-01-2024, 06:09 AM
|
||
New Model Aims to Explain Polygenic Diseases by Connecting Genomic Mutations and Regulatory Networks
by seqadmin
Started by seqadmin, 10-30-2024, 05:31 AM
|
0 responses
21 views
0 likes
|
Last Post
by seqadmin
10-30-2024, 05:31 AM
|
||
Started by seqadmin, 10-24-2024, 06:58 AM
|
0 responses
25 views
0 likes
|
Last Post
by seqadmin
10-24-2024, 06:58 AM
|
||
New AI Model Designs Synthetic DNA Switches for Targeted Gene Expression in Specific Cell Types
by seqadmin
Started by seqadmin, 10-23-2024, 08:43 AM
|
0 responses
54 views
0 likes
|
Last Post
by seqadmin
10-23-2024, 08:43 AM
|
Leave a comment: