SEQanswers - STAR read Mapping

SEQanswers (http://seqanswers.com/forums/index.php)

- RNA Sequencing (http://seqanswers.com/forums/forumdisplay.php?f=26)

- - STAR read Mapping (http://seqanswers.com/forums/showthread.php?t=87776)

NDUFB11

02-19-2019 09:52 AM

STAR read Mapping

Hello everybody,

I kindly need your help again ...

I want to map my reads on the genome using STAR but for some reason I get an error that I don't know how to solve

When I type the command line here is what I get:

Code:

Patrizios-MBP:processing patrizio$ STAR --genomeDir genome/star/ --runThreadN 4 --readFilesIn clipped/SRR359063_1.fastq.zip clipped/ SRR359063_2.fastq.zip --outFileNamePrefix mapping/star/SRR359063. --outSAMattributes All --outSAMtype BAM SortedByCoordinate --readFileCommand gunzip

EXITING: FATAL INPUT ERROR: unrecoginzed parameter name "readFileIn" in input "Command-Line-Initial"
SOLUTION: use correct parameter name (check the manual)

Feb 19 18:49:56 ...... FATAL ERROR, exiting[/B]

I also tried to process unzipped files and remove the command to execute each single input file --readFileCommand, but still I can't go further,

I have a MacOs terminal,

thank you

Patrizio

GenoMax

02-19-2019 09:57 AM

STAR requires a lot of RAM depending on genome size. e.g. human genome will need ~30G of free RAM which you may not have available on your Mac. If that is true, you will not be able to do this on your local mac.

NDUFB11

02-19-2019 10:04 AM

Do you think that this is the reason why I have this error? this is a nematode genome so I thought the process would have been faster than human

I will scan the command, I found an error probably I didn't type correctly the command

later I will also try with topHat, it should be feasible because requires less RAM

thank you

GenoMax

02-19-2019 10:15 AM

Oops. I should have looked at the error carefully. You are correct.

Edit: Found the problem. You have a space between clipped and file name here (clipped/ SRR359063_2.fastq.zip ). Remove that.

NDUFB11

02-19-2019 10:26 AM

It's working :D :D :D !!!! holeee'!

have a look

Code:

Patrizios-MBP:processing patrizio$ STAR --genomeDir genome/star/ --runThreadN 4 --readFilesIn clipped/SRR359063_1.fastq clipped/SRR359063_2.fastq --outFileNamePrefix mapping/star/SRR359063. --outSAMattributes All --outSAMtype BAM SortedByCoordinate 

Feb 19 19:22:44 ..... started STAR run

Feb 19 19:22:44 ..... loading genome

Feb 19 19:22:49 ..... started mapping

Feb 19 19:23:50 ..... finished mapping

Feb 19 19:23:50 ..... started sorting BAM

Feb 19 19:23:54 ..... finished successfully

Patrizios-MBP:processing patrizio$

I can go further with the process !

Thank you GenoMax :)

at some point I will ship a bottle of wine to you :D

-P

GenoMax

02-19-2019 10:34 AM

Looks good!

yoderto

10-22-2019 10:28 AM

Just a quick comment-- STAR does come with RAM limiting arguments: --limitGenomeGenerateRAM, --limitIObufferSize , --genomeSAsparseD

All times are GMT -8. The time now is 12:36 AM.