hyraxAbi: Modules for parsing, generating and manipulating AB1 files.

This is a package candidate release! Here you can preview how this package release will appear once published to the main package index (which can be accomplished via the 'maintain' link below). Please note that once a package has been published to the main package index it cannot be undone! Please consult the package uploading documentation for more information.

[maintain] [Publish]

This library provides functionality for parsing, modifying, writing and generating ABI files

Any AB1 file confirming to the standard at http://www6.appliedbiosystems.com/support/software_community/ABIF_File_Format.pdf should be supported.

This library also support generating a minimal ABI file from a FASTA input sequence.

A basic terminal application that can dump and generate AB1s is included. See https://github.com/hyraxbio/hyraxAbi/blob/master/app/Main.hs

See Examples directory


[Skip to Readme]

Properties

Versions 0.2.3.1
Change log None available
Dependencies base (>=4.7 && <5), binary (>=0.8.5.1 && <0.8.6.0), bytestring (>=0.10.8.2 && <0.10.9.0), directory (>=1.3.0.2 && <1.3.2.0), filepath (>=1.4.1.2 && <1.4.2.0), hscolour (>=1.24.4 && <1.25.0), hyraxAbi, pretty-show (>=1.6.16 && <1.7.0), protolude (>=0.2.2 && <0.2.3), text (>=1.2.3.0 && <1.2.4.0) [details]
License BSD-3-Clause
Copyright 2018 HyraxBio
Author HyraxBio
Maintainer andre@hyraxbio.co.za, andre@andrevdm.com
Category Bioinformatics
Home page https://github.com/hyraxbio/hyraxAbi/#readme
Source repo head: git clone https://github.com/hyraxbio/hyraxAbi
Uploaded by andrevdm at 2018-07-09T15:44:53Z

Modules

[Index]

Downloads

Maintainer's Corner

Package maintainers

For package maintainers and hackage trustees


Readme for hyraxAbi-0.2.3.1

[back to package description]

HyraxBio AB1 parser, writer and generator (beta 0.2)

This project contains

See http://www6.appliedbiosystems.com/support/software_community/ABIF_File_Format.pdf for a high level overview of the AB1 file format.

Building

Build with one of

Terminal app

Run with

Dump AB1

To dump an existing AB1 run

hyraxAbi-exe dump example.ab1

This will output the structure of the AB1 like this

Header { hName = "ABIF" , hVersion = 101 }
Directory
  { dTagName = "tdir"
  , dTagNum = 1
  , dElemTypeCode = 1023
  , dElemTypeDesc = "root"
  , dElemType = ElemRoot
  , dElemSize = 28
  , dElemNum = 13
  , dDataSize = 364
  , dDataOffset = 61980
  , dData = ""
  , dDataDebug = []
  }
