Project

General

Profile

Wiki » History » Revision 4

Revision 3 (Katie Lennard, 09/19/2022 05:04 PM) → Revision 4/26 (Katie Lennard, 09/20/2022 10:34 AM)

# Wiki

# Data location:

The data was transferred from Athena medmicro):

```
/MedMicro/Clinton/CRE Pfizer Feb 2022/CRE study_1A_results_17022022
/MedMicro/Clinton/CRE Pfizer Feb 2022/CRE study_1B_results_21022022
```

to Ilifu:

```
/scratch3/users/katiel/Clinton/CRE_study_August_2022/
```

# Reference #Reference data:

Klebsiella pneumoniae – strain HS11286 (GenBank accession no. CP003200.1) (n=18);
Serratia marcescens – strain KS10 (GenBank accession no. CP027798.1) (n=3);
Escherichia coli – strain ATCC 25922 (GenBank accession no. CP009072.1) (n=1); and
Enterobacter cloacae – strain ATCC 13047 (GenBank accession no. NC_014121.1) (n=1).

```
/scratch3/users/katiel/Clinton/CRE_study_August_2022/ref_genomes
```

# Objectives #Objectives workflow:
![workflow.png]()

# QC: #QC:
11 sample failed QC phred scores before trimming and filtering; none failed after filtering and trimming. Filtering and trimming were executed as follows:

```
nextflow run kviljoen/fastq_QC --reads '/scratch3/users/katiel/Clinton/CRE_study_August_2022/raw/study_1A_B_combined/*_R{1,2}_001.fastq.gz' -profile ilifu
```
QC reports can be found in the 'files' tab

# AMR profiling
The preference from Clinton is to do AMR profiling with the ResFinder DB. I'm getting errors there that I think relate to the header formatting though so in the interim have run with the ARG_annot DB that we used for previous projects as:
```
nextflow run kviljoen/uct-srst2 --reads '/scratch3/users/katiel/Clinton/CRE_study_August_2022/2022-09-19-fastq_QC/bbduk/*_{1,2}.fq' -profile ilifu --gene_db /cbio/users/katie/Nicol/Ps_aerug_srst2_MLST/ARGannot_r3.fasta --outdir /scratch3/users/katiel/Clinton/CRE_study_August_2022/srst2_resFinder/coverage_80_run --min_gene_cov 80
```