CRISPR/Cas9 screen analysis: reads R1 and R2 mixed
Hi! I want to analyse data from a CRISPR/Cas9 screen (control vs. treatment) and I'm using Mageck (https://sourceforge.net/projects/mageck/). The sequencing was performed using Illumina (paired-end).
The problem is that I've noticed that in R1 fastq files the half of the reads containing the sgRNA are R2 (in R2 files there are R1 reads as well). Should I consider these reads in the sgRNA count?.
• 2,629 views
•
link
written
by
Swimming bird •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Problems with zero count values and mapping % in MaGeCK
written by DragosV •
-
Problem with mapping rate while using Mageck to process CRISPR screen data
written by Luwell •
-
Approaches to validate hits from CRISPR-screen in mice
written by GLG •
-
News: Course - Introduction to CRISPR
written by Physalia-courses
-
how to make BBmerge produce same strand fastq file
written by shangguandong1996 •
-
MAGeCK normalization; during count, test or both?
written by sonark •
-
Extract gRNA sequence using cutadapt
written by Swimming bird •
-
Difference in R1 and R2 FASTQ files from RNA seq paired end
written by jaynishshah •
-
Trim Paired-end Fastq Files
written by yuabrahamliu •
-
Illumina paired-end reads R1 and R2 mixed together?
written by lvogel •
How did that happen? Always best to go back and get original data files in cases where you are in doubt.
I don't know because the sequencing was commissioned but I think it could be due to the ligation: https://seekdeep.brown.edu/illumina_paired_info.html
Are you saying that you have short inserts (these being sgRNA) so R1/R2 are likely to overlap (i.e. there is no mixing per se)? You could just use R1 read or look into a tool like
bbmerge.shfrom BBTools that can merge R1/R2 reads. You can then trim adapters and then count the consensus sequence produced (or align to a reference and then count).No, there is a mixing because R1 files contain R2 reads.
If they are truly mixed and you want to separate the reads (and if they have standard Illumina headers then you can do something like):