SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Java Developer eandrade Academic/Non-Profit Jobs 1 08-22-2019 03:59 AM
GATK "Code exception" memento Bioinformatics 2 02-07-2012 12:41 AM
samtools tview produces "Floating point exception" on big file? pmaugeri Bioinformatics 2 10-28-2011 08:06 AM
GATK IndelRealigner Error:Code Exception HongxiangZheng Bioinformatics 1 12-14-2010 07:21 PM
Java for bioinformatics.. vani s kulkarni Bioinformatics 1 07-30-2009 10:19 AM

Reply
 
Thread Tools
Old 02-24-2012, 03:01 PM   #1
mhadidi2002
Member
 
Location: Germany

Join Date: Jun 2011
Posts: 24
Default Scripture Java Exception

Hello,

I have used scripture to do the segmentation task.

The command:

scripture$ java -jar scripture-beta2.jar -alignment lyco.sam -out lyc_scrip.txt -sizeFile spe.sizes.txt -chr SL2.40ch01 -chrSequences S_lycopersicum.fa



It gives me exceptions while running, I have copied the lines, it is repeated, but here I am showing it only for one position (below)
I also have some empty files, which are:
  1. extendedpieces.bed
  2. firstGraph.bed
  3. igv.log
  4. introns.bed

The exception is pasted below.
Your help is highly appreciated.

##################################################

java.lang.NullPointerException
at broad.pda.gene.GeneTools.orientationFromSpliceSites(GeneTools.java:62 8)
at broad.pda.gene.GeneTools.orientationFromSpliceSites(GeneTools.java:65 6)
at broad.pda.seq.segmentation.ContinuousDataAlignmentModel.makeGraphWith Counts(ContinuousDataAlignmentModel.java:1298)
at broad.pda.seq.segmentation.ContinuousDataAlignmentModel.contiguousAcr ossGraph(ContinuousDataAlignmentModel.java:1212)
at broad.pda.seq.segmentation.ContinuousDataAlignmentModel.main(Continuo usDataAlignmentModel.java:2306)
SL2.40ch01 90217921 90217933 SL2.40ch01:90217921-90217933 0

##################################################
mhadidi2002 is offline   Reply With Quote
Old 06-27-2012, 01:08 PM   #2
AlbertOuTw
Junior Member
 
Location: Taipei

Join Date: Aug 2010
Posts: 6
Default

Hi mhadidi2002,

From your command line you use -chr SL2.40ch01 to name a chromosome but as far as my using experience of scripture, a "legal" chromosome name "MUST" start with "chr". So if you name your chromosome chrSL2.40ch01 it might work.

Just make sure the chromosome id in all *.sam/*.bam, sizeFile and chromosomeFile must start with chr. The easiest way is to modify all chromosome ids in the chromosomeFiles and re-built the bowtie index and remapping them again. Or you can use some scripting language to modify them. Good Luck!
AlbertOuTw is offline   Reply With Quote
Reply

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 06:59 AM.


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