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

Hi,

I'm having lot of trouble giving the right syntax for NovoalignCS. I have paired reads of SOLiD data and want to map them to human genome. I'm giving the following syntax and have no luck so far executing the tool.

novoalignCS 
-d h19_index_novo_cs \
-f N12_0056_DefaultLibrary_F3.csfasta N12_0056_DefaultLibrary_F5-BC.csfasta N12_0056_DefaultLibrary_F3.QV.qual N12_0056_DefaultLibrary_F5-BC.QV.qual \
-F CSFASTAnQV \
-o SAM \
-i PE \
-k

The error it's displaying is: "Problem processing command line arguments: Unrecognised option. at N12_0056_DefaultLibrary_F3.QV.qual".

solid
• 2,047 views • link
updated by Ashutosh Pandey Ashutosh Pandey • written by Jordan Jordan

1 answer

Can you give a link to the manual ? I think the csfasta and the qual file for each pair should go together (one after the other). Sth like 

-f N12_0056_DefaultLibrary_F3.csfasta N12_0056_DefaultLibrary_F3.QV.qual N12_0056_DefaultLibrary_F5-BC.csfasta  N12_0056_DefaultLibrary_F5-BC.QV.qual
• 0 views • link
updated by Istvan Albert Istvan Albert • written by Ashutosh Pandey Ashutosh Pandey

Log in to answer this question.

More posts like this
  • Bowtie alignment for PE SOLiD data
    written by nikkihathi • nikkihathi
    Sorry my mistake: solution of simple: bowtie -t -p 4 -C --sam Genome/mouce_database_10/mm10_c -f -1 /SOLiD_data/F3/reads/5500xl_F3.csfasta -2 /SOLiD_data/F5-RNA/reads/5500xl_F5-RNA.csfasta --Q1 /SOLiD_data/F3/reads/5500xl_F3.QV.qual --Q2 /SOLiD_data/F5-RNA/reads/5500xl_F5-RNA.QV.qual output1.sam 2> stdout.err …
  • Need help using Shrimp2 on paired end color-space SOLiD data.
    written by Jordan Jordan
    Hi, I have SOLiD reads which are paried-end (75bp and 35bp) in .csfasta and .QV.qual format. I would like to use Shrimp2 to align them. …
  • Problem With Counting Mapped Reads
    written by njs531 • njs531
    <p>Hi, This is my very first experience analysing RNAseq data. My goal is to do differential analysis between two strains of a bacteria. So far, …
  • Sra-Pileup Question
    written by Sergei Lebedev • Sergei Lebedev
    <p>I'm having trouble extracting a <a href='http://samtools.sourceforge.net/pileup.shtml'>pileup</a> file from an <a href='ftp://ftp-trace.ncbi.nlm.nih.gov/sra/sra-instant/reads/ByExp/sra/SRX/SRX101/SRX101206/SRR364089/SRR364089.sra'>SRA</a> using <code>sra-pileup</code> command. NCBI <a href='http://www.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?view=toolkit_doc&amp;f=sra-pileup'>documentation</a> only covers command line arguments for the …
  • Error Reading Reference Genome Using Bfast
    written by Jordan Jordan
    <p>Hi,</p> <p>I'm getting an error while executing the following step in BFAST.</p> <pre><code> bfast fasta2brg -f ~/Refs/human_hg19.fa "************************************************************" Checking input parameters supplied by the user …
  • Solid Sequencing Protocol And How To Obtain Quality Of Solid Reads?
    written by Jordan Jordan
    <p>I have just started working on NGS few weeks back and I was trying to get the pipeline for SOLiD right.</p> <p>From what I understand …
  • Aat (Analysis And Annotation Of Genomes) - Error
    written by Daniel Fernandez • Daniel Fernandez
    <p>Hi everyone.</p> <p>I'm trying to analysing a mammalian genome by the AAT package (<a href="http://aatpackage.sourceforge.net/">http://aatpackage.sourceforge.net/</a>), the package comes with examples (only I have to run …
  • Tophat Question: No Accepted_Hits.Bam Generated
    written by Rnaer • Rnaer
    <p>I am using tophat to map the SOLiD single-end 50mer reads to S.cerevisae genome. The commands I used are:</p> <pre><code>$ bowtie-build -C reads/sgd_genome.fa reads/sgd_genome $ …