Working with paired reads in BAM format
Project description
888888b. 88888888888
888 "88b 888
888 .88P 888
8888888K. 8888b. 88888b.d88b. 888 888 888 88888b. .d88b. 888d888
888 "Y88b "88b 888 "888 "88b 888 888 888 888 "88b d8P Y8b 888P"
888 888 .d888888 888 888 888 888 888 888 888 888 88888888 888
888 d88P 888 888 888 888 888 888 Y88b 888 888 d88P Y8b. 888
8888888P" "Y888888 888 888 888 888 "Y88888 88888P" "Y8888 888
888 888
Y8b d88P 888
"Y88P" 888
# Overview
Ability to work out the orientation and insert size of a paired read data file
Can estimate relative orientation and gap between pairs of contigs in the bam file (Useful for scaffolding)
# Installation
use pip:
```sh
bamtyper type bamfile.bam
Determining OT for BAM 'bamfile'
Orientation: IN Insert: 301, Stdev: 29
bamtyper links - Parse BAM files and get linking reads
Usage 1:
bamtyper links bamfile1.bam bamfile2.bam -c
contig3 , [ contig2 , 3 , SS , 58 ]
contig2 , [ contig3 , 3 , SS , 58 ] , [ contig1 , 4 , EE , 45 ]
contig1 , [ contig2 , 4 , EE , 45 ]
contig3 0.6206 0.5234
contig2 0.6558 0.0123
contig1 0.6523 0.5634
Where:
contig3 0.6206 0.5234
Reports 0.6206 fragments per base in bamfile1 and 0.5234 in bamfile2
If these were 100bp reads then this would imply
coverages of 62x and 52x respectively
# Administration
Project home page, info on the source tree, documentation, issues and how to contribute, see http://github.com/minillinim/BamTyper
This software is currently unpublished.
Copyright © 2012 Michael Imelfort. See LICENSE.txt for further details.
888 "88b 888
888 .88P 888
8888888K. 8888b. 88888b.d88b. 888 888 888 88888b. .d88b. 888d888
888 "Y88b "88b 888 "888 "88b 888 888 888 888 "88b d8P Y8b 888P"
888 888 .d888888 888 888 888 888 888 888 888 888 88888888 888
888 d88P 888 888 888 888 888 888 Y88b 888 888 d88P Y8b. 888
8888888P" "Y888888 888 888 888 888 "Y88888 88888P" "Y8888 888
888 888
Y8b d88P 888
"Y88P" 888
# Overview
Ability to work out the orientation and insert size of a paired read data file
Can estimate relative orientation and gap between pairs of contigs in the bam file (Useful for scaffolding)
# Installation
use pip:
```sh
Determining OT for BAM 'bamfile'
Orientation: IN Insert: 301, Stdev: 29
bamtyper links - Parse BAM files and get linking reads
Usage 1:
contig3 , [ contig2 , 3 , SS , 58 ]
contig2 , [ contig3 , 3 , SS , 58 ] , [ contig1 , 4 , EE , 45 ]
contig1 , [ contig2 , 4 , EE , 45 ]
contig3 0.6206 0.5234
contig2 0.6558 0.0123
contig1 0.6523 0.5634
Where:
contig3 0.6206 0.5234
Reports 0.6206 fragments per base in bamfile1 and 0.5234 in bamfile2
If these were 100bp reads then this would imply
coverages of 62x and 52x respectively
# Administration
Project home page, info on the source tree, documentation, issues and how to contribute, see http://github.com/minillinim/BamTyper
This software is currently unpublished.
Copyright © 2012 Michael Imelfort. See LICENSE.txt for further details.