Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • yip
    Junior Member
    • Oct 2011
    • 5

    HISAT2: a surprising read shows up in alignment results

    I was spot-checking the output of alignment results from HISAT2 and noticed one read (CATCTTGTTTCACACTAGTCACATTCTTGTTTTAAGTGCCTTTAGTTTTAACAGTTCACTTTTTACAGTACTATTT shown below) that does not exist in my input fastq files. Is this the result of a bug or am I missing anything?

    $ hisat2 -p 40 -x hg19/genome -1 read1.fq.gz -2 read2.fq.gz 2>/dev/null | grep CATCTTGTTTCACACTAGTCACATTCTTGTTTTAAGTGCCTTTAGTTTTAACAGTTCACTTTTTACAGTACTATTT
    HISEQ:287:CCK1AANXX:5:2307:14574:99798 409 chr19 14733987 0 76M = 14733987 0 CATCTTGTTTCACACTAGTCACATTCTTGTTTTAAGTGCCTTTAGTTTTAACAGTTCACTTTTTACAGTACTATTT GGGGGGGGCGGGEGG>FGGGGGGBEGGGD@GEGEGGGGGGGGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGCBCCC AS:i:-5 ZS:i:-5 XN:i:0 XM:i:1 XO:i:0 XG:i:0 NM:i:1 MD:Z:69G6 YT:Z:UP NH:i:2
    $ zcat read1.fq.gz read2.fq.gz|grep CATCTTGTTTCACACTAGTCACATTCTTGTTTTAAGTGCCTTTAGTTTTAACAGTTCACTTTTTACAGTACTATTT
    $ hisat2 -p 40 -x hg19/genome -1 read1.fq.gz -2 read2.fq.gz 2>/dev/null | grep HISEQ:287:CCK1AANXX:5:2307:14574:99798
    HISEQ:287:CCK1AANXX:5:2307:14574:99798 137 chr10 33189363 0 76M = 33189363 0 AAATAGTACTGTAAAAAGTGAACTGTTAAAACTAAAGGCACTTAAAACAAGAATGTGACTAGTGTGAAACAAGATG CCCBCGGGGGGGGGGGGGGGGGGGGFGGGGGGGGGGGGGGGGEGEG@DGGGEBGGGGGGF>GGEGGGCGGGGGGGG AS:i:-5 ZS:i:-5 XN:i:0 XM:i:1 XO:i:0 XG:i:0 NM:i:1 MD:Z:6C69 YT:Z:UP NH:i:2
    HISEQ:287:CCK1AANXX:5:2307:14574:99798 409 chr19 14733987 0 76M = 14733987 0 CATCTTGTTTCACACTAGTCACATTCTTGTTTTAAGTGCCTTTAGTTTTAACAGTTCACTTTTTACAGTACTATTT GGGGGGGGCGGGEGG>FGGGGGGBEGGGD@GEGEGGGGGGGGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGCBCCC AS:i:-5 ZS:i:-5 XN:i:0 XM:i:1 XO:i:0 XG:i:0 NM:i:1 MD:Z:69G6 YT:Z:UP NH:i:2
    HISEQ:287:CCK1AANXX:5:2307:14574:99798 69 chr10 33189363 0 * = 33189363 0 ATAAACATACAGTGGTCTGTTATGGCACTAACTCCCCGTACTCTGCGTTGGTATCAACGCAGAGTACGGGGAGTTA 33<::BFGGGGGGGGDGGGGGGGGGGGGGGGEGGGGGGGGGGGGGGGGFEGGGGGGGGGG@GGEGGGGGGGEGGGG YT:Z:UP

    The other two entries in the sam output with read name of HISEQ:287:CCK1AANXX:5:2307:14574:99798 are from my input fastq files as expected (a read pair).

    I appreciate any comments. Thanks!
  • yip
    Junior Member
    • Oct 2011
    • 5

    #2
    Here's the answer (thanks to Daehwan from the HISAT2 team):
    CATCTTGTTTCACACTAGTCACATTCTTGTTTTAAGTGCCTTTAGTTTTAACAGTTCACTTTTTACAGTACTATTT is the reverse complement of AAATAGTACTGTAAAAAGTGAACTGTTAAAACTAAAGGCACTTAAAACAAGAATGTGACTAGTGTGAAACAAGATG, which is read 1 of HISEQ:287:CCK1AANXX:5:2307:14574:99798.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      New Genomics Tools and Methods Shared at AGBT 2025
      by seqadmin


      This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

      The Headliner
      The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
      03-03-2025, 01:39 PM
    • seqadmin
      Investigating the Gut Microbiome Through Diet and Spatial Biology
      by seqadmin




      The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...
      02-24-2025, 06:31 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 05:03 AM
    0 responses
    16 views
    0 reactions
    Last Post seqadmin  
    Started by seqadmin, 03-19-2025, 07:27 AM
    0 responses
    16 views
    0 reactions
    Last Post seqadmin  
    Started by seqadmin, 03-18-2025, 12:50 PM
    0 responses
    16 views
    0 reactions
    Last Post seqadmin  
    Started by seqadmin, 03-03-2025, 01:15 PM
    0 responses
    185 views
    0 reactions
    Last Post seqadmin  
    Working...