03/21/2018 03:18:00 PM - kneaddata.knead_data - INFO: Running kneaddata v0.7.0 03/21/2018 03:18:00 PM - kneaddata.knead_data - INFO: Output files will be written to: /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main 03/21/2018 03:18:00 PM - kneaddata.knead_data - DEBUG: Running with the following arguments: verbose = False bmtagger_path = None minscore = 50 bowtie2_path = /n/sw/centos6/bowtie2-2.2.1/bowtie2 maxperiod = 500 no_discordant = False serial = True fastqc_start = False bmtagger = False cat_final_output = True log_level = DEBUG log = /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.log max_memory = 500m remove_intermediate_output = True fastqc_path = None output_dir = /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main trf_path = None remove_temp_output = True reference_db = /n/huttenhower_lab/data/kneaddata_databases/Homo_sapiens /n/huttenhower_lab/data/kneaddata_databases/ribosomal_RNA/SILVA_128_LSUParc_SSUParc_ribosomal_RNA /n/huttenhower_lab/data/kneaddata_databases/transcriptome/human_hg38_refMrna input = /n/regal/huttenhower_lab/carze/data/hmp2/workflow/deposition/HMP2/2018-03-20/MTX/HSMA33ME_R1.fastq.gz /n/regal/huttenhower_lab/carze/data/hmp2/workflow/deposition/HMP2/2018-03-20/MTX/HSMA33ME_R2.fastq.gz pi = 10 reorder = False pm = 80 trimmomatic_path = /n/huttenhower_lab/tools/kneaddata/bin_v0.7.0_devel/trimmomatic-0.33.jar store_temp_output = False mismatch = 7 threads = 4 delta = 7 bowtie2_options = --very-sensitive --phred33 bypass_trim = False processes = 1 trimmomatic_quality_scores = -phred33 fastqc_end = False trf = False trimmomatic_options = ILLUMINACLIP:/n/home07/carze/.local/bin/Trimmomatic-0.36/adapters/TruSeq3-PE-2.fa:2:30:10:8:TRUE SLIDINGWINDOW:4:20 MINLEN:50 output_prefix = HSMA33ME match = 2 03/21/2018 03:18:00 PM - kneaddata.utilities - INFO: Decompressing gzipped file ... 03/21/2018 03:21:31 PM - kneaddata.utilities - INFO: Decompressed file created: /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/decompressed_aQoOcn_HSMA33ME_R1.fastq 03/21/2018 03:21:31 PM - kneaddata.utilities - INFO: Decompressing gzipped file ... 03/21/2018 03:35:41 PM - kneaddata.utilities - INFO: Decompressed file created: /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/decompressed_8RnxWR_HSMA33ME_R2.fastq 03/21/2018 03:35:42 PM - kneaddata.utilities - INFO: Reformatting file sequence identifiers ... 03/21/2018 03:57:16 PM - kneaddata.utilities - INFO: Reformatting file sequence identifiers ... 03/21/2018 04:13:46 PM - kneaddata.utilities - INFO: READ COUNT: raw pair1 : Initial number of reads ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/reformatted_identifiersZ6i4CD_decompressed_aQoOcn_HSMA33ME_R1 ): 26787390 03/21/2018 04:15:07 PM - kneaddata.utilities - INFO: READ COUNT: raw pair2 : Initial number of reads ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/reformatted_identifiersjwyu78_decompressed_8RnxWR_HSMA33ME_R2 ): 26787390 03/21/2018 04:15:07 PM - kneaddata.utilities - DEBUG: Checking input file to Trimmomatic : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/reformatted_identifiersZ6i4CD_decompressed_aQoOcn_HSMA33ME_R1 03/21/2018 04:15:07 PM - kneaddata.utilities - DEBUG: Checking input file to Trimmomatic : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/reformatted_identifiersjwyu78_decompressed_8RnxWR_HSMA33ME_R2 03/21/2018 04:15:07 PM - kneaddata.utilities - INFO: Running Trimmomatic ... 03/21/2018 04:15:07 PM - kneaddata.utilities - INFO: Execute command: java -Xmx500m -d64 -jar /n/huttenhower_lab/tools/kneaddata/bin_v0.7.0_devel/trimmomatic-0.33.jar PE -threads 4 -phred33 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/reformatted_identifiersZ6i4CD_decompressed_aQoOcn_HSMA33ME_R1 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/reformatted_identifiersjwyu78_decompressed_8RnxWR_HSMA33ME_R2 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.1.fastq /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.single.1.fastq /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.2.fastq /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.single.2.fastq ILLUMINACLIP:/n/home07/carze/.local/bin/Trimmomatic-0.36/adapters/TruSeq3-PE-2.fa:2:30:10:8:TRUE SLIDINGWINDOW:4:20 MINLEN:50 03/21/2018 04:26:03 PM - kneaddata.utilities - DEBUG: TrimmomaticPE: Started with arguments: -threads 4 -phred33 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/reformatted_identifiersZ6i4CD_decompressed_aQoOcn_HSMA33ME_R1 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/reformatted_identifiersjwyu78_decompressed_8RnxWR_HSMA33ME_R2 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.1.fastq /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.single.1.fastq /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.2.fastq /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.single.2.fastq ILLUMINACLIP:/n/home07/carze/.local/bin/Trimmomatic-0.36/adapters/TruSeq3-PE-2.fa:2:30:10:8:TRUE SLIDINGWINDOW:4:20 MINLEN:50 Using PrefixPair: 'TACACTCTTTCCCTACACGACGCTCTTCCGATCT' and 'GTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT' Using Long Clipping Sequence: 'AGATCGGAAGAGCACACGTCTGAACTCCAGTCAC' Using Long Clipping Sequence: 'TACACTCTTTCCCTACACGACGCTCTTCCGATCT' Using Long Clipping Sequence: 'GTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT' Using Long Clipping Sequence: 'AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGTA' ILLUMINACLIP: Using 1 prefix pairs, 4 forward/reverse sequences, 0 forward only sequences, 0 reverse only sequences Input Read Pairs: 26787390 Both Surviving: 12804335 (47.80%) Forward Only Surviving: 9067880 (33.85%) Reverse Only Surviving: 830516 (3.10%) Dropped: 4084659 (15.25%) TrimmomaticPE: Completed successfully 03/21/2018 04:26:04 PM - kneaddata.utilities - DEBUG: Checking output file from Trimmomatic : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.1.fastq 03/21/2018 04:26:04 PM - kneaddata.utilities - DEBUG: Checking output file from Trimmomatic : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.single.1.fastq 03/21/2018 04:26:04 PM - kneaddata.utilities - DEBUG: Checking output file from Trimmomatic : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.2.fastq 03/21/2018 04:26:04 PM - kneaddata.utilities - DEBUG: Checking output file from Trimmomatic : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.single.2.fastq 03/21/2018 04:26:25 PM - kneaddata.utilities - INFO: READ COUNT: trimmed pair1 : Total reads after trimming ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.1.fastq ): 12804335 03/21/2018 04:27:53 PM - kneaddata.utilities - INFO: READ COUNT: trimmed pair2 : Total reads after trimming ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.2.fastq ): 12804335 03/21/2018 04:28:55 PM - kneaddata.utilities - INFO: READ COUNT: trimmed orphan1 : Total reads after trimming ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.single.1.fastq ): 9067880 03/21/2018 04:28:59 PM - kneaddata.utilities - INFO: READ COUNT: trimmed orphan2 : Total reads after trimming ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.single.2.fastq ): 830516 03/21/2018 04:28:59 PM - kneaddata.run - INFO: Decontaminating ... 03/21/2018 04:28:59 PM - kneaddata.utilities - DEBUG: Checking input file to bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.1.fastq 03/21/2018 04:28:59 PM - kneaddata.utilities - DEBUG: Checking input file to bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.2.fastq 03/21/2018 04:28:59 PM - kneaddata.utilities - DEBUG: Checking input file to bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.single.1.fastq 03/21/2018 04:28:59 PM - kneaddata.utilities - DEBUG: Checking input file to bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.single.2.fastq 03/21/2018 04:28:59 PM - kneaddata.utilities - INFO: Running bowtie2 ... 03/21/2018 04:28:59 PM - kneaddata.utilities - INFO: Execute command: kneaddata_bowtie2_discordant_pairs --bowtie2 /n/sw/centos6/bowtie2-2.2.1/bowtie2 --threads 4 -x /n/huttenhower_lab/data/kneaddata_databases/Homo_sapiens --bowtie2-options "--very-sensitive --phred33" -1 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.1.fastq -2 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.2.fastq --un-pair /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_paired_clean_%.fastq --al-pair /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_paired_contam_%.fastq -U /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.single.1.fastq,/n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.single.2.fastq --un-single /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_unmatched_%_clean.fastq --al-single /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_unmatched_%_contam.fastq -S /dev/null 03/21/2018 06:20:43 PM - kneaddata.utilities - DEBUG: 35507066 reads; of these: 35507066 (100.00%) were unpaired; of these: 18137810 (51.08%) aligned 0 times 5830513 (16.42%) aligned exactly 1 time 11538743 (32.50%) aligned >1 times 48.92% overall alignment rate pair1_aligned : 5574286 pair2_aligned : 5574286 orphan1_unaligned : 3900831 orphan2_unaligned : 719803 orphan2_aligned : 582174 pair2_unaligned : 6758588 pair1_unaligned : 6758588 orphan1_aligned : 5638510 03/21/2018 06:20:43 PM - kneaddata.utilities - DEBUG: Checking output file from bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_paired_clean_1.fastq 03/21/2018 06:20:43 PM - kneaddata.utilities - DEBUG: Checking output file from bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_paired_clean_2.fastq 03/21/2018 06:20:43 PM - kneaddata.utilities - DEBUG: Checking input file to bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.1.fastq 03/21/2018 06:20:43 PM - kneaddata.utilities - DEBUG: Checking input file to bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.2.fastq 03/21/2018 06:20:43 PM - kneaddata.utilities - DEBUG: Checking input file to bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.single.1.fastq 03/21/2018 06:20:43 PM - kneaddata.utilities - DEBUG: Checking input file to bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.single.2.fastq 03/21/2018 06:20:43 PM - kneaddata.utilities - INFO: Running bowtie2 ... 03/21/2018 06:20:43 PM - kneaddata.utilities - INFO: Execute command: kneaddata_bowtie2_discordant_pairs --bowtie2 /n/sw/centos6/bowtie2-2.2.1/bowtie2 --threads 4 -x /n/huttenhower_lab/data/kneaddata_databases/ribosomal_RNA/SILVA_128_LSUParc_SSUParc_ribosomal_RNA --bowtie2-options "--very-sensitive --phred33" -1 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_paired_clean_1.fastq -2 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_paired_clean_2.fastq --un-pair /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_paired_clean_%.fastq --al-pair /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_paired_contam_%.fastq -U /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_unmatched_1_clean.fastq,/n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_unmatched_2_clean.fastq --un-single /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_unmatched_%_clean.fastq --al-single /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_unmatched_%_contam.fastq -S /dev/null 03/21/2018 06:37:42 PM - kneaddata.utilities - DEBUG: 18137810 reads; of these: 18137810 (100.00%) were unpaired; of these: 17667212 (97.41%) aligned 0 times 16748 (0.09%) aligned exactly 1 time 453850 (2.50%) aligned >1 times 2.59% overall alignment rate pair1_aligned : 10410 pair2_aligned : 10410 orphan1_unaligned : 3592386 orphan2_unaligned : 624210 orphan2_aligned : 118463 pair2_unaligned : 6725308 pair1_unaligned : 6725308 orphan1_aligned : 331315 03/21/2018 06:37:44 PM - kneaddata.utilities - DEBUG: Checking output file from bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_paired_clean_1.fastq 03/21/2018 06:37:44 PM - kneaddata.utilities - DEBUG: Checking output file from bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_paired_clean_2.fastq 03/21/2018 06:37:44 PM - kneaddata.utilities - DEBUG: Checking input file to bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.1.fastq 03/21/2018 06:37:44 PM - kneaddata.utilities - DEBUG: Checking input file to bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.2.fastq 03/21/2018 06:37:44 PM - kneaddata.utilities - DEBUG: Checking input file to bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.single.1.fastq 03/21/2018 06:37:44 PM - kneaddata.utilities - DEBUG: Checking input file to bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.trimmed.single.2.fastq 03/21/2018 06:37:44 PM - kneaddata.utilities - INFO: Running bowtie2 ... 03/21/2018 06:37:44 PM - kneaddata.utilities - INFO: Execute command: kneaddata_bowtie2_discordant_pairs --bowtie2 /n/sw/centos6/bowtie2-2.2.1/bowtie2 --threads 4 -x /n/huttenhower_lab/data/kneaddata_databases/transcriptome/human_hg38_refMrna --bowtie2-options "--very-sensitive --phred33" -1 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_paired_clean_1.fastq -2 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_paired_clean_2.fastq --un-pair /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_human_hg38_refMrna_bowtie2_paired_clean_%.fastq --al-pair /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_human_hg38_refMrna_bowtie2_paired_contam_%.fastq -U /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_unmatched_1_clean.fastq,/n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_unmatched_2_clean.fastq --un-single /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_human_hg38_refMrna_bowtie2_unmatched_%_clean.fastq --al-single /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_human_hg38_refMrna_bowtie2_unmatched_%_contam.fastq -S /dev/null 03/21/2018 06:51:20 PM - kneaddata.utilities - DEBUG: 17667212 reads; of these: 17667212 (100.00%) were unpaired; of these: 17311155 (97.98%) aligned 0 times 131567 (0.74%) aligned exactly 1 time 224490 (1.27%) aligned >1 times 2.02% overall alignment rate pair1_aligned : 36291 pair2_aligned : 36291 orphan1_unaligned : 3455309 orphan2_unaligned : 486520 orphan2_aligned : 142044 pair2_unaligned : 6684663 pair1_unaligned : 6684663 orphan1_aligned : 141431 03/21/2018 06:51:20 PM - kneaddata.utilities - DEBUG: Checking output file from bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_human_hg38_refMrna_bowtie2_paired_clean_1.fastq 03/21/2018 06:51:20 PM - kneaddata.utilities - DEBUG: Checking output file from bowtie2 : /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_human_hg38_refMrna_bowtie2_paired_clean_2.fastq 03/21/2018 06:51:29 PM - kneaddata.run - INFO: Total contaminate sequences in file ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_paired_contam_1.fastq ) : 5574286 03/21/2018 06:51:40 PM - kneaddata.run - INFO: Total contaminate sequences in file ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_paired_contam_2.fastq ) : 5574286 03/21/2018 06:51:48 PM - kneaddata.run - INFO: Total contaminate sequences in file ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_unmatched_1_contam.fastq ) : 5638510 03/21/2018 06:51:54 PM - kneaddata.run - INFO: Total contaminate sequences in file ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_unmatched_2_contam.fastq ) : 582174 03/21/2018 06:51:54 PM - kneaddata.run - INFO: Total contaminate sequences in file ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_paired_contam_1.fastq ) : 10410 03/21/2018 06:51:54 PM - kneaddata.run - INFO: Total contaminate sequences in file ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_paired_contam_2.fastq ) : 10410 03/21/2018 06:51:56 PM - kneaddata.run - INFO: Total contaminate sequences in file ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_unmatched_1_contam.fastq ) : 331315 03/21/2018 06:51:56 PM - kneaddata.run - INFO: Total contaminate sequences in file ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_unmatched_2_contam.fastq ) : 118463 03/21/2018 06:51:57 PM - kneaddata.run - INFO: Total contaminate sequences in file ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_human_hg38_refMrna_bowtie2_paired_contam_1.fastq ) : 36291 03/21/2018 06:51:57 PM - kneaddata.run - INFO: Total contaminate sequences in file ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_human_hg38_refMrna_bowtie2_paired_contam_2.fastq ) : 36291 03/21/2018 06:51:57 PM - kneaddata.run - INFO: Total contaminate sequences in file ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_human_hg38_refMrna_bowtie2_unmatched_1_contam.fastq ) : 141431 03/21/2018 06:51:57 PM - kneaddata.run - INFO: Total contaminate sequences in file ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_human_hg38_refMrna_bowtie2_unmatched_2_contam.fastq ) : 142044 03/21/2018 06:52:09 PM - kneaddata.utilities - INFO: READ COUNT: decontaminated Homo_sapiens pair1 : Total reads after removing those found in reference database ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_paired_clean_1.fastq ): 6758588 03/21/2018 06:53:23 PM - kneaddata.utilities - INFO: READ COUNT: decontaminated Homo_sapiens pair2 : Total reads after removing those found in reference database ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_paired_clean_2.fastq ): 6758588 03/21/2018 06:53:55 PM - kneaddata.utilities - INFO: READ COUNT: decontaminated SILVA_128_LSUParc_SSUParc_ribosomal_RNA pair1 : Total reads after removing those found in reference database ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_paired_clean_1.fastq ): 6725308 03/21/2018 06:54:11 PM - kneaddata.utilities - INFO: READ COUNT: decontaminated SILVA_128_LSUParc_SSUParc_ribosomal_RNA pair2 : Total reads after removing those found in reference database ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_paired_clean_2.fastq ): 6725308 03/21/2018 06:55:28 PM - kneaddata.utilities - INFO: READ COUNT: decontaminated human_hg38_refMrna pair1 : Total reads after removing those found in reference database ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_human_hg38_refMrna_bowtie2_paired_clean_1.fastq ): 6684663 03/21/2018 06:55:47 PM - kneaddata.utilities - INFO: READ COUNT: decontaminated human_hg38_refMrna pair2 : Total reads after removing those found in reference database ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_human_hg38_refMrna_bowtie2_paired_clean_2.fastq ): 6684663 03/21/2018 07:01:38 PM - kneaddata.utilities - INFO: READ COUNT: final pair1 : Total reads after merging results from multiple databases ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_paired_1.fastq ): 6684663 03/21/2018 07:01:41 PM - kneaddata.utilities - INFO: READ COUNT: final pair2 : Total reads after merging results from multiple databases ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_paired_2.fastq ): 6684663 03/21/2018 07:02:03 PM - kneaddata.utilities - INFO: READ COUNT: decontaminated Homo_sapiens orphan1 : Total reads after removing those found in reference database ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_unmatched_1_clean.fastq ): 3900831 03/21/2018 07:02:18 PM - kneaddata.utilities - INFO: READ COUNT: decontaminated SILVA_128_LSUParc_SSUParc_ribosomal_RNA orphan1 : Total reads after removing those found in reference database ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_unmatched_1_clean.fastq ): 3592386 03/21/2018 07:03:06 PM - kneaddata.utilities - INFO: READ COUNT: decontaminated human_hg38_refMrna orphan1 : Total reads after removing those found in reference database ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_human_hg38_refMrna_bowtie2_unmatched_1_clean.fastq ): 3455309 03/21/2018 07:03:37 PM - kneaddata.utilities - INFO: READ COUNT: final orphan1 : Total reads after merging results from multiple databases ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_unmatched_1.fastq ): 3446158 03/21/2018 07:03:45 PM - kneaddata.utilities - INFO: READ COUNT: decontaminated Homo_sapiens orphan2 : Total reads after removing those found in reference database ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_Homo_sapiens_bowtie2_unmatched_2_clean.fastq ): 719803 03/21/2018 07:03:45 PM - kneaddata.utilities - INFO: READ COUNT: decontaminated SILVA_128_LSUParc_SSUParc_ribosomal_RNA orphan2 : Total reads after removing those found in reference database ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_unmatched_2_clean.fastq ): 624210 03/21/2018 07:03:46 PM - kneaddata.utilities - INFO: READ COUNT: decontaminated human_hg38_refMrna orphan2 : Total reads after removing those found in reference database ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_human_hg38_refMrna_bowtie2_unmatched_2_clean.fastq ): 486520 03/21/2018 07:03:51 PM - kneaddata.utilities - INFO: READ COUNT: final orphan2 : Total reads after merging results from multiple databases ( /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME_unmatched_2.fastq ): 473180 03/21/2018 07:04:23 PM - kneaddata.knead_data - INFO: Final output file created: /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33ME.fastq