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.
-
-
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. -
Ah, I guess I'll just have to build up my Perl skills. Probably a good idea, anyway. And yes, the consensus sequences are the same size, so it shouldn't be too difficult.
Thanks for your help!Comment
-
Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?
Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating...10-14-2026, 09:55 AM -
Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.
We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing...
Comment