Sanger-Sized Pairwise Alignment To Variants To Translation
0
0
Entering edit mode
12.1 years ago

I'm looking for some implementations of the following pattern (in any language):

  1. One Sanger-sized query sequence is pairwise aligned against a short (<100kbp) reference with some BED/GFF gene calls.
  2. Alignment is parsed - sequence substitutions and indels in the non-genic areas are reported as is.
  3. Mutations in the genic sections are examined for syn/non-syn substitutions, frameshifts.

I imagine this would be a fairly common script so I'd like to see some existing examples and approaches before reinventing the wheel. "Use Bioperl" is not an acceptable answer.

codon alignment parsing • 2.4k views
ADD COMMENT

Login before adding your answer.

Traffic: 1765 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6