Running SeroBA

We are now ready to use SeroBA to determine the serotype of our samples. Move into the data directory where we keep the reads. In this case we have called the directory run_seroba and it is in the directory called data.


In [ ]:
cd data/run_seroba

Have a quick look at the contents of the directory:


In [ ]:
ls -al

As you can see, there are two gzipped fastq files for each sample, one for forward reads and one for reverse reads.

Now run SeroBA using the runSerotyping command:

`seroba runSerotyping database forward_read reverse_read prefix`

Where the options are:

database       path to the database directory
forward_read   forward read file in fastq format
reverse_read   reverse read file in fastq format
prefix         a unique prefix

In [ ]:
seroba runSerotyping ../database/ sample1_1.fq.gz \
    sample1_2.fq.gz sample1

If you are running Seroba on the Sanger farm and instead want to use the central database, you can use:

seroba runSerotyping $SEROBA_DB forward_read reverse_read prefix

Lets have a look at the results in the next section. You can also return to the index or revisit the previous section.