SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Tophat Fusion Post shai_rosenberg RNA Sequencing 0 06-09-2015 07:14 AM
TopHat Fusion post (Error) samar Bioinformatics 5 07-14-2014 09:52 PM
Tophat Fusion Post Processing filtering out confirmed positive fusion mknut Bioinformatics 2 01-09-2014 07:00 AM
Tophat2 with fusion search and tophat-fusion-post problems seqfast Bioinformatics 9 07-30-2013 06:16 PM
MACS: ValueError: invalid literal for int() with base 10: '1:N:0:' feralBiologist Bioinformatics 0 07-13-2013 02:01 AM

Reply
 
Thread Tools
Old 03-17-2017, 07:11 AM   #1
komalsrathi
Member
 
Location: Philadelphia, PA

Join Date: Aug 2013
Posts: 14
Default tophat-fusion-post: ValueError: invalid literal for int() with base 10: 'exonCount'

Hi everyone,

I am running tophat-fusion-post like this:

Code:
tophat-fusion-post -o ./fusion_results -p 8 --num-fusion-reads 1 --num-fusion-pairs 2 --num-fusion-both 5 /mnt/isilon/cbmi/variome/reference/bowtie_indexes/hg38_no_alt/Homo_sapiens/NCBI/GRCh38/Sequence/BowtieIndex/genome
My root folder is tophat-fusion. I have 4 tophat output folders under it: CHP212, SHSY5Y, SKNAS and SKNSH, each of which contain a fusions.out file. I have created symbolic links to refGene.txt, ensGene.txt and blast database (blast) in the same folder.

This is my directory structure where I have run tophat:

Code:
$ tree -L 2 ./tophat-fusion

./
|-- CHP212
|   |-- accepted_hits.bam
|   |-- align_summary.txt
|   |-- deletions.bed
|   |-- fusions.out
|   |-- insertions.bed
|   |-- junctions.bed
|   |-- logs
|   |-- prep_reads.info
|   `-- unmapped.bam
|-- CHP212.sh
|-- CHP212_R1.fastq.gz -> /mnt/isilon/maris_lab/target_nbl_ngs/CellLineRNASeq/rawfiles/cat_fastq/CHP212_R1.fastq.gz
|-- CHP212_R2.fastq.gz -> /mnt/isilon/maris_lab/target_nbl_ngs/CellLineRNASeq/rawfiles/cat_fastq/CHP212_R2.fastq.gz
|-- IMR32.sh
|-- IMR32_R1.fastq.gz -> /mnt/isilon/maris_lab/target_nbl_ngs/CellLineRNASeq/rawfiles/cat_fastq/IMR32_R1.fastq.gz
|-- IMR32_R2.fastq.gz -> /mnt/isilon/maris_lab/target_nbl_ngs/CellLineRNASeq/rawfiles/cat_fastq/IMR32_R2.fastq.gz
|-- SHSY5Y
|   |-- accepted_hits.bam
|   |-- align_summary.txt
|   |-- deletions.bed
|   |-- fusions.out
|   |-- insertions.bed
|   |-- junctions.bed
|   |-- logs
|   |-- prep_reads.info
|   `-- unmapped.bam
|-- SHSY5Y.sh
|-- SHSY5Y_R1.fastq.gz -> /mnt/isilon/maris_lab/target_nbl_ngs/CellLineRNASeq/rawfiles/cat_fastq/SHSY5Y_R1.fastq.gz
|-- SHSY5Y_R2.fastq.gz -> /mnt/isilon/maris_lab/target_nbl_ngs/CellLineRNASeq/rawfiles/cat_fastq/SHSY5Y_R2.fastq.gz
|-- SKNAS
|   |-- accepted_hits.bam
|   |-- align_summary.txt
|   |-- deletions.bed
|   |-- fusions.out
|   |-- insertions.bed
|   |-- junctions.bed
|   |-- logs
|   |-- prep_reads.info
|   `-- unmapped.bam
|-- SKNAS.sh
|-- SKNAS_R1.fastq.gz -> /mnt/isilon/maris_lab/target_nbl_ngs/CellLineRNASeq/rawfiles/cat_fastq/SKNAS_R1.fastq.gz
|-- SKNAS_R2.fastq.gz -> /mnt/isilon/maris_lab/target_nbl_ngs/CellLineRNASeq/rawfiles/cat_fastq/SKNAS_R2.fastq.gz
|-- SKNSH
|   |-- accepted_hits.bam
|   |-- align_summary.txt
|   |-- deletions.bed
|   |-- fusions.out
|   |-- insertions.bed
|   |-- junctions.bed
|   |-- logs
|   |-- prep_reads.info
|   `-- unmapped.bam
|-- SKNSH.sh
|-- SKNSH_R1.fastq.gz -> /mnt/isilon/maris_lab/target_nbl_ngs/CellLineRNASeq/rawfiles/cat_fastq/SKNSH_R1.fastq.gz
|-- SKNSH_R2.fastq.gz -> /mnt/isilon/maris_lab/target_nbl_ngs/CellLineRNASeq/rawfiles/cat_fastq/SKNSH_R2.fastq.gz
|-- blast -> /mnt/isilon/cbmi/variome/reference/blast_db/hg38
|-- ensGene.txt -> /mnt/isilon/cbmi/variome/reference/blast_db/hg38/ensGene.txt
|-- fusion_results
|   |-- fusion_seq.bwtout
|   |-- fusion_seq.fa
|   |-- fusion_seq.map
|   |-- logs
|   `-- tmp
|-- refGene.txt -> /mnt/isilon/cbmi/variome/reference/blast_db/hg38/refGene.txt
`-- tophat-fusion.sh
When I run tophat-fusion-post under this directory, I am getting the following errors:

