Seqanswers Leaderboard Ad



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

  • SpliceMap: all reads less than 50bp after clipping

    I'm using SpliceMap to align my short read data (fastq files with 36 bp pair-end reads). Since it appears that SpliceMap requires reads of length 50 after clipping, all of the reads are being excluded. Does anyone know a way around this?
    I thought about modifying the source code to be fine with reads 36 bp long but it's tricky for me since the the files (e.g. runSpliceMap) don't open as human readable code.
    I'm using version of SpliceMap.

    A snippet of the screen output:

    Preparing the reads!...
    Bases removed from front: 0
    Using as many bases as possible.
    Warning: excluding line temp/read_1_1|2(less than 50bp after clipping)
    Warning: excluding line temp/read_1_2|2(less than 50bp after clipping)
    Warning: excluding line temp/read_1_1|6(less than 50bp after clipping)
    Warning: excluding line temp/read_1_2|6(less than 50bp after clipping)
    Warning: excluding line temp/read_1_1|10(less than 50bp after clipping)
    Warning: excluding line temp/read_1_2|10(less than 50bp after clipping)
    Warning: excluding line temp/read_1_1|14(less than 50bp after clipping)
    Warning: excluding line temp/read_1_2|14(less than 50bp after clipping)
    I'd like to use SpliceMap because it's recommended by IDP-fusion (which is the next step in my analysis) for producing the alignment sam and junction detection bed files.
    Last edited by belgravia; 12-08-2016, 02:42 PM.

  • #2
    I'm Jason, Postdoc in Kin Fai's lab familiar with IDP. Thanks for your question and linking the SeqAnswers post. I'll post up the same response there.

    IDP was written with fairly rigid input requirements but we have several scripts to help open up more options for valid inputs to IDP. We have scripts that allow us to substitute hisat aligned short reads for the splicemap aligned short reads.

    To get the necessary scripts, you can clone our lab's public repository:
    $ git clone

    And then add

    to your python path, i.e.:
    edit ~/.bashrc and add the line
    export PYTHONPATH=$PYTHONPATH:<path to source>/Au-public/iron/pythonlib

    The utilities I refer to are found in

    IDP requires a short read alignment as well as junctions
    of known splices and to add on splices for predicted isoforms and to
    do the quantification.

    We have a script to make the output of hisat alignment compatible with IDP (not yet compatible with
    hisat2, sorry). After hisat alignment you can convert the output sam/bam into a
    splicemap-like sam and a corresponding junction bed file.

    Are the two scripts that can accomplish this. If you have trouble running these please let me know.



    Latest Articles


    • 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
    • seqadmin
      Strategies for Sequencing Challenging Samples
      by seqadmin

      Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
      03-22-2024, 06:39 AM





    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 12:08 PM
    0 responses
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 10:19 PM
    0 responses
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 09:21 AM
    0 responses
    Last Post seqadmin  
    Started by seqadmin, 04-04-2024, 09:00 AM
    0 responses
    Last Post seqadmin