Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- Path to Bowtie 2 specified as: bowtie2
- Output format is BAM (default)
- Alignments will be written out in BAM format. Samtools found here: '/home/enrico16/bin/samtools'
- Reference genome folder provided is /home/enrico16/analysis/00genome/ (absolute path is '/home/enrico16/analysis/00genome/)'
- FastQ format assumed (by default)
- Each Bowtie 2 instance is going to be run with 2 threads. Please monitor performance closely and tune down if needed!
- Library is assumed to be strand-specific (directional), alignments to strands complementary to the original top or bottom strands will be ignored (i.e. not performed!)
- Output will be written into the directory: /home/enrico16/analysis/05alignment/SampleX/
- Using temp directory: /home/enrico16/analysis/05alignment/SampleX
- Temporary files will be written into the directory: /home/enrico16/analysis/05alignment/SampleX/
- Running Bismark Parallel version. Number of parallel instances to be spawned: 2
- Paired-end alignments will be performed
- =======================================
- The provided filenames for paired-end alignments are /home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz and /home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair2.fastq.gz
- Paired-end alignments will be performed
- =======================================
- The provided filenames for paired-end alignments are /home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz and /home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair2.fastq.gz
- Finished subdividing /home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz for PID: 8809 and offset 2
- Using the subset file >/home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.2< as new in-file 1 (instead of >/home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz<)
- Finished subdividing /home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz for PID: 0 and offset 1
- Using the subset file >/home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.1< as new in-file 1 (instead of >/home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz<)
- Finished subdividing /home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair2.fastq.gz for PID: 8809 and offset 2
- Using the subset file >/home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair2.fastq.gz.temp.2< as new in-file 2 (instead of >/home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair2.fastq.gz<)
- Input files are in FastQ format
- Finished subdividing /home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair2.fastq.gz for PID: 0 and offset 1
- Using the subset file >/home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair2.fastq.gz.temp.1< as new in-file 2 (instead of >/home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair2.fastq.gz<)
- Input files are in FastQ format
- Writing a C -> T converted version of the input file SampleX-trimmed-pair1.fastq.gz.temp.2 to /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.2_C_to_T.fastq
- Writing a C -> T converted version of the input file SampleX-trimmed-pair1.fastq.gz.temp.1 to /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.1_C_to_T.fastq
- Created C -> T converted version of the FastQ file SampleX-trimmed-pair1.fastq.gz.temp.2 (45829608 sequences in total)
- Created C -> T converted version of the FastQ file SampleX-trimmed-pair1.fastq.gz.temp.1 (45829609 sequences in total)
- Writing a G -> A converted version of the input file SampleX-trimmed-pair2.fastq.gz.temp.2 to /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair2.fastq.gz.temp.2_G_to_A.fastq
- Writing a G -> A converted version of the input file SampleX-trimmed-pair2.fastq.gz.temp.1 to /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair2.fastq.gz.temp.1_G_to_A.fastq
- Created G -> A converted version of the FastQ file SampleX-trimmed-pair2.fastq.gz.temp.2 (45829608 sequences in total)
- Input files are SampleX-trimmed-pair1.fastq.gz.temp.2_C_to_T.fastq and SampleX-trimmed-pair2.fastq.gz.temp.2_G_to_A.fastq (FastQ)
- Now running 2 instances of Bowtie 2 against the bisulfite genome of /home/enrico16/analysis/00genome/ with the specified options: -q --score-min L,0,-0.2 -p 2 --reorder --ignore-quals --no-mixed --no-discordant --maxins 500
- Now starting a Bowtie 2 paired-end alignment for CTread1GAread2CTgenome (reading in sequences from /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.2_C_to_T.fastq and /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair2.fastq.gz.temp.2_G_to_A.fastq, with the options: -q --score-min L,0,-0.2 -p 2 --reorder --ignore-quals --no-mixed --no-discordant --maxins 500 --norc))
- Created G -> A converted version of the FastQ file SampleX-trimmed-pair2.fastq.gz.temp.1 (45829609 sequences in total)
- Input files are SampleX-trimmed-pair1.fastq.gz.temp.1_C_to_T.fastq and SampleX-trimmed-pair2.fastq.gz.temp.1_G_to_A.fastq (FastQ)
- Now running 2 instances of Bowtie 2 against the bisulfite genome of /home/enrico16/analysis/00genome/ with the specified options: -q --score-min L,0,-0.2 -p 2 --reorder --ignore-quals --no-mixed --no-discordant --maxins 500
- Now starting a Bowtie 2 paired-end alignment for CTread1GAread2CTgenome (reading in sequences from /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.1_C_to_T.fastq and /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair2.fastq.gz.temp.1_G_to_A.fastq, with the options: -q --score-min L,0,-0.2 -p 2 --reorder --ignore-quals --no-mixed --no-discordant --maxins 500 --norc))
- Found first alignment:
- HISEQ-MFG:108:HFC37ADXX:1:1101:1280:1985_1:N:0:TGGTGA/1 99 5_CT_converted 30762821 42 101M = 30762876 156 NTATAAATAAAAAAAAATATATAGGGTAGATAATAATTTAAATATATAATTAGTTTTTATAAATTAATATTTTATGAAGATATGTTAAAGTATTTTTTTAA #1=DDFFFHHHHHJJJIEEGDFFFIGDDFBBGGGIC>FHHGGEEGHGHAHE;FHIII?GHIE:CEHAE@CEFFDD;AAC>CEDCDDCDAC>CDEDEDDDBD AS:i:-1 XN:i:0 XM:i:1 XO:i:0 XG:i:0 NM:i:1 MD:Z:0A100 YS:i:0 YT:Z:CP
- HISEQ-MFG:108:HFC37ADXX:1:1101:1280:1985_2:N:0:TGGTGA/2 147 5_CT_converted 30762876 42 101M = 30762821 -156 TTTATAAATTAATATTTTATGAAGATATGTTAAAGTATTTTTTTAAAATAAAATATATAAAGAGGATGAGAAGTTAATTTATATATTAAAAAAGTATAAAT DCACADDCCEEECEEECC?:?@D@>EHCECEHDCIHGHGFEGGGJIIJIJIJJJJJJJJJIIIHFFBGF@H@HIJIHIJIJJIJJJJIDFGGHFFFFF@B@ AS:i:0 XN:i:0 XM:i:0 XO:i:0 XG:i:0 NM:i:0 MD:Z:101 YS:i:-1 YT:Z:CP
- Now starting a Bowtie 2 paired-end alignment for CTread1GAread2GAgenome (reading in sequences from /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.2_C_to_T.fastq and /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair2.fastq.gz.temp.2_G_to_A.fastq, with the options: -q --score-min L,0,-0.2 -p 2 --reorder --ignore-quals --no-mixed --no-discordant --maxins 500 --nofw))
- Found first alignment:
- HISEQ-MFG:108:HFC37ADXX:1:1101:1102:1988_1:N:0:TGGTGA/1 77 * 0 0 * * 0 0 NTAGTGTAAATATTGGAAAATTAGATTTTTTTTTTATAATAATATTGTATGTTTAAATTTTTAAGATTTTTGTTAAAATTAATATTAGTGAATAAAATTTT #1:ABDFFHHHHHGJIGHIJHHIJHHIJJJJJJIJGIGFGHGIEIJI@DHIGIHJJHHHHHCD>D@EEEDCCCDDDDDEDDCDDFEDEACACCCDDDEDAC YT:Z:UP
- HISEQ-MFG:108:HFC37ADXX:1:1101:1102:1988_2:N:0:TGGTGA/2 141 * 0 0 * * 0 0 ATTAACATTTTTACTATATTAATTCTTCTTATCCAAAATCACAATATATTTTTCCATAATTCAAATATTATATTATATTCATCAATAAAATTTCTCAAACA CCCFFFFFHHHHHJGIJJJJJJJJIGJIJJJJJJJJIEIIJJGGIGIIHJJJJIFHGGJDHHCFCGFGIGIIIJJJIIJIGIJIEHIGHCGHHEHHEDCDD YT:Z:UP
- Now starting a Bowtie 2 paired-end alignment for CTread1GAread2GAgenome (reading in sequences from /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.1_C_to_T.fastq and /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair2.fastq.gz.temp.1_G_to_A.fastq, with the options: -q --score-min L,0,-0.2 -p 2 --reorder --ignore-quals --no-mixed --no-discordant --maxins 500 --nofw))
- Found first alignment:
- HISEQ-MFG:108:HFC37ADXX:1:1101:1280:1985_1:N:0:TGGTGA/1 77 * 0 0 * * 0 0 NTATAAATAAAAAAAAATATATAGGGTAGATAATAATTTAAATATATAATTAGTTTTTATAAATTAATATTTTATGAAGATATGTTAAAGTATTTTTTTAA #1=DDFFFHHHHHJJJIEEGDFFFIGDDFBBGGGIC>FHHGGEEGHGHAHE;FHIII?GHIE:CEHAE@CEFFDD;AAC>CEDCDDCDAC>CDEDEDDDBD YT:Z:UP
- HISEQ-MFG:108:HFC37ADXX:1:1101:1280:1985_2:N:0:TGGTGA/2 141 * 0 0 * * 0 0 ATTTATACTTTTTTAATATATAAATTAACTTCTCATCCTCTTTATATATTTTATTTTAAAAAAATACTTTAACATATCTTCATAAAATATTAATTTATAAA @B@FFFFFHGGFDIJJJJIJJIJIHIJIH@H@FGBFFHIIIJJJJJJJJJIJIJIIJGGGEFGHGHICDHECECHE>@D@?:?CCEEECEEECCDDACACD YT:Z:UP
- >>> Writing bisulfite mapping results to SampleX-trimmed-pair1.fastq.gz.temp.2_bismark_bt2_pe.bam <<<
- Found first alignment:
- HISEQ-MFG:108:HFC37ADXX:1:1101:1102:1988_1:N:0:TGGTGA/1 83 7_GA_converted 16039576 42 101M = 16039304 -373 AAAATTTTATTCACTAATATTAATTTTAACAAAAATCTTAAAAATTTAAACATACAATATTATTATAAAAAAAAAATCTAATTTTCCAATATTTACACTAN CADEDDDCCCACAEDEFDDCDDEDDDDDCCCDEEE@D>DCHHHHHJJHIGIHD@IJIEIGHGFGIGJIJJJJJJIHHJIHHJIHGIJGHHHHHFFDBA:1# AS:i:-1 XN:i:0 XM:i:1 XO:i:0 XG:i:0 NM:i:1 MD:Z:100A0 YS:i:0 YT:Z:CP
- HISEQ-MFG:108:HFC37ADXX:1:1101:1102:1988_2:N:0:TGGTGA/2 163 7_GA_converted 16039304 42 101M = 16039576 373 ATTAACATTTTTACTATATTAATTCTTCTTATCCAAAATCACAATATATTTTTCCATAATTCAAATATTATATTATATTCATCAATAAAATTTCTCAAACA CCCFFFFFHHHHHJGIJJJJJJJJIGJIJJJJJJJJIEIIJJGGIGIIHJJJJIFHGGJDHHCFCGFGIGIIIJJJIIJIGIJIEHIGHCGHHEHHEDCDD AS:i:0 XN:i:0 XM:i:0 XO:i:0 XG:i:0 NM:i:0 MD:Z:101 YS:i:-1 YT:Z:CP
- >>> Writing bisulfite mapping results to SampleX-trimmed-pair1.fastq.gz.temp.1_bismark_bt2_pe.bam <<<
- Current working directory is: /home/enrico16/analysis/02trimming
- Now reading in and storing sequence information of the genome specified in: /home/enrico16/analysis/00genome/
- Current working directory is: /home/enrico16/analysis/02trimming
- Now reading in and storing sequence information of the genome specified in: /home/enrico16/analysis/00genome/
- chr 1 (248956422 bp)
- chr 1 (248956422 bp)
- chr 10 (133797422 bp)
- chr 10 (133797422 bp)
- chr 11 (135086622 bp)
- chr 11 (135086622 bp)
- chr 12 (133275309 bp)
- chr 12 (133275309 bp)
- chr 13 (114364328 bp)
- chr 13 (114364328 bp)
- chr 14 (107043718 bp)
- chr 14 (107043718 bp)
- chr 15 (101991189 bp)
- chr 15 (101991189 bp)
- chr 16 (90338345 bp)
- chr 16 (90338345 bp)
- chr 17 (83257441 bp)
- chr 17 (83257441 bp)
- chr 18 (80373285 bp)
- chr 18 (80373285 bp)
- chr 19 (58617616 bp)
- chr 19 (58617616 bp)
- chr 2 (242193529 bp)
- chr 2 (242193529 bp)
- chr 20 (64444167 bp)
- chr 20 (64444167 bp)
- chr 21 (46709983 bp)
- chr 21 (46709983 bp)
- chr 22 (50818468 bp)
- chr 22 (50818468 bp)
- chr 3 (198295559 bp)
- chr 3 (198295559 bp)
- chr 4 (190214555 bp)
- chr 4 (190214555 bp)
- chr 5 (181538259 bp)
- chr 5 (181538259 bp)
- chr 6 (170805979 bp)
- chr 6 (170805979 bp)
- chr 7 (159345973 bp)
- chr 7 (159345973 bp)
- chr 8 (145138636 bp)
- chr 8 (145138636 bp)
- chr 9 (138394717 bp)
- chr 9 (138394717 bp)
- chr MT (16569 bp)
- chr MT (16569 bp)
- chr X (156040895 bp)
- chr X (156040895 bp)
- chr Y (57227415 bp)
- chr Y (57227415 bp)
- chr KI270728.1 (1872759 bp)
- chr KI270728.1 (1872759 bp)
- chr KI270727.1 (448248 bp)
- chr KI270727.1 (448248 bp)
- chr KI270442.1 (392061 bp)
- chr KI270729.1 (280839 bp)
- chr KI270442.1 (392061 bp)
- chr GL000225.1 (211173 bp)
- chr KI270729.1 (280839 bp)
- chr KI270743.1 (210658 bp)
- chr GL000225.1 (211173 bp)
- chr GL000008.2 (209709 bp)
- chr KI270743.1 (210658 bp)
- chr GL000009.2 (201709 bp)
- chr GL000008.2 (209709 bp)
- chr KI270747.1 (198735 bp)
- chr GL000009.2 (201709 bp)
- chr KI270722.1 (194050 bp)
- chr KI270747.1 (198735 bp)
- chr GL000194.1 (191469 bp)
- chr KI270722.1 (194050 bp)
- chr KI270742.1 (186739 bp)
- chr GL000194.1 (191469 bp)
- chr GL000205.2 (185591 bp)
- chr KI270742.1 (186739 bp)
- chr GL000195.1 (182896 bp)
- chr GL000205.2 (185591 bp)
- chr KI270736.1 (181920 bp)
- chr GL000195.1 (182896 bp)
- chr KI270733.1 (179772 bp)
- chr KI270736.1 (181920 bp)
- chr GL000224.1 (179693 bp)
- chr KI270733.1 (179772 bp)
- chr GL000219.1 (179198 bp)
- chr GL000224.1 (179693 bp)
- chr KI270719.1 (176845 bp)
- chr GL000219.1 (179198 bp)
- chr GL000216.2 (176608 bp)
- chr KI270719.1 (176845 bp)
- chr KI270712.1 (176043 bp)
- chr GL000216.2 (176608 bp)
- chr KI270706.1 (175055 bp)
- chr KI270712.1 (176043 bp)
- chr KI270725.1 (172810 bp)
- chr KI270706.1 (175055 bp)
- chr KI270744.1 (168472 bp)
- chr KI270725.1 (172810 bp)
- chr KI270734.1 (165050 bp)
- chr KI270744.1 (168472 bp)
- chr GL000213.1 (164239 bp)
- chr KI270734.1 (165050 bp)
- chr GL000220.1 (161802 bp)
- chr GL000213.1 (164239 bp)
- chr KI270715.1 (161471 bp)
- chr GL000220.1 (161802 bp)
- chr GL000218.1 (161147 bp)
- chr KI270715.1 (161471 bp)
- chr KI270749.1 (158759 bp)
- chr KI270741.1 (157432 bp)
- chr GL000218.1 (161147 bp)
- chr GL000221.1 (155397 bp)
- chr KI270749.1 (158759 bp)
- chr KI270716.1 (153799 bp)
- chr KI270741.1 (157432 bp)
- chr GL000221.1 (155397 bp)
- chr KI270731.1 (150754 bp)
- chr KI270716.1 (153799 bp)
- chr KI270751.1 (150742 bp)
- chr KI270731.1 (150754 bp)
- chr KI270750.1 (148850 bp)
- chr KI270519.1 (138126 bp)
- chr KI270751.1 (150742 bp)
- chr GL000214.1 (137718 bp)
- chr KI270750.1 (148850 bp)
- chr KI270708.1 (127682 bp)
- chr KI270519.1 (138126 bp)
- chr KI270730.1 (112551 bp)
- chr GL000214.1 (137718 bp)
- chr KI270438.1 (112505 bp)
- chr KI270708.1 (127682 bp)
- chr KI270737.1 (103838 bp)
- chr KI270730.1 (112551 bp)
- chr KI270721.1 (100316 bp)
- chr KI270438.1 (112505 bp)
- chr KI270738.1 (99375 bp)
- chr KI270737.1 (103838 bp)
- chr KI270748.1 (93321 bp)
- chr KI270721.1 (100316 bp)
- chr KI270435.1 (92983 bp)
- chr GL000208.1 (92689 bp)
- chr KI270738.1 (99375 bp)
- chr KI270748.1 (93321 bp)
- chr KI270538.1 (91309 bp)
- chr KI270435.1 (92983 bp)
- chr KI270756.1 (79590 bp)
- chr GL000208.1 (92689 bp)
- chr KI270739.1 (73985 bp)
- chr KI270757.1 (71251 bp)
- chr KI270538.1 (91309 bp)
- chr KI270709.1 (66860 bp)
- chr KI270756.1 (79590 bp)
- chr KI270746.1 (66486 bp)
- chr KI270739.1 (73985 bp)
- chr KI270753.1 (62944 bp)
- chr KI270757.1 (71251 bp)
- chr KI270589.1 (44474 bp)
- chr KI270726.1 (43739 bp)
- chr KI270709.1 (66860 bp)
- chr KI270735.1 (42811 bp)
- chr KI270746.1 (66486 bp)
- chr KI270711.1 (42210 bp)
- chr KI270745.1 (41891 bp)
- chr KI270753.1 (62944 bp)
- chr KI270714.1 (41717 bp)
- chr KI270589.1 (44474 bp)
- chr KI270732.1 (41543 bp)
- chr KI270726.1 (43739 bp)
- chr KI270713.1 (40745 bp)
- chr KI270735.1 (42811 bp)
- chr KI270754.1 (40191 bp)
- chr KI270711.1 (42210 bp)
- chr KI270710.1 (40176 bp)
- chr KI270745.1 (41891 bp)
- chr KI270717.1 (40062 bp)
- chr KI270714.1 (41717 bp)
- chr KI270724.1 (39555 bp)
- chr KI270732.1 (41543 bp)
- chr KI270720.1 (39050 bp)
- chr KI270713.1 (40745 bp)
- chr KI270723.1 (38115 bp)
- chr KI270754.1 (40191 bp)
- chr KI270718.1 (38054 bp)
- chr KI270710.1 (40176 bp)
- chr KI270317.1 (37690 bp)
- chr KI270717.1 (40062 bp)
- chr KI270740.1 (37240 bp)
- chr KI270724.1 (39555 bp)
- chr KI270755.1 (36723 bp)
- chr KI270707.1 (32032 bp)
- chr KI270720.1 (39050 bp)
- chr KI270579.1 (31033 bp)
- chr KI270723.1 (38115 bp)
- chr KI270752.1 (27745 bp)
- chr KI270718.1 (38054 bp)
- chr KI270512.1 (22689 bp)
- chr KI270322.1 (21476 bp)
- chr KI270317.1 (37690 bp)
- chr GL000226.1 (15008 bp)
- chr KI270311.1 (12399 bp)
- chr KI270366.1 (8320 bp)
- chr KI270740.1 (37240 bp)
- chr KI270511.1 (8127 bp)
- chr KI270448.1 (7992 bp)
- chr KI270521.1 (7642 bp)
- chr KI270581.1 (7046 bp)
- chr KI270755.1 (36723 bp)
- chr KI270582.1 (6504 bp)
- chr KI270515.1 (6361 bp)
- chr KI270588.1 (6158 bp)
- chr KI270591.1 (5796 bp)
- chr KI270707.1 (32032 bp)
- chr KI270522.1 (5674 bp)
- chr KI270507.1 (5353 bp)
- chr KI270590.1 (4685 bp)
- chr KI270584.1 (4513 bp)
- chr KI270579.1 (31033 bp)
- chr KI270320.1 (4416 bp)
- chr KI270382.1 (4215 bp)
- chr KI270468.1 (4055 bp)
- chr KI270467.1 (3920 bp)
- chr KI270362.1 (3530 bp)
- chr KI270752.1 (27745 bp)
- chr KI270517.1 (3253 bp)
- chr KI270593.1 (3041 bp)
- chr KI270528.1 (2983 bp)
- chr KI270587.1 (2969 bp)
- chr KI270364.1 (2855 bp)
- chr KI270371.1 (2805 bp)
- chr KI270512.1 (22689 bp)
- chr KI270333.1 (2699 bp)
- chr KI270374.1 (2656 bp)
- chr KI270411.1 (2646 bp)
- chr KI270414.1 (2489 bp)
- chr KI270322.1 (21476 bp)
- chr KI270510.1 (2415 bp)
- chr KI270390.1 (2387 bp)
- chr KI270375.1 (2378 bp)
- chr KI270420.1 (2321 bp)
- chr GL000226.1 (15008 bp)
- chr KI270509.1 (2318 bp)
- chr KI270315.1 (2276 bp)
- chr KI270302.1 (2274 bp)
- chr KI270311.1 (12399 bp)
- chr KI270518.1 (2186 bp)
- chr KI270530.1 (2168 bp)
- chr KI270304.1 (2165 bp)
- chr KI270366.1 (8320 bp)
- chr KI270418.1 (2145 bp)
- chr KI270424.1 (2140 bp)
- chr KI270417.1 (2043 bp)
- chr KI270511.1 (8127 bp)
- chr KI270508.1 (1951 bp)
- chr KI270303.1 (1942 bp)
- chr KI270448.1 (7992 bp)
- chr KI270381.1 (1930 bp)
- chr KI270529.1 (1899 bp)
- chr KI270521.1 (7642 bp)
- chr KI270425.1 (1884 bp)
- chr KI270396.1 (1880 bp)
- chr KI270363.1 (1803 bp)
- chr KI270581.1 (7046 bp)
- chr KI270386.1 (1788 bp)
- chr KI270465.1 (1774 bp)
- chr KI270582.1 (6504 bp)
- chr KI270383.1 (1750 bp)
- chr KI270384.1 (1658 bp)
- chr KI270515.1 (6361 bp)
- chr KI270330.1 (1652 bp)
- chr KI270372.1 (1650 bp)
- chr KI270548.1 (1599 bp)
- chr KI270588.1 (6158 bp)
- chr KI270580.1 (1553 bp)
- chr KI270387.1 (1537 bp)
- chr KI270591.1 (5796 bp)
- chr KI270391.1 (1484 bp)
- chr KI270305.1 (1472 bp)
- chr KI270373.1 (1451 bp)
- chr KI270522.1 (5674 bp)
- chr KI270422.1 (1445 bp)
- chr KI270316.1 (1444 bp)
- chr KI270507.1 (5353 bp)
- chr KI270340.1 (1428 bp)
- chr KI270338.1 (1428 bp)
- chr KI270590.1 (4685 bp)
- chr KI270583.1 (1400 bp)
- chr KI270334.1 (1368 bp)
- chr KI270584.1 (4513 bp)
- chr KI270429.1 (1361 bp)
- chr KI270393.1 (1308 bp)
- chr KI270320.1 (4416 bp)
- chr KI270516.1 (1300 bp)
- chr KI270389.1 (1298 bp)
- chr KI270382.1 (4215 bp)
- chr KI270466.1 (1233 bp)
- chr KI270388.1 (1216 bp)
- chr KI270468.1 (4055 bp)
- chr KI270544.1 (1202 bp)
- chr KI270310.1 (1201 bp)
- chr KI270467.1 (3920 bp)
- chr KI270412.1 (1179 bp)
- chr KI270395.1 (1143 bp)
- chr KI270362.1 (3530 bp)
- chr KI270376.1 (1136 bp)
- chr KI270337.1 (1121 bp)
- chr KI270517.1 (3253 bp)
- chr KI270335.1 (1048 bp)
- chr KI270593.1 (3041 bp)
- chr KI270378.1 (1048 bp)
- chr KI270528.1 (2983 bp)
- chr KI270379.1 (1045 bp)
- chr KI270329.1 (1040 bp)
- chr KI270587.1 (2969 bp)
- chr KI270419.1 (1029 bp)
- chr KI270336.1 (1026 bp)
- chr KI270312.1 (998 bp)
- chr KI270364.1 (2855 bp)
- chr KI270539.1 (993 bp)
- chr KI270371.1 (2805 bp)
- chr KI270385.1 (990 bp)
- chr KI270423.1 (981 bp)
- chr KI270333.1 (2699 bp)
- chr KI270392.1 (971 bp)
- chr KI270394.1 (970 bp)
- chr KI270374.1 (2656 bp)
- chr KI270411.1 (2646 bp)
- chr KI270414.1 (2489 bp)
- chr KI270510.1 (2415 bp)
- chr KI270390.1 (2387 bp)
- chr KI270375.1 (2378 bp)
- chr KI270420.1 (2321 bp)
- chr KI270509.1 (2318 bp)
- chr KI270315.1 (2276 bp)
- chr KI270302.1 (2274 bp)
- chr KI270518.1 (2186 bp)
- chr KI270530.1 (2168 bp)
- chr KI270304.1 (2165 bp)
- chr KI270418.1 (2145 bp)
- chr KI270424.1 (2140 bp)
- chr KI270417.1 (2043 bp)
- chr KI270508.1 (1951 bp)
- chr KI270303.1 (1942 bp)
- chr KI270381.1 (1930 bp)
- chr KI270529.1 (1899 bp)
- chr KI270425.1 (1884 bp)
- chr KI270396.1 (1880 bp)
- chr KI270363.1 (1803 bp)
- chr KI270386.1 (1788 bp)
- chr KI270465.1 (1774 bp)
- chr KI270383.1 (1750 bp)
- chr KI270384.1 (1658 bp)
- chr KI270330.1 (1652 bp)
- chr KI270372.1 (1650 bp)
- chr KI270548.1 (1599 bp)
- chr KI270580.1 (1553 bp)
- chr KI270387.1 (1537 bp)
- chr KI270391.1 (1484 bp)
- chr KI270305.1 (1472 bp)
- chr KI270373.1 (1451 bp)
- chr KI270422.1 (1445 bp)
- chr KI270316.1 (1444 bp)
- chr KI270340.1 (1428 bp)
- chr KI270338.1 (1428 bp)
- chr KI270583.1 (1400 bp)
- chr KI270334.1 (1368 bp)
- chr KI270429.1 (1361 bp)
- chr KI270393.1 (1308 bp)
- chr KI270516.1 (1300 bp)
- chr KI270389.1 (1298 bp)
- chr KI270466.1 (1233 bp)
- chr KI270388.1 (1216 bp)
- chr KI270544.1 (1202 bp)
- chr KI270310.1 (1201 bp)
- chr KI270412.1 (1179 bp)
- chr KI270395.1 (1143 bp)
- chr KI270376.1 (1136 bp)
- chr KI270337.1 (1121 bp)
- chr KI270335.1 (1048 bp)
- chr KI270378.1 (1048 bp)
- chr KI270379.1 (1045 bp)
- chr KI270329.1 (1040 bp)
- chr KI270419.1 (1029 bp)
- chr KI270336.1 (1026 bp)
- chr KI270312.1 (998 bp)
- chr KI270539.1 (993 bp)
- chr KI270385.1 (990 bp)
- chr KI270423.1 (981 bp)
- chr KI270392.1 (971 bp)
- chr KI270394.1 (970 bp)
- Reading in the sequence files /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.1 and /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair2.fastq.gz.temp.1
- Reading in the sequence files /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.2 and /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair2.fastq.gz.temp.2
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:1:1103:7037:86167_1:N:0:TGGTGA KI270467.1 3819
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:1:1104:11259:2224_1:N:0:TGGTGA KI270337.1 2
- Processed 1000000 sequence pairs so far
- Processed 1000000 sequence pairs so far
- Processed 2000000 sequence pairs so far
- Processed 2000000 sequence pairs so far
- Processed 3000000 sequence pairs so far
- Processed 3000000 sequence pairs so far
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:1:1111:12667:101169_1:N:0:TGGTGA KI270738.1 1
- Processed 4000000 sequence pairs so far
- Processed 4000000 sequence pairs so far
- Processed 5000000 sequence pairs so far
- Processed 5000000 sequence pairs so far
- Processed 6000000 sequence pairs so far
- Processed 6000000 sequence pairs so far
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:1:1201:8205:62303_1:N:0:TGGTGA KI270584.1 4414
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:1:1204:2112:18281_1:N:0:TGGTGA KI270709.1 2
- Processed 7000000 sequence pairs so far
- Processed 7000000 sequence pairs so far
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:1:1205:2516:26717_1:N:0:TGGTGA KI270466.1 1133
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:1:1206:6007:67535_1:N:0:TGGTGA MT 1
- Processed 8000000 sequence pairs so far
- Processed 8000000 sequence pairs so far
- Processed 9000000 sequence pairs so far
- Processed 9000000 sequence pairs so far
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:1:1211:3180:52060_1:N:0:TGGTGA KI270709.1 2
- Processed 10000000 sequence pairs so far
- Processed 10000000 sequence pairs so far
- Processed 11000000 sequence pairs so far
- Processed 11000000 sequence pairs so far
- Processed 12000000 sequence pairs so far
- Processed 12000000 sequence pairs so far
- Processed 13000000 sequence pairs so far
- Processed 13000000 sequence pairs so far
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:1:2107:2805:52591_1:N:0:TGGTGA KI270738.1 1
- Processed 14000000 sequence pairs so far
- Processed 14000000 sequence pairs so far
- Processed 15000000 sequence pairs so far
- Processed 15000000 sequence pairs so far
- Processed 16000000 sequence pairs so far
- Processed 16000000 sequence pairs so far
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:1:2115:14930:32016_1:N:0:TGGTGA KI270738.1 1
- Processed 17000000 sequence pairs so far
- Processed 17000000 sequence pairs so far
- Processed 18000000 sequence pairs so far
- Processed 18000000 sequence pairs so far
- Processed 19000000 sequence pairs so far
- Processed 19000000 sequence pairs so far
- Processed 20000000 sequence pairs so far
- Processed 20000000 sequence pairs so far
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:1:2209:17422:84192_1:N:0:TGGTGA KI270735.1 2
- Processed 21000000 sequence pairs so far
- Processed 21000000 sequence pairs so far
- Processed 22000000 sequence pairs so far
- Processed 22000000 sequence pairs so far
- Processed 23000000 sequence pairs so far
- Processed 23000000 sequence pairs so far
- Processed 24000000 sequence pairs so far
- Processed 24000000 sequence pairs so far
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:2:1106:3802:49999_1:N:0:TGGTGA MT 2
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:2:1106:12497:100353_1:N:0:TGGTGA KI270736.1 181820
- Processed 25000000 sequence pairs so far
- Processed 25000000 sequence pairs so far
- Processed 26000000 sequence pairs so far
- Processed 26000000 sequence pairs so far
- Processed 27000000 sequence pairs so far
- Processed 27000000 sequence pairs so far
- Processed 28000000 sequence pairs so far
- Processed 28000000 sequence pairs so far
- Processed 29000000 sequence pairs so far
- Processed 29000000 sequence pairs so far
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:2:1204:9213:71254_1:N:0:TGGTGA KI270709.1 2
- Processed 30000000 sequence pairs so far
- Processed 30000000 sequence pairs so far
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:2:1205:19968:55797_1:N:0:TGGTGA KI270467.1 3818
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:2:1207:11346:89023_1:N:0:TGGTGA KI270589.1 2
- Processed 31000000 sequence pairs so far
- Processed 31000000 sequence pairs so far
- Processed 32000000 sequence pairs so far
- Processed 32000000 sequence pairs so far
- Processed 33000000 sequence pairs so far
- Processed 33000000 sequence pairs so far
- Processed 34000000 sequence pairs so far
- Processed 34000000 sequence pairs so far
- Processed 35000000 sequence pairs so far
- Processed 35000000 sequence pairs so far
- Processed 36000000 sequence pairs so far
- Processed 36000000 sequence pairs so far
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:2:2107:12713:58820_1:N:0:TGGTGA KI270467.1 3819
- Processed 37000000 sequence pairs so far
- Processed 37000000 sequence pairs so far
- Processed 38000000 sequence pairs so far
- Processed 38000000 sequence pairs so far
- Processed 39000000 sequence pairs so far
- Processed 39000000 sequence pairs so far
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:2:2115:17081:9319_1:N:0:TGGTGA KI270466.1 1133
- Processed 40000000 sequence pairs so far
- Processed 40000000 sequence pairs so far
- Processed 41000000 sequence pairs so far
- Processed 41000000 sequence pairs so far
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:2:2205:17864:60244_1:N:0:TGGTGA KI270726.1 1
- Processed 42000000 sequence pairs so far
- Processed 42000000 sequence pairs so far
- Processed 43000000 sequence pairs so far
- Processed 43000000 sequence pairs so far
- Processed 44000000 sequence pairs so far
- Processed 44000000 sequence pairs so far
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:2:2215:20794:36465_1:N:0:TGGTGA KI270735.1 2
- Processed 45000000 sequence pairs so far
- Processed 45000000 sequence pairs so far
- Chromosomal sequence could not be extracted for HISEQ-MFG:108:HFC37ADXX:2:2216:8707:3858_1:N:0:TGGTGA KI270467.1 3818
- 45829609 reads; of these:
- 45829609 (100.00%) were paired; of these:
- 26578709 (57.99%) aligned concordantly 0 times
- 16759687 (36.57%) aligned concordantly exactly 1 time
- 2491213 (5.44%) aligned concordantly >1 times
- 42.01% overall alignment rate
- 45829609 reads; of these:
- 45829609 (100.00%) were paired; of these:
- 26622152 (58.09%) aligned concordantly 0 times
- 16684357 (36.41%) aligned concordantly exactly 1 time
- 2523100 (5.51%) aligned concordantly >1 times
- 41.91% overall alignment rate
- Processed 45829609 sequences in total
- Successfully deleted the temporary files /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.1_C_to_T.fastq and /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair2.fastq.gz.temp.1_G_to_A.fastq
- Final Alignment report
- ======================
- Sequence pairs analysed in total: 45829609
- Number of paired-end alignments with a unique best hit: 33708979
- Mapping efficiency: 73.6%
- Sequence pairs with no alignments under any condition: 9908734
- Sequence pairs did not map uniquely: 2211896
- Sequence pairs which were discarded because genomic sequence could not be extracted: 10
- Number of sequence pairs with unique best (first) alignment came from the bowtie output:
- CT/GA/CT: 16895340 ((converted) top strand)
- GA/CT/CT: 0 (complementary to (converted) top strand)
- GA/CT/GA: 0 (complementary to (converted) bottom strand)
- CT/GA/GA: 16813629 ((converted) bottom strand)
- Number of alignments to (merely theoretical) complementary strands being rejected in total: 0
- Final Cytosine Methylation Report
- =================================
- Total number of C's analysed: 1031115512
- Total methylated C's in CpG context: 34462143
- Total methylated C's in CHG context: 4048566
- Total methylated C's in CHH context: 14008321
- Total methylated C's in Unknown context: 61
- Total unmethylated C's in CpG context: 7964298
- Total unmethylated C's in CHG context: 208418423
- Total unmethylated C's in CHH context: 762213761
- Total unmethylated C's in Unknown context: 1110
- C methylated in CpG context: 81.2%
- C methylated in CHG context: 1.9%
- C methylated in CHH context: 1.8%
- C methylated in unknown context (CN or CHN): 5.2%
- 45829608 reads; of these:
- 45829608 (100.00%) were paired; of these:
- 26572015 (57.98%) aligned concordantly 0 times
- 16764984 (36.58%) aligned concordantly exactly 1 time
- 2492609 (5.44%) aligned concordantly >1 times
- 42.02% overall alignment rate
- 45829608 reads; of these:
- 45829608 (100.00%) were paired; of these:
- 26622711 (58.09%) aligned concordantly 0 times
- 16681809 (36.40%) aligned concordantly exactly 1 time
- 2525088 (5.51%) aligned concordantly >1 times
- 41.91% overall alignment rate
- Processed 45829608 sequences in total
- Successfully deleted the temporary files /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.2_C_to_T.fastq and /home/enrico16/analysis/05alignment/SampleX/SampleX-trimmed-pair2.fastq.gz.temp.2_G_to_A.fastq
- Final Alignment report
- ======================
- Sequence pairs analysed in total: 45829608
- Number of paired-end alignments with a unique best hit: 33711074
- Mapping efficiency: 73.6%
- Sequence pairs with no alignments under any condition: 9904043
- Sequence pairs did not map uniquely: 2214491
- Sequence pairs which were discarded because genomic sequence could not be extracted: 11
- Number of sequence pairs with unique best (first) alignment came from the bowtie output:
- CT/GA/CT: 16891215 ((converted) top strand)
- GA/CT/CT: 0 (complementary to (converted) top strand)
- GA/CT/GA: 0 (complementary to (converted) bottom strand)
- CT/GA/GA: 16819848 ((converted) bottom strand)
- Number of alignments to (merely theoretical) complementary strands being rejected in total: 0
- Final Cytosine Methylation Report
- =================================
- Total number of C's analysed: 1030909873
- Total methylated C's in CpG context: 34455854
- Total methylated C's in CHG context: 4048303
- Total methylated C's in CHH context: 13992717
- Total methylated C's in Unknown context: 64
- Total unmethylated C's in CpG context: 7965039
- Total unmethylated C's in CHG context: 208377228
- Total unmethylated C's in CHH context: 762070732
- Total unmethylated C's in Unknown context: 1100
- C methylated in CpG context: 81.2%
- C methylated in CHG context: 1.9%
- C methylated in CHH context: 1.8%
- C methylated in unknown context (CN or CHN): 5.5%
- Now waiting for all child processes to complete
- Right, cleaning up now...
- Deleting temporary sequence files...
- /home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.1 Failed to delete temporary FastQ file /home/enrico16/analysis/05alignment/SampleX//home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.1: No such file or directory
- /home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair2.fastq.gz.temp.1 Failed to delete temporary FastQ file /home/enrico16/analysis/05alignment/SampleX//home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair2.fastq.gz.temp.1: No such file or directory
- /home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.2 Failed to delete temporary FastQ file /home/enrico16/analysis/05alignment/SampleX//home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.2: No such file or directory
- /home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair2.fastq.gz.temp.2 Failed to delete temporary FastQ file /home/enrico16/analysis/05alignment/SampleX//home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair2.fastq.gz.temp.2: No such file or directory
- Now merging BAM files /home/enrico16/analysis/05alignment/SampleX//home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.1_bismark_bt2_pe.bam /home/enrico16/analysis/05alignment/SampleX//home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.2_bismark_bt2_pe.bam into >>> /home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz_bismark_bt2_pe.bam <<<
- Merging from file >> /home/enrico16/analysis/05alignment/SampleX//home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.1_bismark_bt2_pe.bam <<
- open: No such file or directory
- [main_samview] fail to open "/home/enrico16/analysis/05alignment/SampleX//home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.1_bismark_bt2_pe.bam" for reading.
- Failed to close filehandle
- Merging from file >> /home/enrico16/analysis/05alignment/SampleX//home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.2_bismark_bt2_pe.bam <<
- open: No such file or directory
- [main_samview] fail to open "/home/enrico16/analysis/05alignment/SampleX//home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.2_bismark_bt2_pe.bam" for reading.
- Failed to close filehandle
- Failed to close output filehandle
- Deleting temporary BAM files...
- /home/enrico16/analysis/05alignment/SampleX//home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.1_bismark_bt2_pe.bam Failed to delete temporary BAM file /home/enrico16/analysis/05alignment/SampleX//home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.1_bismark_bt2_pe.bam: No such file or directory
- /home/enrico16/analysis/05alignment/SampleX//home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.2_bismark_bt2_pe.bam Failed to delete temporary BAM file /home/enrico16/analysis/05alignment/SampleX//home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz.temp.2_bismark_bt2_pe.bam: No such file or directory
- Failed to write to /home/enrico16/analysis/02trimming/SampleX/SampleX-trimmed-pair1.fastq.gz_bismark_bt2_PE_report.txt: No such file or directory
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement