<div dir="ltr"><div>I don't think we have anything like this (yet). Are efforts like VCF (variant call format) related but separate in your mind?<br></div><div><br></div><div>Peter<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Oct 31, 2023 at 7:31 PM David Merberg <<a href="mailto:merbergd@gmail.com">merbergd@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg3959824987974020465"><div><div style="font-family:Helvetica,Arial;font-size:13px">Hello biopython world,</div><div style="font-family:Helvetica,Arial;font-size:13px"><br></div><div style="font-family:Helvetica,Arial;font-size:13px">For my last job, I wrote some python code to categorize and describe sequence changes of many types. I used biopython to handle sequences and some basic functions like  IO and translation, but I did not find a module for reading variants/mutants and applying them to sequences.</div><div style="font-family:Helvetica,Arial;font-size:13px"><br></div><div style="font-family:Helvetica,Arial;font-size:13px">Some cases are trivial, but some are not. For example, a small deletion in the nucleotide sequence may have no effect on the amino acid corresponding to the position of the affected codon, but will affect downstream amino acids. Protein changes caused by deletions or insertions of 3, 6, 9 . . . nucleotides can also be tricky to calculate.</div><div style="font-family:Helvetica,Arial;font-size:13px"><br></div><div style="font-family:Helvetica,Arial;font-size:13px">My question is whether there is a biopython module to read variants in a standard format (see for example <a href="http://varnomen.hgvs.org/" target="_blank">http://varnomen.hgvs.org/</a>)? Along with the variant objects there could be a set of methods to operate on mutated sequences. Does the community think that this would be useful if it does not already exist?</div><div style="font-family:Helvetica,Arial;font-size:13px"><br></div><div style="font-family:Helvetica,Arial;font-size:13px">I implemented many functions for these sorts of operations, but I realized soon afterwards that there are probably better ways to do much of it. I always wanted to redo the work, but never had time. Now I have time, but am not at that job. If it would be useful to the community, I may be able to take it on as a contribution to biopython.</div><div style="font-family:Helvetica,Arial;font-size:13px"><br></div><div style="font-family:Helvetica,Arial;font-size:13px">A caveat is that I don’t have experience contributing to multi-developer projects. I try to write clean, well documented code and I’m familiar with the basics of git. So, it’s OK if you’d prefer that I start with something smaller (like unit tests or documentation). Just let me know.</div><div style="font-family:Helvetica,Arial;font-size:13px"><br></div><div style="font-family:Helvetica,Arial;font-size:13px">Dave Merberg</div><br><div class="gmail_signature"></div></div>
_______________________________________________<br>
Biopython mailing list  -  <a href="mailto:Biopython@biopython.org" target="_blank">Biopython@biopython.org</a><br>
<a href="https://mailman.open-bio.org/mailman/listinfo/biopython" rel="noreferrer" target="_blank">https://mailman.open-bio.org/mailman/listinfo/biopython</a><br>
</div></blockquote></div>