This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Is There An Existing Data Format Suitable For Holding Dna Sequence Variations In A Population?

If I have multiple genomic DNA sequences for a small protein, and want to represent this variation in a file, are there any existing formats to do this?

I could save all the sequences one per line and re-calculate the variation information each time, but this is a waste of computational resources.

Before I create YASF (yet-another-sequence-format) I was wondering if anyone knew of an existing one?

It should ideally be able to represent A|C|G|T with an optional 3 or 4 floats for the relative abundance of each. I wouldn't want to store the floats if there was no variation at any particular point in the sequence.

If it could handle gaps/insertions that would be useful too!

[Edit] I should have mentioned that the proteins in question are going to be antibodies, and so the data will consist of large numbers of different sequences based on similar VDJ recombination building-blocks with somatic hypermutation providing a vast number of similar, yet different, DNA sequences.

format
• 3,106 views • link
updated by Casey Bergman Casey Bergman • written by Nickengland Nickengland

3 answers

Perhaps MAF http://www.bioperl.org/wiki/MAF_multiple_alignment_format is applicable? VCF (variant call format) is also useful for this type of information. Some multiple alignment software outputs might be useful. Position Specific Weight Matrices might be applicable. It all depends a bit on what you want to do with the data and from where the data are coming.

• 1 views • link
written by Sean Davis Sean Davis

I think VCF is going to be what you want to use here

• 0 views • link
written by Chris Miller Chris Miller

I think MAF is going to be more appropriate than VCF, as the data I will be using will have ambigious/unknown positional information in the chromosome, as random mucleotides are inserted or deleated during B-cell development.

• 0 views • link
written by Nickengland Nickengland

If you need to create a new format, use a RDF file with an already defined ontology (e.g: http://variationontology.org/ ).

You can then query your data with SPARQL , transform the data using XSLT, share...

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

How about the Genome Variation Format: http://www.sequenceontology.org/resources/gvf.html

• 1 views • link
written by Casey Bergman Casey Bergman

Log in to answer this question.

More posts like this