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 human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...-
Channel: Articles
02-24-2025, 06:31 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 03-03-2025, 01:15 PM
|
0 responses
169 views
0 likes
|
Last Post
by seqadmin
03-03-2025, 01:15 PM
|
||
Started by seqadmin, 02-28-2025, 12:58 PM
|
0 responses
256 views
0 likes
|
Last Post
by seqadmin
02-28-2025, 12:58 PM
|
||
Started by seqadmin, 02-24-2025, 02:48 PM
|
0 responses
636 views
0 likes
|
Last Post
by seqadmin
02-24-2025, 02:48 PM
|
||
Started by seqadmin, 02-21-2025, 02:46 PM
|
0 responses
265 views
0 likes
|
Last Post
by seqadmin
02-21-2025, 02:46 PM
|
Comment