/opt/biosoftware/tophat/tophat/tophat -N 5 --read-edit-dist 5 -p 20 -o IRTestis02FC_1 /home/rneme/PopTranscriptomes/mapping/index/mm10 pSN7640224_16800_IRTestis02FC_1_1_sequence.fq.gz,pSN7640224_16800_IRTestis02FC_1_2_sequence.fq.gz
#>prep_reads:
/opt/biosoftware/tophat/tophat-2.0.8b.Linux_x86_64/prep_reads --min-anchor 8 --splice-mismatches 0 --min-report-intron 50 --max-report-intron 500000 --min-isoform-fraction 0.15 --output-dir IRTestis02FC_1/ --max-multihits 20 --max-seg-multihits 40 --segment-length 25 --segment-mismatches 2 --min-closure-exon 100 --min-closure-intron 50 --max-closure-intron 5000 --min-coverage-intron 50 --max-coverage-intron 20000 --min-segment-intron 50 --max-segment-intron 500000 --read-mismatches 5 --read-gap-length 2 --read-edit-dist 5 --read-realign-edit-dist 6 --max-insertion-length 3 --max-deletion-length 3 -z gzip -p20 --no-closure-search --no-coverage-search --no-microexon-search --fastq --aux-outfile=IRTestis02FC_1/prep_reads.info --index-outfile=IRTestis02FC_1/tmp/left_kept_reads.bam.index --sam-header=IRTestis02FC_1/tmp/mm10_genome.bwt.samheader.sam --outfile=IRTestis02FC_1/tmp/left_kept_reads.bam pSN7640224_16800_IRTestis02FC_1_1_sequence.fq.gz,pSN7640224_16800_IRTestis02FC_1_2_sequence.fq.gz
#>map_start:
/opt/biosoftware/tophat/tophat-2.0.8b.Linux_x86_64/bam2fastx --all --fastq IRTestis02FC_1/tmp/left_kept_reads.bam|/opt/biosoftware/bowtie/bowtie2/bowtie2-align -q -k 20 -D 15 -R 2 -N 0 -L 20 -i S,1,1.25 --gbar 4 --mp 6,2 --np 1 --rdg 5,3 --rfg 5,3 --score-min C,-32,0 -p 20 --sam-no-hd -x /home/rneme/PopTranscriptomes/mapping/index/mm10 -|/opt/biosoftware/tophat/tophat-2.0.8b.Linux_x86_64/fix_map_ordering --bowtie2-min-score 30 --read-mismatches 5 --read-gap-length 2 --read-edit-dist 5 --read-realign-edit-dist 6 --index-outfile IRTestis02FC_1/tmp/left_kept_reads.mapped.bam.index --sam-header IRTestis02FC_1/tmp/mm10_genome.bwt.samheader.sam - IRTestis02FC_1/tmp/left_kept_reads.mapped.bam IRTestis02FC_1/tmp/left_kept_reads_unmapped.bam
#>map_segments:
gzip -cd< IRTestis02FC_1/tmp/left_kept_reads_seg1.fq.z|/opt/biosoftware/bowtie/bowtie2/bowtie2-align -q -k 41 -N 1 -L 20 -p 20 --sam-no-hd -x /home/rneme/PopTranscriptomes/mapping/index/mm10 -|/opt/biosoftware/tophat/tophat-2.0.8b.Linux_x86_64/fix_map_ordering --bowtie2-min-score 10 --read-mismatches 5 --read-gap-length 2 --read-edit-dist 5 --read-realign-edit-dist 6 --index-outfile IRTestis02FC_1/tmp/left_kept_reads_seg1.bam.index --sam-header IRTestis02FC_1/tmp/mm10_genome.bwt.samheader.sam - IRTestis02FC_1/tmp/left_kept_reads_seg1.bam IRTestis02FC_1/tmp/left_kept_reads_seg1_unmapped.bam
gzip -cd< IRTestis02FC_1/tmp/left_kept_reads_seg2.fq.z|/opt/biosoftware/bowtie/bowtie2/bowtie2-align -q -k 41 -N 1 -L 20 -p 20 --sam-no-hd -x /home/rneme/PopTranscriptomes/mapping/index/mm10 -|/opt/biosoftware/tophat/tophat-2.0.8b.Linux_x86_64/fix_map_ordering --bowtie2-min-score 10 --read-mismatches 5 --read-gap-length 2 --read-edit-dist 5 --read-realign-edit-dist 6 --index-outfile IRTestis02FC_1/tmp/left_kept_reads_seg2.bam.index --sam-header IRTestis02FC_1/tmp/mm10_genome.bwt.samheader.sam - IRTestis02FC_1/tmp/left_kept_reads_seg2.bam IRTestis02FC_1/tmp/left_kept_reads_seg2_unmapped.bam
gzip -cd< IRTestis02FC_1/tmp/left_kept_reads_seg3.fq.z|/opt/biosoftware/bowtie/bowtie2/bowtie2-align -q -k 41 -N 1 -L 20 -p 20 --sam-no-hd -x /home/rneme/PopTranscriptomes/mapping/index/mm10 -|/opt/biosoftware/tophat/tophat-2.0.8b.Linux_x86_64/fix_map_ordering --bowtie2-min-score 10 --read-mismatches 5 --read-gap-length 2 --read-edit-dist 5 --read-realign-edit-dist 6 --index-outfile IRTestis02FC_1/tmp/left_kept_reads_seg3.bam.index --sam-header IRTestis02FC_1/tmp/mm10_genome.bwt.samheader.sam - IRTestis02FC_1/tmp/left_kept_reads_seg3.bam IRTestis02FC_1/tmp/left_kept_reads_seg3_unmapped.bam
gzip -cd< IRTestis02FC_1/tmp/left_kept_reads_seg4.fq.z|/opt/biosoftware/bowtie/bowtie2/bowtie2-align -q -k 41 -N 1 -L 20 -p 20 --sam-no-hd -x /home/rneme/PopTranscriptomes/mapping/index/mm10 -|/opt/biosoftware/tophat/tophat-2.0.8b.Linux_x86_64/fix_map_ordering --bowtie2-min-score 10 --read-mismatches 5 --read-gap-length 2 --read-edit-dist 5 --read-realign-edit-dist 6 --index-outfile IRTestis02FC_1/tmp/left_kept_reads_seg4.bam.index --sam-header IRTestis02FC_1/tmp/mm10_genome.bwt.samheader.sam - IRTestis02FC_1/tmp/left_kept_reads_seg4.bam IRTestis02FC_1/tmp/left_kept_reads_seg4_unmapped.bam
#>find_juncs:
/opt/biosoftware/tophat/tophat-2.0.8b.Linux_x86_64/segment_juncs --min-anchor 8 --splice-mismatches 0 --min-report-intron 50 --max-report-intron 500000 --min-isoform-fraction 0.15 --output-dir IRTestis02FC_1/ --max-multihits 20 --max-seg-multihits 40 --segment-length 25 --segment-mismatches 2 --min-closure-exon 100 --min-closure-intron 50 --max-closure-intron 5000 --min-coverage-intron 50 --max-coverage-intron 20000 --min-segment-intron 50 --max-segment-intron 500000 --read-mismatches 5 --read-gap-length 2 --read-edit-dist 5 --read-realign-edit-dist 6 --max-insertion-length 3 --max-deletion-length 3 -z gzip -p20 --no-closure-search --no-coverage-search --no-microexon-search --sam-header IRTestis02FC_1/tmp/mm10_genome.bwt.samheader.sam --ium-reads IRTestis02FC_1/tmp/left_kept_reads_seg1_unmapped.bam,IRTestis02FC_1/tmp/left_kept_reads_seg2_unmapped.bam,IRTestis02FC_1/tmp/left_kept_reads_seg3_unmapped.bam,IRTestis02FC_1/tmp/left_kept_reads_seg4_unmapped.bam /home/rneme/PopTranscriptomes/mapping/index/mm10.fa IRTestis02FC_1/tmp/segment.juncs IRTestis02FC_1/tmp/segment.insertions IRTestis02FC_1/tmp/segment.deletions IRTestis02FC_1/tmp/segment.fusions IRTestis02FC_1/tmp/left_kept_reads.bam IRTestis02FC_1/tmp/left_kept_reads.mapped.bam IRTestis02FC_1/tmp/left_kept_reads_seg1.bam,IRTestis02FC_1/tmp/left_kept_reads_seg2.bam,IRTestis02FC_1/tmp/left_kept_reads_seg3.bam,IRTestis02FC_1/tmp/left_kept_reads_seg4.bam
#>juncs_db:
/opt/biosoftware/tophat/tophat-2.0.8b.Linux_x86_64/juncs_db 3 26 IRTestis02FC_1/tmp/segment.juncs /dev/null /dev/null /dev/null /home/rneme/PopTranscriptomes/mapping/index/mm10.fa > IRTestis02FC_1/tmp/segment_juncs.fa
/opt/biosoftware/bowtie/bowtie2/bowtie2-build IRTestis02FC_1/tmp/segment_juncs.fa IRTestis02FC_1/tmp/segment_juncs
#>map2juncs:
gzip -cd< IRTestis02FC_1/tmp/left_kept_reads_seg1.fq.z|/opt/biosoftware/bowtie/bowtie2/bowtie2-align -q -k 41 -N 1 -L 20 -p 20 --sam-no-hd -x IRTestis02FC_1/tmp/segment_juncs -|/opt/biosoftware/tophat/tophat-2.0.8b.Linux_x86_64/fix_map_ordering --bowtie2-min-score 10 --read-mismatches 5 --read-gap-length 2 --read-edit-dist 5 --read-realign-edit-dist 6 --index-outfile IRTestis02FC_1/tmp/left_kept_reads_seg1.to_spliced.bam.index --sam-header IRTestis02FC_1/tmp/segment_juncs.bwt.samheader.sam - IRTestis02FC_1/tmp/left_kept_reads_seg1.to_spliced.bam
gzip -cd< IRTestis02FC_1/tmp/left_kept_reads_seg2.fq.z|/opt/biosoftware/bowtie/bowtie2/bowtie2-align -q -k 41 -N 1 -L 20 -p 20 --sam-no-hd -x IRTestis02FC_1/tmp/segment_juncs -|/opt/biosoftware/tophat/tophat-2.0.8b.Linux_x86_64/fix_map_ordering --bowtie2-min-score 10 --read-mismatches 5 --read-gap-length 2 --read-edit-dist 5 --read-realign-edit-dist 6 --index-outfile IRTestis02FC_1/tmp/left_kept_reads_seg2.to_spliced.bam.index --sam-header IRTestis02FC_1/tmp/segment_juncs.bwt.samheader.sam - IRTestis02FC_1/tmp/left_kept_reads_seg2.to_spliced.bam
gzip -cd< IRTestis02FC_1/tmp/left_kept_reads_seg3.fq.z|/opt/biosoftware/bowtie/bowtie2/bowtie2-align -q -k 41 -N 1 -L 20 -p 20 --sam-no-hd -x IRTestis02FC_1/tmp/segment_juncs -|/opt/biosoftware/tophat/tophat-2.0.8b.Linux_x86_64/fix_map_ordering --bowtie2-min-score 10 --read-mismatches 5 --read-gap-length 2 --read-edit-dist 5 --read-realign-edit-dist 6 --index-outfile IRTestis02FC_1/tmp/left_kept_reads_seg3.to_spliced.bam.index --sam-header IRTestis02FC_1/tmp/segment_juncs.bwt.samheader.sam - IRTestis02FC_1/tmp/left_kept_reads_seg3.to_spliced.bam
gzip -cd< IRTestis02FC_1/tmp/left_kept_reads_seg4.fq.z|/opt/biosoftware/bowtie/bowtie2/bowtie2-align -q -k 41 -N 1 -L 20 -p 20 --sam-no-hd -x IRTestis02FC_1/tmp/segment_juncs -|/opt/biosoftware/tophat/tophat-2.0.8b.Linux_x86_64/fix_map_ordering --bowtie2-min-score 10 --read-mismatches 5 --read-gap-length 2 --read-edit-dist 5 --read-realign-edit-dist 6 --index-outfile IRTestis02FC_1/tmp/left_kept_reads_seg4.to_spliced.bam.index --sam-header IRTestis02FC_1/tmp/segment_juncs.bwt.samheader.sam - IRTestis02FC_1/tmp/left_kept_reads_seg4.to_spliced.bam
/opt/biosoftware/tophat/tophat-2.0.8b.Linux_x86_64/long_spanning_reads --min-anchor 8 --splice-mismatches 0 --min-report-intron 50 --max-report-intron 500000 --min-isoform-fraction 0.15 --output-dir IRTestis02FC_1/ --max-multihits 20 --max-seg-multihits 40 --segment-length 25 --segment-mismatches 2 --min-closure-exon 100 --min-closure-intron 50 --max-closure-intron 5000 --min-coverage-intron 50 --max-coverage-intron 20000 --min-segment-intron 50 --max-segment-intron 500000 --read-mismatches 5 --read-gap-length 2 --read-edit-dist 5 --read-realign-edit-dist 6 --max-insertion-length 3 --max-deletion-length 3 -z gzip -p20 --no-closure-search --no-coverage-search --no-microexon-search --sam-header IRTestis02FC_1/tmp/mm10_genome.bwt.samheader.sam --bowtie2-max-penalty 6 --bowtie2-min-penalty 2 --bowtie2-penalty-for-N 1 --bowtie2-read-gap-open 5 --bowtie2-read-gap-cont 3 --bowtie2-ref-gap-open 5 --bowtie2-ref-gap-cont 3 /home/rneme/PopTranscriptomes/mapping/index/mm10.fa IRTestis02FC_1/tmp/left_kept_reads.bam IRTestis02FC_1/tmp/segment.juncs IRTestis02FC_1/tmp/segment.insertions IRTestis02FC_1/tmp/segment.deletions /dev/null IRTestis02FC_1/tmp/left_kept_reads.candidates.bam IRTestis02FC_1/tmp/left_kept_reads_seg1.bam,IRTestis02FC_1/tmp/left_kept_reads_seg2.bam,IRTestis02FC_1/tmp/left_kept_reads_seg3.bam,IRTestis02FC_1/tmp/left_kept_reads_seg4.bam IRTestis02FC_1/tmp/left_kept_reads_seg1.to_spliced.bam,IRTestis02FC_1/tmp/left_kept_reads_seg2.to_spliced.bam,IRTestis02FC_1/tmp/left_kept_reads_seg3.to_spliced.bam,IRTestis02FC_1/tmp/left_kept_reads_seg4.to_spliced.bam
#>tophat_reports:
/opt/biosoftware/tophat/tophat-2.0.8b.Linux_x86_64/tophat_reports --min-anchor 8 --splice-mismatches 0 --min-report-intron 50 --max-report-intron 500000 --min-isoform-fraction 0.15 --output-dir IRTestis02FC_1/ --max-multihits 20 --max-seg-multihits 40 --segment-length 25 --segment-mismatches 2 --min-closure-exon 100 --min-closure-intron 50 --max-closure-intron 5000 --min-coverage-intron 50 --max-coverage-intron 20000 --min-segment-intron 50 --max-segment-intron 500000 --read-mismatches 5 --read-gap-length 2 --read-edit-dist 5 --read-realign-edit-dist 6 --max-insertion-length 3 --max-deletion-length 3 -z gzip -p20 --no-closure-search --no-coverage-search --no-microexon-search --sam-header IRTestis02FC_1/tmp/mm10_genome.bwt.samheader.sam --report-discordant-pair-alignments --report-mixed-alignments --samtools=/usr/local/bin/samtools --bowtie2-max-penalty 6 --bowtie2-min-penalty 2 --bowtie2-penalty-for-N 1 --bowtie2-read-gap-open 5 --bowtie2-read-gap-cont 3 --bowtie2-ref-gap-open 5 --bowtie2-ref-gap-cont 3 /home/rneme/PopTranscriptomes/mapping/index/mm10.fa IRTestis02FC_1/junctions.bed IRTestis02FC_1/insertions.bed IRTestis02FC_1/deletions.bed IRTestis02FC_1/fusions.out IRTestis02FC_1/tmp/accepted_hits IRTestis02FC_1/tmp/left_kept_reads.mapped.bam,IRTestis02FC_1/tmp/left_kept_reads.candidates IRTestis02FC_1/tmp/left_kept_reads.bam
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits0.bam IRTestis02FC_1/tmp/accepted_hits0_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits1.bam IRTestis02FC_1/tmp/accepted_hits1_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits2.bam IRTestis02FC_1/tmp/accepted_hits2_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits3.bam IRTestis02FC_1/tmp/accepted_hits3_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits4.bam IRTestis02FC_1/tmp/accepted_hits4_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits5.bam IRTestis02FC_1/tmp/accepted_hits5_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits6.bam IRTestis02FC_1/tmp/accepted_hits6_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits7.bam IRTestis02FC_1/tmp/accepted_hits7_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits8.bam IRTestis02FC_1/tmp/accepted_hits8_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits9.bam IRTestis02FC_1/tmp/accepted_hits9_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits10.bam IRTestis02FC_1/tmp/accepted_hits10_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits11.bam IRTestis02FC_1/tmp/accepted_hits11_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits12.bam IRTestis02FC_1/tmp/accepted_hits12_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits13.bam IRTestis02FC_1/tmp/accepted_hits13_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits14.bam IRTestis02FC_1/tmp/accepted_hits14_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits15.bam IRTestis02FC_1/tmp/accepted_hits15_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits16.bam IRTestis02FC_1/tmp/accepted_hits16_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits17.bam IRTestis02FC_1/tmp/accepted_hits17_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits18.bam IRTestis02FC_1/tmp/accepted_hits18_sorted
/usr/local/bin/samtools sort IRTestis02FC_1/tmp/accepted_hits19.bam IRTestis02FC_1/tmp/accepted_hits19_sorted
/usr/local/bin/samtools merge -f -h IRTestis02FC_1/tmp/mm10_genome.bwt.samheader.sam IRTestis02FC_1/accepted_hits.bam IRTestis02FC_1/tmp/accepted_hits0_sorted.bam IRTestis02FC_1/tmp/accepted_hits1_sorted.bam IRTestis02FC_1/tmp/accepted_hits2_sorted.bam IRTestis02FC_1/tmp/accepted_hits3_sorted.bam IRTestis02FC_1/tmp/accepted_hits4_sorted.bam IRTestis02FC_1/tmp/accepted_hits5_sorted.bam IRTestis02FC_1/tmp/accepted_hits6_sorted.bam IRTestis02FC_1/tmp/accepted_hits7_sorted.bam IRTestis02FC_1/tmp/accepted_hits8_sorted.bam IRTestis02FC_1/tmp/accepted_hits9_sorted.bam IRTestis02FC_1/tmp/accepted_hits10_sorted.bam IRTestis02FC_1/tmp/accepted_hits11_sorted.bam IRTestis02FC_1/tmp/accepted_hits12_sorted.bam IRTestis02FC_1/tmp/accepted_hits13_sorted.bam IRTestis02FC_1/tmp/accepted_hits14_sorted.bam IRTestis02FC_1/tmp/accepted_hits15_sorted.bam IRTestis02FC_1/tmp/accepted_hits16_sorted.bam IRTestis02FC_1/tmp/accepted_hits17_sorted.bam IRTestis02FC_1/tmp/accepted_hits18_sorted.bam IRTestis02FC_1/tmp/accepted_hits19_sorted.bam
/opt/biosoftware/tophat/tophat-2.0.8b.Linux_x86_64/bam_merge -Q --sam-header IRTestis02FC_1/tmp/mm10_genome.bwt.samheader.sam IRTestis02FC_1/unmapped.bam IRTestis02FC_1/tmp/unmapped_left_0.bam IRTestis02FC_1/tmp/unmapped_left_1.bam IRTestis02FC_1/tmp/unmapped_left_2.bam IRTestis02FC_1/tmp/unmapped_left_3.bam IRTestis02FC_1/tmp/unmapped_left_4.bam IRTestis02FC_1/tmp/unmapped_left_5.bam IRTestis02FC_1/tmp/unmapped_left_6.bam IRTestis02FC_1/tmp/unmapped_left_7.bam IRTestis02FC_1/tmp/unmapped_left_8.bam IRTestis02FC_1/tmp/unmapped_left_9.bam IRTestis02FC_1/tmp/unmapped_left_10.bam IRTestis02FC_1/tmp/unmapped_left_11.bam IRTestis02FC_1/tmp/unmapped_left_12.bam IRTestis02FC_1/tmp/unmapped_left_13.bam IRTestis02FC_1/tmp/unmapped_left_14.bam IRTestis02FC_1/tmp/unmapped_left_15.bam IRTestis02FC_1/tmp/unmapped_left_16.bam IRTestis02FC_1/tmp/unmapped_left_17.bam IRTestis02FC_1/tmp/unmapped_left_18.bam IRTestis02FC_1/tmp/unmapped_left_19.bam
#>alldone: