Originally posted by kmcarr
View Post
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Originally posted by chenyao View PostAre you sure? If it is , then there is no difference between "--initial-read-mismatches" and "--bowtie-n". Why they set two options?
You do NOT pass a number with the --bowtie-n option. By default TopHat runs Bowtie allowing 2 mismatches, in either the full read (-v) or seed (-n) mode. If you wish to changes this for either mode you set it with the "--initial-read-mismatches" option but as stated above this parameter can only accept numbers from 0-3.
And I can't allow the mismatch more than 3? This is unfair for the long read.
If you believe that your reads may contain a larger number of mismatches this when you should employ the --bowtie-n option. Understand that when bowtie is running in -n mode it is only concerned with the number of mismatches in the initial 5' seed.
You should really carefully study the manuals for TopHat and Bowtie to fully understand the run parameters, their limitations and implications for your experiment.
Leave a comment:
-
Originally posted by kmcarr View PostWell again, even for the -v option (default without specifying -n) the maximum allowable read mismatches is 3.
Leave a comment:
-
Originally posted by chenyao View PostNo, I did not use bowtie-n when I got this error. Also, I got the same error even I set the --initial-read-mismatches to 4.
Leave a comment:
-
Originally posted by kmcarr View Postchenyao,
You asked in an earlier thread about using the "--bowtie-n" option for TopHat, did you use that option when you got this error? The maximum number of mismatches allowed for -n is 3 so your --initial-read-mismatches 10 would produce an error when the combination of "-n 10" is passed by TopHat to Bowtie. It is possible that TopHat is not trapping this error and then when it tries to open one of the files which should have been produced by Bowtie (but wasn't because Bowtie failed) you get that error.
Leave a comment:
-
chenyao,
You asked in an earlier thread about using the "--bowtie-n" option for TopHat, did you use that option when you got this error? The maximum number of mismatches allowed for -n is 3 so your --initial-read-mismatches 10 would produce an error when the combination of "-n 10" is passed by TopHat to Bowtie. It is possible that TopHat is not trapping this error and then when it tries to open one of the files which should have been produced by Bowtie (but wasn't because Bowtie failed) you get that error.
Leave a comment:
-
Tophat with --initial-read-mismatches command Error:
Dear all,
I use tophat with --initial-read-mismatches 10
However I have the error:
gzip: stdout: Broken pipe
error:could not open pipe gzip:left_kept.fa.z
----
I thought it was the gzip error, so I use -z0 to have unzip file, but I have another error:
self.file=open(filename)
IOError: [Errno 2] No such file or directory: '/left_kept_reads_missing.fq'
------------
I use solid pair-end reads. Is it the reason ?Tags: None
Latest Articles
Collapse
-
by seqadmin
Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...-
Channel: Articles
09-23-2024, 06:35 AM -
-
by seqadmin
During the COVID-19 pandemic, scientists observed that while some individuals experienced severe illness when infected with SARS-CoV-2, others were barely affected. These disparities left researchers and clinicians wondering what causes the wide variations in response to viral infections and what role genetics plays.
Jean-Laurent Casanova, M.D., Ph.D., Professor at Rockefeller University, is a leading expert in this crossover between genetics and infectious...-
Channel: Articles
09-09-2024, 10:59 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 10-02-2024, 04:51 AM
|
0 responses
13 views
0 likes
|
Last Post
by seqadmin
10-02-2024, 04:51 AM
|
||
Started by seqadmin, 10-01-2024, 07:10 AM
|
0 responses
21 views
0 likes
|
Last Post
by seqadmin
10-01-2024, 07:10 AM
|
||
Started by seqadmin, 09-30-2024, 08:33 AM
|
0 responses
25 views
0 likes
|
Last Post
by seqadmin
09-30-2024, 08:33 AM
|
||
Started by seqadmin, 09-26-2024, 12:57 PM
|
0 responses
18 views
0 likes
|
Last Post
by seqadmin
09-26-2024, 12:57 PM
|
Leave a comment: