Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • Olha
    Junior Member
    • Sep 2015
    • 6

    RseQC_read_distribution.py_error

    Dear, SEQanswer community,

    I try to get read distribution statistics after TopHat alignment using RseqQC. Specifically, I am interested in percentage of reads, that aligned to different regions (like CDS exons, UTR exons, intrond and intergenic regions). I used the following command line
    Code:
    /home/olha/RSeQC-2.6.4/build/scripts-2.7/read_distribution.py -i /home/olha/A15.bam -r /home/olha/hg19_RefSeq.bed
    However, I got the error message:
    Code:
    processing /home/olha/hg19_RefSeq.bed ...Traceback (most recent call last):
      File "/home/olha/RSeQC-2.6.4/build/scripts-2.7/read_distribution.py", line 303, in <module>
        main()
      File "/home/olha/RSeQC-2.6.4/build/scripts-2.7/read_distribution.py", line 192, in main
        intergenic_down1kb_base,intergenic_down5kb_base,intergenic_down10kb_base) = process_gene_model(options.ref_gene_model)
      File "/home/olha/RSeQC-2.6.4/build/scripts-2.7/read_distribution.py", line 76, in process_gene_model
        utr_3 = obj.getUTR(utr=3)
      File "/usr/local/lib/python2.7/dist-packages/RSeQC-2.3.7-py2.7-linux-x86_64.egg/qcmodule/BED.py", line 425, in getUTR
        txStart=int(fields[1])
    ValueError: invalid literal for int() with base 10: 'NM_032291'
    Does it something wrong with my bam or bed file? I used hg19 for aligning to the reference genome. Also I tried different bam files and come up with the same mistake.

    Thank You very much for help!

Latest Articles

Collapse

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 06-05-2026, 10:09 AM
0 responses
14 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 06-04-2026, 08:59 AM
0 responses
24 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 06-02-2026, 12:03 PM
0 responses
31 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 06-02-2026, 11:40 AM
0 responses
23 views
0 reactions
Last Post SEQadmin2  
Working...