Samview from bam file aligned to cDNA reference
Hello,
So I have a bam file from alignment to cDNA reference from Ensembl Hg38. Suppose transcript ID is T1. The length of the transcript is 15000. I want to extract the reads that aligned to T1 at position between 200-400. I think this is possible for the ordinary bam file which is aligned to human genome reference. I remember just specify chr1:pos1-pos2. How can I extract reads that correspond to a particular transcript ID and base position.
• 2,315 views
•
link
updated
by
Biostar
•
written
by
bharata1803
•
written
by
bharata1803
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Beagle: Index -l out of bounds for length 2
written by Tom •
-
How to Extract Specific Region on Bam file
written by santos48 •
-
No Sequencing Data at Low Positions
written by ccnn •
-
bam file coordinate system change
written by Gregor Rot
-
Consensus of multiple BED files
written by bruce.moran
-
Biological difference for selecting alignment reference to DNA vs cDNA
written by bharata1803
-
How to convert genomic position to transcript position?
written by biolab
-
Counting reads supporting a specific splice junction from bam file
written by Nicolas Rosewick
-
How To Get The Size Of Sv With Respect To Breakdancer Output
written by skm770 •
-
How To Interpret Breakdancer Output File ?
written by madkitty
if you know the position of T1 on the genome, the process is called "translation" ...