Hi,
I would like to get coverage per a set of genomic ranges, with a little complication that I need coverage over T, G, A, and C provided separately. My idea to do it is first add feature names to mpileup file (with bedtools intersect), and then do something akin to R dyplyr::summarize(). But maybe there is bash alternative for that?
The second step of my struggle is to count specific mismatches based on the mpileup code. This I thought of doing in R, because I know how, but maybe someone could help me to get started with awk on that. A one-liner to count the number of occurrences of "g" in column 5 (see below) and print this number instead in the same column would help to get me started.
Of course, if there is a more efficient way to accomplish the task - let me know (I am sure there is)!
slam_500_spike 390 A 46 .,..,.,.g..,,g,.,,,,.g.,,,,,,,g,,....,,,,,,,. FFFF:FFFFFFFFFJJFFFJFJJJJFJJJJFFJJJJJJJFFJJJJJ
0 answers
No answers yet.
Log in to answer this question.
-
Separate letters in the same column into different rows
written by genomes_and_MGEs •
-
Annotating reads with the name of the feature to which they map
written by lechu •
-
counting the number of reads in FASTQ files
written by Ric
-
Extract allele count from gnomAD using R
written by Nicolas Rosewick
-
metagenomic assembly coverage, for multiple samples
written by willnotburn •
-
Calculating coverage from ref and alt columns?
written by PsFeathers •
-
Find patterns in DNA sequence
written by kindlychung •
-
Questions regarding extracting ref or non-ref allele counts in bam files.
written by mangfu100
-
Counting Read Support At A Particular Position Of A Bam File
written by Shaun Jackman
-
Samtools Mpileup Directional Coverage
written by User 1874 •
