ELF>@@H\$Ld$ALl$Hl$HLt$L|$HHL0IDw H@HT$`D$hG DGot=H\$xH$L$L$L$L$HĨtH5DD$1HHDD$tHIDAH5HcHIAHT$M4(L$D$LHO47HT$L$18Z@H@I,/Hl$8KDCHt$8A1EEE1>5-$t H9|C9D)D)Љ1)ЃA1E)AI‰D$$D$hSAAE)D+t$$t$()9}5)Ɖ47t$$D+t$$))Ɖ40A)ED+\$$D\$(C LcHcILC0LD$HcSLHHHc{HD$@KLD$~_H5LT$@1AAHcAB HcD)AEA)DA:9CHc|$(D9|$$I}YHcD$$ATD$$1H5LHcЃATH҃HcAA9A:Hc{LT$Hc{ILT$HD$0LL$`KHT$@Mt$(HD$LL,$ALT$f{DcK1Lt$(1Lt$8L|$0H\$8.t?t:t A9HcE DAuA9D+T$$BD9}xLcCt@u9uJtGDE8EFDǐCAED\=E<8DHEEHDE;E0Lt$(H\$8s~"1HcЃHLA*@9ƈD$LS~(1HcI 6*T5@K9H5IZHLT$LT$HLT$H|$@LT$LH|$0L1E)1҉A)AI҉T$$+D$$D$(au@|$Lt3@t.f;HpH|$LL$1HT$1|$L@ń@t1ub@K1@@;SHcʃI4 @HcI4 1(DDA@E9|@(@;{|BZ1Hc{LT$Hc{If{IDKLT$KLT$X1LT$8E1Lt$(H\$PIcA4tDt@t;tu>AD9@뿃uA>9}UEˉD+\$$L$8E1ADt9}-HcEtAuAE9uHL$0D4L$8HcTAW9}*HcD\Ft E8ECރ9E|HcTATS9|FHcD\Ft8E8ECރ9E} HcAEA8AFӃT9|H5IZHLT$LT$Lt$(H\$PLT$XS~/1Hcй@A4HT:@8v@(N@ 9CLLT$LKLT$1AD$$ED$(*@B1@AWAVAUATUSH8T$,WHG0DD$,T$ӋWO HcۃHHcL4H,HEL,(ID$,D$,D_|$-D$(D$E1D$$D$ AD$_tmthyGEDD$IL$9L$AAAuE IcDL DHHʄDA< uD@1A )D)щʍ D^B BHcɉ HHDA<9A@AB :D)EA)Dt:t5HcAT v'D$9tt"º!CƒD$ D$$A9\EEDD$E~&T$(0B HcAL A9uT$(E}E~(IcÀ<tB1BHcҀ<t+A9ECC@AD$(2\$ T$ *T$1f(*9^Yu*T$$*L$,YX*D$(^Xf.w8fW_f(f(\f(^Qf.z!uYX,H8[]A\A]A^A_f( $ $]D$(D$$@AWIAVAUATUSHHt$(T$XdH%(H$1L$\OG HW0D$lHcHH HHT$@HL$fAHEG T$XHL$pL|$H1E1D$<HL$PAgfffEǃT$0ZthtcƒtEƄ HL$HD$<A;D$<0 HcD$HЄA>D@DA )D)ˉ1aB AHAHcHt$HȄA, >A@@AB *D)EA)Dt9tuBHc҃ AwAOAG IW0HcHcH H A9GD$ EEAG;D$ 7GEDHH{HH HHH HHH HHH HHH HpHsDD$H HCDD$HQ)w11)~3H|$PȃA))0H҉u؅H;HcHD8H;CHCHHy3HCHct4p@4HH uHCHt=0H HQH;SVHC^HCEHHT$(IcD,EDt$ 1E2fHCD,0HCA9H ~~BD=HD,EH3HNH;KrHH{HH HHH HHH HHH HHH HpHsHT$H3HCHT$HNkEDt$ EDD$l1DD$lD$ EAG;D$ I}L|$HH HQH;SHC0HCH HH{HH HHH HHH HHH HHH HpHsDD$H HCDD$HQRH5LAG HHtEHDL$lAD9"HIO0HH81HLLD$lH5iLH5LAG HBHtUL`H1LHLiL;+LsHIO0HMMH81HLH LCH5ZLDH5LLH0HH{HH HHH HHH HHH HHH HrHsH;HCHcŃDp-HH{HH HHH HHH HHH HHH HrHsDD$H;HCDD$HcŃDp-H HQH;SsjHC0HCHHH{HH HHH HHH HHH HHH HrHsH HCHQHH{HH HHH HHH HHH HHH HpHsDD$H HCDD$HQ>MLsH11HcI9[|I4<; $ > $ >   I : ;  : ;I8 : ; : ; I8 I !I/ &I< & : ; I8  : ;  : ; I : ;  : ; I 8  : ; I 8 : ;I : ; I.? : ; ' I 4 : ; I.? : ;' I : ;I.? : ; ' I@: ; I : ; I!4: ; I"4: ; I #4: ; I$4: ; I%4: ; I& U' U(4: ;I)4: ;I*.? : ;' I@ +: ;I,4: ; I -.: ; ' I .: ; I/.? : ; ' @01RUX Y 11241 34141RUX Y 5 61X Y 7 8.? : ;' I@9: ;I:4: ;I;4: ;I <4: ;I =4: ;I>1X Y?1X Y@.? : ; ' @ A4: ; I? < B4: ;I? < C4: ; I?  r8intii1-o b# # # # # #( #0 #8 #@ #H #P #X "#` $#h &b#p *b#t ,p#x 0F# 1T# 2# 6# ?{# H# I# J# K# L-# Nb# P#   # # b#v  o   ( &[ 'b )i 1? 4Mb d j s b#fd b# # # # b# s#$ b#4 b#8 b#< b#@ #H #P #X #` !#h $#p $#x b  %k )fpr *fpw +X $ %b# &# ' # ' #x ,# 0b# 1b# 2# 3#( 4#0 5b#8 6b#< 7b#@ 8#H 9#P :6HW5 X# Y5# Z;# [# [# [#( \-#0 \-#8 ]#@5^ tid#pos#bin5#5#5#5# 5# # # # #L8X # b# b#$ b#( X#0* 'u{^ibam b#x# #A $l-#m-#s# ]d #e #bw!b#"/111d1-`b``gggb < bE bref b!kb!ib!bwb!xb!yb!ybb!yeb!xbb!xeb"bw#b$;%cE "]!bqX!zqX$X& !opb!lb'%sX%rX%qX!seqX%bqX#@& (X(X')opb)lb')opb)lb*,bw +b,+ref,b bref b!seqX$X$;%cE !ib!xb!yb!mmb,qbw$,lenbw#b$b!t'!jb!lb!opb'!zb!c1b!c2b-Lb* .cLb.sL* %bufN0 %lOb%xOb$ @ -@bd .c@b.s@* /`bref b b,seqXw";w%cE !ib!xb!yb!ub%str* $ ,nm &!j$b!l$b!op$b0@ 81Z 1Q 0@ 51Z 1Q 0 51 1 '2 3 3 4@ P1Z 1Q '!z'b!c1(b!c2(b0 .1 1 '2 3 3 4@ P1Z 1Q 4@ .1Z 1Q 0 D'1 1 '2 3 3 4@ P1Z 1Q &!jHb!lHb!opHb'!zKb!c1Lb!c2Lb5$X5$[X6`11&@$gX'#kb6s117!q{X81b91b915)c3b:3b)tid3b)ret3b;len3b<3bw<3bw<3bw <3bw,:3b<3bw:3b)fp4<4w;fai5w)ref6<6T<6T=b7>L11>H1 >I11 > QX10 1% 5|)qub> P10 1% > O10 1% >T11>U11>V?11>Wf11>X11>Y11>Z11>[11>\)11> RP10 1% >bw11?p11@wbref bAAAA;AIb ?D B4C0  A/]@ samtools/usr/include/bits/usr/include/usr/lib/gcc/x86_64-linux-gnu/4.4.5/includebam_md.c.pysam.cstring3.hkstring.hstdio2.hstdlib.hstddef.htypes.hstdio.hlibio.hfaidx.hstdint.hknetfile.hbgzf.hbam.hsam.hkprobaln.hpysam.hgetopt.h +rH>ZL.Y8L'p@p8D8xTNxJR~NtIK'x7L5v:>v~~<XL.5.KX5.KL6XW=/=H>>2/=YWY%B?XAJ%xfA<YeK_o>Kuc+X+,$8CXlc<$o<;=+r<X=$X=J( t.Jvfs;uWK1WYI0gWYI;=I z81X"t\x~KSIy7O3a͔b gY-er9JY-e)uYYxt r+uYㅈ/0eWYgE`&<xXgV}Ɂ_Sk=zJf&JL4Zz. /=^^zXgY+-e:Bz<qJX'-e-&XAAiJl<XK^J!<h<.gJwX  rf=I6 (KBmJWX*fVJ/fQX/g q VX= u%g*B)XG<8<=I6#pOKI=+-YeM/3J-e*/;!z+ XIf$=I6 փDBiJOKnJ<=nJg $OKnJf=nJ/ J=I6w'XH<%+OKI==I6B>XB<fV>#X<f[v#'#Y=I6=I6pKI==I6I=I6(XJ" .tt~<~ .~ XJ ~X~X&~$*8\*a xS+g~X JK'sY:ZZ]ST Zo/K:=Kju~XR6XZ=~X}"}X1~BQZQNMMDRGEqreuNhbSC:n:Adbhu -[bam_fillmd1] different NM for read '%s': %d -> %d [bam_fillmd1] different MD for read '%s': '%s' -> '%s' [bam_fillmd] unrecognized option '-%c' Usage: samtools fillmd [-eubrS] Options: -e change identical bases to '=' -u uncompressed BAM output (for piping) -b compressed BAM output -S the input is SAM with header -A modify the quality string -r compute the BQ tag (without -A) or cap baseQ by BAQ (with -A) -E extended BAQ for better sensitivity but lower specificity [bam_fillmd] input SAM does not have header. Abort! [bam_fillmd] fail to find sequence '%s' in the reference. ?Zd;_@V-?+w+_ w[U[eSeUS-GSSZSqvSSS_ S[T[}]}T]]] _ ]AQAu\\\\ _ \[QiQU#0U07QQrUQ&3X39PdqXUX`QPU P.CPCLQTP-PXP)P7|PQP_nUnQU0mQmoUoQPQ)Q1[Q[{U{QUUPF Q Q0PHLPP4 9 P[XIUX\X\XXY\Y 9 X[mV^{VP-S-GPGSVPSP%)PVVZUZqSqxUxSUSUSU$S$5P5[U[SSUUP[YgnYYYY[[d>[Hn[[[[/ 9 [[ZaZ.PpPZZZ9 F P[_\______ F _^^#ZZ`QPQtPPP%4P@PP 4 PF L PQ _ P[^ ^^^^^!G^[_ ^[`PiPPP[UirUwUHXUUUU`iX X03P9PPPP07T<TTT0Q7RQ Q%MQ]oQQ4Y;Y|YY` j U` j Tp r wr t wt v wv x w x y w(y z w0z ~ w8~ ' wp Up l T ' Tp Q V ' V ] ' ]  ^o r ^g ^ ^ ' ^P l P [ Q [ P l [ ' [ l Z Z@ l Q Q Q @ Qg l Q QH l ar c c a" ' a o So r Pr S P S P 4 S4 P Pg l P P J XQ l X X r Px P P P  P& P Pg l P P Q \ Q Q Q @ Qg l Q QQ r P P& t P P& P Pg l P P U  U* l U U0 2 w2 7 w7 9 w9 ; w ; < w(< = w0= D w8D 5w0 z Uz __0__Jf_sx_}__0_5_0 R T Q$QwQ(QHkQ \0\\QQQ)@QJf\fQ\E\EHQQQ/5\0 u Q0 m REMPPjQQRQQ&5Q/Q X h_hm]_X[0[0_[P=[t_tPJ_Jf[fs_x}_}PX_[[&5[0_0@P@rX_ ^TYT^0Y0^Y^1J^JfYf^^/5^ T?TqTTTiTTiTT?TT1T TT&-TE~TTiTT_PX4=X=DPD^XP&5P V VVV?VV&EV V/5V XXHXT_X_^0vXX$X0rXX1^1JXfX^XX^E^EXXrXr^X^/5^ 6PSPwP0PP@JP PP VwVV V0VV0VVP2V;V1VV&EVVPV/5VRQQTFQR R Q$QQQQwQ(QHkQ03QQQQfxQ}QEHQQQ$6PP0PTPwPP03P ZZ GZXZZ{ZZwZZ)ZZZ&5ZEZZrZZRP2R;BRRPRUQQUKQ07UUVVV[PkV V0VVVJfVVV&EVVPV/5V9RDFQkQTQERRTS0SS3SJfSX0XXJfXXX&5XIP!PP PJWP$QOQQQ!QQQQ&5QHP{P P V0VkVVJfVVV&EV/5VT_ V$V$,P,V&EV/5Vj VVPEV5V;@P@BwBKwKPwPUw UYw(Y\w0\ew8e w@U$])] ]J e ] ]@T&^)^ ^J ^ ^PP)wPP@P-PaS V)!S!V9S J VJ S V SS(P(S J S S Q SkuP  P ' PX"\) \[ V)V VJ V V(__ J _e _&^^ J ^ ^;QPQ^PP U T Qv< bam_prob_realn_coreK bam_prob_realn bam_cap_mapQd bam_fillmd1_core`bam_fillmdbam_fillmd1Rbam_nt16_nt4_table, pH9 F   0`q -  o 0  X  0 J0@ (1?\_f?0?0 0 00`n] e 00? 00? 0?Bei`cG]0E}0E}}0OR?D@&/&/_IO_read_ptrmemcpysize_ttamr_shortbuftamwmode_wisizetamFilebaq_flagtypepath_IO_buf_baseknetFile_slong long unsigned intdata__srcGNU C 4.4.5ctrl_fdthresbam_header_t__tamFile_tresponseint16_tlong long intsigned charleftqual_filenobam_nt16_table_IO_read_endconfstatelong int_flagsis_diff_cur_columndoublen_cigarl_text_lock_old_offsetbam1_tsamfile_tcache_sizebam_prob_realn_corem_datais_readyold_nm_ibam_prob_realntarget_nameportbam1_core_tstrcatis_realnstdinpasv_portuncompressed_blockunsigned intclip_qmax_responsefprintffaidx_t__streamlong unsigned int_IO_write_ptrBGZF_sbuftarget_lenshort unsigned intoptargsize_cmdfile_descriptor_IO_save_baseseek_offsetstrcpyno_reconnectold_md_flags2kpa_par_told_nmis_sam_instdoutretrhttp_hostcigarblock_lengthflt_flagl_qnamebamFile__faidx_t_offsetmtidrghtoptindl_qseq_IO_write_endkpa_par_defblock_offsetcompressed_block_size_IO_lock_t_IO_FILE__kstring_t__off_t__nptrrg2lib_modefloat_IO_marker_poserroropen_mode_markersclip_lpasv_ipatoiint32_tunsigned charbam_fillmd1_IO_buf_endextend_baqshort int_chainknetFilekputcflagcompressed_block_vtable_offsetdata_lenbam_nt16_nt4_tableFILEcorekputwdictcharblock_addresstextn_targetsbam_fillmd1_corehostis_uncompressedfpoutpysamerrapply_baq_next__off64_toffset_IO_read_basebam_fillmd_IO_save_endhash__fmtbam_cap_mapQ__pad1__pad2__pad3__pad4__pad5_unused2uncompressed_block_size__destowned_filemode_ruint8_t_IO_backup_basekstring_tmax_nmargccapQmposint64_t__lencompress_levell_auxargv_IO_write_base/data/martha/PythonPackages/pysam-0.7.4cachen_textis_bam_outheaderuint32_tsamtools/bam_md.c.pysam.cfile_sizeGCC: (Ubuntu/Linaro 4.4.4-14ubuntu5) 4.4.5zRx  _ M^K@ 0TBBB B(A0A8DpW4BCBB B(A0A8GX8BGBCB CB(CA0BA8ED.symtab.strtab.shstrtab.rela.text.data.bss.debug_abbrev.rela.debug_info.rela.debug_line.rodata.str1.1.rodata.str1.8.rela.rodata.rodata.cst8.debug_loc.rela.debug_pubnames.rela.debug_aranges.debug_ranges.debug_str.comment.note.GNU-stack.rela.eh_frame @ ȟ & !,0!10!D$v?x$UZ;P`a2B)p20B$Ex E F6|p?}00o}p0߄0Ë,P0 `f     !&+0 5 : ? 8E K Q pW !] 'c i %o u {  H x   0 #  _ #.=BN` ]p jns0 @!*3;EQZ bam_md.c.pysam.c.LC0.LC1.LC2.LC3.LC4.LC6.LC7.LC8.LC9.LC10.LC11.LC13.LC12.LC15.LC27.LC14.LC17.LC18.LC19.LC20.LC21.LC22.LC23.LC24.LC25.LC16.LC26.LC28bam_prob_realn_core_GLOBAL_OFFSET_TABLE_kpa_par_defbam_aux_getcallocmemcpybam_nt16_nt4_tablebam_nt16_tablekpa_glocalbam_aux_appendfreebam_aux_delbam_prob_realnbam_cap_mapQlogsqrtbam_fillmd1_corereallocbam_aux2ipysamerr__fprintf_chkbam_aux_drop_other__ctype_toupper_loc__stack_chk_failbam_fillmdgetoptoptargstrtol__strcat_chkoptindsamopenfai_loadsamwritesamreadfai_fetchfai_destroysamclosebam_fillmd1. 233j4|54 6 43 6: 7448`x9:::::0;449: : 7D h >|    ?   4v 7A::wA 7sAAA 3(B@ CKZDe;q93 CD; 9)13<EAAnAAF1GE I CDs JK JKP ZLo MN!N MO4M@ZPeQ:R0="L")L5 C<#KDU$aDk%wD&D'D(D)D*D+D : :# :- S5 T? TR ,\ Lq Cx - D C . D Gf 0 @    p! ) . ; VB GI P W ^ xl (q | U    1   f a h   W( a7 F  U d "s    h @    b c   / ? O _ p w  O {    L E Y   ?  9+ 6 I P ol dx M _ !   2  0    - ; I W Re #     G  "  0 > L 3Z h v   v Q        [ ' B   I  S | k 4 6    V. < J _ tj <|   * M  { % ^ j  |   A = M     :  & 1 > I Q _ Y j 9y ~   1 8 o    e   E m! -+ J 9 Q _ d n s  R @ C "   p  y ! 0 D> M Y ` a j z " E  p ' h  <  r    -  p+ 8 F Ke +p )z n  @ L  c    5 A f m 0 u 5}  B #   #  -    % / A B S ` !n |#w P?    $$ %  P % &3 r'A (J R 0g y Z) ) ?     * + }  ( P9 F,F ,T J-Y f -t t. .|  x/  /=^   0  @ 8!)A&IEZ 0b tn@v ~ 0 * `1 \ 1 d2 @ 2 3 64    &" #1 < 4A /P D[ 4j $5o 6 5 / p6H>^Yaw 522OOee%{D{Lks .J 6e Un ] |   U   %6 H6  k6   < ( uE S `6i M   $( ,048<@"D&HL.P2T6X:\>`BdFhJlNpRtVxZ|^bDj<0vz,$(   D` Xp 0 @