Stringtie --merge with PE and SE files
0
0
Entering edit mode
5.4 years ago

I'm trying to create a merged gtf from TCGA bam files with Stringtie --merge. But some of the files use SE reads and some - PE reads. Are those bam file compatible and can be used together to generate a new gtf? I tried to find this info in the manual and online but couldn't.

Thank you!

assembly bam • 1.2k views
ADD COMMENT
0
Entering edit mode

How is it that some are single- while others paired-end? I had assumed that TCGA data was almost exclusively paired end. Have you tried to run the command?

ADD REPLY
0
Entering edit mode

Most of the data there is PE but there are some samples that are SE. Some of them are even from the same individuals, it's weird. I'm trying to compare tumor to normal tissue from the same individual and there are few of them with PE and SE reads samples as different biological samples. Read lengths also vary so I don't even know if it would make sense to use samples from the same person with different read lengths or PE and SE but I was wondering if I can combine all of them in one gtf. I haven't tried running it but it will probably work, would it be correct to combine SE and PE bams though?

ADD REPLY
0
Entering edit mode

I see. Most of them should be PE, though, right? Can you achieve a large sample n with just the PE samples? I processed both SE and PE RNA-seq samples in the past with no apparent technical / batch effects.

ADD REPLY

Login before adding your answer.

Traffic: 2040 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6