bamurai 0.2.2

A package for splitting reads BAM files into smaller fragments.

Bamurai

A Python package for splitting reads BAM files into smaller fragments.

Description

Bamurai is a command-line tool that splits reads in BAM files to a rough target length and outputs a FASTQ file. This is useful for testing the performance of bioinformatics tools on long-reads data at different read lengths.

Currently there are 3 main features of Bamurai:

1. Splitting reads in a BAM file to a target length
2. Splitting reads in a BAM file to a target number of pieces per read
3. Getting statistics from a BAM or FASTQ(.gz) file

The split command splits reads into a target length, each read will be split into fragments as close to the target length as possible. Reads shorter than the target length will not be split.

The divide command splits reads into a target number of pieces, each read will be split into the number of pieces specified. A further minimum length can be specified to ensure that reads are not split if the resultant fragments are less than the minimum length.

The stats command will output the following information by default:

Statistics for input.bam:
  Total reads: 8160
  Average read length: 30638
  Throughput: 250006998
  N50: 82547

It can be used with the --tsv argument to output the statistics in a tab-separated format for computational analysis.

file_name       total_reads     avg_read_len    throughput      n50
input.bam      8160    30638   250006998       82547

Installation

pip install git+https://github.com/Shians/Bamurai.git

Usage

To get help on the command-line interface and list available commands

bamurai --help

To get help on a specific command

bamurai <command> --help

Splitting reads to target size

To split a file into 10,000 bp reads

bamurai split input.bam --target-length 10000 --output output.fastq

To create a gzipped output file

bamurai split input.bam --target-length 10000 | gzip > output.fastq.gz

Dividing reads into a target number of pieces

To divide reads into 2 pieces

bamurai divide input.bam --num_fragments 2 --output output.fastq

To divide reads into 2 pieces unless resultant fragments are less than 1000 bp

bamurai divide input.bam --num_fragments 2 --min_length 1000 --output output.fastq

Getting statistics from a BAM or FASTQ file

To get stats from a BAM file

bamurai stats input.bam

To get stats from a FASTQ file or Gzipped FASTQ file

bamurai stats input.fastq
bamurai stats input.fastq.gz