Snap aligner to Igv?
0
0
Entering edit mode
20 months ago
newbie • 0

Hi, I try to use Snap aligner (http://snap.cs.berkeley.edu/) and view the result in Igv for some viruses. Here is the Snap output:

After that, I use samtools to index and sort the .sam file for .bam and .bam.bai. (They are not 0 bytes.) I added the links to Igv (js) and this is what I see:

The files I used are working in Bowtie for example, so they cannot be the problem. What am I missing here?

snap igv aligner • 564 views
0
Entering edit mode

Just to confirm. You first sorted and then indexed BAM files, correct?

0
Entering edit mode

Yes. I followed the exact same routine as with Bowtie.

0
Entering edit mode

do you see any read with samtools view in the region you're looking with IGV ?

0
Entering edit mode

Yes, but like this: MN908947.3_Here_is_the_ref_genome_name... This could be the problem? In other bam file which worked fine there is no ref genome name there.

0
Entering edit mode

Yes, but like this: MN908947.3_Here_is_the_ref_genome_name... This could be the problem

I don't understand your answer, but try to write the correct contig name in IGV....

0
Entering edit mode

Snap made my contig like: MN908947.3_Ecoli and in Igv there is MN908947.3 so they dont match. maybe thats why this is Snap beta.

0
Entering edit mode

Why would SNAP change your contig names? Did it do that during the index creation?