Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • where is real problem, tophat2 code or my fastq files

    Dear FOlks,
    I am so new, an early stage researcher.

    I am using TopHat2 to map the reads, I guess, I am fulfilling all the requirements, my code is

    /usr/local/bin/tophat2 -p 8 -G ~/path/to/Homo_sapiens.GRCh37.72.gtf -o
    ~/path/to/Human_mapping_iPS_s7_rep1
    --splice-mismatches 1 --max-multihits 30 --microexon-search --fusion-search
    ~/path/to/bowtie2_index/hg19
    ~/path/to/myfile.fastq

    I am submitting on grid engine cluster with qsub -l h_vmem=50G [above_script]
    this is showing error as:
    """""TopHat requires all reads be either FASTQ or FASTA. Mixing formats is not supported"""

    I am bit frustrated because my fastq files look fine to me as shown in code

    @SOLEXA-GA05_00009_SRi_AD_MS_BN_VW:7:1:2364:933#ATGAGCA
    NGGCCTTCCCACATTCTTTACACTCATAGGTTTTCTCACCAGTGTGAGTTCTCTTGTGCACAATAAGGTAAGAGCC
    +SOLEXA-GA05_00009_SRi_AD_MS_BN_VW:7:1:2364:933#ATGAGCA
    !454478347;09977778<655476;69;8588380745<75;57495945158::=677976:7674:64763-

    Please help???????
    Manvendra Singh
    PhD student,
    Mobile DNA, MDC, Berlin

  • #2
    My Guess is that there is a problem somewhere in the fastq file, although the first 4 lines look good! A few checks that you can do are:

    1. there are 4 lines for each sequence read in the file, i.e. total count of lines in the fastq file is 4x the total number of reads
    2. every first line of each record starts with @ and every third line starts with +
    3. the length of the quality sequence is the same as the length of the sequence read for every record

    Comment


    • #3
      Looking at the sequence identifiers I wonder if this is old data from a GAII machine. It is then likely in the older illumina (1.3) Fastq format. If that is the case then you may need to add the relevant options for tophat to take that into account.

      From TopHat manual

      --solexa-quals Use the Solexa scale for quality values in FASTQ files.
      --solexa1.3-quals As of the Illumina GA pipeline version 1.3, quality scores are encoded in Phred-scaled base-64. Use this option for FASTQ files from pipeline 1.3 or later.

      Comment


      • #4
        Thanks guys,
        My Problem is figured out. There was a problem with my fastq file
        Manvendra Singh
        PhD student,
        Mobile DNA, MDC, Berlin

        Comment


        • #5
          Originally posted by manvendra7 View Post
          Dear FOlks,
          I am so new, an early stage researcher.

          I am using TopHat2 to map the reads, I guess, I am fulfilling all the requirements, my code is

          /usr/local/bin/tophat2 -p 8 -G ~/path/to/Homo_sapiens.GRCh37.72.gtf -o
          ~/path/to/Human_mapping_iPS_s7_rep1
          --splice-mismatches 1 --max-multihits 30 --microexon-search --fusion-search
          ~/path/to/bowtie2_index/hg19
          ~/path/to/myfile.fastq

          I am submitting on grid engine cluster with qsub -l h_vmem=50G [above_script]
          this is showing error as:
          """""TopHat requires all reads be either FASTQ or FASTA. Mixing formats is not supported"""

          I am bit frustrated because my fastq files look fine to me as shown in code

          @SOLEXA-GA05_00009_SRi_AD_MS_BN_VW:7:1:2364:933#ATGAGCA
          NGGCCTTCCCACATTCTTTACACTCATAGGTTTTCTCACCAGTGTGAGTTCTCTTGTGCACAATAAGGTAAGAGCC
          +SOLEXA-GA05_00009_SRi_AD_MS_BN_VW:7:1:2364:933#ATGAGCA
          !454478347;09977778<655476;69;8588380745<75;57495945158::=677976:7674:64763-

          Please help???????
          For all what I know is there are diff verions of fastq format. you better have to check of the software requirements.

          FASTQ: three main versions, illumina 1.3+, 1.5+ and 1.8+

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Genetic Variation in Immunogenetics and Antibody Diversity
            by seqadmin



            The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
            11-06-2024, 07:24 PM
          • seqadmin
            Choosing Between NGS and qPCR
            by seqadmin



            Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
            10-18-2024, 07:11 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 11-08-2024, 11:09 AM
          0 responses
          230 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 11-08-2024, 06:13 AM
          0 responses
          166 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 11-01-2024, 06:09 AM
          0 responses
          80 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 10-30-2024, 05:31 AM
          0 responses
          27 views
          0 likes
          Last Post seqadmin  
          Working...
          X