viewing read map data from Subread/Rsubread

[akh22]

Hi all,

I've been trying to view mapping BAM outputs from Rsubread by using IGV and IGB but these viewers don't appear to like BAM files created by Rsubread. Can anybody recommend other map viewers or any pointers to to accomplish this ?

Thanks in advance.

[shi]

Reads in BAM output from Subreads are in the same order as in your input FASTQ file. You need to sort them by chromosomal locations before feeding them to IGV or IGB.

[Brian Bushnell]

Reads in BAM output from Subreads are in the same order as in your input FASTQ file. You need to sort them by chromosomal locations before feeding them to IGV or IGB.

More precisely, you need to sort (by coordinate, as stated) AND index them. You can do both with Samtools or Picard.

[akh22]

I did sort and index BAM files but still IGB viewer can't open it.

[Brian Bushnell]

I have never used IGB, just IGV.

Anyway...

It would help if you described the command lines you used, and posted the header of the bam file (if it's reasonably short), and the specific error messages you get.

[GenoMax]

Are you getting a specific error or is nothing happening when you try to open the file?

I do not know about IGB but for IGV you may need to adjust the memory assigned to IGV at startup (or use the appropriate web launcher).