Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • 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 ?

  • #2
    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.

    Comment


    • #3
      Originally posted by kmcarr View Post
      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.
      No, 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.

      Comment


      • #4
        Originally posted by chenyao View Post
        No, 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.
        Well again, even for the -v option (default without specifying -n) the maximum allowable read mismatches is 3.

        Comment


        • #5
          Originally posted by kmcarr View Post
          Well again, even for the -v option (default without specifying -n) the maximum allowable read mismatches is 3.
          Are you sure? If it is , then there is no difference between "--initial-read-mismatches" and "--bowtie-n". Why they set two options? And I can't allow the mismatch more than 3? This is unfair for the long read.

          Comment


          • #6
            Originally posted by chenyao View Post
            Are you sure? If it is , then there is no difference between "--initial-read-mismatches" and "--bowtie-n". Why they set two options?
            No, they are not the same thing. These two TopHat parameters control the parameters which are passed to Bowtie; "--bowtie-n" runs bowtie with the -n option which specifies counting mismatches only in the initial seed match. The alternative mode is -v which counts mismatches over the entire read length. By default TopHat calls Bowtie with the -v option.

            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.
            Life isn't fair .

            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.

            Comment


            • #7
              Originally posted by kmcarr View Post
              No, they are not the same thing. These two TopHat parameters control the parameters which are passed to Bowtie; "--bowtie-n" runs bowtie with the -n option which specifies counting mismatches only in the initial seed match. The alternative mode is -v which counts mismatches over the entire read length. By default TopHat calls Bowtie with the -v option.

              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.



              Life isn't fair .

              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.
              Thank you so much. So --bowtie-n do not need the number, that's my mistake.

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Latest Developments in Precision Medicine
                by seqadmin



                Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

                Somatic Genomics
                “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
                05-24-2024, 01:16 PM
              • seqadmin
                Recent Advances in Sequencing Analysis Tools
                by seqadmin


                The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
                05-06-2024, 07:48 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 05-24-2024, 07:15 AM
              0 responses
              198 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 05-23-2024, 10:28 AM
              0 responses
              220 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 05-23-2024, 07:35 AM
              0 responses
              229 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 05-22-2024, 02:06 PM
              0 responses
              13 views
              0 likes
              Last Post seqadmin  
              Working...
              X