Re: [Samtools-help] About BAM->CRAM->BAM

2016-05-20 Thread jkb
On 2016-05-18 18:35, Lin, Chih-Hsu wrote: > Thanks for your reply. I used UCSC hg19 from Illumina iGenomes. Is the value > 8192 abnormally large for NM tag? This is a BAM from RNA-seq, so there is an > intron of 8192 nt. How does NM compute? Does NM consider CIGAR "N" into > account? Thanks

Re: [Samtools-help] About BAM->CRAM->BAM

2016-05-18 Thread Lin, Chih-Hsu
rom: James Bonfield Sent: Wednesday, May 18, 2016 3:25:17 AM To: Lin, Chih-Hsu Cc: samtools-help@lists.sourceforge.net Subject: Re: [Samtools-help] About BAM->CRAM->BAM On Tue, May 17, 2016 at 09:17:20PM +, Lin, Chih-Hsu wrote: > During the conversion, BAM->CRAM->BAM, using samto

Re: [Samtools-help] About BAM->CRAM->BAM

2016-05-18 Thread James Bonfield
On Tue, May 17, 2016 at 09:17:20PM +, Lin, Chih-Hsu wrote: > During the conversion, BAM->CRAM->BAM, using samtools 1.3.1, I found the NM > tags were changed. Does anyone have solution to that? Can you give us a concrete example please; an alignment record along with the relevant @SQ line so w