Get identical mapped nucleotide length from bowtie2 generated SAM files
I used bowtie2 to map my reads to a reference genome and get the SAM files returned. I am wondering if there is a way to get the length of the identical mapped nucleotide length. Basically I want to get a hamming distance or identical mapping vs total read length ratio from the SAM files.
1 answer
find out the answer in another post: How To Find Out The Mismatchs Of An Alignment Entry In The Sam File?
Log in to answer this question.
More posts like this
-
warning[vg::giraffe]: Cannot cluster reads with a fragment distance smaller than read distance
written by cassiwatt •
-
vg call coverage/suppport
written by Danny •
-
Calculating Average Coverage or Read Depth for a Sequence (WES)
written by Matt01 •
-
Pair ended samples cannot map to the genome
written by SayHey •
-
Estimate edit distance to a reference for a subset of bases in reads (bam files)
written by AQ7 •
-
how to get total read count, mapped read count and unmapped read count from multiple sam or bam fi…
written by Bioinfonext
-
bbsplit not using all reads in library
written by ghd1998 •
-
ratio mRNA:pre-mRNA in total RNA-seq data
written by JindrichK •
-
strand handling of of reference sequence for index
written by yifangt86 •
-
samtools mapping problem
written by BSP •
