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 SRR5264057 ––outFilterMismatchNoverReadLmax 0.04 ––outFilterMismatchNmax 999 ––outFilterMultimapNmax 1 ––alignEndsType Extend5pOfRead1 --outReadsUnmapped Fastx --quantMode GeneCounts --readFilesCommand zcat --readFilesIn SRR5264057_1.fastq.gz SRR5264057_2.fastq.gz ##### Initial USER parameters from Command Line: outFileNamePrefix SRR5264057 ###### 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 SRR5264057 ~RE-DEFINED outReadsUnmapped Fastx ~RE-DEFINED quantMode GeneCounts ~RE-DEFINED readFilesCommand zcat ~RE-DEFINED readFilesIn SRR5264057_1.fastq.gz SRR5264057_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 SRR5264057_1.fastq.gz SRR5264057_2.fastq.gz readFilesCommand zcat outFileNamePrefix SRR5264057 outReadsUnmapped Fastx outSAMtype BAM SortedByCoordinate quantMode GeneCounts ------------------------------- ##### Final effective command line: STAR --runThreadN 10 --genomeDir /no_backup/jvalcarcel/emancini/anieto/02_INDEX/GRCm38 --readFilesIn SRR5264057_1.fastq.gz SRR5264057_2.fastq.gz --readFilesCommand zcat --outFileNamePrefix SRR5264057 --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 SRR5264057_1.fastq.gz SRR5264057_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 SRR5264057 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 SRR5264057_1.fastq.gz -rw-r--r-- 1 emancini CRG_Lab_Juan_Valcarcel 872311919 ene 29 14:07 SRR5264057_1.fastq.gz readsCommandsFile: exec > "SRR5264057_STARtmp/tmp.fifo.read1" echo FILE 0 zcat "SRR5264057_1.fastq.gz" Input read files for mate 2, from input string SRR5264057_2.fastq.gz -rw-r--r-- 1 emancini CRG_Lab_Juan_Valcarcel 901696076 ene 29 14:07 SRR5264057_2.fastq.gz readsCommandsFile: exec > "SRR5264057_STARtmp/tmp.fifo.read2" echo FILE 0 zcat "SRR5264057_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:39:16 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:39:58 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: SRR5264057_STARtmp//Unmapped.out.mate1.thread0 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate2.thread0 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate1.thread1 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate2.thread1 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate1.thread2 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate2.thread2 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate1.thread3 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate2.thread3 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate1.thread4 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate2.thread4 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate1.thread5 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate2.thread5 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate1.thread6 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate2.thread6 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate1.thread7 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate2.thread7 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate1.thread8 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate2.thread8 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate1.thread9 ... ok Opening the file: SRR5264057_STARtmp//Unmapped.out.mate2.thread9 ... ok Created thread # 1 Created thread # 2 Created thread # 3 Created thread # 4 Created thread # 5 Starting to map file # 0 mate 1: SRR5264057_1.fastq.gz mate 2: SRR5264057_2.fastq.gz Created thread # 6 Created thread # 7 Created thread # 8 Created thread # 9 BAM sorting: 255460 mapped reads BAM sorting bins genomic start loci: 1 0 80828333 2 0 128814630 3 1 58044008 4 1 79127657 5 1 109011957 6 2 27337720 7 2 78351068 8 2 108593098 9 3 34990725 10 3 104188084 11 4 54593569 12 5 10478896 13 5 85375374 14 5 99779125 15 6 42955850 16 7 6650150 17 7 35835115 18 7 45569036 19 8 23846047 20 8 75001364 21 9 9984894 22 10 24613735 23 10 82900481 24 10 161037504 25 11 24666667 26 11 71543773 27 11 130681394 28 11 179178551 29 12 55477123 30 12 97861386 31 13 32505277 32 13 103565078 33 13 134473552 34 14 30063247 35 14 108575217 36 14 135858260 37 15 18465568 38 15 86432508 39 15 140746957 40 16 45083564 41 16 114705753 42 17 22577355 43 17 85539929 44 18 21505709 45 18 62373626 46 18 108944108 47 19 8078 48 20 93253979 Thread #2 end of input stream, nextChar=-1 Completed: thread #6 Completed: thread #8 Completed: thread #2 Completed: thread #7 Completed: thread #3 Completed: thread #4 Completed: thread #0 Completed: thread #9 Completed: thread #5 Completed: thread #1 Joined thread # 1 Joined thread # 2 Joined thread # 3 Joined thread # 4 Joined thread # 5 Joined thread # 6 Joined thread # 7 Joined thread # 8 Joined thread # 9 Feb 09 16:43:09 ..... started sorting BAM Max memory needed for sorting = 99482292 ALL DONE!