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

HI,

I am using multiple alignment of non coding sequence to select conserved block. All the sequence has different length and when i run command : Gblocks result.pir -t=d -e=-gb1 -b4=5 -d=y, it is not accepting my data and showing---

Sequences with different number of positions. Execution terminated.

Is there any way to process data in Gblock?

Thanks in advance.

alignment
• 3,493 views • link
written by baurumon • baurumon

i tried it but did not work.

I converted my .fasta file to .pir with EMBL sequence conversion tool. https://www.ebi.ac.uk/Tools/sfc/ . Sequences look like as below-

>D1;3amexicanus
, 119 bases
 AGAAGCCAGC TCTGAAGAAA GTGCGCATGT CTGCTGACGC TATGCTGCAA
 GCTCTGCTCG GCTCCAAACA CAAGGTGTCA ATGGACCTGA GGGCCAACCT
 CAAACAAGTC AAGAAGGAG*

>D1;4amexicanus
, 127 bases
 AGAGGCATGA AGGTCATTGA GAACAGGGCA CTGAAGGATG AGGAGAAGAT
 GGAGCTGCAA GAGATCCAGC TCAAAGAGGC CAAACACATT GCTGAGGAAG
 CTGATCGCAA ATATGAGGAA GTGAGTT*

>D1;5amexicanus
, 519 bases
 GGTCTAATGT TCATTCCTTG TGTTTCTTGG CCCAAACAAA TCTCTTCTGC
 TTGTTGCCTT TCCTAAGCAG TGGATTCCTA GCAGCTATTT GATCATGAAG
 GCCTGATTTG CGCCTTAGTC TCCTCTTAAC AGTTGTTCTA GAGATGTGTC
 TGCTGCTAGA ACTCTGTGTG GCATTCAACT GCTCTCTGAT CTgagctgct
 gttaacttgc ggtttctgag gctggtgact CGGATGTACT TATCCTGAGA
 AGCAGCAGAG gtgactcttg gtcttccttt cctgggttgG TCTTCATGTG
 TGCCAGTTTC GCTGTAGCGC TTGATGGTTT TTGTGACTCC ACTTGGGgac
 acatttaaag tttttgcaat tttctggact gactgacctt catttcttaa
 agtaatggTG GCCACTCATT TTTCTTTAGT TAGCTGAttg gttcttgcca
 taatatgaat tttaacagtt gtCCAATAAG GCTGTCGGGT GTGTAATAAC
 CAGACTTCTG CGCAACACA*
• 0 views • link
updated by GenoMax GenoMax • written by baurumon • baurumon

1 answer

You need to input an alignment in .pir format not plain fasta sequences converted to .pir format for this tool.

• 0 views • link
written by GenoMax GenoMax

Log in to answer this question.

More posts like this