m

src/lib.rs

@@ -308,11 +308,6 @@ impl SequencesGraph {
                 consensus_sequence.push(base);
             }
         }
-
-        
-        // println!("RES {:?}", String::from_utf8_lossy(&consensus_sequence));
-        // println!("last = {:?}", consensus_sequence.last().unwrap().to_owned());
-        // println!("last = {:?}", consensus_sequence.last().unwrap().to_owned() == b'N');
         consensus_sequence
     }
 
@@ -452,7 +447,7 @@ impl NeoContig {
                         }
                     }
                     
-                    if let (Some(qname), Some(flag), Some(rname), Some(pos), Some(cigar), Some(sequence)) = (qname, flag, rname, pos, cigar, sequence) {
+                    if let (Some(qname), Some(flag), Some(rname), Some(pos), Some(cigar), Some(_)) = (qname, flag, rname, pos, cigar, sequence) {
                         alignments.push(Sam::new(qname, flag, rname, pos, cigar, self.contig.clone(), &fa));
                     }
                 }
@@ -680,7 +675,7 @@ fn matched_range (cigar: &str, flag: &i32, matched_seq: &mut str) -> (Range<usiz
     (range_ref, range_query, ref_cigar_string)
 }
 
-fn revcomp(dna: &str) -> String{
+pub fn revcomp(dna: &str) -> String{
     let mut rdna: String = String::with_capacity(dna.len()); 
     for c in dna.chars().rev() {
         rdna.push(switch_base(c));