Using the 27-primates UCSC Multiz Alignment · Biostar

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

Using the 27-primates UCSC Multiz Alignment

Hello,

I am attempting to use the 27-primates UCSC multiz alignment data so that I can find an orthogonal sequence (across different mammals) to a human reference sequence (whether an ORF, a transcript, etc). I downloaded all the data in the maf folders here https://hgdownload.soe.ucsc.edu/goldenPath/hg38/multiz30way/maf/, but am very lost on how to extract anything meaningful out of it. Specifically, even with its description, I am unclear on what the alignments and maf folders are, and how to read the files within them. If anyone could send a resources/explanations, it would be greatly appreciated. Thanks!

multiz ucsc conservation alignment

• 1,358 views • link

written 3.3 years ago by Daniel •

There are several maf related utility programs included in Kent utilities at UCSC. You can find their descriptions in this file (there are other programs in there as well so look for maf programs). I linked linux version but there are macOS binaries available as well.

Here is an example of how to extract regions: https://groups.google.com/a/soe.ucsc.edu/g/genome/c/87HvXW137SI

• 0 views • link

written 3.3 years ago by GenoMax

Thank you so much. The documentation seems a bit ambiguous, but I'll try my best and will post some questions here if I have them!

• 0 views • link

written 3.3 years ago by Daniel •

0 answers

No answers yet.

Log in to answer this question.

More posts like this

Phased haplotype using WGS data from Nebula Genomics and GATK
written 1.2 years ago by biostars •

I am a technically-competent non-bioinformatician with a set of three WGS sequences for a grandmother, son, and daughter that were provided by Nebula Genomics (for …
Converting genome build from hg38 to hg19 using liftOver via bigsnpr - error fetching file
written 2.0 years ago by dam4l

I am trying to use LiftOver via bigsnpr to change genome coordinates from hg38 to hg19. I have tried the following: ```r hg19 <- snp_modifyBuild( …
Guidance on Reference-Based Assembly and Variant Analysis of Exome Sequencing Data
written 2.0 years ago by Rohan •

Hello, I am working with paired-end whole genome sequencing (WGS) data focused on the exome. My goal is to assemble these reads based on the …
Can I liftover hs37d5 to hg38 using hg19ToHg38 chain file?
written 3.1 years ago by ConvolutedGenome •

Hello, I am currently working on analyzing the somatic mutation data from the PCAWG project. In the main PCAWG paper (https://www.nature.com/articles/s41586-020-1969-6), they mentioned aligning their …
Funcotator reference file error in GATK4
written 4.4 years ago by KE •

I am using a function called Funcotator ([here][1]) in GATK4. The code I am trying to run is: gatk Funcotator \ --variant chr21.vcf \ --reference …
How to extract MAF block regions that have full alignment?
written 4.6 years ago by nattzy94 •

Hi, I am trying to assess the evolutionary conservation of putative small protein coding regions. I am following this [paper][1] which does a similar thing. …
ROAST-Multiz inconsistent row size
written 6.1 years ago by Diana MoSa •

Hello, I am having issues whith ROAST in order to perform Multiz aligment of 37 genomes. I performed pairwise alignments with lastz for all the …
Assemble bacterial .fastq files and find differences (SPAdes followed by MUMmer)
written 6.9 years ago by LRStar

This is somewhat of a follow-up from a previous post (https://www.biostars.org/p/395813/). I was recently given 4 paired-end .fastq files (where each read has about **150 …
Job: Bioinformatics Research Consultant needed for constructing TFBS and regulatory networks from t…
written 9.5 years ago by tfhahn •

If you are interested in teaching me any of this but you don't have the time to read everything I wrote please email me at …