If the the two duplicated regions are similar in sequence (90% identity) and very closely localized to each other (100bp distant), how likely do they show identical sequences due to assembly error? Thanks very much!
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Originally posted by maimaiti2008 View PostIf the the two duplicated regions are similar in sequence (90% identity) and very closely localized to each other (100bp distant), how likely do they show identical sequences due to assembly error? Thanks very much!
It would probably depend on the assembler you use, and on the length of your reads. Different assemblers give different results. If you have long reads or paired reads that can span the duplicated regions, then the assemblers are less likely to make an 'error'.
-
To me, 90% identity doesn't sound like too much for de novo assemblers to tell the difference.
If you're using a k-mer based approach, you have to have an exact k-1 identity to join neighboring k-mers. So, if you have 1/10 bases different and you're assembling at k=35, your kmers covering these genes will have, on average, 3.5 nucleotide differences. So unless you have regions of near 100% identity, where the assembler might think the more divergent region between them is just a bubble, you should be fine. And even then, paired end reads or longer kmers would likely take that problem away entirely.
Comment
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist on Modified Bases...-
Channel: Articles
Yesterday, 07:01 AM -
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
55 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
51 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
45 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
||
Started by seqadmin, 04-04-2024, 09:00 AM
|
0 responses
55 views
0 likes
|
Last Post
by seqadmin
04-04-2024, 09:00 AM
|
Comment