Code:
[Fri Mar 17 15:09:18 2017] Beginning TopHat-Fusion post-processing run (v2.1.0)
-----------------------------------------------
[Fri Mar 17 15:09:18 2017] Extracting 23-mer around fusions and mapping them using Bowtie
[Fri Mar 17 15:09:30 2017] Filtering fusions
Traceback (most recent call last):
  File "/home/rathik/tools/miniconda3/envs/fusion-env/bin/tophat-fusion-post", line 2924, in <module>
    sys.exit(main())
  File "/home/rathik/tools/miniconda3/envs/fusion-env/bin/tophat-fusion-post", line 2895, in main
    filter_fusion(bwt_idx_prefix, params)
  File "/home/rathik/tools/miniconda3/envs/fusion-env/bin/tophat-fusion-post", line 965, in filter_fusion
    ensGene_list = read_genes("ensGene.txt")
  File "/home/rathik/tools/miniconda3/envs/fusion-env/bin/tophat-fusion-post", line 917, in read_genes
    num_exons = int(line[7])
ValueError: invalid literal for int() with base 10: 'exonCount'
komalsrathi is offline   Reply With Quote
Old 06-14-2019, 12:35 AM   #2
Enraico
Junior Member
 
Location: Venice

Join Date: Feb 2012
Posts: 5
Default

remove header line from the ensGene.txt file
Code:
mv ensGene.txt ensGene.txt.bk
grep -v "^#" ensGene.txt.bk > ensGene.txt
Enraico is offline   Reply With Quote
Old 04-28-2020, 11:07 PM   #3
warrenfelsh
Junior Member
 
Location: dubai

Join Date: Apr 2020
Posts: 1
Default

The error message invalid literal for int() with base 10 would seem to indicate that you are passing a string that's not an integer to the int() function . In other words it's either empty, or has a character in it other than a digit.

You can solve this error by using Python isdigit() method to check whether the value is number or not. The returns True if all the characters are digits, otherwise False .

if val.isdigit():

The other way to overcome this issue is to wrap your code inside a Python try...except block to handle this error.

Python2.x and Python3.x

Sometimes the difference between Python2.x and Python3.x that leads to this ValueError: invalid literal for int() with base 10 .

With Python2.x , int(str(3/2)) gives you "1". With Python3.x , the same gives you ("1.5"): ValueError: invalid literal for int() with base 10: "1.5".
warrenfelsh is offline   Reply With Quote
Reply

Tags
tophat-fusion-post, tophat2.1.0

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 09:08 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO