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
-
NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification, detection...Today, 06:46 AM -
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
Comment