I have been attempting to use CASH to analyze a dataset for alternative splicing, but I am having difficulty getting the program to even read in my gtf.
[al2025@premise isoform]$ java -jar -Xmx900m /pylon2/mc3bg6p/al2025/storage/cash_v2.1.0-alpha1/cash.jar \ > --Case:stress /mnt/lustre/macmaneslab/al2025/isoform/STAR_mappings/blu7-x_female_gonad_stress_STAR_anno_Aligned.sortedByCoord.out.bam \ > --Control:control /mnt/lustre/macmaneslab/al2025/isoform/STAR_mappings/L-G107_male_gonad_control_STAR_anno_Aligned.sortedByCoord.out.bam \ > --GTF /mnt/lustre/macmaneslab/al2025/isoform/genome_compare/Columba_livia_1.0.gtf \ > --Output /mnt/lustre/macmaneslab/al2025/isoform/CASH_out cash version 2.1.0-alpha1 2017-09-19 15:58 [runtime:0ms] [com.novelbio.nbcgui.controlseq.CtrlSplicing:main]-[INFO] Start reading gtf file /mnt/lustre/macmaneslab/al2025/isoform/genome_compare/Columba_livia_1.0.gtf java.lang.NullPointerException at com.novelbio.analysis.seq.genome.gffOperate.GffHashGeneAbs.getMapName2Detail(GffHashGeneAbs.java:152) at com.novelbio.listOperate.ListHashSearch.ReadGffarray(ListHashSearch.java:247) at com.novelbio.analysis.seq.genome.gffOperate.GffHashGeneAbs.ReadGffarray(GffHashGeneAbs.java:50) at com.novelbio.analysis.seq.genome.gffOperate.GffHashGene.read(GffHashGene.java:187) at com.novelbio.analysis.seq.genome.gffOperate.GffHashGene.<init>(GffHashGene.java:90) at com.novelbio.nbcgui.controlseq.CtrlSplicing.mainCmd(CtrlSplicing.java:355) at com.novelbio.nbcgui.controlseq.CtrlSplicing.main(CtrlSplicing.java:67) Exception in thread "main" com.novelbio.analysis.seq.genome.ExceptionNbcGFF: GffFile Formate Error:/mnt/lustre/macmaneslab/al2025/isoform/genome_compare/Columba_livia_1.0.gtf at com.novelbio.listOperate.ListHashSearch.ReadGffarray(ListHashSearch.java:255) at com.novelbio.analysis.seq.genome.gffOperate.GffHashGeneAbs.ReadGffarray(GffHashGeneAbs.java:50) at com.novelbio.analysis.seq.genome.gffOperate.GffHashGene.read(GffHashGene.java:187) at com.novelbio.analysis.seq.genome.gffOperate.GffHashGene.<init>(GffHashGene.java:90) at com.novelbio.nbcgui.controlseq.CtrlSplicing.mainCmd(CtrlSplicing.java:355) at com.novelbio.nbcgui.controlseq.CtrlSplicing.main(CtrlSplicing.java:67) Caused by: java.lang.NullPointerException at com.novelbio.analysis.seq.genome.gffOperate.GffHashGeneAbs.getMapName2Detail(GffHashGeneAbs.java:152) at com.novelbio.listOperate.ListHashSearch.ReadGffarray(ListHashSearch.java:247) ... 5 more
my gtf file is formatted as follows
NC_013978.1 RefSeq exon 1 69 . + . transcript_id "rna19564"; NC_013978.1 RefSeq exon 70 1042 . + . transcript_id "rna19565"; NC_013978.1 RefSeq exon 1043 1114 . + . transcript_id "rna19566"; NC_013978.1 RefSeq exon 1115 2700 . + . transcript_id "rna19567"; NC_013978.1 RefSeq exon 2701 2774 . + . transcript_id "rna19568"; NC_013978.1 RefSeq CDS 2786 3751 . + 0 transcript_id "gene15879"; gene_id "gene15879"; gene_name "ND1"; NC_013978.1 RefSeq exon 3769 3839 . + . transcript_id "rna19569"; NC_013978.1 RefSeq exon 3845 3915 . - . transcript_id "rna19570"; NC_013978.1 RefSeq exon 3915 3983 . + . transcript_id "rna19571"; NC_013978.1 RefSeq CDS 3984 5023 . + 0 transcript_id "gene15880"; gene_id "gene15880"; gene_name "ND2"; NC_013978.1 RefSeq exon 5024 5094 . + . transcript_id "rna19572"; NC_013978.1 RefSeq exon 5096 5164 . - . transcript_id "rna19573"; NC_013978.1 RefSeq exon 5167 5239 . - . transcript_id "rna19574"; NC_013978.1 RefSeq exon 5242 5308 . - . transcript_id "rna19575"; NC_013978.1 RefSeq exon 5308 5379 . - . transcript_id "rna19576";
Any insight into what I am doing wrong?
Log in to post a comment.
I have been attempting to use CASH to analyze a dataset for alternative splicing, but I am having difficulty getting the program to even read in my gtf.
my gtf file is formatted as follows
Any insight into what I am doing wrong?