Hi all, I am using HISAT2 and trying to write unmapped reads from paired-end data and the flag --un is resulting in empty files. I know they are not all mapping so does anyone have any recommendations to write those reads to a file?
-
-
Hi I had the same problem, you need to use --un-conc for paired end reads (write paired-end reads that fail to align concordantly to file(s) at <path>). There is an explanation of the relevant SAM flags in the help manual for Hisat2. -
Thanks for the reply! I was hoping there was a means of writing out only those mates that did not map anywhere rather than all reads that did not map concordantly, but it is a good place to start I suppose.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