Header Leaderboard Ad


unmapped reads - bowtie



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
      A Brief Overview and Common Challenges in Single-cell Sequencing Analysis
      by seqadmin

      ​​​​​​The introduction of single-cell sequencing has advanced the ability to study cell-to-cell heterogeneity. Its use has improved our understanding of somatic mutations1, cell lineages2, cellular diversity and regulation3, and development in multicellular organisms4. Single-cell sequencing encompasses hundreds of techniques with different approaches to studying the genomes, transcriptomes, epigenomes, and other omics of individual cells. The analysis of single-cell sequencing data i...

      01-24-2023, 01:19 PM
    • seqadmin
      Introduction to Single-Cell Sequencing
      by seqadmin
      Single-cell sequencing is a technique used to investigate the genome, transcriptome, epigenome, and other omics of individual cells using high-throughput sequencing. This technology has provided many scientific breakthroughs and continues to be applied across many fields, including microbiology, oncology, immunology, neurobiology, precision medicine, and stem cell research.

      The advancement of single-cell sequencing began in 2009 when Tang et al. investigated the single-cell transcriptomes
      01-09-2023, 03:10 PM