output unmapped file in STAR alignment
Hi All. I am using STAR alignment to remove rRNA reads, and collect unmapped reads for further alignment. My command line is: STAR --genomeDir ./genome/index/rRNA/ --readFilesIn RA1_R1.fastq RA1_R2.fastq --clip3pNbases 12 --clip5pNbases 50 --outFileNamePrefix rRNA1 --outReadsUnmapped Fastx
However, the output unmapped file mix read1 and read2 together, not in a separated R1.fastq and R2.fastq. Does anyone know how to set parameters? Thanks a lot.
• 3,881 views
•
link
written
by
linjc.xmu •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Aligning fastq files using STAR
written by kcarey •
-
For loop for aligning multiple paired end fastq files
written by kcarey •
-
for paired-end unmapped reads after star mapping
written by tvibhaps •
-
star --readFilesIn error when using unmapped.out.mate1/2 files
written by kokoko •
-
RNAseq STAR mapping with clipping
written by asta.tranholm •
-
snakemake errors and threads
written by Assa Yeroslaviz
-
Alignment loop using STAR
written by lamia_203
-
cufflinks large bundles warning
written by Sam
-
How can I only took the mapped ones from STAR ?
written by rekren •
-
STAR aligner fails to produce unmapped file
written by dec986
Which version of STAR are you using?
STAR_2.5.3a, thanks.
Maybe this is fixed in the latest release (2.6.1a).