aaChanges

classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

aaChanges

Scooter Willis-2
Trying to work with aaChanges under Human Genome Variation using an Amazon EC2 cloudman launched instance.

The output is generated by Filter Pileup and has HG19 associated for the Database/build. 

The gene set annotations imported from USCS also have HG19.

I get the following errors when selecting input data set and appears to be the same problem reported 2/2011(http://user.list.galaxyproject.org/Problems-using-aaChanges-tool-td4134510.html)

Also trying to sort through what looks to be an issue with using pileup for SNP calling in the 10 column format. Column 3 is the reference base and column 4 is the consensus base. I have numerous examples where reference base = consensus base with a breakout of the number of A T C G represented as a count. For aaChanges I need an input line for each non-zero count of A T C G from the pileup report to test which of those SNPs would actually result in an amino acid change. Doesn't seem to be an easy thing to do in Galaxy.

Any recommendations on how others are doing SNP calling from exome data and detecting amino acid changes? Does anyone have a published workflow they would recommend or could share?



___________________________________________________________
Please keep all replies on the list by using "reply all"
in your mail client.  To manage your subscriptions to this
and other Galaxy lists, please use the interface at:

  http://lists.bx.psu.edu/