BWA error mapping RNA-seq data to human reference

[Natasha Arora]

Dear ngless Team,

I have RNA-seq data from human body sites, and I would like to remove the human reads (keeping all microorganism sequences). Managed to install ngless on a macOS Big Sur after Luis pointed me out to the location for the file (https://github.com/ngless-toolkit/ngless/releases/tag/v1.2.0), and after changing permissions to get it to run. The code I am running is very simple, and am doing so on a test fastq file containing only a small number of reads. Please note that in the code I had to specify I was using ngless 1.1 in order to get the reference hg19 to download (did not work if I specified I was using ngless 1.2 as it could not access the website to download it).

$ngless "1.1"

$input = fastq('test.fq')

$mapped_hg19 = map(input, reference='hg19')

The error I encounter is with bwa mem and I wonder if it has to do with not having paired end reads. Maybe I need to specify I have single read sequencing somewhere, but how? 

"Failed command

Executable used:: /Users/natashaarora/Documents/Programmes/NGLess-1/bin/../share/ngless/bin/ngless-1.2.0-bwa

Command line was::

mem -t 1 -K 100000000 /Users/natashaarora/Documents/Programmes/NGLess-1/bin/../share/ngless/data/References/hg19/Sequence/BWAIndex/reference.fa-bwa-0.7.17.gz -p -

Error code was 1.

Stdout: <captured by inner process>

Stderr: [M::bwa_idx_load_from_disk] read 0 ALT contigs

[E::main_mem] fail to open file `-p'."

I apologise in advance if I have missed something very basic. Thank you very much for your help.

Natasha

[Luis Pedro Coelho]

Dear Natasha,

Indeed, I don't see any reason why this should not work. In fact, it works for me (on Linux), as it should.

I CC Shaojun who did the MacOS X building. @Shaojun: can you test Natasha's issue on Mac and see if you see the same issue.

Thanks

Luis

Luis Pedro Coelho | Fudan University | http://luispedro.org

https://orcid.org/0000-0002-9280-7885

--

You received this message because you are subscribed to the Google Groups "NGLess" group.

To unsubscribe from this group and stop receiving emails from it, send an email to ngless+un...@googlegroups.com.

To view this discussion on the web visit https://groups.google.com/d/msgid/ngless/47ce50b4-6ffa-4b27-af23-f1d3c8f5beadn%40googlegroups.com.