Does anyone know if BWA is handling masked sequences properly? I'm having a problem where I find some reads aligned to genome gaps. I have tried converting the N's to X's, but still encountering the same problem.
-
-
This is a known issue. Ambiguous bases are converted to random sequences. If you use bwa-short, the XN tag tells you how many "N" in the alignment. If you use bwasw, you do not get this information, but you may find such hits by checking you reference genome.
The chance of a good alignment arising from random hits is vanishingly small. Just ignore those hits. -
thanks
Thanks!
I just want to add that the N's in my sequence are not for ambigious bases, they represent a sequence gap of unknown size.. So basically its a long stretch [ATGC]+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN[ATGC]+Last edited by rand; 11-09-2009, 08:08 AM.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