How to count positively and negatively aligned reads in a bam file?

I have a bam file and I would like to count the number of plus and minus strand aligned tothe reference sequence. How can I get this data from a bam file? Thanks for your help in advance.

bam

• 1,173 views • link

written 8.4 years ago by MAPK

Hello MAPK!

Questions similar to yours can already be found at:

Forward Strand Or Reverse Strand

We have closed your question to allow us to keep similar content in the same thread.

If you disagree with this please tell us why in a reply below. We'll be happy to talk about it.

Cheers!

PS: Answers in the other thread should help you.

• 1 views • link

written 8.4 years ago by Ram

0 answers

No answers yet.

More posts like this

Assembly of DNA from .bam file
written 6.9 years ago by Triple Nipple •

From a .bam file of RNA sequence reads that are aligned to a mouse reference genome, how can I get a new bam file that …
Extraction of batch subject sequence from Blast result.
written 8.2 years ago by kabir.deb •

``` ~/ncbi-blast-2.5.0+/bin/makeblastdb -in sub.fna -dbtype nucl -out sub_DB -parse_seqids ~/ncbi-blast-2.5.0+/bin/tblastn -query que.fasta -db sub_DB -outfmt 6 -out result_0.01_out -evalue 0.01 -num_threads 24 ~/ncbi-blast-2.5.0+/bin/blastdbcmd -db sub_DB …
What cigar to use to count reads with last base mismatch in Rsamtools?
written 8.4 years ago by MAPK

I have been struggling with this problem for long time and would appreciate if somebody answer some questions I have below. I need to count …
How to use scanbamparam and CIGAR from Rsamtools to count specific reads?
written 8.4 years ago by MAPK

Hi, I have used default bowtie option which allows for three mismatches to get bam file aligned with my reference sequence. I want to know …
How do we find the number of sequences with terminal mismatches?
written 8.5 years ago by MAPK

I have created a bowtie aligned .sam file with this command: `bowtie -S -p 18 /media/owner/ref.fa test1.fastq test1.sam` which I then converted to .bam file …
How can I plot alignment between fastq reads and reference sequence?
written 8.5 years ago by MAPK

I have a fastq file and a fasta file with one sequence. I want to plot a figure to show all the reads in fastq …
BigWig file of strand-specific RNA-seq data
written 8.7 years ago by yuabrahamliu •

Hi all, I began to process strand-specific RNA-seq data from dUTR method library these days. I added the flag `--rna-strandness RF` to the default `hisat2` …
Is there any reason to analyze plus/minus strands separately in bisulfite seq/methylation data?
written 10.3 years ago by biohack92

I have bisulfite sequencing data that consists of read count data. The data is split based on strand orientation (plus & minus). What sort of …
Program to extract sequences corresponding to the miRNA alignment with flanking sequences from both…
written 12.2 years ago by shaanbdu •

I am working on miRNA identification in plants and for this I initially ran a blast with EST sequences as queries against the mine database. …
Blast Nucleotide Vs. Nucleotide, Plus-Minus Results Only
written 14.0 years ago by Burlappsack •

<p>Hi Folks, <br/> I am interested in running a blastn search to search for complementary sequence (potential biological function) between two groups of sequences, (Query …