Question: Samtools mpileup output miss some bases.
0
gravatar for hxlei613
18 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 • 419 views
ADD COMMENTlink modified 18 months ago • written 18 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: 904 users visited in the last hour