Seqanswers Leaderboard Ad



No announcement yet.
  • 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

  • #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.


    Latest Articles


    • seqadmin
      Latest Developments in Precision Medicine
      by seqadmin

      Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

      Somatic Genomics
      “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
      05-24-2024, 01:16 PM
    • seqadmin
      Recent Advances in Sequencing Analysis Tools
      by seqadmin

      The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
      05-06-2024, 07:48 AM





    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 01:32 PM
    0 responses
    Last Post seqadmin  
    Started by seqadmin, 05-24-2024, 07:15 AM
    0 responses
    Last Post seqadmin  
    Started by seqadmin, 05-23-2024, 10:28 AM
    0 responses
    Last Post seqadmin  
    Started by seqadmin, 05-23-2024, 07:35 AM
    0 responses
    Last Post seqadmin