Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • HTseq-count errors

    I put 10 .bam files (from Tophat) through the same samtools & htseq-count commands, and 3/10 didn't work (each with their own error) in htseq.

    $ samtools view -h -o sample.sam sample.bam

    $ htseq-count --stranded no -i gene_id sample.sam genes.gtf > sample.gene.counts

    The errors are:

    1) ...596437 GFF lines processed.
    Error occured when reading beginning of SAM/BAM file.
    [Exception type: StopIteration, raised in count.py:84]

    This one particularly confuses me since the program didn't choke on any of the other sam files that were prepared the same way.

    2) ...20000000 SAM alignment records processed.
    Error occured when processing SAM input (line 20080682 of file TH2/DEprep/BR2_N2In.sam):
    ('SAM line does not contain at least 11 tab-delimited fields.', 'line 20080682 of file TH2/DEprep/BR2_N2In.sam')
    [Exception type: ValueError, raised in _HTSeq.pyx:1276]

    3) ...35089676 SAM alignments processed.
    [Errno 5] Input/output error
    [Exception type: IOError, raised in count.py:212]


    Any suggestions for how I should fix these problems?

  • #2
    I wonder if you're having hard drive corruption issues.

    In any case, make sure you're using the most recent version of htseq-count. Also, check line 20080682 of the sam file mentioned and ensure it's not malformed.

    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
    8 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 06:07 PM
    0 responses
    8 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-22-2024, 10:03 AM
    0 responses
    49 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-21-2024, 07:32 AM
    0 responses
    66 views
    0 likes
    Last Post seqadmin  
    Working...
    X