...
| Code Block |
|---|
#!/bin/bash -l
#PBS -N nano-Q
#PBS -l walltime=24:00:00
#PBS -l mem=16gb
#PBS -l ncpus=8
cd $PBS_O_WORKDIR
################################################################################################################################
# USER DEFINE VARIABLES
################################################################################################################################
SAMPLEID=NC483
REFNAME=NC001477_RefGenome
LIST='/work/phylo/OxfordNanopore/nextflow/assembly/data/NC483/NC483_NC001477_ref_list.txt'
REF='/work/phylo/OxfordNanopore/nextflow/assembly/data/NC483/NC483_NC001477_reference_sequence.fasta'
ONT='/work/phylo/OxfordNanopore/nextflow/assembly/data/NC483/NC483_FAU10290_pass_barcode96_0cf303ee.fastq'
################################################################################################################################
#activate conda environment containing tools for analysis
conda activate nanoQ2
#generic mapping reads
minimap2 -a $REF $ONT > ${SAMPLEID}_aln.sam
#mapping noisy reads
#minimap2 -ax $REF $ONT > ${SAMPLEID}_aln2.sam
#Samtools
samtools view -bt ${LIST} -o ${SAMPLEID}_aln.bam ${SAMPLEID}_aln.sam
samtools sort -T /tmp/aln.sorted -o ${SAMPLEID}_aln.sorted.bam ${SAMPLEID}_aln.bam
samtools index ${SAMPLEID}_aln.sorted.bam
#run nano-Q
python /work/phylo/OxfordNanopore/nextflow/tools/git/Nano-Q/nano-q.py -b ${SAMPLEID}_aln.sorted.bam -c 1 -l 9000 -nr 1 -q 5 -j 10 |
creating a ref_list.txt file
| Code Block |
|---|
#use a tool provided by emboss
infoseq $REF -only -name -length | sed 1d > ${REF}_list.txt |
submit the job to the scheduler
...