Question: GESS - exon skipping scanner
0
gravatar for hinkel2
23 months ago by
hinkel210
hinkel210 wrote:

Hi!

I'm trying to detect exon skipping events from my RNAseq data, following the described workflow. First, I aligned my reads to reference genome using tophat2. Output bam file was sorted using samtools and used for further analysis by GESS as recommended:

python GESS.py -i myfile_sorted.bam -g Arabidopsis_thaliana.TAIR10.31.dna.genome.fa -o outputfolder -n outputdata -l 5

whereafter GESS starts running

ARGUMENTS LIST: name = outputdata RNA-seq bamfile = /path/to/myfile_sorted.bam Params l = 5 INFO @ Thu, 20 Apr 2017 18:13:26: Initializing... INFO @ Thu, 20 Apr 2017 18:13:26: Constructing SpliceGraph: Intron gap linking...

However, after 5 to 10 minutes GESS stops running showing following error message:

Traceback (most recent call last): File "GESS.py", line 122, in <module> main() File "GESS.py", line 60, in main global_data = GESS_graphBuild(options,global_data) File "/path/to/gess/impl/SpliceSiteGraphProc.py", line 145, in GESS_graphBuild chr_ssgraph,chr_depthdict = ssGraphProc.process() File "/path/to/gess/impl/SpliceSiteGraphProc.py", line 53, in process self.chromSplicGraphMap = self.updateToDiGraph() File "/path/to/gess/impl/SpliceSiteGraphProc.py", line 103, in updateToDiGraph refseq = SeqIO.parse(genomedir + chrid + '.fa','fasta').next() File "/usr/lib/python2.7/dist-packages/Bio/SeqIO/__init__.py", line 583, in parse with as_handle(handle, mode) as fp: File "/usr/lib/python2.7/contextlib.py", line 17, in __enter__ return self.gen.next() File "/usr/lib/python2.7/dist-packages/Bio/File.py", line 90, in as_handle with open(handleish, mode, **kwargs) as fp: IOError: [Errno 2] No such file or directory: './Pt.fa'

I have absolutely no idea what's going wrong here. Especially I dont know where the './Pt.fa' comes from. Is there someone who faced this problem before and/or could help me? :)

Thanks in advance!

ADD COMMENTlink modified 20 months ago by adglink0 • written 23 months ago by hinkel210
0
gravatar for adglink
20 months ago by
adglink0
adglink0 wrote:

I think it's because some where in your bam file you have reads aligned to a chromosome referred as "Pt". So GESS is expecting there to be Pt.fa file.

ADD COMMENTlink written 20 months ago by adglink0
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: 1086 users visited in the last hour