cytobandread - Read cytogenetic banding information

Syntax

Input Arguments

File	String specifying a file containing cytogenetic G-banding data, such as an NCBI ideogram text file or a UCSC Genome Browser cytoband text file.

Output Arguments

CytoStruct

Structure containing cytogenetic G-banding data in the following fields: ChromLabels BandStartBPs BandEndBPs BandLabels GieStains

Description

CytoStruct = cytobandread(File) reads File, which is a string specifying a file containing cytogenetic G-banding data, and returns CytoStruct, which is a structure containing the following fields.

Field	Description
ChromLabels	Cell array containing the chromosome label (number or letter) on which each band is located.
BandStartBPs	Column vector containing the number of the base pair at the start of each band.
BandEndBPs	Column vector containing the number of the base pair at the end of each band.
BandLabels	Cell array containing the FISH label of each band, for example, p32.3.
GieStains	Cell array containing the Giemsa staining result for each band. Possible stain results depend on the species. For example, for Homo sapiens, the possibilities are: gneg gpos25 gpos50 gpos75 gpos100 acen stalk gvar

Tip   You can download files containing cytogenetic G-banding data from the NCBI or UCSC Genome Browser ftp site. For example, you can download the cytogenetic banding data for Homo sapiens from: ftp://ftp.ncbi.nlm.nih.gov/genomes/H_sapiens/mapview/ideogram.gz or ftp://hgdownload.cse.ucsc.edu/goldenPath/hg18/database/cytoBandIdeo.txt.gz

Examples

Read the cytogenetic banding information for Homo sapiens into a structure.

hs_cytobands = cytobandread('hs_cytoBand.txt')

hs_cytobands = 

     ChromLabels: {862x1 cell}
    BandStartBPs: [862x1 int32]
      BandEndBPs: [862x1 int32]
      BandLabels: {862x1 cell}
       GieStains: {862x1 cell}

See Also

chromosomeplot

© 1984-2012- The MathWorks, Inc.    -   Site Help   -   Patents   -   Trademarks   -   Privacy Policy   -   Preventing Piracy   -   RSS