03/21/2018 03:18:22 PM - kneaddata.knead_data - INFO: Running kneaddata v0.7.0 03/21/2018 03:18:22 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:22 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/HSMA33MG.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/HSMA33MG_R1.fastq.gz /n/regal/huttenhower_lab/carze/data/hmp2/workflow/deposition/HMP2/2018-03-20/MTX/HSMA33MG_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 = HSMA33MG match = 2 03/21/2018 03:18:22 PM - kneaddata.utilities - INFO: Decompressing gzipped file ... 03/21/2018 03:19:25 PM - kneaddata.utilities - INFO: Decompressed file created: /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/decompressed_wlVumu_HSMA33MG_R1.fastq 03/21/2018 03:19:25 PM - kneaddata.utilities - INFO: Decompressing gzipped file ... 03/21/2018 03:22:46 PM - kneaddata.utilities - INFO: Decompressed file created: /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/decompressed_jKadqr_HSMA33MG_R2.fastq 03/21/2018 03:22:47 PM - kneaddata.utilities - INFO: Reformatting file sequence identifiers ... 03/21/2018 03:25:33 PM - kneaddata.utilities - INFO: Reformatting file sequence identifiers ... 03/21/2018 03:31:56 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_identifiersgL3Tm8_decompressed_jKadqr_HSMA33MG_R2 ): 27348378 03/21/2018 03:33:44 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_identifiersl3Cu7G_decompressed_wlVumu_HSMA33MG_R1 ): 27348378 03/21/2018 03:33:44 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_identifiersl3Cu7G_decompressed_wlVumu_HSMA33MG_R1 03/21/2018 03:33:44 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_identifiersgL3Tm8_decompressed_jKadqr_HSMA33MG_R2 03/21/2018 03:33:44 PM - kneaddata.utilities - INFO: Running Trimmomatic ... 03/21/2018 03:33:44 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_identifiersl3Cu7G_decompressed_wlVumu_HSMA33MG_R1 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/reformatted_identifiersgL3Tm8_decompressed_jKadqr_HSMA33MG_R2 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG.trimmed.1.fastq /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG.trimmed.single.1.fastq /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG.trimmed.2.fastq /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG.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 03:45:49 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_identifiersl3Cu7G_decompressed_wlVumu_HSMA33MG_R1 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/reformatted_identifiersgL3Tm8_decompressed_jKadqr_HSMA33MG_R2 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG.trimmed.1.fastq /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG.trimmed.single.1.fastq /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG.trimmed.2.fastq /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG.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: 27348378 Both Surviving: 12472024 (45.60%) Forward Only Surviving: 8950914 (32.73%) Reverse Only Surviving: 1162685 (4.25%) Dropped: 4762755 (17.42%) TrimmomaticPE: Completed successfully 03/21/2018 03:45:49 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/HSMA33MG.trimmed.1.fastq 03/21/2018 03:45:49 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/HSMA33MG.trimmed.single.1.fastq 03/21/2018 03:45:49 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/HSMA33MG.trimmed.2.fastq 03/21/2018 03:45:49 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/HSMA33MG.trimmed.single.2.fastq 03/21/2018 03:46:20 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/HSMA33MG.trimmed.1.fastq ): 12472024 03/21/2018 03:46:48 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/HSMA33MG.trimmed.2.fastq ): 12472024 03/21/2018 03:49:06 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/HSMA33MG.trimmed.single.1.fastq ): 8950914 03/21/2018 03:49:08 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/HSMA33MG.trimmed.single.2.fastq ): 1162685 03/21/2018 03:49:08 PM - kneaddata.run - INFO: Decontaminating ... 03/21/2018 03:49:13 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/HSMA33MG.trimmed.1.fastq 03/21/2018 03:49:13 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/HSMA33MG.trimmed.2.fastq 03/21/2018 03:49:13 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/HSMA33MG.trimmed.single.1.fastq 03/21/2018 03:49:13 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/HSMA33MG.trimmed.single.2.fastq 03/21/2018 03:49:13 PM - kneaddata.utilities - INFO: Running bowtie2 ... 03/21/2018 03:49:13 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/HSMA33MG.trimmed.1.fastq -2 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG.trimmed.2.fastq --un-pair /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG_Homo_sapiens_bowtie2_paired_clean_%.fastq --al-pair /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG_Homo_sapiens_bowtie2_paired_contam_%.fastq -U /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG.trimmed.single.1.fastq,/n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG.trimmed.single.2.fastq --un-single /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG_Homo_sapiens_bowtie2_unmatched_%_clean.fastq --al-single /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG_Homo_sapiens_bowtie2_unmatched_%_contam.fastq -S /dev/null 03/21/2018 05:53:26 PM - kneaddata.utilities - DEBUG: 35057647 reads; of these: 35057647 (100.00%) were unpaired; of these: 18446990 (52.62%) aligned 0 times 4538779 (12.95%) aligned exactly 1 time 12071878 (34.43%) aligned >1 times 47.38% overall alignment rate pair1_aligned : 5289145 pair2_aligned : 5289145 orphan1_unaligned : 4301929 orphan2_unaligned : 758819 orphan2_aligned : 893624 pair2_unaligned : 6693121 pair1_unaligned : 6693121 orphan1_aligned : 5138743 03/21/2018 05:53:26 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/HSMA33MG_Homo_sapiens_bowtie2_paired_clean_1.fastq 03/21/2018 05:53:26 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/HSMA33MG_Homo_sapiens_bowtie2_paired_clean_2.fastq 03/21/2018 05:53:26 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/HSMA33MG.trimmed.1.fastq 03/21/2018 05:53:26 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/HSMA33MG.trimmed.2.fastq 03/21/2018 05:53:26 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/HSMA33MG.trimmed.single.1.fastq 03/21/2018 05:53:26 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/HSMA33MG.trimmed.single.2.fastq 03/21/2018 05:53:26 PM - kneaddata.utilities - INFO: Running bowtie2 ... 03/21/2018 05:53:26 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/HSMA33MG_Homo_sapiens_bowtie2_paired_clean_1.fastq -2 /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG_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/HSMA33MG_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/HSMA33MG_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/HSMA33MG_Homo_sapiens_bowtie2_unmatched_1_clean.fastq,/n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG_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/HSMA33MG_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/HSMA33MG_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_unmatched_%_contam.fastq -S /dev/null 03/21/2018 06:09:08 PM - kneaddata.utilities - DEBUG: 18446990 reads; of these: 18446990 (100.00%) were unpaired; of these: 17989074 (97.52%) aligned 0 times 22608 (0.12%) aligned exactly 1 time 435308 (2.36%) aligned >1 times 2.48% overall alignment rate pair1_aligned : 14569 pair2_aligned : 14569 orphan1_unaligned : 4079348 orphan2_unaligned : 606060 orphan2_aligned : 179478 pair2_unaligned : 6651833 pair1_unaligned : 6651833 orphan1_aligned : 249300 03/21/2018 06:09:08 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/HSMA33MG_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_paired_clean_1.fastq 03/21/2018 06:09:08 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/HSMA33MG_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_paired_clean_2.fastq 03/21/2018 06:09:08 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/HSMA33MG.trimmed.1.fastq 03/21/2018 06:09:08 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/HSMA33MG.trimmed.2.fastq 03/21/2018 06:09:08 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/HSMA33MG.trimmed.single.1.fastq 03/21/2018 06:09:08 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/HSMA33MG.trimmed.single.2.fastq 03/21/2018 06:09:08 PM - kneaddata.utilities - INFO: Running bowtie2 ... 03/21/2018 06:09:08 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/HSMA33MG_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/HSMA33MG_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/HSMA33MG_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/HSMA33MG_human_hg38_refMrna_bowtie2_paired_contam_%.fastq -U /n/regal/huttenhower_lab/carze/data/hmp2/workflow/processing/HMP2/2018-03-20/MTX/kneaddata/main/HSMA33MG_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/HSMA33MG_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/HSMA33MG_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/HSMA33MG_human_hg38_refMrna_bowtie2_unmatched_%_contam.fastq -S /dev/null 03/21/2018 06:22:06 PM - kneaddata.utilities - DEBUG: 17989074 reads; of these: 17989074 (100.00%) were unpaired; of these: 17655718 (98.15%) aligned 0 times 110804 (0.62%) aligned exactly 1 time 222552 (1.24%) aligned >1 times 1.85% overall alignment rate pair1_aligned : 27513 pair2_aligned : 27513 orphan1_unaligned : 3973565 orphan2_unaligned : 440497 orphan2_aligned : 169055 pair2_unaligned : 6620828 pair1_unaligned : 6620828 orphan1_aligned : 109275 03/21/2018 06:22:06 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/HSMA33MG_human_hg38_refMrna_bowtie2_paired_clean_1.fastq 03/21/2018 06:22:06 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/HSMA33MG_human_hg38_refMrna_bowtie2_paired_clean_2.fastq 03/21/2018 06:23:01 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/HSMA33MG_Homo_sapiens_bowtie2_paired_contam_1.fastq ) : 5289145 03/21/2018 06:23:07 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/HSMA33MG_Homo_sapiens_bowtie2_paired_contam_2.fastq ) : 5289145 03/21/2018 06:23:24 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/HSMA33MG_Homo_sapiens_bowtie2_unmatched_1_contam.fastq ) : 5138743 03/21/2018 06:23:25 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/HSMA33MG_Homo_sapiens_bowtie2_unmatched_2_contam.fastq ) : 893624 03/21/2018 06:23:25 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/HSMA33MG_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_paired_contam_1.fastq ) : 14569 03/21/2018 06:23:26 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/HSMA33MG_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_paired_contam_2.fastq ) : 14569 03/21/2018 06:23:27 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/HSMA33MG_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_unmatched_1_contam.fastq ) : 249300 03/21/2018 06:23:28 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/HSMA33MG_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_unmatched_2_contam.fastq ) : 179478 03/21/2018 06:23:28 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/HSMA33MG_human_hg38_refMrna_bowtie2_paired_contam_1.fastq ) : 27513 03/21/2018 06:23:28 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/HSMA33MG_human_hg38_refMrna_bowtie2_paired_contam_2.fastq ) : 27513 03/21/2018 06:23:28 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/HSMA33MG_human_hg38_refMrna_bowtie2_unmatched_1_contam.fastq ) : 109275 03/21/2018 06:23:28 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/HSMA33MG_human_hg38_refMrna_bowtie2_unmatched_2_contam.fastq ) : 169055 03/21/2018 06:23:40 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/HSMA33MG_Homo_sapiens_bowtie2_paired_clean_1.fastq ): 6693121 03/21/2018 06:24:25 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/HSMA33MG_Homo_sapiens_bowtie2_paired_clean_2.fastq ): 6693121 03/21/2018 06:25:59 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/HSMA33MG_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_paired_clean_1.fastq ): 6651833 03/21/2018 06:26: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/HSMA33MG_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_paired_clean_2.fastq ): 6651833 03/21/2018 06:26:24 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/HSMA33MG_human_hg38_refMrna_bowtie2_paired_clean_1.fastq ): 6620828 03/21/2018 06:26:33 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/HSMA33MG_human_hg38_refMrna_bowtie2_paired_clean_2.fastq ): 6620828 03/21/2018 06:29:44 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/HSMA33MG_paired_1.fastq ): 6620828 03/21/2018 06:29:48 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/HSMA33MG_paired_2.fastq ): 6620828 03/21/2018 06:30:02 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/HSMA33MG_Homo_sapiens_bowtie2_unmatched_1_clean.fastq ): 4301929 03/21/2018 06:30:09 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/HSMA33MG_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_unmatched_1_clean.fastq ): 4079348 03/21/2018 06:30:23 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/HSMA33MG_human_hg38_refMrna_bowtie2_unmatched_1_clean.fastq ): 3973565 03/21/2018 06:31:00 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/HSMA33MG_unmatched_1.fastq ): 3965073 03/21/2018 06:31:03 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/HSMA33MG_Homo_sapiens_bowtie2_unmatched_2_clean.fastq ): 758819 03/21/2018 06:31:18 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/HSMA33MG_SILVA_128_LSUParc_SSUParc_ribosomal_RNA_bowtie2_unmatched_2_clean.fastq ): 606060 03/21/2018 06:31:25 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/HSMA33MG_human_hg38_refMrna_bowtie2_unmatched_2_clean.fastq ): 440497 03/21/2018 06:31:31 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/HSMA33MG_unmatched_2.fastq ): 425383 03/21/2018 06:31:54 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/HSMA33MG.fastq