Extract primary STAR mapped reads to specific chromosomes
Hi,
I have RNA-seq reads which I mapped with STAR and I'd like to extract from the bam file the primary alignments that map to chr1, chr2, and chr3.
I know the second part is achieved by:
samtools view -b -L <my_chrs_bed_file> <my_bam_file>
So I'm mostly interested in the first part, and how to pipe the two parts.
Thanks
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Samtools mapped and unmapped flags, what is the difference between the flags?
written by DNAngel
-
Samtools flagstat
written by rakszewska •
-
Retrieving multiple regions (from a list) in a bam file and outputing as a table
written by dodausp
-
How to download Grch38p7 database associated DBSNP151
written by clear.choi •
-
How can I get this sorted by ascending order? (when the number and character is together)
written by woojoy14 •
-
extract region from the bam file
written by 9521ljh •
-
Cufflinks with large BAM file
written by qudrat
-
how to randomly select intervals of different lengths (not records) from a BED file
written by jing.mengrabbit •
-
samtools view by chromosome region chr:1000-1001000 creates "unknown reference name" error
written by ariel.balter
-
Tutorial: How To Separate Illumina Based Strand Specific Rna-Seq Alignments By Strand
written by Istvan Albert
I don't follow, what's the 'first' part you want to achieve? You have one question: how to extract reads aligning to chr1 chr2 etc, and you have the solution.