@@ -34,9 +34,9 @@ expand the data to a directory called ``ocean-short``.
3434This is a toy dataset. It is based on real data, but the samples were trimmed
3535so that they contains only 250k paired-end reads.
3636
37- The dataset is organized in classical MOCAT style. Ngless does not require this
38- structure, but this tutorial also demonstrates how to upgrade from your
39- existing MOCAT-based projects. ::
37+ The dataset is organized in so that each directory contains a sample with
38+ multiple fastq files). NGLess does not require this structure, but it is
39+ convenient ::
4040
4141 $ find
4242 ./SAMEA2621229.sampled
@@ -63,10 +63,9 @@ The rest of this tutorial is an explanation of the steps in this script.
6363
6464To run ngless, we need write a script. We start with a few imports::
6565
66- ngless "0.7"
67- import "parallel" version "0.6"
68- import "mocat" version "0.0"
69- import "omrgc" version "0.0"
66+ ngless "1.4"
67+ import "parallel" version "1.0"
68+ import "omrgc" version "1.0"
7069
7170These will all be used in the tutorial.
7271
@@ -94,7 +93,7 @@ sample.
9493
9594First, we load the data (the FastQ files)::
9695
97- input = load_mocat_sample (sample)
96+ input = load_fastq_directory (sample)
9897
9998And, now, we preprocess the data::
10099
@@ -153,15 +152,14 @@ Full script
153152
154153Here is the full script::
155154
156- ngless "0.8 "
155+ ngless "1.4 "
157156 import "parallel" version "0.0"
158- import "mocat" version "0.0"
159- import "omrgc" version "0.0"
157+ import "omrgc" version "1.0"
160158
161159
162160 samples = readlines('tara.demo.short')
163161 sample = lock1(samples)
164- input = load_mocat_sample (sample)
162+ input = load_fastq_directory (sample)
165163
166164 input = preprocess(input, keep_singles=False) using |read|:
167165 read = substrim(read, min_quality=25)
0 commit comments