STAR version=2.7.1a STAR compilation time,server,dir=vie jun 21 18:03:03 CEST 2019 ant-login8.linux.crg.es:/software/mi/el7.2/software/STAR-2.7.1a/source ##### DEFAULT parameters: versionGenome 2.7.1a parametersFiles - sysShell - runMode alignReads runThreadN 1 runDirPerm User_RWX runRNGseed 777 genomeDir ./GenomeDir/ genomeLoad NoSharedMemory genomeFastaFiles - genomeChainFiles - genomeSAindexNbases 14 genomeChrBinNbits 18 genomeSAsparseD 1 genomeSuffixLengthMax 18446744073709551615 genomeFileSizes 0 genomeConsensusFile - readFilesType Fastx readFilesIn Read1 Read2 readFilesPrefix - readFilesCommand - readMatesLengthsIn NotEqual readMapNumber 18446744073709551615 readNameSeparator / inputBAMfile - bamRemoveDuplicatesType - bamRemoveDuplicatesMate2basesN 0 limitGenomeGenerateRAM 31000000000 limitIObufferSize 150000000 limitOutSAMoneReadBytes 100000 limitOutSJcollapsed 1000000 limitOutSJoneRead 1000 limitBAMsortRAM 0 limitSjdbInsertNsj 1000000 limitNreadsSoft 18446744073709551615 outTmpDir - outTmpKeep None outStd Log outReadsUnmapped None outQSconversionAdd 0 outMultimapperOrder Old_2.4 outSAMtype SAM outSAMmode Full outSAMstrandField None outSAMattributes Standard outSAMunmapped None outSAMorder Paired outSAMprimaryFlag OneBestScore outSAMreadID Standard outSAMmapqUnique 255 outSAMflagOR 0 outSAMflagAND 65535 outSAMattrRGline - outSAMheaderHD - outSAMheaderPG - outSAMheaderCommentFile - outBAMcompression 1 outBAMsortingThreadN 0 outBAMsortingBinsN 50 outSAMfilter None outSAMmultNmax 18446744073709551615 outSAMattrIHstart 1 outSAMtlen 1 outSJfilterReads All outSJfilterCountUniqueMin 3 1 1 1 outSJfilterCountTotalMin 3 1 1 1 outSJfilterOverhangMin 30 12 12 12 outSJfilterDistToOtherSJmin 10 0 5 10 outSJfilterIntronMaxVsReadN 50000 100000 200000 outWigType None outWigStrand Stranded outWigReferencesPrefix - outWigNorm RPM outFilterType Normal outFilterMultimapNmax 10 outFilterMultimapScoreRange 1 outFilterScoreMin 0 outFilterScoreMinOverLread 0.66 outFilterMatchNmin 0 outFilterMatchNminOverLread 0.66 outFilterMismatchNmax 10 outFilterMismatchNoverLmax 0.3 outFilterMismatchNoverReadLmax 1 outFilterIntronMotifs None outFilterIntronStrands RemoveInconsistentStrands clip5pNbases 0 clip3pNbases 0 clip3pAfterAdapterNbases 0 clip3pAdapterSeq - clip3pAdapterMMp 0.1 winBinNbits 16 winAnchorDistNbins 9 winFlankNbins 4 winAnchorMultimapNmax 50 winReadCoverageRelativeMin 0.5 winReadCoverageBasesMin 0 scoreGap 0 scoreGapNoncan -8 scoreGapGCAG -4 scoreGapATAC -8 scoreStitchSJshift 1 scoreGenomicLengthLog2scale -0.25 scoreDelBase -2 scoreDelOpen -2 scoreInsOpen -2 scoreInsBase -2 seedSearchLmax 0 seedSearchStartLmax 50 seedSearchStartLmaxOverLread 1 seedPerReadNmax 1000 seedPerWindowNmax 50 seedNoneLociPerWindow 10 seedMultimapNmax 10000 seedSplitMin 12 alignIntronMin 21 alignIntronMax 0 alignMatesGapMax 0 alignTranscriptsPerReadNmax 10000 alignSJoverhangMin 5 alignSJDBoverhangMin 3 alignSJstitchMismatchNmax 0 -1 0 0 alignSplicedMateMapLmin 0 alignSplicedMateMapLminOverLmate 0.66 alignWindowsPerReadNmax 10000 alignTranscriptsPerWindowNmax 100 alignEndsType Local alignSoftClipAtReferenceEnds Yes alignEndsProtrude 0 ConcordantPair alignInsertionFlush None peOverlapNbasesMin 0 peOverlapMMp 0.01 chimSegmentMin 0 chimScoreMin 0 chimScoreDropMax 20 chimScoreSeparation 10 chimScoreJunctionNonGTAG -1 chimMainSegmentMultNmax 10 chimJunctionOverhangMin 20 chimOutType Junctions chimFilter banGenomicN chimSegmentReadGapMax 0 chimMultimapNmax 0 chimMultimapScoreRange 1 chimNonchimScoreDropMin 20 chimOutJunctionFormat 0 sjdbFileChrStartEnd - sjdbGTFfile - sjdbGTFchrPrefix - sjdbGTFfeatureExon exon sjdbGTFtagExonParentTranscript transcript_id sjdbGTFtagExonParentGene gene_id sjdbGTFtagExonParentGeneName gene_name sjdbGTFtagExonParentGeneType gene_type gene_biotype sjdbOverhang 100 sjdbScore 2 sjdbInsertSave Basic varVCFfile - waspOutputMode None quantMode - quantTranscriptomeBAMcompression 1 quantTranscriptomeBan IndelSoftclipSingleend twopass1readsN 18446744073709551615 twopassMode None soloType None soloCBstart 1 soloUMIstart 17 soloCBlen 16 soloUMIlen 10 soloBarcodeReadLength 1 soloCBwhitelist - soloStrand Forward soloOutFileNames Solo.out/ genes.tsv barcodes.tsv matrix.mtx matrixSJ.mtx matrixGeneFull.mtx soloFeatures Gene soloUMIdedup 1MM_All ##### Command Line: STAR --genomeDir /no_backup/jvalcarcel/emancini/anieto/02_INDEX/GRCm38 ––outSJfilterReads Unique --outSAMtype BAM SortedByCoordinate --runThreadN 10 --outFileNamePrefix SRR5264070 ––outFilterMismatchNoverReadLmax 0.04 ––outFilterMismatchNmax 999 ––outFilterMultimapNmax 1 ––alignEndsType Extend5pOfRead1 --outReadsUnmapped Fastx --quantMode GeneCounts --readFilesCommand zcat --readFilesIn SRR5264070_1.fastq.gz SRR5264070_2.fastq.gz ##### Initial USER parameters from Command Line: outFileNamePrefix SRR5264070 ###### All USER parameters from Command Line: genomeDir /no_backup/jvalcarcel/emancini/anieto/02_INDEX/GRCm38 ~RE-DEFINED outSAMtype BAM SortedByCoordinate ~RE-DEFINED runThreadN 10 ~RE-DEFINED outFileNamePrefix SRR5264070 ~RE-DEFINED outReadsUnmapped Fastx ~RE-DEFINED quantMode GeneCounts ~RE-DEFINED readFilesCommand zcat ~RE-DEFINED readFilesIn SRR5264070_1.fastq.gz SRR5264070_2.fastq.gz ~RE-DEFINED ##### Finished reading parameters from all sources ##### Final user re-defined parameters-----------------: runThreadN 10 genomeDir /no_backup/jvalcarcel/emancini/anieto/02_INDEX/GRCm38 readFilesIn SRR5264070_1.fastq.gz SRR5264070_2.fastq.gz readFilesCommand zcat outFileNamePrefix SRR5264070 outReadsUnmapped Fastx outSAMtype BAM SortedByCoordinate quantMode GeneCounts ------------------------------- ##### Final effective command line: STAR --runThreadN 10 --genomeDir /no_backup/jvalcarcel/emancini/anieto/02_INDEX/GRCm38 --readFilesIn SRR5264070_1.fastq.gz SRR5264070_2.fastq.gz --readFilesCommand zcat --outFileNamePrefix SRR5264070 --outReadsUnmapped Fastx --outSAMtype BAM SortedByCoordinate --quantMode GeneCounts ##### Final parameters after user input--------------------------------: versionGenome 2.7.1a parametersFiles - sysShell - runMode alignReads runThreadN 10 runDirPerm User_RWX runRNGseed 777 genomeDir /no_backup/jvalcarcel/emancini/anieto/02_INDEX/GRCm38 genomeLoad NoSharedMemory genomeFastaFiles - genomeChainFiles - genomeSAindexNbases 14 genomeChrBinNbits 18 genomeSAsparseD 1 genomeSuffixLengthMax 18446744073709551615 genomeFileSizes 0 genomeConsensusFile - readFilesType Fastx readFilesIn SRR5264070_1.fastq.gz SRR5264070_2.fastq.gz readFilesPrefix - readFilesCommand zcat readMatesLengthsIn NotEqual readMapNumber 18446744073709551615 readNameSeparator / inputBAMfile - bamRemoveDuplicatesType - bamRemoveDuplicatesMate2basesN 0 limitGenomeGenerateRAM 31000000000 limitIObufferSize 150000000 limitOutSAMoneReadBytes 100000 limitOutSJcollapsed 1000000 limitOutSJoneRead 1000 limitBAMsortRAM 0 limitSjdbInsertNsj 1000000 limitNreadsSoft 18446744073709551615 outFileNamePrefix SRR5264070 outTmpDir - outTmpKeep None outStd Log outReadsUnmapped Fastx outQSconversionAdd 0 outMultimapperOrder Old_2.4 outSAMtype BAM SortedByCoordinate outSAMmode Full outSAMstrandField None outSAMattributes Standard outSAMunmapped None outSAMorder Paired outSAMprimaryFlag OneBestScore outSAMreadID Standard outSAMmapqUnique 255 outSAMflagOR 0 outSAMflagAND 65535 outSAMattrRGline - outSAMheaderHD - outSAMheaderPG - outSAMheaderCommentFile - outBAMcompression 1 outBAMsortingThreadN 0 outBAMsortingBinsN 50 outSAMfilter None outSAMmultNmax 18446744073709551615 outSAMattrIHstart 1 outSAMtlen 1 outSJfilterReads All outSJfilterCountUniqueMin 3 1 1 1 outSJfilterCountTotalMin 3 1 1 1 outSJfilterOverhangMin 30 12 12 12 outSJfilterDistToOtherSJmin 10 0 5 10 outSJfilterIntronMaxVsReadN 50000 100000 200000 outWigType None outWigStrand Stranded outWigReferencesPrefix - outWigNorm RPM outFilterType Normal outFilterMultimapNmax 10 outFilterMultimapScoreRange 1 outFilterScoreMin 0 outFilterScoreMinOverLread 0.66 outFilterMatchNmin 0 outFilterMatchNminOverLread 0.66 outFilterMismatchNmax 10 outFilterMismatchNoverLmax 0.3 outFilterMismatchNoverReadLmax 1 outFilterIntronMotifs None outFilterIntronStrands RemoveInconsistentStrands clip5pNbases 0 clip3pNbases 0 clip3pAfterAdapterNbases 0 clip3pAdapterSeq - clip3pAdapterMMp 0.1 winBinNbits 16 winAnchorDistNbins 9 winFlankNbins 4 winAnchorMultimapNmax 50 winReadCoverageRelativeMin 0.5 winReadCoverageBasesMin 0 scoreGap 0 scoreGapNoncan -8 scoreGapGCAG -4 scoreGapATAC -8 scoreStitchSJshift 1 scoreGenomicLengthLog2scale -0.25 scoreDelBase -2 scoreDelOpen -2 scoreInsOpen -2 scoreInsBase -2 seedSearchLmax 0 seedSearchStartLmax 50 seedSearchStartLmaxOverLread 1 seedPerReadNmax 1000 seedPerWindowNmax 50 seedNoneLociPerWindow 10 seedMultimapNmax 10000 seedSplitMin 12 alignIntronMin 21 alignIntronMax 0 alignMatesGapMax 0 alignTranscriptsPerReadNmax 10000 alignSJoverhangMin 5 alignSJDBoverhangMin 3 alignSJstitchMismatchNmax 0 -1 0 0 alignSplicedMateMapLmin 0 alignSplicedMateMapLminOverLmate 0.66 alignWindowsPerReadNmax 10000 alignTranscriptsPerWindowNmax 100 alignEndsType Local alignSoftClipAtReferenceEnds Yes alignEndsProtrude 0 ConcordantPair alignInsertionFlush None peOverlapNbasesMin 0 peOverlapMMp 0.01 chimSegmentMin 0 chimScoreMin 0 chimScoreDropMax 20 chimScoreSeparation 10 chimScoreJunctionNonGTAG -1 chimMainSegmentMultNmax 10 chimJunctionOverhangMin 20 chimOutType Junctions chimFilter banGenomicN chimSegmentReadGapMax 0 chimMultimapNmax 0 chimMultimapScoreRange 1 chimNonchimScoreDropMin 20 chimOutJunctionFormat 0 sjdbFileChrStartEnd - sjdbGTFfile - sjdbGTFchrPrefix - sjdbGTFfeatureExon exon sjdbGTFtagExonParentTranscript transcript_id sjdbGTFtagExonParentGene gene_id sjdbGTFtagExonParentGeneName gene_name sjdbGTFtagExonParentGeneType gene_type gene_biotype sjdbOverhang 100 sjdbScore 2 sjdbInsertSave Basic varVCFfile - waspOutputMode None quantMode GeneCounts quantTranscriptomeBAMcompression 1 quantTranscriptomeBan IndelSoftclipSingleend twopass1readsN 18446744073709551615 twopassMode None soloType None soloCBstart 1 soloUMIstart 17 soloCBlen 16 soloUMIlen 10 soloBarcodeReadLength 1 soloCBwhitelist - soloStrand Forward soloOutFileNames Solo.out/ genes.tsv barcodes.tsv matrix.mtx matrixSJ.mtx matrixGeneFull.mtx soloFeatures Gene soloUMIdedup 1MM_All ---------------------------------------- Input read files for mate 1, from input string SRR5264070_1.fastq.gz -rw-r--r-- 1 emancini CRG_Lab_Juan_Valcarcel 1672000667 ene 29 14:17 SRR5264070_1.fastq.gz readsCommandsFile: exec > "SRR5264070_STARtmp/tmp.fifo.read1" echo FILE 0 zcat "SRR5264070_1.fastq.gz" Input read files for mate 2, from input string SRR5264070_2.fastq.gz -rw-r--r-- 1 emancini CRG_Lab_Juan_Valcarcel 1725385212 ene 29 14:17 SRR5264070_2.fastq.gz readsCommandsFile: exec > "SRR5264070_STARtmp/tmp.fifo.read2" echo FILE 0 zcat "SRR5264070_2.fastq.gz" WARNING: --limitBAMsortRAM=0, will use genome size as RAM limit for BAM sorting Finished loading and checking parameters Reading genome generation parameters: ### STAR --runMode genomeGenerate --runThreadN 10 --genomeDir GRCm38 --genomeFastaFiles genome.fa --outFileNamePrefix GRCm38 --sjdbGTFfile genes.gtf --sjdbOverhang 99 ### GstrandBit=32 versionGenome 2.7.1a ~RE-DEFINED genomeFastaFiles genome.fa ~RE-DEFINED genomeSAindexNbases 14 ~RE-DEFINED genomeChrBinNbits 18 ~RE-DEFINED genomeSAsparseD 1 ~RE-DEFINED sjdbOverhang 99 ~RE-DEFINED sjdbFileChrStartEnd - ~RE-DEFINED sjdbGTFfile genes.gtf ~RE-DEFINED sjdbGTFchrPrefix - ~RE-DEFINED sjdbGTFfeatureExon exon ~RE-DEFINED sjdbGTFtagExonParentTranscripttranscript_id ~RE-DEFINED sjdbGTFtagExonParentGene gene_id ~RE-DEFINED sjdbInsertSave Basic ~RE-DEFINED genomeFileSizes 2776514406 22235024810 ~RE-DEFINED Genome version is compatible with current STAR Number of real (reference) chromosomes= 22 1 10 130694993 0 2 11 122082543 130809856 3 12 120129022 252968960 4 13 120421639 373293056 5 14 124902244 493879296 6 15 104043685 618921984 7 16 98207768 722993152 8 17 94987271 821297152 9 18 90702639 916455424 10 19 61431566 1007419392 11 1 195471971 1069023232 12 2 182113224 1264582656 13 3 160039680 1446772736 14 4 156508116 1606942720 15 5 151834684 1763704832 16 6 149736546 1915748352 17 7 145441459 2065694720 18 8 129401213 2211184640 19 9 124595110 2340683776 20 MT 16299 2465464320 21 X 171031299 2465726464 22 Y 91744698 2636906496 --sjdbOverhang = 99 taken from the generated genome Started loading the genome: Wed Feb 9 16:48:26 2022 Genome: size given as a parameter = 2776514406 SA: size given as a parameter = 22235024810 SAindex: size given as a parameter = 1 Read from SAindex: pGe.gSAindexNbases=14 nSAi=357913940 nGenome=2776514406; nSAbyte=22235024810 GstrandBit=32 SA number of indices=5390309044 Shared memory is not used for genomes. Allocated a private copy of the genome. Genome file size: 2776514406 bytes; state: good=1 eof=0 fail=0 bad=0 Loading Genome ... done! state: good=1 eof=0 fail=0 bad=0; loaded 2776514406 bytes SA file size: 22235024810 bytes; state: good=1 eof=0 fail=0 bad=0 Loading SA ... done! state: good=1 eof=0 fail=0 bad=0; loaded 22235024810 bytes Loading SAindex ... done: 1565873619 bytes Finished loading the genome: Wed Feb 9 16:48:54 2022 Processing splice junctions database sjdbN=240490, pGe.sjdbOverhang=99 alignIntronMax=alignMatesGapMax=0, the max intron size will be approximately determined by (2^winBinNbits)*winAnchorDistNbins=589824 Opening the file: SRR5264070_STARtmp//Unmapped.out.mate1.thread0 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate2.thread0 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate1.thread1 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate2.thread1 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate1.thread2 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate2.thread2 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate1.thread3 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate2.thread3 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate1.thread4 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate2.thread4 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate1.thread5 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate2.thread5 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate1.thread6 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate2.thread6 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate1.thread7 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate2.thread7 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate1.thread8 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate2.thread8 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate1.thread9 ... ok Opening the file: SRR5264070_STARtmp//Unmapped.out.mate2.thread9 ... ok Created thread # 1 Created thread # 2 Created thread # 3 Created thread # 4 Created thread # 5 Created thread # 6 Created thread # 7 Created thread # 8 Starting to map file # 0 mate 1: SRR5264070_1.fastq.gz mate 2: SRR5264070_2.fastq.gz Created thread # 9 BAM sorting: 254526 mapped reads BAM sorting bins genomic start loci: 1 0 80826685 2 1 3291209 3 1 58933349 4 1 80467634 5 1 110513187 6 2 31335547 7 2 89637811 8 3 13232318 9 3 77759468 10 4 24471544 11 4 77200953 12 5 76233439 13 5 98950536 14 6 20692445 15 6 96066868 16 7 33638684 17 7 39846913 18 7 39847643 19 7 57080793 20 8 35582732 21 9 5133930 22 9 46300738 23 10 42710243 24 10 123597976 25 10 180803122 26 11 34776263 27 11 105697334 28 11 155670724 29 12 19249461 30 12 89467608 31 12 146732193 32 13 82290930 33 13 129849014 34 14 5782659 35 14 88963113 36 14 129899389 37 15 3201427 38 15 72436412 39 15 125139719 40 16 30527242 41 16 105558737 42 17 15519970 43 17 84596788 44 18 15315955 45 18 57954210 46 18 108837699 47 19 7598 48 20 74271679 Thread #7 end of input stream, nextChar=-1 Completed: thread #2 Completed: thread #1 Completed: thread #6 Completed: thread #5 Completed: thread #7 Completed: thread #3 Completed: thread #4 Completed: thread #8 Completed: thread #0 Joined thread # 1 Joined thread # 2 Joined thread # 3 Joined thread # 4 Joined thread # 5 Joined thread # 6 Joined thread # 7 Joined thread # 8 Completed: thread #9 Joined thread # 9 Feb 09 16:54:13 ..... started sorting BAM Max memory needed for sorting = 177675365 ALL DONE!