Question: STAR aligner fails to produce unmapped file
0
gravatar for dec986
2.4 years ago by
dec986190
United States
dec986190 wrote:

Hello,

I'm trying to generate an unmapped bam file for further processing using STAR alignment.

but I've tried numerous iterations of this command

/usr/bin/time STAR --outSAMtype BAM SortedByCoordinate --outSAMattributes All --genomeDir /home/con/GENE_DATA/mm10 --readFilesIn TN_R2_cat.fastq
  /usr/bin/time STAR --outSAMtype BAM SortedByCoordinate --outSAMattributes All --genomeDir /home/con/GENE_DATA/mm10 --readFilesIn TN_R2_cat.fastq 2> STAR.err &
  cat STAR.err 
  echo "/usr/bin/time STAR --outSAMtype BAM SortedByCoordinate --outSAMattributes All --genomeDir /home/con/GENE_DATA/mm10 --readFilesIn TN_R2_cat.fastq 2> STAR.err &" >> notes.txt 
  /usr/bin/time STAR --outSAMtype BAM SortedByCoordinate --outSAMattributes All --genomeDir /home/con/GENE_DATA/mm10 --readFilesIn TN_R2_cat.fastq --outSAMunmapped
  /usr/bin/time STAR --outSAMtype BAM SortedByCoordinate --outSAMattributes All --genomeDir /home/con/GENE_DATA/mm10 --readFilesIn TN_R2_cat.fastq --outSAMunmapped='unmapped'
  /usr/bin/time STAR --outSAMtype BAM SortedByCoordinate --outSAMattributes All --genomeDir /home/con/GENE_DATA/mm10 --readFilesIn TN_R2_cat.fastq --outSAMunmapped='unmapped'
  /usr/bin/time STAR --outSAMtype BAM SortedByCoordinate --outSAMattributes All --genomeDir /home/con/GENE_DATA/mm10 --outReadsUnmapped --readFilesIn TN_R2_cat.fastq 
  /usr/bin/time STAR --outSAMtype BAM SortedByCoordinate --outSAMattributes All --genomeDir /home/con/GENE_DATA/mm10 --outReadsUnmapped='fastx' --readFilesIn TN_R2_cat.fastq 
  /usr/bin/time STAR --outSAMtype BAM SortedByCoordinate --outSAMattributes All --genomeDir /home/con/GENE_DATA/mm10 --outReadsUnmapped='fastx' --runThreadN 2 --readFilesIn TN_R2_cat.fastq 
  /usr/bin/time STAR --outSAMtype BAM SortedByCoordinate --outSAMattributes All --genomeDir /home/con/GENE_DATA/mm10 --outReadsUnmapped='fastx' --readFilesIn TN_R2_cat.fastq 
  /usr/bin/time STAR --outSAMtype BAM SortedByCoordinate --outSAMattributes All --genomeDir /home/con/GENE_DATA/mm10 --outReadsUnmapped='fastx' --readFilesIn TN_R2_cat.fastq 
  /usr/bin/time STAR --outSAMtype BAM SortedByCoordinate --outSAMattributes All --genomeDir /home/con/GENE_DATA/mm10 --outReadsUnmapped --readFilesIn TN_R2_cat.fastq 
  /usr/bin/time STAR --outSAMtype BAM SortedByCoordinate --outSAMattributes All --genomeDir /home/con/GENE_DATA/mm10 --outReadsUnmapped='unmapped' --readFilesIn TN_R2_cat.fastq 
  STAR   --genomeDir /home/con/GENE_DATA/mm10   --readFilesIn TN_R2_cat.fastq      --outReadsUnmapped unmapped   --outSAMtype BAM   SortedByCoordinate      --outSAMattrib
  STAR --genomeDir /home/con/GENE_DATA/mm10 --readFilesIn TN_R2_cat.fastq --outReadsUnmapped unmapped SortedByCoordinate  --outSAMattributes All
  STAR --genomeDir /home/con/GENE_DATA/mm10 --readFilesIn TN_R2_cat.fastq --outReadsUnmapped unmapped SortedByCoordinate  --outSAMattributes All
  STAR   --genomeDir /home/con/GENE_DATA/mm10   --readFilesIn TN_R2_cat.fastq --outSAMattributes All --outSAMunmapped=unmapped
  STAR   --genomeDir /home/con/GENE_DATA/mm10   --readFilesIn TN_R2_cat.fastq --outSAMattributes All --outReadsUnmapped=fast

but I cannot get the output Unmapped.out.mate1/2 that the manual promises How can I get STAR to produce the unmapped file? What am I getting wrong from the manual?

star • 2.0k views
ADD COMMENTlink modified 2.4 years ago by WouterDeCoster38k • written 2.4 years ago by dec986190
4
gravatar for WouterDeCoster
2.4 years ago by
Belgium
WouterDeCoster38k wrote:

The flag to set is just --outReadsUnmapped Fastx.

ADD COMMENTlink written 2.4 years ago by WouterDeCoster38k

thanks WouterDeCoster!

ADD REPLYlink written 2.4 years ago by dec986190
1
gravatar for Manvendra Singh
2.4 years ago by
Manvendra Singh2.0k
Berlin, Germany
Manvendra Singh2.0k wrote:

try not to use SortedByCoordinate, get the output in sam format

try using this argument in your command and lets see if you get

--outReadsUnmapped Fastx
ADD COMMENTlink written 2.4 years ago by Manvendra Singh2.0k
1

--outSAMtype BAM SortedByCoordinate works together with --outReadsUnmapped Fastx

ADD REPLYlink written 2.4 years ago by WouterDeCoster38k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1239 users visited in the last hour