Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • gringer
    replied
    What for your experience are the best parameters for short read mapping?
    Anything that gets you quicker to re-analysing your target sample with a non-SOLiD system. You shouldn't be mapping SOLiD reads at all, especially not 35bp ones; interpreting the results is just too confusing (even for the people writing code to do the mapping).

    Leave a comment:


  • colindaven
    replied
    I recently found another non-commercial aligner for SOLiD data - it is called subread. I have only tried for Illumina data but it's very good there.

    Again, I wouldn't use tophat for such short read data.

    Leave a comment:


  • Ilarius
    replied
    I did use tophat2 with bowtie1, anyway if I let the default -g option (i think it is 20) it takes ages to finish.

    As regards shrimp I didn't find any working link to download it!

    Leave a comment:


  • colindaven
    replied
    In general, in comparison to Illumina PE reads don't expect high mapping percentages. From memory I think 50-70% was pretty good, as the data were not extensively filtered (on the 5500xl machine at least).

    You probably don't want to use Tophat for such short reads as I don't think it will be able to find anchors/seeds successfully.

    My best experiences with SOLiD were with the commercial program NovoalignCS. You might be able to get a free trial.

    Otherwise, Shrimp2 was also not too bad.

    Why not use just bowtie1 instead of tophat2 with bowtie1 ?

    Leave a comment:


  • Ilarius
    started a topic mapping 35bp solid reads with tophat unsuccessful

    mapping 35bp solid reads with tophat unsuccessful

    Hi, I recently downloaded some solid short reads (35 bp), single end, and convert them in cfasta with ahi-dump.

    I found that there are very few mappers that support solid read mapping in color space and finally i decided to use tophat using bowtie1 (which supported color space mapping).
    Unfortunately I less less than 30 % of the reads mapping.
    My previous experience has always been with 100 bp paired end reads, so I am not completely sure if I am using the right parameters for short read mapping:

    tophat -G $annotation --segment-length 17 --segment-mismatches 1 -g 1 --bowtie1 -p 12 -o $OUTPUT_FOLDER --coverage-search --color --quals $genome $READ $QUAL

    maybe is the segment-mismatch and maximum multi-hits too restrictive? (-g is 20 by default). What for your experience are the best parameters for short read mapping?

    Thanks
    Ilario

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 on Modified Bases...
    Today, 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
37 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-10-2024, 10:19 PM
0 responses
41 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-10-2024, 09:21 AM
0 responses
35 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-04-2024, 09:00 AM
0 responses
54 views
0 likes
Last Post seqadmin  
Working...
X