How To Figure Out The Base On The Special Positions
Assume I have RNASeq data and a list of chromosome loci. Can I use some tools to get an output like this:
chr 2032121 A:20 T:350 C:300 G:55
chr 12098311 A:190 T:0 C:20 G:760
I want to calculate in each loci, how many A/T/C/G in the RNASeq reads. I am using samtools mpileup -l FILE -f hg19.fa in1.bam. But seems that can only identify the base on the reference genome. Thanks for help.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Per Sample Mpileup Output From Merged Bam File
written by thomasvangurp •
-
How To Make Samtools Mpileup Produce Genotype Information At Every Base?
written by komal.rathi
-
Vcf2Fq Not Calling Heterozygous Sites
written by Justin
-
Samtools Mpileup Defined Regions
written by Vova Naumov
-
Allele Depth (~ Dp4) Values For All (Including Multi-Allelic Alt) Possible Alleles Reported?
written by Rm
-
Is There A Tool To Extract The Reference Genomic Sequence From A Given Coordinate?
written by bibb77
-
Problems In Running Samtools Mpileup/Depth For Specific Regions
written by C Shao •
-
Why Does Vcf Aggregate Indels Like This?
written by Jeremy Leipzig
-
Samtools Mpileup With Refrence Addition
written by Tiffani •
-
Why Does Mpileup Skip My Mutation ?
written by Pierre Lindenbaum
Pileup format gives you information about variations, it is shown here, I think this discussion is also highly relevant How to get a summary for a special chromosome position in BAM file
I wrote a script to parse that info from mpileup output a while back: http://blog.nextgenetics.net/?e=56
However, I am not sure that it works 100% correctly. So I would spot check the output.
Thanks Damian. I have a question here: I see if I have 2000 positions, but mpileup will only calculate about 1200 positions, is there a parameter that I can figure out on the all positions. Or does that mean the left 800 positions are in the sequencing gaps?