This is a test version of Biostars. For the public version, visit https://www.biostars.org.
loaded a sorted bam file in IGV but I can't see anything

Hi,

I have genomic SRR file (DNA-seq) for a pseudomonas genome. I used a file containing corresponding cDNA sequence and nucleotide sequence of my gene of interest to build the index by using the salmon tool. Then, used samtools to create SRR.bam file 

salmon quant -p 10 -i  SRR_index -l A -1 SRR_1 .fastq -2 SRR_2.fastq -z -o SRR_qt | samtools view -bs > SRR.bam
samtools sort SRR.bam SRR.sorted.bam
samtools index SRR.sorted.bam

Then load the .fasta files ( the index.fasta.fai was in the same directory), *.gff files of the to the IGV. I also loaded .sorted.bam file in IGV but I can't see anything.

Please give me a solution.

assembly bam igv
• 4,606 views • link
updated by Biostar Biostar • written by saadleeshehreen saadleeshehreen

Hello,

but I can't see anything

A screenshot would be useful to understand what you see and what not. Try to zoom in, as IGV doesn't show a coverage track or read information if the region of the current view is to large.

fin swimmer

• 1 views • link
written by finswimmer finswimmer

Hello again,

there would be grey bars (those would be the aligned reads) if you:

  1. have a valid alignment file
  2. zoom in enough

I don't know salmon. But after a quick look at the manual it doesn't look like the output is an alignment file, isn't it? I guess you have to choose this before with another program like bwa or bbmapor ...

fin swimmer

• 1 views • link
written by finswimmer finswimmer

Zoom in to the maximum perhaps

• 1 views • link
written by rse rse

Thanks. I use BWA and now the alignment look like that.. I didn't split the window. But I get two alignment near the position of my gene of interest. What does it mean?

Untitled_2

• 1 views • link
written by saadleeshehreen saadleeshehreen

1 answer

Salmon is a pseudoaligner so clearly, the output you are getting does not have mapping information. Try --writeMappings argument to include the alignment information.

Ref: https://salmon.readthedocs.io/en/latest/salmon.html

• 1 views • link
written by Arup Ghosh Arup Ghosh

To add information , i never be able to read the BAM (convert from the SAM output ) with IGV with this option.

• 1 views • link
written by Titus Titus

Log in to answer this question.

More posts like this