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: 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: Igv

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
ADD COMMENT
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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

ADD REPLY
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.

ADD REPLY
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....

ADD REPLY
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.

ADD REPLY
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 1115 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