Extracting the mapped sequence of a certain gene
1
0
Entering edit mode
9 months ago
randy.kang • 0

Hello, guys.

First, I feel embarrassed to ask this seemingly fundamental question for all of you and I am very unfamiliar with FASTQ/BAM handling..... Recently, we are trying to design a probe for RNA scope (or maybe possibly antibody as well in the future) that we found in the scRNA sequencing study recently. Is there anyway that we extract the sequence of the gene (probably partial or intron-region genes) which was recognized as our gene of interest?

Thanks a tone in advance guys!!

gene scRNA-seq BAM • 321 views
ADD COMMENT
1
Entering edit mode
9 months ago
Shred ▴ 960

Easy:

samtools view -L your_gene.bed your_alignment.bam

Where your_gene.bed is the coordinates of your gene as a BED and your_alignment.bam contains aligned sequences.

ADD COMMENT
0
Entering edit mode

Thank you so much!!!!

ADD REPLY

Login before adding your answer.

Traffic: 1387 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6