Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • chayan
    Member
    • Nov 2012
    • 52

    bash loop to count multiple words from multiple corresponding files

    I have generated multiple read mapping files (.sam) generated with bowtie2 "--no-unal" parameter. As the number of files are huge (7 reference genomes mapped against >6000 reads) i want a loop to count the % (or total number) mapped reads for each file.

    HTSeq could be an option but don't know how to loop it for multiple files.
    Any help would be much appreciated.
  • GenoMax
    Senior Member
    • Feb 2008
    • 7142

    #2
    You should consider using featureCounts. You may not need a loop. You can feed all SAM files to featureCounts and generate a read count matrix directly.

    Code:
    featureCounts <options> file1.sam file2.sam file3.sam etc.

    Comment

    Latest Articles

    Collapse

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, 06-09-2026, 11:58 AM
    0 responses
    27 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-05-2026, 10:09 AM
    0 responses
    34 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-04-2026, 08:59 AM
    0 responses
    40 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-02-2026, 12:03 PM
    0 responses
    62 views
    0 reactions
    Last Post SEQadmin2  
    Working...