tseemann/nullarbor

Error - report not written and empty samples

開放

#262 建立於 2020年3月12日

 (1 則留言) (0 個反應) (1 位負責人)Perl (39 個分叉)github user discovery
help wanted

倉庫指標

星標
 (144 顆星)
PR 合併指標
 (PR 指標待抓取)

描述

Hi, We are running Nullarbor on a mixture of samples we sequenced ourselves, and some downloaded from ENA. When we run, we don't get a final report. Looking more closely, it looks like, for the samples we downlaoded from ENA, the assembled contigs (fasta) and snps (vcf) are empty (fasta all gaps and vcf is just the header). The nohup.out says there has been some error with java - I've compressed it and added it below.

I ran the offending samples through a workflow I use for mapping/snp calling (bwa mem, samtools, freebayes) and I found that there is a problem with bwa mem throwing an error due to orphan reads in the split paired read files generated by using the ena toolkit 'fastq-dump' command. When I download the .fastq files as a single interleaved file seems to work fine. Not sure if this is contributing to the problem or not.

The command we ran was: Command 1: nullarbor2.pl --name ancientA --ref /home/ubuntu/volume_sdb/Anthrax/AncientA/CZC5_NZ_AP018443.1.fasta --input ./ancientA.tab --outdir ancientA Command 2: nohup nice make -j 2 -C /home/ubuntu/volume_sdb/Nullarbor/AncientA/ancientA &

Any assistance would be greatly appreciated, thank you so much!

nohup.out.gz

貢獻者指南