sorted BAM output

[ullo...@googlemail.com]

Dear all,

I would like to have a sorted bam from mapping to a custom reference. It seems that I'm using the wrong syntax. Could somebody please guide me? Using ngless v1.2 samtools module v 1.0

Best,

Ulrike

SARScov2_mapped = map(non_human_reads, fafile='allSARScov2variants210221.fna',allbest=true)

write(mapstats(SARScov2_mapped),ofile=RESULTS</> 'SARScov2.stats210221.txt')
counts = count(SARScov2_mapped,
            features=['seqname'],
            normalization={raw})

collect(counts,
        current=current,
        allneeded=samples,
        ofile='SARScov2.raw.210221.txt')
sorted = samtools_sort(SARScov2_mapped)
write(sorted, ofile=RESULTS<\> 'SARScov2.sorted.bam')
write(SARScov2_mapped, ofile=RESULTS<\> 'SARScov2.fastq.gz')

Exiting after fatal error while loading and running script
Script Error (there is likely an error in your script)
"human.ngl" (line 77, column 29):
unexpected '\\'
expecting "#", "//", "/*", "{", "'", "\"", digit, "0x", "_", letter, "!=", "==", "</>", "<=", "<", ">=", ">", "+", "*", " ", tab, ";", "\r", "\n" or end of input

[Luis Pedro Coelho]

(I'm just reading on my phone, while on the subway, so could be wrong, but)

The operator is </>, inspired on Unix-like path separators.

Hth
Luis

--
Luis Pedro Coelho

[ullo...@googlemail.com]

True, typo.

Thanks

[Luis Pedro Coelho]

I just committed a slight improvement so that you'd get a much better error message in the future:

https://github.com/ngless-toolkit/ngless/commit/6f63c99f37b9254cb10dd73bbb46c05ac1f7264a

HTH,

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/b120c95a-d88b-490b-a2a9-08f8e77de3b1n%40googlegroups.com.