This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Viewing the sai file generated using bwa

Hi. I am working on RNA-Seq analysis which was partially done by a researcher who left the team. The only mapping files that I have is the sai files but not the bam files. May I know if there is a way to view the content of the sai files? Thanks in advance.

bwa mapping
• 2,832 views • link
updated by noahaus • noahaus • written by Louis Kok • Louis Kok

Do you still have the fastq files?

• 1 views • link
written by h.mon h.mon

Yes. I have the fastq but do not know what the reference it was mapped to.

• 1 views • link
written by Louis Kok • Louis Kok

Isn't it better to start all over? If you don't know the genome used, I wonder what else you don't know. Do you know if the reads are good quality? Do you know if adapters have been trimmed? Were the reads mapped to a reference genome or transcriptome?

These sai files are output from BWA, which is not an appropriate program to map RNAseq to a reference genome.

• 1 views • link
written by h.mon h.mon

are you sure that files sai files? They are index files of sam files. See if you have .sam files along with .sai. If you have .sai files only, you can not proceed further.

• 1 views • link
written by cpad0112 cpad0112

1 answer

I'd checkout the BWA command samse. If you have transcript data such as fastq, then you can save the output of the command as a SAM which can be later on made into a BAM.

• 1 views • link
written by noahaus • noahaus

Log in to answer this question.

More posts like this