I have an interleaved Nexus file, were each block correspond to a gene in the alignment.
How could I split that file into separate alignments, one for each of those genes.
I have tried using the AlignIO module from biopython, but it just load the whole alignment as one big sequence for each taxa without any clear way to iterate through the genes (blocks) or interleaves of the alignment.
Do you know if there is a way to iterate over the blocks in the nexus file with Biopython? I could not find anything related on the Biopython documentation.
Thanks!
•
written
by
rafa.rios.50 •
1 answer
The easiest way to do this is to export the interleaved file as a sequential one. Lots of tools will do this (PAUP*, EMBOSS, etc.).
I don't have an example to test, but read.nexus.data() in ape in R might be able to read it in. If so, it's easy to write it out:
library(ape)
a <- read.nexus.data("yourfile.nex")
write.nexus.data(a, "yourfile.sequential.nex", interleaved=FALSE)
•
written
by
Brice Sarver
Log in to answer this question.
-
seurat: Average expression function and find markers DE fold changes do not add up
written by rasing02 •
-
Concatenating multiple nexus file in a supermatrix using Biopython
written by CUCIO199 •
-
map both interleaved and orphaned reads to assembly for each sample (with bowtie2)
written by willnotburn •
-
Ensembl Multiple Genome alignments do not yield the same as pairwise alignments
written by except.for.me.and.my.monkey •
-
MrBayes nexus file execute problem
written by Mehmet
-
samtools mpileup Aborted (Core dumped)
written by rafa.rios.50 •
-
Converting Mugsy Output (Maf) To Other Formats, Problems With +/- Strand And Overlapping Blocks
written by tmortimer •
-
Alignio Gives 'Assertionerror' When Reading Emboss Alignment Files
written by a1ultima
-
Biopython Indentation Error
written by priyasshah •
-
How Can I Automate Blast Of >1 Sequence And Output The Top Hits For Each Input?
written by bow
