Question: using IGV to visualize RNA seq alignment to transcriptome
0
gravatar for nikelle.petrillo
2.5 years ago by
Providence College, Providence, RI
nikelle.petrillo90 wrote:

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

igv rna seq rna alignment • 2.2k views
ADD COMMENTlink written 2.5 years ago by nikelle.petrillo90

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 REPLYlink written 2.5 years ago by genomax59k

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 REPLYlink written 2.5 years ago by nikelle.petrillo90

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 REPLYlink modified 2.5 years ago • written 2.5 years ago by igor7.1k

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 REPLYlink modified 2.1 years ago • written 2.1 years ago by ando.kelli40
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: 1099 users visited in the last hour