Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • unmapped reads - bowtie

    Hi everybody,

    I am running bowtie to analyse my paired-end Illumina data sets. I have two files one from wild type and one from a mutation we are working with.

    to run bowtie i used this command:
    bowtie -a --best --strata -m 10 -n 2 -l 22 -q --un dilptotal.unmapped -t -p 6 --chunkmbs 256 --max dilptotal.maxHits -S d_melanogaster_fb5_32 -1 dilp_s2_1_sequence.fq -2 dilp_s2_2_sequence.fq dilptotal.sam

    and to run it with the option of trimmed ends I used this command:
    bowtie -a --best --strata -m 10 -n 2 -l 22 -q --un dilptotal_trimmed.unmapped -t -p 6 -3 11 -5 11 --chunkmbs 256 --max dilptotal_trimmed.maxHits -S d_melanogaster_fb5_32 -1 dilp_s2_1_sequence.fq -2 dilp_s2_2_sequence.fq dilptotal_trimmed.sam

    Beside the fact that only later, I read that the options 'best' and 'strata' are not working with paired-ends , I run it as it is and got the following results:
    for the option w.o trimming:
    Read number 24781053 (57.81%)
    Unmapped 17831504 (41.60%)
    Max. hit 255358 (0.60%)
    for the trimmed command:
    read number 28,575,776 (66.66%)
    Unmapped 14,003,690 (32.67%)
    Max. hit 288,449 (0.67%)

    As one can see here I have over two third of my data mapped when I am trimming both ends (while only 57% with the complete reads). If I understand it correctly I can assume that the ends are adapter residues which were needed to be cut out (I run quality control using fastqc software).

    My question here regards the unmapped reads. I can see that there are almost a third of my data which wasn't mapped to the genome. Is there a way of telling why this reads where not mapped.
    These are the reason i could think of for a read to not be mapped:
    1. quality of read is too bad
    2. repetitive elements - more than 10 possible hits (with this preferences).
    3. reads were not found on the reference genome
    4. ???
    Can you please tell me what other possible reasons there can be for not being able to map a read?

    Is there a way of telling what reads are falling into which group?

    What abut the reads which were mapped. As I am working with RNA-Seq data I am interested to know if I have some rRNA residues in my reads.
    Is there a way of finding out how much of this two third mapped reads are really mapped to mRNA and not to some rRNA (which was left over after cleaning/purification of the mRNA) or even to some intergenic region or interons?

    Thanks for any help
    Assa

  • #2
    there is an option to dump the unaligned reads in bowtie.
    you may check to see whether the quality of those reads are bad or not.
    also, if you are aligning mRNA, you shouldnt be aligning against the reference genome. since a transcript would be genome sequence minus introns and stuff like that.
    since bowtie cant align gaps, all those mRNA reads will show as unmapped.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Quality Control Essentials for Next-Generation Sequencing Workflows
      by seqadmin




      Like all molecular biology applications, next-generation sequencing (NGS) workflows require diligent quality control (QC) measures to ensure accurate and reproducible results. Proper QC begins at nucleic acid extraction and continues all the way through to data analysis. This article outlines the key QC steps in an NGS workflow, along with the commonly used tools and techniques.

      Nucleic Acid Quality Control
      Preparing for NGS starts with isolating the...
      02-10-2025, 01:58 PM
    • seqadmin
      An Introduction to the Technologies Transforming Precision Medicine
      by seqadmin


      In recent years, precision medicine has become a major focus for researchers and healthcare professionals. This approach offers personalized treatment and wellness plans by utilizing insights from each person's unique biology and lifestyle to deliver more effective care. Its advancement relies on innovative technologies that enable a deeper understanding of individual variability. In a joint documentary with our colleagues at Biocompare, we examined the foundational principles of precision...
      01-27-2025, 07:46 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 02-07-2025, 09:30 AM
    0 responses
    54 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 02-05-2025, 10:34 AM
    0 responses
    88 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 02-03-2025, 09:07 AM
    0 responses
    70 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 01-31-2025, 08:31 AM
    0 responses
    44 views
    0 likes
    Last Post seqadmin  
    Working...
    X