|
|
@@ -400,7 +400,7 @@ impl NeoContig {
|
|
|
"POS", "MAPQ", "CIGAR", "RNEXT", "PNEXT", "TLEN", "SEQ",
|
|
|
"QUAL", "ALN"];
|
|
|
|
|
|
- let stdout = cmd!("echo", "-e", format!("\">sequence\n{}\"", self.contig))
|
|
|
+ let stdout = cmd!("echo", "-e", format!("\">sequence_{}\n{}\"", self.name.to_string(), self.contig))
|
|
|
.pipe(cmd!("bwa", "mem", fasta_ref_path, "-"))
|
|
|
.stdout_capture()
|
|
|
.stderr_capture()
|
|
|
@@ -408,9 +408,9 @@ impl NeoContig {
|
|
|
|
|
|
let mut alignments: Vec<Sam> = Vec::new();
|
|
|
stdout.split("\n").for_each(|line|{
|
|
|
- println!("{:?}", line);
|
|
|
- if line.len() > 9 {
|
|
|
- if &line[0..9] == "sequence" {
|
|
|
+ // println!("{:?}", line);
|
|
|
+ if line.len() > 8 {
|
|
|
+ if &line[0..9] == "sequence_" {
|
|
|
|
|
|
let mut qname: Option<String> = None;
|
|
|
let mut flag: Option<i32> = None;
|
