Filter R1 and its matching R2 based on kmer sequence
Hello!
I am working with RNA-seq data, from paired end sequencing (10x precisely, and sequenced on a NextSeq platform).
Basically I have a list of sequences (15 bp long) which are present in only a subset of my R1 reads. What I would like to do is to filter the R2 reads which match this criteria. Do you know how I could achieve this?
Many thanks for your help.
1 answer
You should be able to use bbduk.sh from BBMap suite in filter mode. There is a guide available here.
Log in to answer this question.
More posts like this
-
How to Analyze 10x Genomics Data with Full-Length R1 and R2 (100bp Each)?
written by zbidav •
-
How to run TRUST4 for BCR/TCR detection using 10X data?
written by O.rka
-
Extract the true single-cell RNA sequencing reads for running SAHMI
written by 573704669 •
-
Some questions on RAD paired-end sequencing
written by john •
-
STARsolo mapping for paired-ended 3' 10X reads
written by tvibhaps •
-
Help Understanding PairedEnd read merging
written by m.radz •
-
Paired end data with R1 demultiplexed but not R2
written by singlegenome •
-
high duplication levels in R1 reads of 10x Genomics samples
written by Assa Yeroslaviz
-
When paired end reads doesn't join, can i consider them single-end in metagenomic analysis?
written by fernandalpcosta •
-
Tool to append Illumina Read 1 to Read 2 for downstream demultiplexing
written by achamess
Please give a representative example of input and expected output. R1 in 10X are barcodes and UMIs which are not present in R2. It is not clear what you aim to do.