Hi all,
I have RNASeq data for 20 different samples (human cell types and tissues). 10 of these samples are represented by single biological replicate, while, the other 10 samples are represented by >=2 biological replicates. In total I have 40 RNASeq fastq files mapped independently.
Now, the purpose is to compare the samples (cell types) against each other (not between the biological replicates of each sample). So, my question is: how can I deal with the biological replicates that I mapped? I have some options:
1) Take one (random) replicate per sample
2) Generate read count of mapped reads for all replicates per sample independently and then average the read counts
3) Merge the mapped reads (using bedtools merge) of the replicates
Any other suggestion?
I have RNASeq data for 20 different samples (human cell types and tissues). 10 of these samples are represented by single biological replicate, while, the other 10 samples are represented by >=2 biological replicates. In total I have 40 RNASeq fastq files mapped independently.
Now, the purpose is to compare the samples (cell types) against each other (not between the biological replicates of each sample). So, my question is: how can I deal with the biological replicates that I mapped? I have some options:
1) Take one (random) replicate per sample
2) Generate read count of mapped reads for all replicates per sample independently and then average the read counts
3) Merge the mapped reads (using bedtools merge) of the replicates
Any other suggestion?
Comment