Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • #31
    Originally posted by TEFA View Post
    Hi EGrassi,

    Did you solve this issue??????
    Actually no, right now I'm trying to avoid using '-' as an htseq-count argument without using it in a pipeline with samtools view, I'll let you know what happens...

    Comment


    • #32
      Originally posted by EGrassi View Post
      Woa, I used samtools -n and htseq-count gave me all counts of zero...that's strange, I will check with IGV or similar tools but cuffdiff gave me fpkg values with the same gtf and bam file...
      Was this by any chance a TopHat output? Somebody alerted me recently to a bug that causes mishandling of SAM lines that contain the "NH" flag, and that seems to appear in newer tophat output. I'll try to fix this ASAP.

      Simon

      Comment


      • #33
        Originally posted by Simon Anders View Post
        Was this by any chance a TopHat output? Somebody alerted me recently to a bug that causes mishandling of SAM lines that contain the "NH" flag, and that seems to appear in newer tophat output. I'll try to fix this ASAP.
        Tophat 2.0 output, yes!

        Thank you,
        E.

        Comment


        • #34
          Yep, I can confirm you that removing the NH flags from a subset of my sam files I get some counts different from 0. This could be a temporary workaroud.

          Comment


          • #35
            Ok, on the whole data it falls back to all 0:
            samtools view S1.sorted.bam.bam | sed 's/NM:i:\d//g' | htseq-count --stranded=no - /rogue/bioinfotree/prj/ewing-rnaseq/local/share/data/Homo_sapiens/UCSC/hg19/Annotation/Genes/genes.gtf > counts_S1_htseqq

            I'm really a little surprised by the amount of errors and other glitches in rna seq software (this is a generic rant, not focused on htseq!)...it's a really complicated and new area, that's true, but right now using any of the existing tool seems a bet and to tell the truth I never feel confident on the results, as long as with every version results changes, and when something runs without errors I never know if that means that it worked or that an exception just had not screamed enough :/

            Comment


            • #36
              The bug with the NH flags, which was introduced during some refactoring in the version of a week ago or so is now fixed in version 0.5.3p8. Sorry that the recent changes turned out to be so volatile.

              Comment


              • #37
                Thanks, 0.5.3p9 has produced counts different from zero. Let's see what deSeq thinks of this data!

                Comment


                • #38
                  Originally posted by Simon Anders View Post
                  Sorry, it seems we made some mix-up between version 0.5.3p4 and 0.5.3p5. Essentially, p5 undid some fixes in p4, including the one for "*" qualities. Now, there is version 0.5.3p6, which should clean up this mess. Please let me know if you still have problems.
                  Hello, I have performed the mapping of miRNA reads to the mirbase precursor file by SHRiMP. But I'm facing a probem in HTSeq while taking the read counts from the output.sam file.

                  The error which I'm trying to solve is;

                  Error occured when processing SAM input (line 2305 of file HBI_10.sam): ("'seq' and 'qualstr' do not have the same length.", 'line 2305 of file HBI_10.sam') [Exception type: ValueError, raised in _HTSeq.pyx:808]

                  Has anyone of you come across this error? What might be the solution?

                  Comment


                  • #39
                    Sushant - what does line 2305 of your SAM file look like? Unless the quality entry is the special value '*' only then it is a likely an unrelated issue.

                    Comment


                    • #40
                      Originally posted by maubp View Post
                      Sushant - what does line 2305 of your SAM file look like? Unless the quality entry is the special value '*' only then it is a likely an unrelated issue.
                      Actually, I did see that SAM file. The problem is that my sequence string & its respective sequence quality string is not of equal length.
                      However, I converted my fastq files to fasta and then performed the alignment, which solved this problem.

                      But the mapping percentage I'm getting with SHRiMP is too low(2.77%), for the mirna aligning to mature.
                      What parameters should I consider?

                      Comment

                      Latest Articles

                      Collapse

                      • seqadmin
                        Strategies for Sequencing Challenging Samples
                        by seqadmin


                        Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                        03-22-2024, 06:39 AM
                      • seqadmin
                        Techniques and Challenges in Conservation Genomics
                        by seqadmin



                        The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                        Avian Conservation
                        Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                        03-08-2024, 10:41 AM

                      ad_right_rmr

                      Collapse

                      News

                      Collapse

                      Topics Statistics Last Post
                      Started by seqadmin, Yesterday, 06:37 PM
                      0 responses
                      10 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, Yesterday, 06:07 PM
                      0 responses
                      9 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 03-22-2024, 10:03 AM
                      0 responses
                      51 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 03-21-2024, 07:32 AM
                      0 responses
                      67 views
                      0 likes
                      Last Post seqadmin  
                      Working...
                      X