dalexander / DumbView

A PacBio genome browser for command line nerds

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

DumbView: a PacBio genome browser for command-line nerds

Screenshot

Requirements: pbcore, GenomicConsensus

Install: python setup.py install

Use:

  • View alignments to a window of the reference: dumbview aligned_reads.cmp.h5 [-r reference.fasta] -w<contigName>:<start>-<end>
  • View alignments around all variant calls in a variants.gff file: dumbview variants.gff
  • To enable "unaligned" mode (useful for eyeballing insertion variants) use -u

Notes:

  • DumbView expects a FASTA index (.fai) file in the same directory as the reference FASTA. Make it with samtools faidx reference.fasta.
  • Coordinates are 1-based, in agreement with the "GFF" convention and most other genome browsers.
  • The consensus track shown at the bottom presently uses the Quiver "NoQVs" model, so it may not be the same consensus as you get from running Quiver. This will be fixed soon.

This is an unsupported tool and the code is horrendous. But I find it useful.

David Alexander, 2013

About

A PacBio genome browser for command line nerds


Languages

Language:Python 83.2%Language:Perl 12.1%Language:Perl 6 4.2%Language:Makefile 0.5%