wip/nanolyse 1.2.0.20200915 NetBSD

Remove reads mapping to the lambda phage genome from a fastq file. This script uses Heng Li's minimap2 (https://github.com/lh3/minimap2) and his mappy Python binding (https://pypi.python.org/pypi/mappy).

### Warning

If (some of) the reads of your genome of interest are sufficiently similar to the lambda genome those reads will be lost.

### Examples

gunzip -c reads.fastq.gz | NanoLyse | gzip > reads_without_lambda.fastq.gz

In combination with NanoFilt (https://github.com/wdecoster/nanofilt):

gunzip -c reads.fastq.gz | NanoLyse | NanoFilt -q 12 \ | gzip > filtered_reads_without_lambda.fastq.gz

Using a different genome to filter on (rather than lambda phage):

gunzip -c reads.fastq.gz | NanoLyse --reference mygenome.fa.gz \ | gzip > reads_without_mygenome.fastq.gz