using IGV to visualize RNA seq alignment to transcriptome
0
0
Entering edit mode
6.5 years ago

Hello,

I created a de novo transcriptome assembly using Trinity. I used Trinity's align_and_estimate_abundance.pl to align my reads using bowtie. I now want to visualize my reads and I have heard a lot about using IGV, however, do I need to sort and index my bowtie.bam file before loading it into IGV?

I tried samtools sort bowtie.bam and I am not sure what is supposed to happen to my .bam file but here is a picture of what my terminal is doing and it does not look right

If anyone can shed some light on this that would be much appreciated, Thank you, Nikelle

rna seq igv rna alignment • 5.1k views
ADD COMMENT
0
Entering edit mode

Yes you do need to sort and index your alignment file before loading it into IGV. If you are using the new samtools (v.1.3) you need to provide a filename for the sorted version by using the -o option with samtools sort. Then samtools index the sorted file.

ADD REPLY
0
Entering edit mode

Thanks Genomax.

I followed your instructions. When I went to upload the file into IGV, there was an error saying:

ERROR: Loading SAM/BAM index files are not supported: /Users/nikellepetrillo/Documents/Senior/RichardsonResearch/Research Meeting/bowtiesorted.bai Load the SAM or BAM file directly.

I am not sure what the problem is

ADD REPLY
0
Entering edit mode

As the error suggests, load the BAM file directly. It looks like you are loading a BAI file, not BAM file. However, the loaded BAM file should have a corresponding BAI file in the same directory.

ADD REPLY
0
Entering edit mode

Hi Nikelle,

I'm not sure if you figured this out. But here's my input...

The extension if you index the BAM files using samtools is usually .bam.bai. I've found heaps of resources online saying that you need to upload a coordinate-sorted, and indexed BAM file. However, I've found that IGV will only take a BAM file without the .bai extension. So I use samtools to sort the BAM file by coordinates, then upload that one to IGV.

samtools sort <filename.bam> <filename_sorted>

During the process samtools will create a heap of intermediate files, then join them all together to give a single file at the end.

ADD REPLY

Login before adding your answer.

Traffic: 1684 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6