This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Error : samtools tview

Hi: I have a sorted and merged BAM file, and I would like to check the alignment, but the result is wrong. So, where is the problem? Thanks a lot.

The wrong result：

$ samtools tview KPGP-00001_picard_sort.bam /DataBase/Human/hg19/hg19.fa
1         11        21        31        41        51        61        71        
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
samtools visualization alignment
• 2,620 views • link
written by 1191106787 • 1191106787

samtools tview is not the most powerful genome browser available, I would recommend IGV, IGB or Tablet (among others) to get a richer view of alignments.

How is the result wrong? You mean the long stretch of Ns? Try going to some position - press g then write down the position, something like chr1:156000.

See some further discussion here.

• 0 views • link
written by h.mon h.mon

all the result was only 71bp and END, but obvious human genome greater, so where is the other sequence?

• 0 views • link
written by 1191106787 • 1191106787

1 answer

I found it was because the wrongly samtools installed conda, and I installed it with source code and got the right result. I have submitted this issue at https://github.com/conda/conda/issues/6191 .

• 0 views • link
written by 1191106787 • 1191106787

Log in to answer this question.

More posts like this