Hello friendly people,
I have a question about how to properly and efficiently deal with repeat regions within coding sequence when mapping RNA-Seq reads.
I have illumina 29bp reads (after trimming off a multiplex index), and am doing RNA-Seq with tissue from Arabidopsis thaliana, accession Col-0. I am currently using bowtie/tophat to perform my initial mapping.
To avoid signal from one gene leaking across to the signal from another, unrelated gene that happens to share one or more 29-mers (roughly, in actuality, 29 with as many as 2 mismatches), I throw away reads that map to multiple locations.
That is all fine and dandy, except that sometimes the repeats come from within the same gene. Transmembrane proteins, for example, are famous (in my book) for having long stretches of exact repeats. Throwing these away is undesirable, since throwing away clean signal always increases noise.
I can easily map my reads to all possible locations, and then write a perl script to sift through them afterwards, keeping multiple-mappers whose mapping locations are all within the same gene. (My organism is very, very well annotated.) But this seems like a basic function. Does anyone know of an option/switch in bowtie that allows what I am looking for (that I just can't find)?
Cheers!
~Rachel
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
Latest Articles
Collapse
-
by seqadmin
Innovations in next-generation sequencing technologies and techniques are driving more precise and comprehensive exploration of complex biological systems. Current advancements include improved accessibility for long-read sequencing and significant progress in single-cell and 3D genomics. This article explores some of the most impactful developments in the field over the past year.
Long-Read Sequencing
Long-read sequencing has...-
Channel: Articles
Yesterday, 01:49 PM -
-
by seqadmin
The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben MartÃnez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...-
Channel: Articles
11-06-2024, 07:24 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 09:29 AM
|
0 responses
14 views
0 likes
|
Last Post
by seqadmin
Yesterday, 09:29 AM
|
||
Started by seqadmin, Yesterday, 09:06 AM
|
0 responses
11 views
0 likes
|
Last Post
by seqadmin
Yesterday, 09:06 AM
|
||
Started by seqadmin, Yesterday, 08:03 AM
|
0 responses
11 views
0 likes
|
Last Post
by seqadmin
Yesterday, 08:03 AM
|
||
Started by seqadmin, 11-22-2024, 07:36 AM
|
0 responses
65 views
0 likes
|
Last Post
by seqadmin
11-22-2024, 07:36 AM
|