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

now output in the exact format that works with sequenom software. 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@1798 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
jmaguire 2009-10-08 20:06:27 +00:00
parent 98e3a0bf1a
commit 8fdb8922b8
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
java/src/org/broadinstitute/sting/playground/gatk/walkers/PickSequenomProbes.java
View File

@ -80,9 +80,9 @@ public class PickSequenomProbes extends RefWalker<String, String>
if (snp != null)
{
String snp_string = new String(leading_bases) + new String("[" + refBase + "/" + snp.getAlternativeBaseForSNP() + "]") + new String(trailing_bases);
String fasta_string = new String(">" + context.getLocation().toString() + "_" + ref.getWindow().toString());
return fasta_string + "\n" + snp_string + "\n";
String assay_sequence = new String(leading_bases) + new String("[" + refBase + "/" + snp.getAlternativeBaseForSNP() + "]") + new String(trailing_bases);
String assay_id = new String(context.getLocation().toString() + "_" + ref.getWindow().toString()).replace(':', '_');
return assay_id + "\t" + assay_sequence + "\n";
}
else
{