Question: Samtools mpileup output miss some bases.
0
gravatar for hxlei613
22 months ago by
hxlei61380
Bangladesh
hxlei61380 wrote:

Alignment files are sorted by Picard. Then I use Samtools mpileup function to create a pileup file of mapped data prepared for DNA methylation–level calculation. The output seems good. But I find it miss some base at start. What about the first five bases ? Why does it start from the sixth base of chrM ?

This is my command : samtools mpileup -f <genome.fa> <read1_val_1.sort.bam> > <read1_val_1.pileup>

Thank you very much !

=========================================================================

chrM 6 C 0

chrM 7 A 0

chrM 8 G 0

chrM 9 G 0

chrM 10 T 4 .... U=RR

chrM 11 C 4 .... _:fc

chrM 12 T 4 .... dBfg

chrM 13 A 7 T...^K.^!.^). dBec033

chrM 14 T 11 .......^K.^K.^K.^K. cFgh0U3C@C@

chrM 15 C 11 ....TT...., XFVc03C<c@:< p="">

chrM 16 A 11 .....C....,-1c `F\03OY]W8

chrM 17 C 11 ....TAA.... iFek033igii

......

=========================================================================

mpileup samtools • 533 views
ADD COMMENTlink modified 22 months ago • written 22 months ago by hxlei61380
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1586 users visited in the last hour