This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Trinity - Type Of Reads: ( Cfa, Cfq, Fa, Or Fq )

Hi there!

I accidentaly ran Trinity with cfq parameter instead of fq which I believe is what I wanted - fastq with quality scores. I would really like to know what cfq stands for as I tried to google it with no success :/

from Trinity help: --seqType <string> :type of reads: ( cfa, cfq, fa, or fq )

Thanks for suggestions.

trinity
• 4,443 views • link
updated by Farhat Farhat • written by Biomonika (Noolean) Biomonika (Noolean)

1 answer

It stands for colorspace fastq, the output from a SOLiD sequencer. For a brief primer check http://marketing.appliedbiosystems.com/images/Product_Microsites/Solid_Knowledge_MS/pdf/CSHL_Fu.pdf

• 0 views • link
written by Farhat Farhat

Thank you. I had Illumina data, but all with really high quality scores. Can I expect that results will be similiar if I ran it with cfq instead of fq? (Will compare both runs though)

• 0 views • link
written by Biomonika (Noolean) Biomonika (Noolean)

No, it will be completely nonsensical most likely.

• 1 views • link
written by Farhat Farhat

Log in to answer this question.

More posts like this
  • Vcf File From Mpileup ... What Next
    written by Dataminer Dataminer
    <p>Hi,</p> <p>I have done variant calling on my file using <a href='http://samtools.sourceforge.net/'>samtools</a> mpile up and I have converetd my file from bcf to vcf. </p> …
  • How To Find The Most Frequent Allel Snp In Vcf File
    written by Chen Sun Chen Sun
    <p>When I use <code>samtools mpileup</code>, there are usually several allel SNP in ALT column according to REF column. just as follows:</p> <pre><code>##fileformat=VCFv4.1 ##samtoolsVersion=0.1.19-44428cd ##reference=file:///gpfs/home/cxs1031/backup/ref/chr22.fa ##contig=&lt;ID=chr22,length=51304566&gt; …
  • Tophat And Multiple-Mapped Reads
    written by Biomonika (Noolean) Biomonika (Noolean)
    <p>I would like to ask about update on this topic:</p> <p><a href='http://seqanswers.com/forums/showthread.php?t=24270'>http://seqanswers.com/forums/showthread.php?t=24270</a></p> <pre><code>-g/--max-multihits &lt;int&gt; Instructs TopHat to allow up to this many alignments to the …
  • Possible Bug In Fastx Toolkit
    written by Abhi Abhi
    <p>I am getting the following error with fastx_quality_stats with version 0.0.13.2. May be this has to do with the changing Illumina scores. Now they are …
  • Error Loading Vcf File In Ucsc Genome Browser
    written by win win
    <p>Hi all, I generated a vcf file using samtools and the wanted to load it into UCSC Genome Browser as suggested at this link</p> <p><a …
  • Fastq Format Confusion
    written by Biomonika (Noolean) Biomonika (Noolean)
    <p>Hi!</p> <p>In the following article <a href='http://file:///home/misa/Sta%C5%BEen%C3%A9/NGS_workshop_1_Intro.pdf'>The Sanger FASTQ file format for sequences with quality scores, and the Solexa/Illumina FASTQ variants.</a> I found this table:</p> …
  • Reporting More Variants With Samtools
    written by Biomonika (Noolean) Biomonika (Noolean)
    <p>I am using samtools mpileup to call variants from my bam file. However, there are situation where for example there is G in reference sequence …
  • Error While Running Samse (Bwa)
    written by Higherdefender • Higherdefender
    <p>Hi,</p> <p>I am trying to align short sequences of chromosomeX w.r.t the reference gnome using bwa.</p> <p>After successfully running aln function when I ran samse …
  • Genomic Alignment And Snp/Indel Calling - My First Ever "Pipeline"
    written by Travis Travis
    <p>So I have finally generated some reads and run it through what I guess could be called a very rudimentary 'pipeline'. I generated a million …
  • Using Seqan To Build A Suffix Tree For Millions Of Reads
    written by Weaklearner • Weaklearner
    <p>I am trying to build a suffix tree for millions of reads with Seqan. Since the number of sequences is so big, I need some …