This document seeks to lay out my process in poking at the DNAsequencing results of a series of Pseudomonas aeruginosa PA14 and PAK strains.
If I understand Dr. Lee and co.’s goal, they wish to ensure that these strains are still reasonably close to the associated reference strains. I therefore am running my default trimming/mapping/variant search methods.
I have a single command that can run all of these commands at the same time, but I have been actively breaking my tools recently; so I decided to run them one at a time with the assumption that something would not work (but everything did work on the first try, so that was nice).
I downloaded the .zip archive file using the link in Dr. Lee’s email. I did not save it though, so if we need to download the data again, we will have to go to him. I created my usual work directory ‘preprocessing/’ within this tree and moved it there. I unzipped it and moved each pair of reads to a directory which follows Dr. Lee’s desired naming convention.
I then created the directories: ‘reference/’ and ‘sample_sheets/’. The sample_sheets remained empty for a while, but I immediately downloaded the full genbank flat file for the Pseudomonas PAK strain from NCBI, found here:
https://www.ncbi.nlm.nih.gov/nuccore/NZ_CP020659
Note, that when downloading, one must hit the ‘customize view’ button on the right and ensure that the entire sequence and all annotations are included. Then hit the ‘send to’ button and send it to a file. This file I copied to reference/paeruginosa_pak.gb.
Given the full PAK genbank file, I converted it to the expected fasta/gff file for mapping:
cd reference
cyoa --method gb2gff --input paeruginosa_pak.gb
This command created a series of fasta and gff files which provide the coordinates for the various annotations (genes/cds/rRNA/intercds) and sequence for the genome, CDS nucleotides, and amino acids. I then copied the genome/gff files to my global reference directory and prepared it for usage by my favorite mapper:
cd ~/libraries/genome
cyoa --method indexhisat --species paeruginosa_pak
Now all of the pieces are in place for me to play. Each of the following steps was performed twice, once for the PA14 samples, once for the PAK samples. The only difference in the invocations was due to the fact that the PAK annotations provide different tags. E.g. I used the ‘Alias’ tag for PA14 and the ‘locus_tag’ tag for PAK. As a result I am only going to write down in this document the PA14 invocations and assume the reader can figure out the difference.
I have a couple of trimming methods, in this instance I just used the default and will operate under the assumption that it is sufficient until I see otherwise.
cd preprocessing
start=$(pwd)
for i in $(/bin/ls -d PA14*); do
cd $i
cyoa --method trim --input $(/bin/ls *.fastq.gz | tr '\n' ':' | sed 's/:$//g')
cd $start
done
The above command line invocation produced a series of trimming jobs which when examined look like this (I am only showing examples from PA14_exoUTY, and am leaving off the beginning and end).
## This is a portion of file:
## preprocessing/PA14_exoUTY/scripts/01trim_7_UTY_S138_R1_001.sh
module add trimomatic
mkdir -p outputs/01trimomatic
## Note that trimomatic prints all output and errors to STDERR, so send both to output
trimmomatic PE \
-threads 1 \
-phred33 \
\
7_UTY_S138_R1_001.fastq.gz 7_UTY_S138_R2_001.fastq.gz \
7_UTY_S138_R1_001-trimmed_paired.fastq 7_UTY_S138_R1_001-trimmed_unpaired.fastq \
7_UTY_S138_R2_001-trimmed_paired.fastq 7_UTY_S138_R2_001-trimmed_unpaired.fastq \
ILLUMINACLIP:/fs/cbcb-software/RedHat-8-x86_64/local/cyoa/202302/prefix/lib/perl5/auto/share/dist/Bio-Adventure/genome/adapters.fa:2:20:10:2:keepBothReads \
SLIDINGWINDOW:4:20 MINLEN:50 1>outputs/01trimomatic/7_UTY_S138_R1_001-trimomatic.stdout \
2>outputs/01trimomatic/7_UTY_S138_R1_001-trimomatic.stderr
excepted=$( { grep "Exception" "outputs/01trimomatic/7_UTY_S138_R1_001-trimomatic.stdout" || test $? = 1; } )
One thing I did not include in the above: upon completion, the script aggressively compresses the trimmed output and symbolically links it to r1_trimmed.fastq.xz and r2_trimmed.fastq.xz. Thus any following steps can use the same input name (r1_trimmed.fastq.xz:r2_trimmed.fastq.xz).
My default mappers run the actual alignment, convert it to a compressed/indexed bam, and count it against the reference genome. In this context, the counting is a little silly, but does have the potential to help find duplications and such.
cd preprocessing
start=$(pwd)
for i in $(/bin/ls -d PA14*); do
cd $i
cyoa --method hisat --input r1_trimmed.fastq.xz:r2_trimmed.fastq.xz \
--species paeruginosa_pa14 --gff_type gene --gff_tag Alias
cd $start
done
## Here is what I ran for PAK
cd preprocessing
start=$(pwd)
for i in $(/bin/ls -d PAK*); do
cd $i
cyoa --method hisat --input r1_trimmed.fastq.xz:r2_trimmed.fastq.xz \
--species paeruginosa_pak --gff_type gene --gff_tag locus_tag
cd $start
done
Similarly, I am just putting the meaty part.
module add hisat2 samtools htseq bamtools
mkdir -p outputs/40hisat2_paeruginosa_pa14
hisat2 -x ${HOME}/libraries/genome/indexes/paeruginosa_pa14 \
-p 8 \
-q -1 <(less /home/trey/sshfs/scratch/atb/dnaseq/paeruginosa_strains_202304/preprocessing/PA14_exoUTY/r1_trimmed.fastq.xz) -2 <(less /home/trey/sshfs/scratch/atb/dnaseq/paeruginosa_strains_202304/preprocessing/PA14_exoUTY/r2_trimmed.fastq.xz) \
--phred33 \
--un outputs/40hisat2_paeruginosa_pa14/unaldis_paeruginosa_pa14_genome.fastq \
--al outputs/40hisat2_paeruginosa_pa14/aldis_paeruginosa_pa14_genome.fastq \
--un-conc outputs/40hisat2_paeruginosa_pa14/unalcon_paeruginosa_pa14_genome.fastq \
--al-conc outputs/40hisat2_paeruginosa_pa14/alcon_paeruginosa_pa14_genome.fastq \
-S outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.sam \
2>outputs/40hisat2_paeruginosa_pa14/hisat2_paeruginosa_pa14_genome_PA14_exoUTY.stderr \
1>outputs/40hisat2_paeruginosa_pa14/hisat2_paeruginosa_pa14_genome_PA14_exoUTY.stdout
The above cyoa invocation also creates this script. It is a little long because it does some checks and creates a couple of filtered versions of the output.
module add samtools bamtools
echo "Starting samtools"
if [[ -f "outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam" && -f "outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.sam" ]]; then
echo "Both the bam and sam files exist, rerunning."
elif [[ -f "outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam" ]]; then
echo "The output file exists, quitting."
exit 0
elif [[ ! -f "outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.sam" ]]; then
echo "Could not find the samtools input file."
exit 1
fi
## If a previous sort file exists due to running out of memory,
## then we need to get rid of them first.
## hg38_100_genome-sorted.bam.tmp.0000.bam
if [[ -f "outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam.tmp.000.bam" ]]; then
rm -f outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam.tmp.*.bam
fi
samtools view -u -t ${HOME}/libraries/genome/paeruginosa_pa14.fasta \
-S outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.sam -o outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam \
2>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam_samtools.stderr \
1>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam_samtools.stdout
echo "First samtools command finished with $?"
samtools sort -l 9 outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam \
-o outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-sorted.bam \
2>>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam_samtools.stderr \
1>>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam_samtools.stdout
rm outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam
rm outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.sam
mv outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-sorted.bam outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam
samtools index outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam \
2>>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam_samtools.stderr \
1>>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam_samtools.stdout
echo "Second samtools command finished with $?"
bamtools stats -in outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam \
2>>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam_samtools.stats 1>&2
echo "Bamtools finished with $?"
## The following will fail if this is single-ended.
samtools view -b -f 2 \
-o outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-paired.bam \
\
outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam 2>>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam_samtools.stderr \
1>>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam_samtools.stdout
samtools index outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-paired.bam \
2>>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam_samtools.stderr \
1>>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam_samtools.stdout
bamtools stats -in outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-paired.bam \
2>>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam_samtools.stats 1>&2
bamtools filter -tag XM:0 \
-in outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam \
-out outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-sorted_nomismatch.bam \
2>>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam_samtools.stats 1>&2
echo "bamtools filter finished with: $?"
samtools index \
\
outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-sorted_nomismatch.bam 2>>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam_samtools.stderr \
1>>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam_samtools.stdout
echo "final samtools index finished with: $?"
Note that this step is not really useful for a dnaseq dataset in most instances. I also have the default orientation set to reverse because most of the samples off our sequencer are reversed; but that is likely not true for this dataset. If it turns out we actually care about these counts, I may need to come back and rerun these.
module add htseq
htseq-count \
-q -f bam \
-s reverse -a 0 \
--type all --idattr Alias \
\
outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-paired.bam \
/home/trey/libraries/genome/paeruginosa_pa14.gff 2>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-paired_sreverse_all_Alias.stderr \
1>outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-paired_sreverse_all_Alias.count
xz -f -9e outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-paired_sreverse_all_Alias.count
I tend to like to use freebayes for this. It is a little conservative, but I think it seems to work quite well. I can also use mpileup and snippy. freebayes and mpileup are setup to feed a post-processing script which I think is kind of fun and will be decribed momentarily.
cd preprocessing
start=$(pwd)
for i in $(/bin/ls -d PA14*); do
cd $i
cyoa --method freebayes \
--input outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-paired.bam \
--species paeruginosa_pa14 --gff_type gene --gff_tag Alias --intron 0
cd $start
done
Unlike hisat, I include the conversion to the binary/compressed/indexed format with the invocation of the variant search. I also include the duplicate search functionality from gatk.
module add gatk freebayes libgsl libhts samtools bcftools vcftools
mkdir -p outputs/50freebayes_paeruginosa_pa14
gatk MarkDuplicates \
-I outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome.bam \
-O outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14_genome_deduplicated.bam \
-M outputs/50freebayes_paeruginosa_pa14/deduplication_stats.txt --REMOVE_DUPLICATES true --COMPRESSION_LEVEL 9 \
2>outputs/50freebayes_paeruginosa_pa14/deduplication.stderr \
1>outputs/50freebayes_paeruginosa_pa14/deduplication.stdout
echo "Finished gatk deduplication." >> outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.stdout
samtools index outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14_genome_deduplicated.bam
echo "Finished samtools index." >> outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.stdout
freebayes -f /home/trey/libraries/genome/paeruginosa_pa14.fasta \
-v outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.vcf \
\
outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14_genome_deduplicated.bam 1>>outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.stdout \
2>>outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.stderr
echo "Finished freebayes." >> outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.stdout
bcftools convert outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.vcf \
-Ob -o outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.bcf \
2>>outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.stderr \
1>>outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.stdout
echo "Finished bcftools convert." >> outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.stdout
bcftools index outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.bcf \
2>>outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.stderr \
1>>outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.stdout
echo "Finished bcftools index." >> outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.stdout
rm outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.vcf
The result from the above freebayes script is a bcf containing the high-quality observed variants. The cyoa invocation also creates the following script, which will require a bit of explanation.
use Bio::Adventure::SNP;
my $result = $h->Bio::Adventure::SNP::SNP_Ratio_Worker(
'outputs/50freebayes_paeruginosa_pa14/paeruginosa_pa14.bcf',
input => 'paeruginosa_pa14',
species => 'freebayes',
vcf_method => '5',
vcf_cutoff => '0.8',
vcf_minpct => 'Alias',
gff_tag => 'gene',
gff_type => 'outputs/50freebayes_paeruginosa_pa14',
output_dir => 'outputs/50freebayes_paeruginosa_pa14/all_tags.txt',
output => 'outputs/50freebayes_paeruginosa_pa14/count.txt',
output_count => 'outputs/50freebayes_paeruginosa_pa14/modified.fasta',
output_genome => 'outputs/50freebayes_paeruginosa_pa14/variants_by_gene.txt',
output_by_gene => 'outputs/50freebayes_paeruginosa_pa14/pkm.txt',
output_pkm => );
The function ‘SNP_Ratio_Worker()’ reads the reference genome, the set of variants, and the genome annotations in order to create a new copy of the genome (modified.fasta) which should be equivalent to the input reads. It also rewrites the bcf data into a matrix which is easier to play with in R/python (all_tags.txt). Finally, it uses the annotation information to explicitly show the amino acid substitions observed in every ORF (variants_by_gene.txt). In theory it should also give a rpkm-esque copy of the variants observed / ORF, but I turned that off because it doesn’t seem very useful and it is a little tricky to get right.
In order to play further with the data, I will need a sample sheet. So I will start out by creating a blank one in excel (libreoffice) which contains only the samplenames in the same format as my directories in preprocessing/.
Once completed, I can use it as the input for my hpgltools package and it should extract the interesting information from the preprocessing logs and fill out the sample sheet accordingly. Lets see if it works!
Here is the before:
::kable(extract_metadata("sample_sheets/all_samples.xlsx")) knitr
## Did not find the condition column in the sample sheet.
## Filling it in as undefined.
## Did not find the batch column in the sample sheet.
## Filling it in as undefined.
sampleid | condition | batch | |
---|---|---|---|
PA14_exoUTY | PA14_exoUTY | undefined | undefined |
PA14_JC | PA14_JC | undefined | undefined |
PA14_lux | PA14_lux | undefined | undefined |
PA14_NBH | PA14_NBH | undefined | undefined |
PA14_pscD_A5 | PA14_pscD_A5 | undefined | undefined |
PA14_pscdD_E4 | PA14_pscdD_E4 | undefined | undefined |
PA14_xcp | PA14_xcp | undefined | undefined |
PA14_xcp_pscD | PA14_xcp_pscD | undefined | undefined |
PAK | PAK | undefined | undefined |
PAK_pscC | PAK_pscC | undefined | undefined |
PAK_xcp | PAK_xcp | undefined | undefined |
PAK_xcp_pscC | PAK_xcp_pscC | undefined | undefined |
Like I said, not much going on. Lets see what it looks like after I run the gatherer on it… (Note, I have been meaning to change this to drop the unused columns, but not yet).
<- gather_preprocessing_metadata("sample_sheets/all_samples.xlsx") modified
## Did not find the condition column in the sample sheet.
## Filling it in as undefined.
## Did not find the batch column in the sample sheet.
## Filling it in as undefined.
## Starting assembly_fasta_nt.
## Example filename: preprocessing/PA14_exoUTY/unicycler_assembly.fasta.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for: preprocessing/PA14_exoUTY/unicycler_assembly.fasta.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for: preprocessing/PA14_JC/unicycler_assembly.fasta.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for: preprocessing/PA14_lux/unicycler_assembly.fasta.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for: preprocessing/PA14_NBH/unicycler_assembly.fasta.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for: preprocessing/PA14_pscD_A5/unicycler_assembly.fasta.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for: preprocessing/PA14_pscdD_E4/unicycler_assembly.fasta.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for: preprocessing/PA14_xcp/unicycler_assembly.fasta.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for: preprocessing/PA14_xcp_pscD/unicycler_assembly.fasta.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for: preprocessing/PAK/unicycler_assembly.fasta.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for: preprocessing/PAK_pscC/unicycler_assembly.fasta.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for: preprocessing/PAK_xcp/unicycler_assembly.fasta.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for: preprocessing/PAK_xcp_pscC/unicycler_assembly.fasta.
## Starting assembly_genbank_annotated.
## Example filename: preprocessing/PA14_exoUTY/outputs/*mergeannot/PA14_exoUTY.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*mergeannot/PA14_exoUTY.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_JC/outputs/*mergeannot/PA14_JC.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_lux/outputs/*mergeannot/PA14_lux.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_NBH/outputs/*mergeannot/PA14_NBH.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*mergeannot/PA14_pscD_A5.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*mergeannot/PA14_pscdD_E4.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp/outputs/*mergeannot/PA14_xcp.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*mergeannot/PA14_xcp_pscD.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for: preprocessing/PAK/outputs/*mergeannot/PAK.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_pscC/outputs/*mergeannot/PAK_pscC.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp/outputs/*mergeannot/PAK_xcp.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*mergeannot/PAK_xcp_pscC.gbk.
## Starting assembly_genbank_stripped.
## Example filename: preprocessing/PA14_exoUTY/outputs/*mergeannot/PA14_exoUTY_stripped.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*mergeannot/PA14_exoUTY_stripped.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_JC/outputs/*mergeannot/PA14_JC_stripped.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_lux/outputs/*mergeannot/PA14_lux_stripped.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_NBH/outputs/*mergeannot/PA14_NBH_stripped.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*mergeannot/PA14_pscD_A5_stripped.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*mergeannot/PA14_pscdD_E4_stripped.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp/outputs/*mergeannot/PA14_xcp_stripped.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*mergeannot/PA14_xcp_pscD_stripped.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK/outputs/*mergeannot/PAK_stripped.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_pscC/outputs/*mergeannot/PAK_pscC_stripped.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp/outputs/*mergeannot/PAK_xcp_stripped.gbk.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*mergeannot/PAK_xcp_pscC_stripped.gbk.
## Starting assembly_cds_amino_acids.
## Example filename: preprocessing/PA14_exoUTY/outputs/*merge_cds_predictions/PA14_exoUTY.faa.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*merge_cds_predictions/PA14_exoUTY.faa.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_JC/outputs/*merge_cds_predictions/PA14_JC.faa.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_lux/outputs/*merge_cds_predictions/PA14_lux.faa.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_NBH/outputs/*merge_cds_predictions/PA14_NBH.faa.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*merge_cds_predictions/PA14_pscD_A5.faa.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*merge_cds_predictions/PA14_pscdD_E4.faa.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp/outputs/*merge_cds_predictions/PA14_xcp.faa.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*merge_cds_predictions/PA14_xcp_pscD.faa.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK/outputs/*merge_cds_predictions/PAK.faa.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_pscC/outputs/*merge_cds_predictions/PAK_pscC.faa.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp/outputs/*merge_cds_predictions/PAK_xcp.faa.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*merge_cds_predictions/PAK_xcp_pscC.faa.
## Starting assembly_cds_nucleotides.
## Example filename: preprocessing/PA14_exoUTY/outputs/*merge_cds_predictions/PA14_exoUTY.ffn.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*merge_cds_predictions/PA14_exoUTY.ffn.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_JC/outputs/*merge_cds_predictions/PA14_JC.ffn.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_lux/outputs/*merge_cds_predictions/PA14_lux.ffn.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_NBH/outputs/*merge_cds_predictions/PA14_NBH.ffn.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*merge_cds_predictions/PA14_pscD_A5.ffn.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*merge_cds_predictions/PA14_pscdD_E4.ffn.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp/outputs/*merge_cds_predictions/PA14_xcp.ffn.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*merge_cds_predictions/PA14_xcp_pscD.ffn.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK/outputs/*merge_cds_predictions/PAK.ffn.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_pscC/outputs/*merge_cds_predictions/PAK_pscC.ffn.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp/outputs/*merge_cds_predictions/PAK_xcp.ffn.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*merge_cds_predictions/PAK_xcp_pscC.ffn.
## Starting assembly_gff.
## Example filename: preprocessing/PA14_exoUTY/outputs/*merge_cds_predictions/PA14_exoUTY.gff.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*merge_cds_predictions/PA14_exoUTY.gff.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_JC/outputs/*merge_cds_predictions/PA14_JC.gff.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_lux/outputs/*merge_cds_predictions/PA14_lux.gff.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_NBH/outputs/*merge_cds_predictions/PA14_NBH.gff.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*merge_cds_predictions/PA14_pscD_A5.gff.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*merge_cds_predictions/PA14_pscdD_E4.gff.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp/outputs/*merge_cds_predictions/PA14_xcp.gff.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*merge_cds_predictions/PA14_xcp_pscD.gff.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK/outputs/*merge_cds_predictions/PAK.gff.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_pscC/outputs/*merge_cds_predictions/PAK_pscC.gff.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp/outputs/*merge_cds_predictions/PAK_xcp.gff.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*merge_cds_predictions/PAK_xcp_pscC.gff.
## Starting assembly_tsv.
## Example filename: preprocessing/PA14_exoUTY/outputs/*merge_cds_predictions/PA14_exoUTY.tsv.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*merge_cds_predictions/PA14_exoUTY.tsv.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_JC/outputs/*merge_cds_predictions/PA14_JC.tsv.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_lux/outputs/*merge_cds_predictions/PA14_lux.tsv.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_NBH/outputs/*merge_cds_predictions/PA14_NBH.tsv.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*merge_cds_predictions/PA14_pscD_A5.tsv.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*merge_cds_predictions/PA14_pscdD_E4.tsv.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp/outputs/*merge_cds_predictions/PA14_xcp.tsv.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*merge_cds_predictions/PA14_xcp_pscD.tsv.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK/outputs/*merge_cds_predictions/PAK.tsv.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_pscC/outputs/*merge_cds_predictions/PAK_pscC.tsv.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp/outputs/*merge_cds_predictions/PAK_xcp.tsv.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*merge_cds_predictions/PAK_xcp_pscC.tsv.
## Starting assembly_xls.
## Example filename: preprocessing/PA14_exoUTY/outputs/*mergeannot/PA14_exoUTY.xlsx.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*mergeannot/PA14_exoUTY.xlsx.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_JC/outputs/*mergeannot/PA14_JC.xlsx.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_lux/outputs/*mergeannot/PA14_lux.xlsx.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_NBH/outputs/*mergeannot/PA14_NBH.xlsx.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*mergeannot/PA14_pscD_A5.xlsx.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*mergeannot/PA14_pscdD_E4.xlsx.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp/outputs/*mergeannot/PA14_xcp.xlsx.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*mergeannot/PA14_xcp_pscD.xlsx.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for: preprocessing/PAK/outputs/*mergeannot/PAK.xlsx.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_pscC/outputs/*mergeannot/PAK_pscC.xlsx.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp/outputs/*mergeannot/PAK_xcp.xlsx.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*mergeannot/PAK_xcp_pscC.xlsx.
## Starting input_r1.
## Example filename: preprocessing/PA14_exoUTY/scripts/*trim_*.sh.
## Starting input_r2.
## Example filename: preprocessing/PA14_exoUTY/scripts/*trim_*.sh.
## Starting trimomatic_input.
## Example filename: preprocessing/PA14_exoUTY/outputs/*trimomatic/*-trimomatic.stderr.
## Starting trimomatic_output.
## Example filename: preprocessing/PA14_exoUTY/outputs/*trimomatic/*-trimomatic.stderr.
## Starting trimomatic_ratio.
## The numerator column is: trimomatic_output.
## The denominator column is: trimomatic_input.
## Starting host_filter_species.
## Example filename: preprocessing/PA14_exoUTY/host_species.txt.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: host_filter_species.
## Starting hisat_genome_single_concordant.
## Example filename: preprocessing/PA14_exoUTY/outputs/*/hisat2_*.stderr.
## Starting hisat_genome_multi_concordant.
## Example filename: preprocessing/PA14_exoUTY/outputs/*/hisat2_*.stderr.
## Starting hisat_genome_single_all.
## Example filename: preprocessing/PA14_exoUTY/outputs/*/hisat2_*.stderr.
## Starting hisat_genome_multi_all.
## Example filename: preprocessing/PA14_exoUTY/outputs/*/hisat2_*.stderr.
## Starting hisat_count_table.
## Example filename: preprocessing/PA14_exoUTY/outputs/*/*_*_genome*.count.xz.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*/*_*_genome*.count.xz.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_pscC/outputs/*/*_*_genome*.count.xz.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for: preprocessing/PAK_xcp/outputs/*/*_*_genome*.count.xz.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*/*_*_genome*.count.xz.
## Starting jellyfish_count_table.
## Example filename: preprocessing/PA14_exoUTY/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_JC/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_lux/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_NBH/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_pscC/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_filename_search(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*jellyfish_*/*_matrix.csv.xz.
## Starting jellyfish_observed.
## Example filename: preprocessing/PA14_exoUTY/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*jellyfish_*/*_matrix.csv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*jellyfish_*/*_matrix.csv.xz.
## Starting kraken_viral_classified.
## Example filename: preprocessing/PA14_exoUTY/outputs/*kraken_viral*/kraken.stderr.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: kraken_viral_classified.
## Starting kraken_viral_unclassified.
## Example filename: preprocessing/PA14_exoUTY/outputs/*kraken_viral*/kraken.stderr.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: kraken_viral_unclassified.
## Starting kraken_first_viral_species.
## Example filename: preprocessing/PA14_exoUTY/outputs/*kraken_viral*/kraken_report.txt.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: kraken_first_viral_species.
## Starting kraken_first_viral_species_reads.
## Example filename: preprocessing/PA14_exoUTY/outputs/*kraken_viral*/kraken_report.txt.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: kraken_first_viral_species_reads.
## Starting kraken_standard_classified.
## Example filename: preprocessing/PA14_exoUTY/outputs/*kraken_standard*/kraken.stderr.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: kraken_standard_classified.
## Starting kraken_standard_unclassified.
## Example filename: preprocessing/PA14_exoUTY/outputs/*kraken_standard*/kraken.stderr.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: kraken_standard_unclassified.
## Starting kraken_first_standard_species.
## Example filename: preprocessing/PA14_exoUTY/outputs/*kraken_standard*/kraken_report.txt.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: kraken_first_standard_species.
## Starting kraken_first_standard_species_reads.
## Example filename: preprocessing/PA14_exoUTY/outputs/*kraken_standard*/kraken_report.txt.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: kraken_first_standard_species_reads.
## Starting possible_host_species.
## Example filename: preprocessing/PA14_exoUTY/outputs/*filter_kraken_host/*.log.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*filter_kraken_host/*.log.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*filter_kraken_host/*.log.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*filter_kraken_host/*.log.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*filter_kraken_host/*.log.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*filter_kraken_host/*.log.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*filter_kraken_host/*.log.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*filter_kraken_host/*.log.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*filter_kraken_host/*.log.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*filter_kraken_host/*.log.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*filter_kraken_host/*.log.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*filter_kraken_host/*.log.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*filter_kraken_host/*.log.
## Starting pernt_mean_coverage.
## Example filename: preprocessing/PA14_exoUTY/outputs/??assembly_coverage_*/base_coverage.tsv.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: pernt_mean_coverage.
## Starting pernt_median_coverage.
## Example filename: preprocessing/PA14_exoUTY/outputs/??assembly_coverage_*/base_coverage.tsv.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: pernt_median_coverage.
## Starting pernt_min_coverage.
## Example filename: preprocessing/PA14_exoUTY/outputs/??assembly_coverage_*/base_coverage.tsv.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: pernt_min_coverage.
## Starting pernt_max_coverage.
## Example filename: preprocessing/PA14_exoUTY/outputs/??assembly_coverage_*/base_coverage.tsv.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: pernt_max_coverage.
## Starting salmon_mapped.
## Example filename: preprocessing/PA14_exoUTY/outputs/*salmon_*/salmon.stderr.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: salmon_mapped.
## Starting shovill_contigs.
## Example filename: preprocessing/PA14_exoUTY/outputs/*shovill_*/shovill.log.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: shovill_contigs.
## Starting shovill_length.
## Example filename: preprocessing/PA14_exoUTY/outputs/*shovill_*/shovill.log.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: shovill_length.
## Starting shovill_estlength.
## Example filename: preprocessing/PA14_exoUTY/outputs/*shovill_*/shovill.log.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: shovill_estlength.
## Starting shovill_minlength.
## Example filename: preprocessing/PA14_exoUTY/outputs/*shovill_*/shovill.log.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: shovill_minlength.
## Starting unicycler_lengths.
## Example filename: preprocessing/PA14_exoUTY/outputs/*unicycler/*final_assembly.fasta.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: unicycler_lengths.
## Starting unicycler_relative_coverage.
## Example filename: preprocessing/PA14_exoUTY/outputs/*unicycler/*final_assembly.fasta.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: unicycler_relative_coverage.
## Starting filtered_relative_coverage.
## Example filename: preprocessing/PA14_exoUTY/outputs/??filter_depth/final_assembly.fasta.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: filtered_relative_coverage.
## Starting phastaf_num_hits.
## Example filename: preprocessing/PA14_exoUTY/outputs/*phastaf_*/phage.bed.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*phastaf_*/phage.bed.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*phastaf_*/phage.bed.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*phastaf_*/phage.bed.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*phastaf_*/phage.bed.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*phastaf_*/phage.bed.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*phastaf_*/phage.bed.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*phastaf_*/phage.bed.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*phastaf_*/phage.bed.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*phastaf_*/phage.bed.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*phastaf_*/phage.bed.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*phastaf_*/phage.bed.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*phastaf_*/phage.bed.
## Starting phageterm_dtr_length.
## Example filename: preprocessing/PA14_exoUTY/outputs/*phageterm_*/phageterm_final_dtr.fasta.
## Warning in dispatch_fasta_lengths(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*phageterm_*/phageterm_final_dtr.fasta.
## Warning in dispatch_fasta_lengths(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_JC/outputs/*phageterm_*/phageterm_final_dtr.fasta.
## Warning in dispatch_fasta_lengths(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_lux/outputs/*phageterm_*/phageterm_final_dtr.fasta.
## Warning in dispatch_fasta_lengths(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_NBH/outputs/*phageterm_*/phageterm_final_dtr.fasta.
## Warning in dispatch_fasta_lengths(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*phageterm_*/phageterm_final_dtr.fasta.
## Warning in dispatch_fasta_lengths(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*phageterm_*/phageterm_final_dtr.fasta.
## Warning in dispatch_fasta_lengths(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp/outputs/*phageterm_*/phageterm_final_dtr.fasta.
## Warning in dispatch_fasta_lengths(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*phageterm_*/phageterm_final_dtr.fasta.
## Warning in dispatch_fasta_lengths(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK/outputs/*phageterm_*/phageterm_final_dtr.fasta.
## Warning in dispatch_fasta_lengths(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_pscC/outputs/*phageterm_*/phageterm_final_dtr.fasta.
## Warning in dispatch_fasta_lengths(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp/outputs/*phageterm_*/phageterm_final_dtr.fasta.
## Warning in dispatch_fasta_lengths(meta, input_file_spec, verbose = verbose, :
## The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*phageterm_*/phageterm_final_dtr.fasta.
## Starting prodigal_positive_strand.
## Example filename: preprocessing/PA14_exoUTY/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*prodigal_*/predicted_cds.gff.
## Starting prodigal_negative_strand.
## Example filename: preprocessing/PA14_exoUTY/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*prodigal_*/predicted_cds.gff.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*prodigal_*/predicted_cds.gff.
## Starting glimmer_positive_strand.
## Example filename: preprocessing/PA14_exoUTY/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*glimmer/glimmer3.predict.
## Starting glimmer_negative_strand.
## Example filename: preprocessing/PA14_exoUTY/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*glimmer/glimmer3.predict.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*glimmer/glimmer3.predict.
## Starting phanotate_positive_strand.
## Example filename: preprocessing/PA14_exoUTY/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*phanotate/*_phanotate.tsv.xz.
## Starting phanotate_negative_strand.
## Example filename: preprocessing/PA14_exoUTY/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*phanotate/*_phanotate.tsv.xz.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*phanotate/*_phanotate.tsv.xz.
## Starting final_gc_content.
## Example filename: preprocessing/PA14_exoUTY/outputs/*prokka/PA14_exoUTY.fna.
## Warning in dispatch_gc(meta, input_file_spec, basedir = basedir, verbose =
## verbose): The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*prokka/PA14_exoUTY.fna.
## Warning in dispatch_gc(meta, input_file_spec, basedir = basedir, verbose =
## verbose): The input file is NA for:
## preprocessing/PA14_JC/outputs/*prokka/PA14_JC.fna.
## Warning in dispatch_gc(meta, input_file_spec, basedir = basedir, verbose =
## verbose): The input file is NA for:
## preprocessing/PA14_lux/outputs/*prokka/PA14_lux.fna.
## Warning in dispatch_gc(meta, input_file_spec, basedir = basedir, verbose =
## verbose): The input file is NA for:
## preprocessing/PA14_NBH/outputs/*prokka/PA14_NBH.fna.
## Warning in dispatch_gc(meta, input_file_spec, basedir = basedir, verbose =
## verbose): The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*prokka/PA14_pscD_A5.fna.
## Warning in dispatch_gc(meta, input_file_spec, basedir = basedir, verbose =
## verbose): The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*prokka/PA14_pscdD_E4.fna.
## Warning in dispatch_gc(meta, input_file_spec, basedir = basedir, verbose =
## verbose): The input file is NA for:
## preprocessing/PA14_xcp/outputs/*prokka/PA14_xcp.fna.
## Warning in dispatch_gc(meta, input_file_spec, basedir = basedir, verbose =
## verbose): The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*prokka/PA14_xcp_pscD.fna.
## Warning in dispatch_gc(meta, input_file_spec, basedir = basedir, verbose =
## verbose): The input file is NA for: preprocessing/PAK/outputs/*prokka/PAK.fna.
## Warning in dispatch_gc(meta, input_file_spec, basedir = basedir, verbose =
## verbose): The input file is NA for:
## preprocessing/PAK_pscC/outputs/*prokka/PAK_pscC.fna.
## Warning in dispatch_gc(meta, input_file_spec, basedir = basedir, verbose =
## verbose): The input file is NA for:
## preprocessing/PAK_xcp/outputs/*prokka/PAK_xcp.fna.
## Warning in dispatch_gc(meta, input_file_spec, basedir = basedir, verbose =
## verbose): The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*prokka/PAK_xcp_pscC.fna.
## Starting interpro_signalp_hits.
## Example filename: preprocessing/PA14_exoUTY/outputs/*interproscan_*/PA14_exoUTY.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*interproscan_*/PA14_exoUTY.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*interproscan_*/PA14_JC.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*interproscan_*/PA14_lux.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*interproscan_*/PA14_NBH.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*interproscan_*/PA14_pscD_A5.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*interproscan_*/PA14_pscdD_E4.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*interproscan_*/PA14_xcp.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*interproscan_*/PA14_xcp_pscD.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*interproscan_*/PAK.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*interproscan_*/PAK_pscC.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*interproscan_*/PAK_xcp.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*interproscan_*/PAK_xcp_pscC.faa.tsv.
## Starting interpro_phobius_hits.
## Example filename: preprocessing/PA14_exoUTY/outputs/*interproscan_*/PA14_exoUTY.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*interproscan_*/PA14_exoUTY.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*interproscan_*/PA14_JC.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*interproscan_*/PA14_lux.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*interproscan_*/PA14_NBH.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*interproscan_*/PA14_pscD_A5.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*interproscan_*/PA14_pscdD_E4.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*interproscan_*/PA14_xcp.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*interproscan_*/PA14_xcp_pscD.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*interproscan_*/PAK.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*interproscan_*/PAK_pscC.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*interproscan_*/PAK_xcp.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*interproscan_*/PAK_xcp_pscC.faa.tsv.
## Starting interpro_pfam_hits.
## Example filename: preprocessing/PA14_exoUTY/outputs/*interproscan_*/PA14_exoUTY.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*interproscan_*/PA14_exoUTY.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*interproscan_*/PA14_JC.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*interproscan_*/PA14_lux.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*interproscan_*/PA14_NBH.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*interproscan_*/PA14_pscD_A5.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*interproscan_*/PA14_pscdD_E4.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*interproscan_*/PA14_xcp.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*interproscan_*/PA14_xcp_pscD.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*interproscan_*/PAK.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*interproscan_*/PAK_pscC.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*interproscan_*/PAK_xcp.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*interproscan_*/PAK_xcp_pscC.faa.tsv.
## Starting interpro_tmhmm_hits.
## Example filename: preprocessing/PA14_exoUTY/outputs/*interproscan_*/PA14_exoUTY.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*interproscan_*/PA14_exoUTY.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*interproscan_*/PA14_JC.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*interproscan_*/PA14_lux.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*interproscan_*/PA14_NBH.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*interproscan_*/PA14_pscD_A5.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*interproscan_*/PA14_pscdD_E4.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*interproscan_*/PA14_xcp.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*interproscan_*/PA14_xcp_pscD.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*interproscan_*/PAK.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*interproscan_*/PAK_pscC.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*interproscan_*/PAK_xcp.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*interproscan_*/PAK_xcp_pscC.faa.tsv.
## Starting interpro_cdd_hits.
## Example filename: preprocessing/PA14_exoUTY/outputs/*interproscan_*/PA14_exoUTY.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*interproscan_*/PA14_exoUTY.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*interproscan_*/PA14_JC.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*interproscan_*/PA14_lux.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*interproscan_*/PA14_NBH.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*interproscan_*/PA14_pscD_A5.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*interproscan_*/PA14_pscdD_E4.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*interproscan_*/PA14_xcp.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*interproscan_*/PA14_xcp_pscD.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*interproscan_*/PAK.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*interproscan_*/PAK_pscC.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*interproscan_*/PAK_xcp.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*interproscan_*/PAK_xcp_pscC.faa.tsv.
## Starting interpro_smart_hits.
## Example filename: preprocessing/PA14_exoUTY/outputs/*interproscan_*/PA14_exoUTY.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*interproscan_*/PA14_exoUTY.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*interproscan_*/PA14_JC.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*interproscan_*/PA14_lux.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*interproscan_*/PA14_NBH.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*interproscan_*/PA14_pscD_A5.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*interproscan_*/PA14_pscdD_E4.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*interproscan_*/PA14_xcp.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*interproscan_*/PA14_xcp_pscD.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*interproscan_*/PAK.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*interproscan_*/PAK_pscC.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*interproscan_*/PAK_xcp.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*interproscan_*/PAK_xcp_pscC.faa.tsv.
## Starting interpro_gene3d_hits.
## Example filename: preprocessing/PA14_exoUTY/outputs/*interproscan_*/PA14_exoUTY.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*interproscan_*/PA14_exoUTY.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*interproscan_*/PA14_JC.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*interproscan_*/PA14_lux.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*interproscan_*/PA14_NBH.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*interproscan_*/PA14_pscD_A5.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*interproscan_*/PA14_pscdD_E4.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*interproscan_*/PA14_xcp.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*interproscan_*/PA14_xcp_pscD.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*interproscan_*/PAK.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*interproscan_*/PAK_pscC.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*interproscan_*/PAK_xcp.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*interproscan_*/PAK_xcp_pscC.faa.tsv.
## Starting interpro_superfamily_hits.
## Example filename: preprocessing/PA14_exoUTY/outputs/*interproscan_*/PA14_exoUTY.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_exoUTY/outputs/*interproscan_*/PA14_exoUTY.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_JC/outputs/*interproscan_*/PA14_JC.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_lux/outputs/*interproscan_*/PA14_lux.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_NBH/outputs/*interproscan_*/PA14_NBH.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscD_A5/outputs/*interproscan_*/PA14_pscD_A5.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_pscdD_E4/outputs/*interproscan_*/PA14_pscdD_E4.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp/outputs/*interproscan_*/PA14_xcp.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PA14_xcp_pscD/outputs/*interproscan_*/PA14_xcp_pscD.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK/outputs/*interproscan_*/PAK.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_pscC/outputs/*interproscan_*/PAK_pscC.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp/outputs/*interproscan_*/PAK_xcp.faa.tsv.
## Warning in dispatch_count_lines(meta, search, input_file_spec, verbose =
## verbose, : The input file is NA for:
## preprocessing/PAK_xcp_pscC/outputs/*interproscan_*/PAK_xcp_pscC.faa.tsv.
## Starting tRNA_hits.
## Example filename: preprocessing/PA14_exoUTY/outputs/*prokka_*/PA14_exoUTY.log.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: tRNA_hits.
## Starting aragorn_tRNAs.
## Example filename: preprocessing/PA14_exoUTY/outputs/*aragorn/aragorn.txt.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: aragorn_tRNAs.
## Starting ictv_taxonomy.
## Example filename: preprocessing/PA14_exoUTY/outputs/*classify_*/*_filtered.tsv.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: ictv_taxonomy.
## Starting ictv_accession.
## Example filename: preprocessing/PA14_exoUTY/outputs/*classify_*/*_filtered.tsv.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: ictv_accession.
## Starting ictv_genus.
## Example filename: preprocessing/PA14_exoUTY/outputs/*classify_*/*_filtered.tsv.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: ictv_genus.
## Starting notes.
## Example filename: preprocessing/PA14_exoUTY/notes.txt.
## The first filename does not exist, assuming this method was not performed.
## Not including new entries for: notes.
## Writing new metadata to: sample_sheets/all_samples_modified.xlsx
## Deleting the file sample_sheets/all_samples_modified.xlsx before writing the tables.
::kable(extract_metadata("sample_sheets/all_samples_modified.xlsx")) knitr
## Did not find the condition column in the sample sheet.
## Filling it in as undefined.
## Did not find the batch column in the sample sheet.
## Filling it in as undefined.
rownames | assemblyfastant | assemblygenbankannotated | assemblygenbankstripped | assemblycdsaminoacids | assemblycdsnucleotides | assemblygff | assemblytsv | assemblyxls | inputr1 | inputr2 | trimomaticinput | trimomaticoutput | trimomaticpercent | hisatgenomesingleconcordant | hisatgenomemulticoncordant | hisatgenomesingleall | hisatgenomemultiall | hisatcounttable | jellyfishcounttable | jellyfishobserved | possiblehostspecies | phastafnumhits | phagetermdtrlength | prodigalpositivestrand | prodigalnegativestrand | glimmerpositivestrand | glimmernegativestrand | phanotatepositivestrand | phanotatenegativestrand | finalgccontent | interprosignalphits | interprophobiushits | interpropfamhits | interprotmhmmhits | interprocddhits | interprosmarthits | interprogene3dhits | interprosuperfamilyhits | condition | batch | |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
PA14_exoUTY | PA14_exoUTY | 0 | 0 | 4706372 | 4350712 | 0.924 | 4280240 | 25722 | 34468 | 25722 | preprocessing/PA14_exoUTY/outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-paired_sreverse_all_Alias.count.xz | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | undefined | undefined | |||||||||
PA14_JC | PA14_JC | 0 | 0 | 5786839 | 5336197 | 0.922 | 5275687 | 32239 | 23521 | 32239 | preprocessing/PA14_JC/outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-paired_sreverse_all_Alias.count.xz | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | undefined | undefined | |||||||||
PA14_lux | PA14_lux | 0 | 0 | 6622570 | 6099776 | 0.921 | 6028065 | 36827 | 24917 | 36827 | preprocessing/PA14_lux/outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-paired_sreverse_all_Alias.count.xz | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | undefined | undefined | |||||||||
PA14_NBH | PA14_NBH | 0 | 0 | 5151127 | 4581433 | 0.889 | 4516421 | 26915 | 32560 | 26915 | preprocessing/PA14_NBH/outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-paired_sreverse_all_Alias.count.xz | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | undefined | undefined | |||||||||
PA14_pscD_A5 | PA14_pscD_A5 | 0 | 0 | 5898210 | 5417082 | 0.918 | 5359077 | 32852 | 21806 | 32852 | preprocessing/PA14_pscD_A5/outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-paired_sreverse_all_Alias.count.xz | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | undefined | undefined | |||||||||
PA14_pscdD_E4 | PA14_pscdD_E4 | 0 | 0 | 0 | 0 | NaN | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | undefined | undefined | ||||||||||
PA14_xcp | PA14_xcp | 0 | 0 | 5683132 | 5214875 | 0.918 | 5134054 | 30352 | 39479 | 30352 | preprocessing/PA14_xcp/outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-paired_sreverse_all_Alias.count.xz | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | undefined | undefined | |||||||||
PA14_xcp_pscD | PA14_xcp_pscD | 0 | 0 | 2026150 | 1514509 | 0.747 | 1471930 | 10238 | 27094 | 10238 | preprocessing/PA14_xcp_pscD/outputs/40hisat2_paeruginosa_pa14/paeruginosa_pa14_genome-paired_sreverse_all_Alias.count.xz | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | undefined | undefined | |||||||||
PAK | PAK | 0 | 0 | 4779558 | 4318745 | 0.904 | 4049183 | 1093 | 179265 | 1093 | preprocessing/PAK/outputs/40hisat2_paeruginosa_pak/paeruginosa_pak_genome-paired_sreverse_all_locus_tag.count.xz | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | undefined | undefined | |||||||||
PAK_pscC | PAK_pscC | 0 | 0 | 5734960 | 5271470 | 0.919 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | undefined | undefined | ||||||||||
PAK_xcp | PAK_xcp | 0 | 0 | 4843414 | 4443669 | 0.917 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | undefined | undefined | ||||||||||
PAK_xcp_pscC | PAK_xcp_pscC | 0 | 0 | 5195158 | 4611474 | 0.888 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | undefined | undefined |
Lets go one small step further. I have a series of modified genomes as well as the reference. We can do a quickie tree of them: First I will copy each modified genome to the tree/ directory and rename them to the sampleID.
start=$(pwd)
mkdir tree
cd preprocessing
for i in $(/bin/ls -d PA*); do
cp $i/outputs/50*/modified.fasta ${start}/tree/${i}.fasta
done
cp ~/libraries/genome/paeruginosa_pa14.fa ${start}/tree
cp ~/libraries/genome/paeruginosa_pak.fa ${start}/tree
Oh, it turns out that at the time of this writing, I forgot to run 3 samples, so this section will need to be redone. But I can at least run it for the samples that I didn’t forget.
<- genomic_sequence_phylo("tree", root = "paeruginosa_pa14") funkytown
## Reading tree/PA14_exoUTY.fasta
## Reading tree/PA14_JC.fasta
## Reading tree/PA14_lux.fasta
## Reading tree/PA14_NBH.fasta
## Reading tree/PA14_pscD_A5.fasta
## Reading tree/PA14_pscD_E4.fasta
## Reading tree/PA14_xcp_pscD.fasta
## Reading tree/paeruginosa_pa14.fasta
## Reading tree/paeruginosa_pak.fasta
## Reading tree/PAK.fasta
plot(funkytown$phy)
::pander(sessionInfo()) pander
R version 4.2.0 (2022-04-22)
Platform: x86_64-pc-linux-gnu (64-bit)
locale: LC_CTYPE=en_US.UTF-8, LC_NUMERIC=C, LC_TIME=en_US.UTF-8, LC_COLLATE=en_US.UTF-8, LC_MONETARY=en_US.UTF-8, LC_MESSAGES=en_US.UTF-8, LC_PAPER=en_US.UTF-8, LC_NAME=C, LC_ADDRESS=C, LC_TELEPHONE=C, LC_MEASUREMENT=en_US.UTF-8 and LC_IDENTIFICATION=C
attached base packages: stats4, stats, graphics, grDevices, utils, datasets, methods and base
other attached packages: hpgltools(v.1.0), testthat(v.3.1.7), reticulate(v.1.28), SummarizedExperiment(v.1.28.0), GenomicRanges(v.1.50.2), GenomeInfoDb(v.1.34.9), IRanges(v.2.32.0), S4Vectors(v.0.36.2), MatrixGenerics(v.1.10.0), matrixStats(v.0.63.0), Biobase(v.2.58.0) and BiocGenerics(v.0.44.0)
loaded via a namespace (and not attached): rappdirs(v.0.3.3), rtracklayer(v.1.58.0), tidyr(v.1.3.0), ggplot2(v.3.4.2), clusterGeneration(v.1.3.7), bit64(v.4.0.5), knitr(v.1.42), DelayedArray(v.0.24.0), data.table(v.1.14.8), KEGGREST(v.1.38.0), RCurl(v.1.98-1.12), doParallel(v.1.0.17), generics(v.0.1.3), GenomicFeatures(v.1.50.4), callr(v.3.7.3), RhpcBLASctl(v.0.23-42), cowplot(v.1.1.1), usethis(v.2.1.6), RSQLite(v.2.3.1), shadowtext(v.0.1.2), bit(v.4.0.5), enrichplot(v.1.18.3), xml2(v.1.3.3), httpuv(v.1.6.9), viridis(v.0.6.2), xfun(v.0.38), hms(v.1.1.3), jquerylib(v.0.1.4), evaluate(v.0.20), promises(v.1.2.0.1), fansi(v.1.0.4), restfulr(v.0.0.15), progress(v.1.2.2), caTools(v.1.18.2), dbplyr(v.2.3.2), igraph(v.1.4.1), DBI(v.1.1.3), htmlwidgets(v.1.6.2), purrr(v.1.0.1), ellipsis(v.0.3.2), dplyr(v.1.1.1), backports(v.1.4.1), annotate(v.1.76.0), aod(v.1.3.2), biomaRt(v.2.54.1), vctrs(v.0.6.1), remotes(v.2.4.2), cachem(v.1.0.7), withr(v.2.5.0), ggforce(v.0.4.1), HDO.db(v.0.99.1), GenomicAlignments(v.1.34.1), treeio(v.1.22.0), prettyunits(v.1.1.1), kmer(v.1.1.2), DOSE(v.3.24.2), ape(v.5.7-1), lazyeval(v.0.2.2), crayon(v.1.5.2), genefilter(v.1.80.3), edgeR(v.3.40.2), pkgconfig(v.2.0.3), tweenr(v.2.0.2), nlme(v.3.1-162), pkgload(v.1.3.2), devtools(v.2.4.5), rlang(v.1.1.0), lifecycle(v.1.0.3), miniUI(v.0.1.1.1), downloader(v.0.4), filelock(v.1.0.2), BiocFileCache(v.2.6.1), rprojroot(v.2.0.3), polyclip(v.1.10-4), graph(v.1.76.0), Matrix(v.1.5-4), aplot(v.0.1.10), boot(v.1.3-28.1), processx(v.3.8.0), png(v.0.1-8), viridisLite(v.0.4.1), rjson(v.0.2.21), bitops(v.1.0-7), gson(v.0.1.0), KernSmooth(v.2.23-20), pander(v.0.6.5), Biostrings(v.2.66.0), blob(v.1.2.4), phylogram(v.2.1.0), stringr(v.1.5.0), qvalue(v.2.30.0), remaCor(v.0.0.11), gridGraphics(v.0.5-1), scales(v.1.2.1), memoise(v.2.0.1), GSEABase(v.1.60.0), magrittr(v.2.0.3), plyr(v.1.8.8), gplots(v.3.1.3), zlibbioc(v.1.44.0), compiler(v.4.2.0), scatterpie(v.0.1.8), BiocIO(v.1.8.0), RColorBrewer(v.1.1-3), lme4(v.1.1-32), Rsamtools(v.2.14.0), cli(v.3.6.1), XVector(v.0.38.0), urlchecker(v.1.0.1), patchwork(v.1.1.2), ps(v.1.7.4), MASS(v.7.3-58.3), mgcv(v.1.8-41), tidyselect(v.1.2.0), stringi(v.1.7.12), highr(v.0.10), yaml(v.2.3.7), GOSemSim(v.2.24.0), locfit(v.1.5-9.7), ggrepel(v.0.9.3), grid(v.4.2.0), sass(v.0.4.5), fastmatch(v.1.1-3), tools(v.4.2.0), parallel(v.4.2.0), rstudioapi(v.0.14), foreach(v.1.5.2), gridExtra(v.2.3), farver(v.2.1.1), ggraph(v.2.1.0), digest(v.0.6.31), shiny(v.1.7.4), Rcpp(v.1.0.10), broom(v.1.0.4), later(v.1.3.0), httr(v.1.4.5), AnnotationDbi(v.1.60.2), Rdpack(v.2.4), colorspace(v.2.1-0), brio(v.1.1.3), XML(v.3.99-0.14), fs(v.1.6.1), splines(v.4.2.0), yulab.utils(v.0.0.6), PROPER(v.1.30.0), tidytree(v.0.4.2), graphlayouts(v.0.8.4), ggplotify(v.0.1.0), plotly(v.4.10.1), sessioninfo(v.1.2.2), xtable(v.1.8-4), jsonlite(v.1.8.4), nloptr(v.2.0.3), ggtree(v.3.6.2), tidygraph(v.1.2.3), ggfun(v.0.0.9), R6(v.2.5.1), RUnit(v.0.4.32), profvis(v.0.3.7), pillar(v.1.9.0), htmltools(v.0.5.5), mime(v.0.12), glue(v.1.6.2), fastmap(v.1.1.1), minqa(v.1.2.5), clusterProfiler(v.4.6.2), BiocParallel(v.1.32.6), codetools(v.0.2-19), fgsea(v.1.24.0), pkgbuild(v.1.4.0), mvtnorm(v.1.1-3), utf8(v.1.2.3), lattice(v.0.20-45), bslib(v.0.4.2), tibble(v.3.2.1), sva(v.3.46.0), pbkrtest(v.0.5.2), curl(v.5.0.0), gtools(v.3.9.4), zip(v.2.2.2), GO.db(v.3.16.0), openxlsx(v.4.2.5.2), survival(v.3.5-5), limma(v.3.54.2), rmarkdown(v.2.21), desc(v.1.4.2), munsell(v.0.5.0), GenomeInfoDbData(v.1.2.9), iterators(v.1.0.14), variancePartition(v.1.28.9), reshape2(v.1.4.4), gtable(v.0.3.3) and rbibutils(v.2.2.13)
message("This is hpgltools commit: ", get_git_commit())
## If you wish to reproduce this exact build of hpgltools, invoke the following:
## > git clone http://github.com/abelew/hpgltools.git
## > git reset 367a6c7f6eb2f7364d1dc34f49e28f6a538747ec
## This is hpgltools commit: Mon Apr 17 16:38:05 2023 -0400: 367a6c7f6eb2f7364d1dc34f49e28f6a538747ec
<- paste0(gsub(pattern = "\\.Rmd", replace = "", x = rmd_file), "-v", ver, ".rda.xz")
this_save message("Saving to ", this_save)
## Saving to index-v20230420.rda.xz
<- sm(saveme(filename = this_save)) tmp