Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Hisat2: Differing amount of reads in mates

    Hello all,

    I am aligning paired-end reads using Hisat2. Unfortunately, some of the fastq files were corrupted in transfer. I was able to recover them using a gzip recovery protocol, but was left with about half the data. I have previously used Tophat2 to align these "recovered" fastq files, but am getting this error when trying to align with Hisat2:
    Code:
    Error, fewer reads in file specified with -2 than in file specified with -1
    How is Tophat2 dealing with this differently than Hisat2?

    To fix this, I am trying repair.sh from the bbtools package to keep the reads that do have pairs in both files and output singletons to a seperate file and then try using all three with Hisat2.

    Although I can't seem to find reference to this error anywhere else, any advice on how I should deal with this?

  • #2
    Don't use paired end reads with unmatched pairs since that can lead to erroneous discordant alignments. repair.sh is the way to go.

    Comment


    • #3
      Thanks Genomax. I would assume in the Tophat2 run with these "recovered" files there were erroneous discordant alignments that would be discarded and not affect overall alignment.

      I was successfully able to run Hisat2 following repair using the new paired-end files, including the singleton file. A high percentage of the singletons mapped uniquely - surely these are not erroneous alignments?

      Comment


      • #4
        Firstly, since the underlying issue is a data corruption issue, I would strongly recommend you re-download the corrupted data. As it is, your results will be not be reproducable from the original data.

        Originally posted by ronaldrcutler View Post
        A high percentage of the singletons mapped uniquely - surely these are not erroneous alignments?
        If only one of the .1.fq.gz/.2.fq.gz pair was corrupted, then there will be a large number of singleton reads from the file that was successfully copied. You would expect the unique mapping rate of these singleton reads to be only slightly less than the the unique mapping rate for the paired end reads. The difference between the two will be due to the aligner being able to use the partner read to disambiguate the mapping location for the pair end reads, but not the singletons.

        TLDR: that behaviour is expected; they're probably correct; redownload the correct data before continuing

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Techniques and Challenges in Conservation Genomics
          by seqadmin



          The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

          Avian Conservation
          Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
          03-08-2024, 10:41 AM
        • seqadmin
          The Impact of AI in Genomic Medicine
          by seqadmin



          Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
          02-26-2024, 02:07 PM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 03-14-2024, 06:13 AM
        0 responses
        32 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-08-2024, 08:03 AM
        0 responses
        71 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-07-2024, 08:13 AM
        0 responses
        80 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-06-2024, 09:51 AM
        0 responses
        68 views
        0 likes
        Last Post seqadmin  
        Working...
        X