Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • question about Solid data analysis using bowtie?


    hello ,everyone
    I am trying to run bowtie with the solid sequence data flower9_1_F3.csfasta and flower9_1_F3_QV.qual
    I built the bowtie-build -C index successfully and moved the generated index files to the Bowtie indexes subdirectory. I run the command:
    bowtie -n 3 --best --strata -a -e 150 -p 4 ../cDNA_reference/TAIR9_cdna -C -f flower9_1_F3.csfasta -Q flower9_1_F3_QV.qual --un fllower9-unmap > flower9-map &
    and the output i got is an error showing on console

    [wanglei@mu01 data]$ [wanglei@mu01 data]$ Too few quality values for read: 425_2031_2008_F3
    are you sure this is a FASTQ-int file?
    terminate called after throwing an instance of 'int'
    Too few quality values for read: 425_2042_1847_F3
    are you sure this is a FASTQ-int file?
    terminate called after throwing an instance of 'int'
    Too few quality values for read: 1277_2042_2013_F3
    are you sure this is a FASTQ-int file?
    Could you please help me to sort out this issue?
    thanks
    wanglei

  • #2
    Please only post to one forum. Cross-posting is like spam, unacceptable.

    Comment


    • #3
      ok,i have deleted it, thanks!

      Comment


      • #4
        Truncated Qual file

        I have had this error when SAET failed and a truncated qual file was produced.

        Comment


        • #5
          Hey guys,
          How did you solve this problem ? Could you please help me with this ?
          I am also getting this error, when using tophat mapping on colorspace file.

          This is my command: [Works fine until half-way]

          tophat --output-dir /results --color -p 8 --no-coverage-search /bowtie_alignment_reference/mm9_c F3.csfasta F3.QV.qual


          I get this error:
          Mapping left_kept_reads_seg1 to genome mm9_c with Bowtie (1/3)
          Mapping left_kept_reads_seg2 to genome mm9_c with Bowtie (2/3)
          Mapping left_kept_reads_seg3 to genome mm9_c with Bowtie (3/3)

          [2012-10-20 22:11:02] Mapping right_kept_reads to genome mm9_c with Bowtie
          [FAILED]
          Error running bowtie:
          Too few quality values for read: 39 I
          are you sure this is a FASTQ-int file?
          terminate called after throwing an instance of 'int'


          regards
          Chirag

          Comment


          • #6
            Did anyone of you solve this problem?
            I have exactly the same problem now.

            Mapping left_kept_reads to genome GRCh37_gatk_colorspace with Bowtie
            Mapping left_kept_reads_seg1 to genome GRCh37_gatk_colorspace with Bowtie (1/2)
            Mapping left_kept_reads_seg2 to genome GRCh37_gatk_colorspace with Bowtie (2/2)
            Mapping right_kept_reads to genome GRCh37_gatk_colorspace with Bowtie
            [FAILED]
            Error running bowtie:
            Too few quality values for read: 28 I
            are you sure this is a FASTQ-int file?
            terminate called after throwing an instance of 'int'

            All 4 files are exactly 200000 lines (so 100000 reads). Mapping the left reads works, but when mapping the right reads it fails...

            [EDIT]
            For anyone who has this error too, check whether you have added the --quals option, also check the command on order of files (first all csfasta files than the qual files). This solved my error...
            Last edited by Jetse; 03-26-2013, 08:02 AM. Reason: solved

            Comment

            Latest Articles

            Collapse

            • seqadmin
              Essential Discoveries and Tools in Epitranscriptomics
              by seqadmin




              The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
              04-22-2024, 07:01 AM
            • seqadmin
              Current Approaches to Protein Sequencing
              by seqadmin


              Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
              04-04-2024, 04:25 PM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 04-11-2024, 12:08 PM
            0 responses
            59 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 04-10-2024, 10:19 PM
            0 responses
            57 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 04-10-2024, 09:21 AM
            0 responses
            53 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 04-04-2024, 09:00 AM
            0 responses
            56 views
            0 likes
            Last Post seqadmin  
            Working...
            X