Hi,
I have bacterial RNA-seq data from paired-end reads (2x75). I'm interested in sense/antisense differential expression.
I would like to know if makes sense to merge R1 and R2 into one read (if overlaps) or maybe work only with 75 bp reads (R2 ones).
Thanks, Bernardo
1 answer
Why bother? Presuming that you used a stranded protocol when making the libraries, you can just align things as pairs and tell whichever program you use to do counting to observe strandedness of pairs. Both htseq-count and featureCounts can do that.
To get at the antisense transcription, one could just run the counters again and just tell them to use the other read as the strand reference (e.g., running htseq-count -s yes ... and then htseq-count -s reverse ...).
Any decent program that does counting will count pairs as a single unit, rather than individually, so yes it is the same.
Log in to answer this question.
-
Racon command line with paired-ends Illumina reads
written by A_heath •
-
Extract the overlapping paired-end reads
written by ATRX
-
Should I merge the Paired-end RNA-seq reads for BLAT for intron identification
written by praasu •
-
Trim Paired-end Fastq Files
written by yuabrahamliu •
-
RNA seq workflow
written by bioinfo •
-
Tool to append Illumina Read 1 to Read 2 for downstream demultiplexing
written by achamess
-
Assembly: collapsed reads vs overlapping paired-end reads?
written by palu •
-
How to download a large amount of NGS metadata for studies?
written by RafaelMP
-
RNA-seq reads mapping to paralogs with high MAPQ. How to ensure this?
written by biotech
-
Ncbi Standards For Genome Annotation After Contig Reduction
written by biotech
If your experiment was conducted so that the fragment size allowed pairs to overlap, merging them would make a lot of sense. So, finding out your fragment size would be the first thing..
200 bp but maybe shorter. The thing is that I need to convert paired-end to single-end data anyway in order to quantify antisense transcription too.