This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How does IGB provide strand specific output

I want to know on what basis IGB separates the reads into forward and reverse strand and is there anyway to generate strand specific bed file output for a set of coordinates.

rna-seq

In what format is your input ?

I wish to know when i give sorted .bam file as input file. How does IGB separate the reads into + and - strand? can i use it to look for antisense transcription?

1 answer

How does IGB separate the reads into + and - strand?

The alignment file (bam or sam) countains flags that tell you if the read is mapped on the reverse or forward strand. More info on this related post.

Can I use it to look for antisense transcription?

Short answer, yes. For instance if you use featureCounts to quantify the number of reads per genes, you can play with the -s option to count reads mapping in anti-sense of genes.

Is there anyway to generate strand specific bed file output for a set of coordinates ?

Yes, with bedtools bamtobed for instance.

How to generate 12 column bed file using bedtools or coverageBed?

This question is completely unrelated to the original one. Please submit as a new post with appropriate title. See this post for guidance.

Log in to answer this question.