This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Vmatch'S Mkvtree Gives Illegal Character

I am trying to contruct an index with Vmatch's mkvtree. If I use two input files together, I get an error:

../software/vmatch.distribution/mkvtree -db data/EST.gb data/mRNA.gb -dna -indexname index/napus -pl -v -allout
>data/EST.gb
reading file "data/EST.gb"
>data/mRNA.gb
reading file "data/mRNA.gb"
../software/vmatch.distribution/mkvtree: Illegal character 'L' in file "data/EST.gb" line 2

But if I use either of them separately, it works fine:

../software/vmatch.distribution/mkvtree -db data/EST.gb  -dna -indexname index/napus -pl -v -allout

or 

../software/vmatch.distribution/mkvtree -db data/mRNA.gb -dna -indexname index/napus -pl -v -allout

The files are at http://www.putlocker.com/file/7C776C1CAFD1FA75 and a copy at http://www.filedropper.com/rawdatatar . The EST.gb file is large, so I needed a Linux machine of about 15GB memory to do make the index for it.

alignment
• 2,531 views • link
written by hbw hbw

Try concatenating the files and constructing an index from the one combined file. Maybe mkvtree does not like your invocation with two file names and is throwing an error, or there may be something wrong with your file format.

• 0 views • link
written by SES SES

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • Database Error When Running Blastn
    written by Anima Mundi Anima Mundi
    <p>Hello, as initially stated <a href='http://www.biostars.org/p/92157/' title='here'>here</a>, I am experiencing problems when I try to run the blastn command <code>blastn -query myquery -db mydatabase.fasta.nhr</code>.</p> <p>What …
  • Galaxy "Stitch Maf Blocks" Produced Null Bases
    written by lam • lam
    <p>problem: "stitch MAF blocks" produced NULL bases </p> <p>steps: From <a href='http://genome.ucsc.edu/cgi-bin/das/hg19/dna?segment=chrY:23310464,23331995'>http://genome.ucsc.edu/cgi-bin/das/hg19/dna?segment=chrY:23310464,23331995</a></p> <p>we got the following bases:</p> <pre><code>agtttacccaactcacaggtatttgaggatacaaacccatgactgggcttggctttaaaaattcctatctgaaatttctagt.......... </code></pre> <p>We want to do multiple-species alignment …
  • Inconsistent Cap3 Assembly
    written by hbw hbw
    <p>This question is cross-posted at <a href='http://seqanswers.com/forums/showthread.php?t=31895'>SeqAnswers</a> with the alignments. I have 5 related EST sequences and I want to assemble contigs. I am using …
  • Good Poly-A Filtering Rules Or Tools
    written by hbw hbw
    <p>I am aligning a large number of ESTs. It seems poly-A tails show in many different ways. In addition to occurring at the very end, …
  • Nucleotide Determination In Samtools Pileup Function
    written by coryprzy • coryprzy
    <p>I'm working with coding in c/c++ using the <a href='http://samtools.sourceforge.net/'>samtools</a> suite.</p> <p>Starting with the example 'calDepth.c' included in the <a href='http://samtools.sourceforge.net/'>samtools</a> package, I've tried to …
  • Do Different Versions Of Bwa Create Differently Named Index Files?
    written by mathiasf52 • mathiasf52
    <p>Hello,</p> <p>I installed <a href='http://bio-bwa.sourceforge.net/'>bwa</a>-0.6.1 and index human genome like </p> <pre><code>bwa index -a bwtsw hg19.fa </code></pre> <p>it generate five files (.amb, .ann, .bwt, .pac, …
  • Plink/Seq Denovo Option
    written by Caddymob Caddymob
    <p>Was looking at the <a href='http://atgu.mgh.harvard.edu/plinkseq/index.shtml'>PLINK/SEQ</a> package today for pulling out denovo variants and munging around in my VCFs. The manual really has nothing on …
  • Blast Gives Cryptic Errors
    written by hbw hbw
    <p>I have a list of proteins in fasta format (say goodProteins.fasta). I first make a compatible database using NCBI's formatdb v 2.2.18:</p> <pre><code>formatdb -i goodProteins.fasta …
  • How To Interpret Funky Vcf Output?
    written by Dan D Dan D
    <p>A customer with twelve RNASEQ mouse samples sequenced on an Illumina HiSeq 2000 wanted to compare the expression numbers for a single gene. This is …
  • Phylogenetic Tree With Dna, Binary And Multi Data - Raxml Update
    written by Pavid Pavid
    <p>Hey!</p> <p>I've made a <a href="http://biostar.stackexchange.com/questions/10406/supertree-with-matrix-representation-with-parsimony">question</a> about supertrees. <br/> My goal is to integrate data from SNPs (DNA), MIRUs (multi -range around 1 to 25) …