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

Useful utility function to print out the Allele as a String since toString prints out * for refs. It was annoying to keep seeing new String(Allele.getBases()). 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@3989 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
ebanks 2010-08-09 20:35:56 +00:00
parent d0ecb8875a
commit 0eeb659aa3
3 changed files with 6 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
java/src/org/broadinstitute/sting/gatk/contexts/variantcontext/VariantContextUtils.java
View File

@ -598,7 +598,7 @@ public class VariantContextUtils {
for (Allele a : inputVC.getAlleles()) {
// get bases for current allele and create a new one with trimmed bases
String newBases = new String(new byte[]{refByte}) + new String(a.getBases());
String newBases = new String(new byte[]{refByte}) + a.getBaseString();
alleles.add(Allele.create(newBases,a.isReference()));
}
@ -610,7 +610,7 @@ public class VariantContextUtils {
List<Allele> newGenotypeAlleles = new ArrayList<Allele>();
for (Allele a : inAlleles) {
if (a.isCalled()) {
String newBases = new String(new byte[]{refByte}) + new String(a.getBases());
String newBases = new String(new byte[]{refByte}) + a.getBaseString();
newGenotypeAlleles.add(Allele.create(newBases,a.isReference()));
}
else {

3
java/src/org/broadinstitute/sting/gatk/walkers/sequenom/PickSequenomProbes.java
View File

@ -129,7 +129,6 @@ public class PickSequenomProbes extends RodWalker<String, String> {
RODRecordList snpList = snpMaskIterator.seekForward(GenomeLocParser.createGenomeLoc(contig,offset-200,offset+200));
if ( snpList != null && snpList.size() != 0 ) {
Iterator<GATKFeature> snpsInWindow = snpList.iterator();
int i = 0;
while ( snpsInWindow.hasNext() ) {
GenomeLoc snp = snpsInWindow.next().getLocation();
// we don't really want to mask out multi-base indels
@ -157,7 +156,7 @@ public class PickSequenomProbes extends RodWalker<String, String> {
else if ( vc.isInsertion() )
assay_sequence = leading_bases + "[-/" + vc.getAlternateAllele(0).toString() + "]" + (char)ref.getBase() + trailing_bases;
else if ( vc.isDeletion() )
assay_sequence = leading_bases + "[" + new String(vc.getReference().getBases()) + "/-]" + trailing_bases.substring(vc.getReference().length()-1);
assay_sequence = leading_bases + "[" + vc.getReference().getBaseString() + "/-]" + trailing_bases.substring(vc.getReference().length()-1);
else
return "";

6
java/src/org/broadinstitute/sting/utils/genotype/vcf/VCFWriter.java
View File

@ -140,7 +140,7 @@ public class VCFWriter {
// REF
alleleMap.put(vc.getReference(), "0");
String refString = new String(vc.getReference().getBases());
String refString = vc.getReference().getBaseString();
mWriter.write(refString);
mWriter.write(VCFConstants.FIELD_SEPARATOR);
@ -148,13 +148,13 @@ public class VCFWriter {
if ( vc.isVariant() ) {
Allele altAllele = vc.getAlternateAllele(0);
alleleMap.put(altAllele, "1");
String alt = new String(altAllele.getBases());
String alt = altAllele.getBaseString();
mWriter.write(alt);
for (int i = 1; i < vc.getAlternateAlleles().size(); i++) {
altAllele = vc.getAlternateAllele(i);
alleleMap.put(altAllele, String.valueOf(i+1));
alt = new String(altAllele.getBases());
alt = altAllele.getBaseString();
mWriter.write(",");
mWriter.write(alt);
}