I uninstalled and reinstalled kneaddataa, and redownloaded the data, but still see the same issue. Commands run and log file are below.
bash$ cat kneaddataOutputPairedEnd/seq1_kneaddata.log
09/18/2020 10:08:37 AM - kneaddata.knead_data - INFO: Running kneaddata v0.7.10
09/18/2020 10:08:37 AM - kneaddata.knead_data - INFO: Output files will be written to: /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd
09/18/2020 10:08:37 AM - kneaddata.knead_data - DEBUG: Running with the following arguments:
verbose = False
bypass_trf = False
bmtagger_path = None
minscore = 50
bowtie2_path = /gfs/apps/bio/bowtie2-2.3.0/bowtie2
maxperiod = 500
no_discordant = False
serial = False
fastqc_start = False
bmtagger = False
cat_final_output = False
log_level = DEBUG
log = /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.log
sequencer_source = NexteraPE
max_memory = 500m
remove_intermediate_output = False
fastqc_path = None
output_dir = /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd
trf_path = /gfs/archive/jostins/microbiome/software/trf/trf
remove_temp_output = True
reference_db = /gfs/archive/jostins/microbiome/kneadata/demo/input/demo_db
input = /gfs/archive/jostins/microbiome/kneadata/demo/input/seq1.fastq /gfs/archive/jostins/microbiome/kneadata/demo/input/seq2.fastq
pi = 10
reorder = False
pm = 80
trimmomatic_path = /gfs/apps/bio/trimmomatic-0.35/trimmomatic-0.35.jar
store_temp_output = False
mismatch = 7
threads = 1
delta = 7
bowtie2_options = --very-sensitive --phred33
bypass_trim = False
processes = 1
trimmomatic_quality_scores = -phred33
fastqc_end = False
trimmomatic_options = None
output_prefix = seq1_kneaddata
match = 2
09/18/2020 10:08:37 AM - kneaddata.utilities - INFO: READ COUNT: raw pair1 : Initial number of reads ( /gfs/archive/jostins/microbiome/kneadata/demo/input/seq1.fastq ): 42473
09/18/2020 10:08:37 AM - kneaddata.utilities - INFO: READ COUNT: raw pair2 : Initial number of reads ( /gfs/archive/jostins/microbiome/kneadata/demo/input/seq2.fastq ): 42473
09/18/2020 10:08:37 AM - kneaddata.utilities - DEBUG: Checking input file to Trimmomatic : /gfs/archive/jostins/microbiome/kneadata/demo/input/seq1.fastq
09/18/2020 10:08:37 AM - kneaddata.utilities - DEBUG: Checking input file to Trimmomatic : /gfs/archive/jostins/microbiome/kneadata/demo/input/seq2.fastq
09/18/2020 10:08:37 AM - kneaddata.utilities - INFO: Running Trimmomatic …
09/18/2020 10:08:37 AM - kneaddata.utilities - INFO: Execute command: java -Xmx500m -jar /gfs/apps/bio/trimmomatic-0.35/trimmomatic-0.35.jar PE -threads 1 -phred33 /gfs/archive/jostins/microbiome/kneadata/demo/input/seq1.fastq /gfs/archive/jostins/microbiome/kneadata/demo/input/seq2.fastq /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.1.fastq /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.single.1.fastq /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.2.fastq /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.single.2.fastq ILLUMINACLIP:/gfs/devel/ljostins/python-venv-2.7.13/lib/python2.7/site-packages/kneaddata/adapters/NexteraPE-PE.fa:2:30:10:8:TRUE SLIDINGWINDOW:4:20 MINLEN:87
09/18/2020 10:08:38 AM - kneaddata.utilities - DEBUG: TrimmomaticPE: Started with arguments:
-threads 1 -phred33 /gfs/archive/jostins/microbiome/kneadata/demo/input/seq1.fastq /gfs/archive/jostins/microbiome/kneadata/demo/input/seq2.fastq /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.1.fastq /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.single.1.fastq /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.2.fastq /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.single.2.fastq ILLUMINACLIP:/gfs/devel/ljostins/python-venv-2.7.13/lib/python2.7/site-packages/kneaddata/adapters/NexteraPE-PE.fa:2:30:10:8:TRUE SLIDINGWINDOW:4:20 MINLEN:87
Using PrefixPair: ‘AGATGTGTATAAGAGACAG’ and ‘AGATGTGTATAAGAGACAG’
Using Long Clipping Sequence: ‘GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG’
Using Long Clipping Sequence: ‘TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG’
Using Long Clipping Sequence: ‘CTGTCTCTTATACACATCTCCGAGCCCACGAGAC’
Using Long Clipping Sequence: ‘CTGTCTCTTATACACATCTGACGCTGCCGACGA’
ILLUMINACLIP: Using 1 prefix pairs, 4 forward/reverse sequences, 0 forward only sequences, 0 reverse only sequences
Input Read Pairs: 42473 Both Surviving: 35341 (83.21%) Forward Only Surviving: 5385 (12.68%) Reverse Only Surviving: 847 (1.99%) Dropped: 900 (2.12%)
TrimmomaticPE: Completed successfully
09/18/2020 10:08:38 AM - kneaddata.utilities - DEBUG: Checking output file from Trimmomatic : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.1.fastq
09/18/2020 10:08:38 AM - kneaddata.utilities - DEBUG: Checking output file from Trimmomatic : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.single.1.fastq
09/18/2020 10:08:38 AM - kneaddata.utilities - DEBUG: Checking output file from Trimmomatic : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.2.fastq
09/18/2020 10:08:38 AM - kneaddata.utilities - DEBUG: Checking output file from Trimmomatic : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.single.2.fastq
09/18/2020 10:08:39 AM - kneaddata.utilities - INFO: READ COUNT: trimmed pair1 : Total reads after trimming ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.1.fastq ): 35341
09/18/2020 10:08:39 AM - kneaddata.utilities - INFO: READ COUNT: trimmed pair2 : Total reads after trimming ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.2.fastq ): 35341
09/18/2020 10:08:39 AM - kneaddata.utilities - INFO: READ COUNT: trimmed orphan1 : Total reads after trimming ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.single.1.fastq ): 5385
09/18/2020 10:08:39 AM - kneaddata.utilities - INFO: READ COUNT: trimmed orphan2 : Total reads after trimming ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.single.2.fastq ): 847
09/18/2020 10:08:40 AM - kneaddata.utilities - DEBUG: Checking input file to trf : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.1.fasta
09/18/2020 10:08:40 AM - kneaddata.utilities - INFO: Running trf …
09/18/2020 10:08:40 AM - kneaddata.utilities - INFO: Execute command: /gfs/archive/jostins/microbiome/software/trf/trf /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.1.fasta 2 7 7 80 10 50 500 -h -ngs
09/18/2020 10:08:43 AM - kneaddata.utilities - DEBUG: 0
09/18/2020 10:08:43 AM - kneaddata.utilities - DEBUG: Checking output file from trf : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.1.fasta.trf.parameters.2.7.7.80.10.50.500.dat
09/18/2020 10:08:43 AM - kneaddata.utilities - DEBUG: Checking input file to trf : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.2.fasta
09/18/2020 10:08:43 AM - kneaddata.utilities - INFO: Running trf …
09/18/2020 10:08:43 AM - kneaddata.utilities - INFO: Execute command: /gfs/archive/jostins/microbiome/software/trf/trf /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.2.fasta 2 7 7 80 10 50 500 -h -ngs
09/18/2020 10:08:46 AM - kneaddata.utilities - DEBUG: 0
09/18/2020 10:08:46 AM - kneaddata.utilities - DEBUG: Checking output file from trf : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.2.fasta.trf.parameters.2.7.7.80.10.50.500.dat
09/18/2020 10:08:47 AM - kneaddata.run - INFO: Total number of sequences with repeats removed from file ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.1.fastq ): 0
09/18/2020 10:08:47 AM - kneaddata.run - INFO: Total number of sequences with repeats removed from file ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.2.fastq ): 0
09/18/2020 10:08:47 AM - kneaddata.utilities - DEBUG: Checking input file to trf : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.single.1.fasta
09/18/2020 10:08:47 AM - kneaddata.utilities - INFO: Running trf …
09/18/2020 10:08:47 AM - kneaddata.utilities - INFO: Execute command: /gfs/archive/jostins/microbiome/software/trf/trf /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.single.1.fasta 2 7 7 80 10 50 500 -h -ngs
09/18/2020 10:08:47 AM - kneaddata.utilities - DEBUG: 0
09/18/2020 10:08:47 AM - kneaddata.utilities - DEBUG: Checking output file from trf : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.single.1.fasta.trf.parameters.2.7.7.80.10.50.500.dat
09/18/2020 10:08:47 AM - kneaddata.run - INFO: Total number of sequences with repeats removed from file ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.single.1.fastq ): 0
09/18/2020 10:08:47 AM - kneaddata.utilities - DEBUG: Checking input file to trf : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.single.2.fasta
09/18/2020 10:08:47 AM - kneaddata.utilities - INFO: Running trf …
09/18/2020 10:08:47 AM - kneaddata.utilities - INFO: Execute command: /gfs/archive/jostins/microbiome/software/trf/trf /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.single.2.fasta 2 7 7 80 10 50 500 -h -ngs
09/18/2020 10:08:48 AM - kneaddata.utilities - DEBUG: 0
09/18/2020 10:08:48 AM - kneaddata.utilities - DEBUG: Checking output file from trf : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.single.2.fasta.trf.parameters.2.7.7.80.10.50.500.dat
09/18/2020 10:08:48 AM - kneaddata.run - INFO: Total number of sequences with repeats removed from file ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.trimmed.single.2.fastq ): 1
09/18/2020 10:08:48 AM - kneaddata.run - INFO: Decontaminating …
09/18/2020 10:08:48 AM - kneaddata.utilities - DEBUG: Checking input file to bowtie2 : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.repeats.removed.1.fastq
09/18/2020 10:08:48 AM - kneaddata.utilities - DEBUG: Checking input file to bowtie2 : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.repeats.removed.2.fastq
09/18/2020 10:08:48 AM - kneaddata.utilities - DEBUG: Checking input file to bowtie2 : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.repeats.removed.unmatched.1.fastq
09/18/2020 10:08:48 AM - kneaddata.utilities - DEBUG: Checking input file to bowtie2 : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.repeats.removed.unmatched.2.fastq
09/18/2020 10:08:48 AM - kneaddata.utilities - INFO: Running bowtie2 …
09/18/2020 10:08:48 AM - kneaddata.utilities - INFO: Execute command: kneaddata_bowtie2_discordant_pairs --bowtie2 /gfs/apps/bio/bowtie2-2.3.0/bowtie2 --threads 1 -x /gfs/archive/jostins/microbiome/kneadata/demo/input/demo_db --bowtie2-options “–very-sensitive --phred33” -1 /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.repeats.removed.1.fastq -2 /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.repeats.removed.2.fastq --un-pair /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_paired_clean_%.fastq --al-pair /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_paired_contam_%.fastq -U /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.repeats.removed.unmatched.1.fastq,/gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata.repeats.removed.unmatched.2.fastq --un-single /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_unmatched_%clean.fastq --al-single /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_unmatched%_contam.fastq -S /dev/null
09/18/2020 10:08:49 AM - kneaddata.utilities - DEBUG: 35341 reads; of these:
35341 (100.00%) were unpaired; of these:
35341 (100.00%) aligned 0 times
0 (0.00%) aligned exactly 1 time
0 (0.00%) aligned >1 times
0.00% overall alignment rate
pair1_aligned : 0
pair2_aligned : 0
orphan1_unaligned : 35341
orphan2_unaligned : 0
orphan2_aligned : 0
pair2_unaligned : 0
pair1_unaligned : 0
orphan1_aligned : 0
09/18/2020 10:08:49 AM - kneaddata.utilities - DEBUG: Checking output file from bowtie2 : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_paired_clean_1.fastq
09/18/2020 10:08:49 AM - kneaddata.utilities - DEBUG: Checking output file from bowtie2 : /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_paired_clean_2.fastq
09/18/2020 10:08:49 AM - kneaddata.run - INFO: Total contaminate sequences in file ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_paired_contam_1.fastq ) : 0
09/18/2020 10:08:49 AM - kneaddata.run - INFO: Total contaminate sequences in file ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_paired_contam_2.fastq ) : 0
09/18/2020 10:08:49 AM - kneaddata.run - INFO: Total contaminate sequences in file ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_unmatched_1_contam.fastq ) : 0
09/18/2020 10:08:49 AM - kneaddata.run - INFO: Total contaminate sequences in file ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_unmatched_2_contam.fastq ) : 0
09/18/2020 10:08:49 AM - kneaddata.utilities - INFO: READ COUNT: decontaminated demo_db pair1 : Total reads after removing those found in reference database ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_paired_clean_1.fastq ): 0
09/18/2020 10:08:49 AM - kneaddata.utilities - INFO: READ COUNT: decontaminated demo_db pair2 : Total reads after removing those found in reference database ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_paired_clean_2.fastq ): 0
09/18/2020 10:08:49 AM - kneaddata.utilities - INFO: READ COUNT: final pair1 : Total reads after merging results from multiple databases ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_paired_1.fastq ): 0
09/18/2020 10:08:49 AM - kneaddata.utilities - INFO: READ COUNT: final pair2 : Total reads after merging results from multiple databases ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_paired_2.fastq ): 0
09/18/2020 10:08:49 AM - kneaddata.utilities - WARNING: Unable to remove file: /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_paired_clean_1.fastq
09/18/2020 10:08:49 AM - kneaddata.utilities - WARNING: Unable to remove file: /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_paired_clean_2.fastq
09/18/2020 10:08:49 AM - kneaddata.utilities - INFO: READ COUNT: decontaminated demo_db orphan1 : Total reads after removing those found in reference database ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_unmatched_1_clean.fastq ): 35341
09/18/2020 10:08:49 AM - kneaddata.utilities - INFO: READ COUNT: final orphan1 : Total reads after merging results from multiple databases ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_unmatched_1.fastq ): 35341
09/18/2020 10:08:49 AM - kneaddata.utilities - WARNING: Unable to remove file: /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_unmatched_1_clean.fastq
09/18/2020 10:08:49 AM - kneaddata.utilities - INFO: READ COUNT: decontaminated demo_db orphan2 : Total reads after removing those found in reference database ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_unmatched_2_clean.fastq ): 0
09/18/2020 10:08:49 AM - kneaddata.utilities - INFO: READ COUNT: final orphan2 : Total reads after merging results from multiple databases ( /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_unmatched_2.fastq ): 0
09/18/2020 10:08:49 AM - kneaddata.utilities - WARNING: Unable to remove file: /gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_demo_db_bowtie2_unmatched_2_clean.fastq
09/18/2020 10:08:49 AM - kneaddata.knead_data - INFO:
Final output files created:
/gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_paired_1.fastq
/gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_paired_2.fastq
/gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_unmatched_1.fastq
/gfs/archive/jostins/microbiome/kneadata/demo/kneaddataOutputPairedEnd/seq1_kneaddata_unmatched_2.fastq