EI-CoreBioinformatics / portcullis

Splice junction analysis and filtering from BAM files
https://ei-corebioinformatics.github.io/portcullis/
GNU General Public License v3.0
38 stars 9 forks source link

Why would max_min_anc be larger than half of read length? #66

Open ritetsu opened 1 month ago

ritetsu commented 1 month ago

Hello! After reading the metrics info, I assume that max_min_anc represents the maximum value of all shorter side of junction reads, which means it should be no more than half of the read length (otherwise it won't be the shorter anchor of a junction read). However, I got max_min_anc=41, which is >38 while dealing with 76bp read length bam file. Does mismatch or something else have influence on max_min_anc value? Or did I misunderstood the definition of max_min_anc?

Thanks, LZ