Error : samtools tview

Question

[1191106787]

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

[h.mon]

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.

[1191106787]

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

[Pierre Lindenbaum]

hu ? kamoulox ?

[Devon Ryan]

XRef: https://github.com/samtools/samtools/issues/728

Answer 1

[1191106787]

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 .