Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Mapping using STAR and Bowtie

    Hi all,

    I want to perform an alignment of short reads against a reference genome (human genome from Ensembl). For that I've used STAR, and Bowtie (not Bowtie2). I've some questions about that.

    1- I've used STAR by setting the more important parameters as follows:

    --outFilterScoreMinOverLread 0.4 --outFilterMatchNminOverLread 0.4 --outFilterMatchNmin 15

    Here is the final log file content generated by STAR:

    Number of input reads | 5845995
    Average input read length | 27
    UNIQUE READS:
    Uniquely mapped reads number | 446299
    Uniquely mapped reads % | 7.63%
    Average mapped length | 25.10
    Number of splices: Total | 4043
    Number of splices: Annotated (sjdb) | 622
    Number of splices: GT/AG | 3930
    Number of splices: GC/AG | 108
    Number of splices: AT/AC | 5
    Number of splices: Non-canonical | 0
    Mismatch rate per base, % | 1.74%
    Deletion rate per base | 0.08%
    Deletion average length | 1.36
    Insertion rate per base | 0.04%
    Insertion average length | 1.07
    MULTI-MAPPING READS:
    Number of reads mapped to multiple loci | 4346149
    % of reads mapped to multiple loci | 74.34%
    Number of reads mapped to too many loci | 963965
    % of reads mapped to too many loci | 16.49%
    UNMAPPED READS:
    % of reads unmapped: too many mismatches | 0.00%
    % of reads unmapped: too short | 0.31%
    % of reads unmapped: other | 1.22%
    CHIMERIC READS:
    Number of chimeric reads | 1503
    % of chimeric reads | 0.03%

    According to this, I got a very low rate of uniquely mapped reads (I used the parameter --outFilterMultimapNmax 10).

    How can I improve the uniquely mapped reads rate ?

    2 - I've also used Bowtie to map the same sample to the same reference genome. Bowtie didn't output the number of uniquely mapped reads and the number of multimapped reads as STAR did. I checked the log file generated by Bowtie:

    # reads processed: 5845995
    # reads with at least one reported alignment: 4206928 (71.96%)
    # reads that failed to align: 1639067 (28.04%)
    Reported 98304317 alignments to 1 output stream(s)

    I want to output the number of uniquely mapped reads and multimapped reads from the sam file generated by bowtie. So I have done the following:

    samtools view -Sb myfile.sam > myfile.bam
    samtools view -F 4 myfile.bam | grep -v "XS:" | wc -l

    But it didn't work for me, and I got a wrong number of uniquely mapped reads. I know that those commands work for files generated by bowtie2.

    Do you have any solution able to retrieve the number of uniquely mapped, and multimapped reads from the sam file generated by bowtie ?

    Thanks a lot

  • #2
    Cross-posted: https://www.biostars.org/p/309209/

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Techniques and Challenges in Conservation Genomics
      by seqadmin



      The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

      Avian Conservation
      Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
      03-08-2024, 10:41 AM
    • seqadmin
      The Impact of AI in Genomic Medicine
      by seqadmin



      Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
      02-26-2024, 02:07 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 03-14-2024, 06:13 AM
    0 responses
    33 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-08-2024, 08:03 AM
    0 responses
    72 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-07-2024, 08:13 AM
    0 responses
    81 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-06-2024, 09:51 AM
    0 responses
    68 views
    0 likes
    Last Post seqadmin  
    Working...
    X