This is a test version of Biostars. For the public version, visit https://www.biostars.org.
454 Sequencing With Bwa Mapping

Well, I indexed my ref.fasta of hg19 (using bwa index -a bwtsw database.fasta) , and mapped my fastq file using bwa, using ( bwa bwasw database.fasta long_read.fastq > aln.sam )... converted my sam to bam.... using samtools and then I tried to visualize in UCSC custom track and IGV...and appeared nothing like this screenshot : http://www.4freeimagehost.com/show.php?i=776af632ebca.png

so why is it happening ??? thanks

bwa igv samtools visualization
• 4,283 views • link
updated by 2184687-1231-83- 2184687-1231-83- • written by Gabriel Wajnberg Gabriel Wajnberg

2 answers

This is actually very counterintuitive (and burned me as well the first time I used IGV) - if you are zoomed out too far then IGV displays nothing ... there should be a more prominent indicator that the content is not shown (you have that small text there that asks you to zoom in)

• 1 views • link
written by Istvan Albert Istvan Albert

Just to confirm, if you do a "samtools view mybam.bam | head" and look at the coordinates for the first 10 lines (third and fourth column), then zoom in to those coordinates in IGV? Is it still blank?

• 0 views • link
written by 2184687-1231-83- 2184687-1231-83-

Log in to answer this question.

More posts like this