This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Parse allele database

Dear Biostars users.

I have this variant format.

ex)

CHROM POS   REF   ALT
1             150    CAC  CAAC

Can I this format change following format using python?

CHROM POS   REF   ALT
1             150      C        CA
1             151      A        AA
python
• 1,656 views • link
updated by Ram Ram • written by oghzzang • oghzzang

Please use the formatting bar (especially the code option) to present your post better. You can use backticks for inline code (`text` becomes text), or select a chunk of text and use the highlighted button to format it as a code block. I've done it for you this time.
code_formatting

• 0 views • link
written by Ram Ram

Thanks. From now, I'll use this button. :)

• 1 views • link
written by oghzzang • oghzzang

1 answer

What you're asking for is called "left aligning normalization". It represents variants in the most parsimonious notation and is one of the best practices I've encountered and continue to use all the time.

If you have the VCF file this data comes from and the reference sequence used in the analysis, you can use either bcftools norm (bcftools) or vt decompose | vt norm (vt) to get to where you need from the VCF file. I'd recommend the latter as it makes tracking changes easier by adding OLD_MULTIALLELIC and OLD_VARIANT INFO fields.

If not, it becomes a much more challenging task because you're going to need to compare the REFERENCE sequence and ALT alleles manually to get to your solution.

• 0 views • link
written by Ram Ram

Log in to answer this question.

More posts like this