Read counts at specific locations in SAM files
I have sam files created using hisat2 of RNAseq reads assembled to a reference of 7 genes (not whole genome). What I need to be able to do is for a given position in the assembly I need too: 1. Count the total # of reads at that position 2. Count of reads of each nucleotide at that position
I have ~30 Rnaseq libraries and ~70 positions I need to assay. This was done in a allohexaploid species so I am looking at positions know to be polymorphic between genomes.
Any ideas of how best to do this?
• 967 views
•
link
written
by
fogarty.mc •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
RNASeq vs TermSeq
written by npb27 •
-
How can I count the number of UMIs and reads at each position after mapping with a reference sequen…
written by naeem40thju •
-
How can I count the barcodes of alligned reads at each position?
written by naeem40thju •
-
no alignment from plasma derived small RNA
written by memento1984_1 •
-
Question: Error HISAT2 output for HTseq-count
written by 11yj3312 •
-
How to find mismatches in an assembly
written by mary •
-
count reads with samtools view
written by schelarina •
-
How to do coverage counting for each genomic position in R ?
written by MAPK
-
how to count un-mapped reads at specific position using samtools?
written by mangfu100
-
Nucletotide distribution, at each position, in a .sam/.bam file ?
written by remi.maglione •
bam-readcount will do this for you. Samtools depth/mpileup would be other options.