So I have an RNA-seq dataset that should have been sequenced as three separate replicates (per condition) that got sequenced as one sample per condition for reasons. The data are paired-end.
This, I suppose, means (in theory) that I would not be able to do any differential gene expression analysis using this data. But I would love for the data not to go to waste.
Is there any way now that I could take the reads files R1 and R2 for each condition and create:
pseudoreplicate1_R1 & pseudoreplicate1_R2
pseudoreplicate2_R1 & pseudoreplicate2_R2
pseudoreplicate3_R1 & pseudoreplicate3_R2
such that each the reads from R1 and R2 are uniformly (evenly?) represented in their subsets pseudoreplicate[1-3]_R[1-2]?
If not, are there any other suggestions for dealing with this scenario? In specific, I must mention that this data is/was intended for a time-series/time-course analysis.
0 answers
No answers yet.
Log in to answer this question.
-
Expected number of raw reads in fq.gz files in PE sequencing
written by Lada •
-
Bulk RNA-seq time course data with unknown sample sizes?
written by Dunois •
-
FastQC report analysis for scRNA seq data
written by aimanbarki •
-
Dropseq, GSE124872 and structure of the reads in the FastQ files
written by tmms •
-
BBDuk quality filtering not producing expected results
written by rebeliscu •
-
TruSeq sequencing output
written by Rogerio Ribeiro •
-
high duplication levels in R1 reads of 10x Genomics samples
written by Assa Yeroslaviz
-
Offline Tools or packages (prefer R and bioconductor packeges) for miRNA-seq time series counts dat…
written by bioinforupesh2009.au •
-
MiSeq: R1 vs. R2 reads show big quality difference - is that normal?
written by mschmid
-
Quality Drop In The Reverse Read Of Paired End Library (Illumina Hiseq 1000)
written by ancient_learner
They will simply become technical/sequencing replicates after split. Nothing is going to change that.