zzh
/
gatk-3.8
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Make transcript2info capable of handling b37/hg19 contigs 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@3915 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
chartl 2010-08-02 17:32:08 +00:00
parent c203e0fb02
commit 2bc69572cb
1 changed files with 6 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
java/src/org/broadinstitute/sting/playground/gatk/walkers/annotator/TranscriptToInfo.java
View File

@ -916,6 +916,7 @@ public class TranscriptToInfo extends RodWalker<TreeMap<String, String>, TreeMap
final char lastChromChar = txChrom.charAt(txChrom.length() -1);
switch( Character.toLowerCase(lastChromChar) ) {
case 'm':
case 't': // for hg19 or b36
key = 0;
break;
case 'x':
@ -925,7 +926,11 @@ public class TranscriptToInfo extends RodWalker<TreeMap<String, String>, TreeMap
key = 24;
break;
default:
key = Integer.parseInt(txChrom.substring(3));
if( txChrom.startsWith("chr") ) {
key = Integer.parseInt(txChrom.substring(3));
} else {
key = Integer.parseInt(txChrom);
}
break;
}