I'm working with a tetraploid soybean species, and I want to determine which DNA-seq reads exist on only one of the two diploid progenitor genomes. I have two consensus fasta files (one for each diploid genome) and want to map reads to them to determine which reads map to only one genome based on mapping scores. I already have a pipeline for accomplishing this, but I'm running into a problem. One consensus file has 30 million more Ns than the other. Reads that should map equally to both species favor one because it has fewer Ns, and therefore a higher mapping score. Is there a way to replace A/C/G/T's in one consensus file with an N where the other consensus file has an N? The consensus sequences were made by mapping RNA-seq and DNA-seq reads to the Glycine max genome, so they should be comparable.
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
I'm not aware of a ready-made tool to do this job but it sounds like something you could do with Perl. I'm sure there are online tutorials you could use otherwise I found this book useful when first starting (http://shop.oreilly.com/product/9780596000806.do)
Are both consensus sequences the same length for each chromosome/scaffold? Doesn't really matter if they aren't, would just make the perl script to do this job easier.
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