Tweet
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
-
The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...11-06-2024, 07:24 PM -
Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...10-18-2024, 07:11 AM
Comment