More and more researchers are now using HLA sequences as part of the genome reference, which means the sequence name will have * and : in the sequence names http://hla.alleles.org/nomenclature/naming.html. While using HLA in the reference the bam file have sequence name (chrom) have : in the name which makes sam.py fail at these location (chrom, i_st, i_end) = i.split(":") .
Will it be possible to update the code to use different delimiter instead of :
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
More and more researchers are now using HLA sequences as part of the genome reference, which means the sequence name will have * and : in the sequence names http://hla.alleles.org/nomenclature/naming.html. While using HLA in the reference the bam file have sequence name (chrom) have : in the name which makes sam.py fail at these location (chrom, i_st, i_end) = i.split(":") .
Will it be possible to update the code to use different delimiter instead of :