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.