How can I only took the mapped ones from STAR ?
Hello,
I wanted to obtain only the "mapped" reads as an output of the STAR. I forgot to delete "--outSAMunmapped Within " and all of my mapped output files also contain all unmapped reads,too. Data size is huge and re-mapping properly wiil take so much time...
How can I fix the BAM files which were supposed to be mapped reads but also include unmapped reads in it ?
#!/bin/bash
mkdir /mnt/data/Toxo_scan/GBR_Male/ToxoMap
while read -r line
do
mkdir ToxoMap/$line
echo $line" -> Running STAR - Toxo now"
STAR --runThreadN 12 --alignIntronMax 1 --outSAMunmapped Within --outSAMtype BAM SortedByCoordinate --genomeDir "/mnt/data/Toxo_scan/toxo_genome" --readFilesIn "/mnt/data/Toxo_scan/GBR_Male/sickle/"$line"_1_clean.fastq" "/mnt/data/Toxo_scan/GBR_Male/sickle/"$line"_2_clean.fastq" --outFileNamePrefix "/mnt/data/Toxo_scan/GBR_Male/ToxoMap/"$line"_" --outReadsUnmapped Fastx
done
• 7,789 views
•
link
updated
by
andrew.j.skelton73
•
written
by
rekren •
•
written
by
rekren •
2 answers
To get your mapped reads:
samtools view -b -F 4 alignment.bam > mapped.bam
To get your unmapped reads:
samtools view -b -f 4 alignment.bam > unmapped.bam
The -f filters based on bitwise samflags. Flag 4 (or 0x0004) is for unmapped reads. -F does the inverse.
• 0 views
•
link
written
by
andrew.j.skelton73
This is what you're looking for:
samtools view -f 0x2 foo.bam > foo.filtered.bam
It will keep only the proper pairs in your alignment (if paired end). Proper pairs: correct orientation and within insert size.
• 0 views
•
link
written
by
Matteo Schiavinato
Log in to answer this question.
More posts like this
-
STAR pair-end gene count discrepancy vs. RSEM
written by anon •
-
unmapped output of the STAR
written by Wakala •
-
trouble with 'too short' reads using STAR. each mate maps fine, but the paired end mapping does not
written by james.alan.gregory •
-
star --readFilesIn error when using unmapped.out.mate1/2 files
written by kokoko •
-
snakemake errors and threads
written by Assa Yeroslaviz
-
STAR alignment error: quantmod not found
written by prabin.dm
-
STAR alignment error: ERROR in input reads
written by prabin.dm
-
Very low unique mapped reads - STAR alignment
written by jsl •
-
STAR alignment error- Segmentation fault (core dumped)
written by EagleEye
-
STAR aligner fails to produce unmapped file
written by dec986
Is the presence of unmapped reads causing issues? If you were going to do counts etc those reads will be ignored.
They were causing the issues because, I was supposed to process only the mapped reads in the workflow. For my case, presence of the unmapped reads in the mapped reads is like having false positive results into the true positive results...
Out of curiosity what was the downstream workflow component that was causing a problem? Most programs should understand reads that are unmapped easily.
The issue was not related with downstream workflow components, it was about my supervisor's satisfaction :) When I came with millions of undesired "unmapped" reads as a result of "mapped" reads, I was in the position of not doing the task well. That's all :)
Aha. No logical solution would suffice in that case.