Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: Filtering paired reads when forward is all scRNA-seq barcode sequence
0
gravatar for connor.driscoll88
8 months ago by
connor.driscoll88 • 0
connor.driscoll88 • 0 wrote:

Hi all,

I have paired-end fastq files from a single-cell RNA seq experiment. The forward files consist of barcode reads (no usable sequence information) while the reverse reads contain transcript sequences. I'm looking for an efficient way to filter out rRNA reads, but it's a little tricky considering I can't align the forward reads. Is there an efficient way of doing this either within an aligner or with a quick downstream script that will result in rRNA-filtered and orphanless paired end (non-interleaved) fastq files?

Thanks.
rna-seq alignment scrna-seq • 317 views
ADD COMMENT • link •
modified 6 months ago by Biostar ♦♦ 20 • written 8 months ago by connor.driscoll88 • 0
1
gravatar for Pierre Lindenbaum
8 months ago by
Pierre Lindenbaum ♦ 122k
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum ♦ 122k wrote:

but it's a little tricky considering I can't align the forward reads.

so R1 won't be mapped but R2 will. And then filter with samtools

ADD COMMENT • link written 8 months ago by Pierre Lindenbaum ♦ 122k

Ah, I see that is quite simple! So the correct samtools flag be 9 (read paired (0x1) & mate unmapped (0x8))?

Edit: Okay, I think it should be 13 (read paired (0x1) & read unmapped (0x4) & mate unmapped (0x8)) since the read 1 is the "read" and read 2 is the "mate." Also, I want the mates which are unmapped.

Edit 2: If anyone is interested, here is what I used to convert sam output to bam, filter only for pairs where reverse read non-mapped, and converted back to two non-interleaved fastqs in a single command:

samtools view -bS input.sam | samtools view -u -f 13 - | samtools fastq -1 input.norRNA_R1.fastq.gz -2 input.norRNA_R2.fastq.gz -
ADD REPLY • link modified 8 months ago • written 8 months ago by connor.driscoll88 • 0
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 812 users visited in the last hour