Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Java error message when using own data #26

Closed
melnel000 opened this issue Oct 18, 2018 · 4 comments
Closed

Java error message when using own data #26

melnel000 opened this issue Oct 18, 2018 · 4 comments

Comments

@melnel000
Copy link

I can run mutation-server locally using the provided test data (input fastq files). However, when I use my own fastq files, I get the following error message (below). Can you offer any assistance on what the problem might be?

Exception in thread "main" java.lang.NullPointerException
at genepi.mut.pileup.PileupToolLocal.run(PileupToolLocal.java:118)
at genepi.base.Tool.start(Tool.java:193)
at genepi.base.Toolbox.start(Toolbox.java:44)
at genepi.mut.Server.main(Server.java:42)

@haansi
Copy link
Collaborator

haansi commented Oct 18, 2018

Hi there, can you please provide the commands you are trying to execute for both files?

@melnel000
Copy link
Author

These are the commands I am using:

Map mitochondrial reads to mitochondrial genome

/opt/exp_soft/bwa-0.7.12/bwa mem /home/melnel000/mtDNA-server/rCRS.fasta
/home/melnel000/mtDNA-server/MTS81_S2_L001_R1_001.fastq_1.gz /home/melnel000/mtDNA-server/MTS81_S2_L001_R2_001.fastq_1.gz
/home/melnel000/mtDNA-server/MTS81.sam

Convert sam to bam file, sort and index

/usr/bin/samtools view -bS /home/melnel000/mtDNA-server/MTS81.sam | samtools sort - /home/melnel000/mtDNA-server/MTS81_sorted
/usr/bin/samtools index /home/melnel000/mtDNA-server/MTS81_sorted.bam

/opt/exp_soft/java/jdk1.8.0_31/bin/java -jar /home/melnel000/mtDNA-server/mutation-server-1.1.9.jar analyse-local --input
/home/melnel000/mtDNA-server/home/melnel000/mtDNA-server/MTS81_sorted.bam --indel --output /home/melnel000/mtDNA-server/MTS81_sorted_results
--reference /home/melnel000/mtDNA-server/rCRS.fasta --level 0.01

When I run the same script with the small_r1.fastq and small_r2.fastq files I get an output file.

Thanks for the assistance

@seppinho
Copy link
Owner

seppinho commented Oct 18, 2018

Can you double check the path. Looks like we can't find a bam file (I think it's a copy/paste error of /home/melnel000)
/home/melnel000/mtDNA-server/home/melnel000/mtDNA-server/MTS81_sorted.bam

@seppinho
Copy link
Owner

New version will include a check for this.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants