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 SRR5264067 ––outFilterMismatchNoverReadLmax 0.04 ––outFilterMismatchNmax 999 ––outFilterMultimapNmax 1 ––alignEndsType Extend5pOfRead1 --outReadsUnmapped Fastx --quantMode GeneCounts --readFilesCommand zcat --readFilesIn SRR5264067_1.fastq.gz SRR5264067_2.fastq.gz ##### Initial USER parameters from Command Line: outFileNamePrefix SRR5264067 ###### 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 SRR5264067 ~RE-DEFINED outReadsUnmapped Fastx ~RE-DEFINED quantMode GeneCounts ~RE-DEFINED readFilesCommand zcat ~RE-DEFINED readFilesIn SRR5264067_1.fastq.gz SRR5264067_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 SRR5264067_1.fastq.gz SRR5264067_2.fastq.gz readFilesCommand zcat outFileNamePrefix SRR5264067 outReadsUnmapped Fastx outSAMtype BAM SortedByCoordinate quantMode GeneCounts ------------------------------- ##### Final effective command line: STAR --runThreadN 10 --genomeDir /no_backup/jvalcarcel/emancini/anieto/02_INDEX/GRCm38 --readFilesIn SRR5264067_1.fastq.gz SRR5264067_2.fastq.gz --readFilesCommand zcat --outFileNamePrefix SRR5264067 --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 SRR5264067_1.fastq.gz SRR5264067_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 SRR5264067 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 SRR5264067_1.fastq.gz -rw-r--r-- 1 emancini CRG_Lab_Juan_Valcarcel 1011407581 ene 29 14:08 SRR5264067_1.fastq.gz readsCommandsFile: exec > "SRR5264067_STARtmp/tmp.fifo.read1" echo FILE 0 zcat "SRR5264067_1.fastq.gz" Input read files for mate 2, from input string SRR5264067_2.fastq.gz -rw-r--r-- 1 emancini CRG_Lab_Juan_Valcarcel 1044084049 ene 29 14:08 SRR5264067_2.fastq.gz readsCommandsFile: exec > "SRR5264067_STARtmp/tmp.fifo.read2" echo FILE 0 zcat "SRR5264067_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:44:08 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:44:25 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: SRR5264067_STARtmp//Unmapped.out.mate1.thread0 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate2.thread0 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate1.thread1 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate2.thread1 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate1.thread2 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate2.thread2 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate1.thread3 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate2.thread3 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate1.thread4 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate2.thread4 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate1.thread5 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate2.thread5 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate1.thread6 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate2.thread6 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate1.thread7 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate2.thread7 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate1.thread8 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate2.thread8 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate1.thread9 ... ok Opening the file: SRR5264067_STARtmp//Unmapped.out.mate2.thread9 ... ok Created thread # 1 Created thread # 2 Created thread # 3 Created thread # 4 Created thread # 5 Created thread # 6 Starting to map file # 0 mate 1: SRR5264067_1.fastq.gz mate 2: SRR5264067_2.fastq.gz Created thread # 7 Created thread # 8 Created thread # 9 BAM sorting: 254232 mapped reads BAM sorting bins genomic start loci: 1 0 80252411 2 0 128882192 3 1 58949129 4 1 80403497 5 1 115270708 6 2 38823834 7 2 83533886 8 2 111671433 9 3 58394489 10 4 18268964 11 4 66307882 12 5 30825736 13 5 89373733 14 6 16882573 15 6 75756704 16 7 26267283 17 7 39846915 18 7 90346105 19 8 61259784 20 9 7427316 21 10 16102736 22 10 55720313 23 10 132465941 24 10 180572431 25 11 34776162 26 11 106693799 27 11 150636348 28 12 8554883 29 12 88218375 30 12 127244291 31 13 59715819 32 13 126435159 33 14 5538017 34 14 89236549 35 14 129899184 36 15 17751893 37 15 85929608 38 15 135726504 39 16 29174572 40 16 98715415 41 17 3513090 42 17 70868600 43 17 123025405 44 18 49502883 45 18 90102372 46 19 1348 47 19 8813 48 20 74300714 Thread #8 end of input stream, nextChar=-1 Completed: thread #7 Completed: thread #3 Completed: thread #8 Completed: thread #1 Completed: thread #6 Completed: thread #4 Completed: thread #2 Completed: thread #0 Joined thread # 1 Joined thread # 2 Joined thread # 3 Joined thread # 4 Completed: thread #5 Joined thread # 5 Joined thread # 6 Joined thread # 7 Joined thread # 8 Completed: thread #9 Joined thread # 9 Feb 09 16:48:25 ..... started sorting BAM Max memory needed for sorting = 100045378 ALL DONE!