@@ -456,7 +456,7 @@ headerIndexBam <- function(bam.list) {
456
456
map.prog <- try(strsplit(header [[1 ]]$ text $ ' @PG' [[1 ]], ' :' )[[1 ]][2 ],
457
457
silent = T )
458
458
if (class(map.prog ) != " try-error" ) {
459
- map.style <- grepl(' tophat|bowtie|bedtools|star' , map.prog ,
459
+ map.style <- grepl(' tophat|bowtie|bedtools|star|hisat ' , map.prog ,
460
460
ignore.case = T )
461
461
if (map.style ){
462
462
v.map.bowtie [i ] <- TRUE
@@ -467,29 +467,7 @@ headerIndexBam <- function(bam.list) {
467
467
v.map.bowtie [i ] <- FALSE
468
468
next
469
469
}
470
- # if(estiMapqStyle(bam.file)){
471
- # warning(sprintf("Aligner for: %s cannot be determined. Style of
472
- # standard SAM mapping score will be used.", bam.file))
473
- # v.map.bowtie[i] <- FALSE
474
- # }else{
475
- # warning(sprintf("Aligner for: %s cannot be determined. Style of
476
- # Bowtie-like SAM mapping score will be used. Would you mind to tell us what
477
- # aligner you are using?", bam.file))
478
- # v.map.bowtie[i] <- TRUE
479
- # }
480
- }
481
- # else {
482
- # cat("\n")
483
- # if(estiMapqStyle(bam.file)){
484
- # warning(sprintf("Aligner for: %s cannot be determined. Style of
485
- # standard SAM mapping score will be used.", bam.file))
486
- # v.map.bowtie[i] <- FALSE
487
- # }else{
488
- # warning(sprintf("Aligner for: %s cannot be determined. Style of
489
- # Bowtie-like SAM mapping score will be used.", bam.file))
490
- # v.map.bowtie[i] <- TRUE
491
- # }
492
- # }
470
+ }
493
471
warning(sprintf(" Aligner for: %s cannot be determined. Style of
494
472
standard SAM mapping score will be used. Would you mind submitting an issue
495
473
report to us on Github? This will benefit people using the same aligner." ,
0 commit comments