[ Directory
    { dTagName = "DATA"
    , dTagNum = 9
    , dElemTypeCode = 4
    , dElemTypeDesc = "short"
    , dElemType = ElemShort
    , dElemSize = 2
    , dElemNum = 7440
    , dDataSize = 14880
    , dDataOffset = 128
    , dData = ""
    , dDataDebug = []
    }
    
.
.
.

DATA {short} tagNum=9 size=2 count=7440 offset=128  []
DATA {short} tagNum=10 size=2 count=7440 offset=15008  []
DATA {short} tagNum=11 size=2 count=7440 offset=29888  []
DATA {short} tagNum=12 size=2 count=7440 offset=44768  []
FWO_ {char} tagNum=1 size=1 count=4 offset=1195463747  ["GATC"]
LANE {short} tagNum=1 size=2 count=1 offset=65536  ["1"]
PBAS {char} tagNum=1 size=1 count=744 offset=59648  ["GGGGGCAACTAAAGGAAGCTCTATTAGATACAGGAGCAGATGATACAGTATTAGAAGAAATGAGTTTGCCAGGAAGATGGAAACCAAAAATGATAGGGGGAATTGGAGGTTTTATCAAAGTAAGACAGTATGATCAGATACTCATAGAAATCTGTGGACATAAAGCTATAGGTACAGTATTAGTAGGACCTACACCTGTCAACATAATTGGAAGAAATCTGTTGACTCAGATTGGTTGCACTTTAAATTTTCCCATTAGCCCTATTGAGACTGTACCAGTAAAATTAAAGCCAGGAATGGATGGCCCAAAAGTTAAACAATGGCCATTGACAGAAGAAAAAATAAAAGCATTAGTAGAAATTTGTACAGAGATGGAAAAGGAAGGGAAAATTTCAAAAATTGGGCCTGAAAATCCATACAATACTCCAGTATTTGCCATAAAGAAAAAAGACAGTACTAAATGGAGAAAATTAGTAGATTTCAGAGAACTTAATAAGAGAACTCAAGACTTCTGGGAAGTTCAATTAGGAATACCACATCCCGCAGGGTTAAAAAAGAAAAAATCAGTAACAGTACTGGATGTGGGTGATGCATATTTTTCAGTTCCCTTAGATGAAGACTTCAGGAAGTATACTGCATTTACCATACCTAGTATAAACAATGAGACACCAGGGATTAGATATCAGTACAATGTGCTTCCACAGGGATGGAAAGGATCACCAGCAATATTCCAAAGTAGCATGA"]
PDMF {pString} tagNum=1 size=1 count=23 offset=60392  ["KB_3500_POP7_BDTv3.mob"]
PDMF {pString} tagNum=2 size=1 count=23 offset=60415  ["KB_3500_POP7_BDTv3.mob"]
PLOC {short} tagNum=1 size=2 count=744 offset=60438  []
S/N% {short} tagNum=1 size=2 count=4 offset=61926  []
SMPL {pString} tagNum=1 size=1 count=10 offset=61934  ["S17-SeqF1"]
CMNT {pString} tagNum=1 size=1 count=1 offset=61944  ["Generated by HyraxBio AB1 generator"]

The data is output twice. The first section is the detail, the second is the summary.

Selected data types have the "debug data" element populated. e.g. the PBAS (FASTA)

Generate minimal AB1s from FASTAs

To create an AB1 run

hyraxAbi-exe gen "./pathContainingFastas" "./pathForOutputAb1s"

This will create an AB1 per input FASTA

Input FASTA format

Each input data should have the following format

> weight
read
> weight
read

Weighted reads

For example

> 0.5
ACG
> 0.3
AAAA
> 1
__AC

Results in the following weighted nucleotide per position

Note that the reads do not need to be the same length.


Example FASTA - single file

eg1.fasta

> 1
ACTG

Here there is a single FASTA with a single read with a weigh of 1 (100%). The chromatogram for this AB1 shows perfect traces for the input ACTG nucleotides


Example FASTA - two FASTA files

eg1.fasta

> 1
ACAG

eg2.fasta

> 1
ACTG

Two input FASTA files both with a weigh of 1. You can see in the second trace that the third nucleotide is a T (the trace is green). Exactly what the base-calling software (phred & recall etc) decide to call the base as depends on your settings and software choices.


Example FASTA - two FASTA files with different weights

eg1.fasta

> 1
ACAG

eg2.fasta

> 0.3
ACTG

Here the second fasta has a weight of 0.3 and you can see the traces are 30% of the height of the top ones.


Example FASTA - single FASTA with a mix

eg1.fasta

> 1
ACAG
> 0.3
ACTG

The single input FASTA has an AT mix at the third nucleotide. The first read has a weight of 1 and the second a weight of 0.3. Notice that the maximum weight is 1, e.g. the first A has the same intensity as the second even though the first one has the reads weighted both 1 and 0.3


Example FASTA - Multiple mixes

eg1.fasta

> 1
ACAG
> 0.3
_GT
> 0.2
_G


Using the modules

For a detailed overview of the code see TODO and the haddock documentation TODO

For now the terminal app (Main.hs) serves as an example and the best starting point to understand the code

E.g. Add a comment to an existing AB1 file

import qualified Hyrax.Abi as Abi
import qualified Hyrax.Abi.Read as Abi
import qualified Hyrax.Abi.Write as Abi

addComment :: IO ()
addComment = do
  abi' <- Abi.readAbi "example.ab1"

  case abi' of
    Left e -> putText $ "error reading ABI: " <> e
    Right abi -> do
      let modified = Abi.addDirectory abi $ Abi.mkComment "new comment"
      Abi.writeAbi "example.modified.ab1" modified

For additional examples see the Examples directory