Skip to main content
Search type Search syntax
Tags [tag]
Exact "words here"
Author user:1234
user:me (yours)
Score score:3 (3+)
score:0 (none)
Answers answers:3 (3+)
answers:0 (none)
isaccepted:yes
hasaccepted:no
inquestion:1234
Views views:250
Code code:"if (foo != bar)"
Sections title:apples
body:"apples oranges"
URL url:"*.example.com"
Saves in:saves
Status closed:yes
duplicate:no
migrated:no
wiki:no
Types is:question
is:answer
Exclude -[tag]
-apples
For more details on advanced search visit our help page
Results tagged with
Search options questions only not deleted user 4001

Questions specific to interacting with and post-processing sequence alignments using the SAMtools package. For questions specifically about formats SAM, BAM or CRAM use tags sam or bam,

3 votes
1 answer
2k views

Running htseq-count over BAM files

I have read the htseq-count documentation as well as samtools and figured that the following command should work: samtools sort -on file.bam | samtools view -h - | htseq-count --type=exon --idattr=gene_id … Just note that the newer versions of htseq-count don't require sorted .bam files and, so following the editing of the .gtf file, all I needed to do was convert it to .sam using samtools view -h and then …
3 votes
2 answers
316 views

Frequency of specific viral sequence in .BAM or .fastq

So, my question is: if I were to use samtools - or any other program - to convert the .bam files to .fastq, is there any possibility that I could lose such information - e.g. viral seqs - that would have …
1 vote
2 answers
99 views

Using preprocessing/alignment functions on the server

I have a series of raw RNA sequencing files: .bam and wish to utilize samtools (or others) to begin the data preprocessing steps defined in many workflows. … : -bash-samtools: command not found When I try to install the package, say samtools, by running sudo apt install samtools, I get an error stating I do not have permission to do such. …