This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Human brain samples (BAM files) differential gene expression

I have been handed 8 sets of RNA-seq data. These were already aligned, and are in BAM format. The alignment was done by an outside company long time ago, and i have no idea what algorithm they used, and not even which reference sequence, all I know is that it's human brain samples.

Any ideas how i could proceed with these files (using cufflinks for example) to compare differential gene expression? i have cufflinks already up and running.

Thanks in advance.

rna-seq
• 1,234 views • link
written by nermze86 • nermze86

Thank you very much, problem solved :) !

• 0 views • link
written by nermze86 • nermze86

If an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted.
Upvote|Bookmark|Accept

• 0 views • link
written by GenoMax GenoMax

1 answer

Hi, you can view the header of the BAM-files with samtools:

samtools view -H your.bam

The header tells you the sequence ID of the genome which was used for the mapping and also the alignment tool that was used.

• 0 views • link
written by Mathias_H • Mathias_H

Adding on this, typically the last lines of the header contain the commands used for alignment. In any case you can always format bam back to fastq with e.g. samtools fastq and realign / quantify against a reference of choice. That is what I would do to be sure it was done properly (not saying the company is incompetent, but simply to be sure that you can confidently stand up for the analysis).

• 0 views • link
written by ATpoint ATpoint

Log in to answer this question.

More posts like this