Question: Subsampling aligned SAM/BAM files for RNA-seq analysis
0
gravatar for mohammedtoufiq91
5 months ago by
mohammedtoufiq9140 wrote:

Hi,

I would like to know if there are any specific tools to sub-sample aligned SAM or BAM files. What type of files could be considered for the sub-sampling sorted or unsorted SAM or BAM files.

Thank you, Toufiq

ADD COMMENTlink modified 5 months ago by WouterDeCoster40k • written 5 months ago by mohammedtoufiq9140
1
gravatar for WouterDeCoster
5 months ago by
Belgium
WouterDeCoster40k wrote:

You did not tell us anything about the aim of your analysis, which is probably important. For expression analysis I'd recommend subsampling the read counts. That is probably more efficient than subsampling alignments. First package that comes to mind is this R package: subseq

ADD COMMENTlink written 5 months ago by WouterDeCoster40k

We have generated a set of RNA-seq samples. These are human PE samples with read length of 150bp. We would like to identify the total number of genes that could be detected by reducing the number of reads to 1/2, 1/3rd and 1/4th. After this, identify the differentially expressed genes between the two conditions i.e. before/after treatment.

ADD REPLYlink written 5 months ago by mohammedtoufiq9140

There is a historical post. Check it out here.

ADD REPLYlink written 5 months ago by ivivek_ngs4.8k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1807 users visited in the last hour