Question: HTseq: Error occured when reading beginning of SAM/BAM file.
0
gravatar for maple964
2.1 years ago by
maple9640
maple9640 wrote:

Dear all,

I tried to use HTseq for a week because I want to use "--nonunique" function. However

Error occured when reading beginning of SAM/BAM file.
  'pysam.csamtools.AlignedRead' object has no attribute 'reference_start'
  [Exception type: AttributeError, raised in _HTSeq.pyx:1371]

I tried to update the newest version of HTseq, and someone suggested to convert bam>sam>bam. I also tried, but it didn't work. (for this I also updated samtools 1.9 version)

I checked this page htseq-count: error when reading beginning of SAM/BAM file., but I couldn't find any solution.

This is my comment line

htseq-count -f bam -s no -m union --nonunique all -q accepted_hits_sam2bam.sorted /usr1/thhuang/index/TAIR10_GFF3_genes_transposons.gtf > Hiseq_count3.count

I am too new in bioinformatics. I really need some help.

Thank you all in advance

rna-seq software error • 796 views
ADD COMMENTlink modified 2.1 years ago by GenoMax95k • written 2.1 years ago by maple9640
1

I wonder if this is applicable in your case: https://github.com/hall-lab/svtyper/issues/37


Side note: you can use -M option with featureCounts to count multi-mapping reads (caution: this is not a good idea for RNAseq data).

ADD REPLYlink modified 2.1 years ago • written 2.1 years ago by GenoMax95k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 956 users visited in the last hour
_