|
122 | 122 | //! let mut read_counter: u32 = 0; |
123 | 123 | //! let mut seq_region: BTreeMap<String, (u32,u32)> = BTreeMap::new(); |
124 | 124 | //! let mut record_vec: Vec<Record> = Vec::new(); |
125 | | -//! loop { |
| 125 | +//! loop { |
126 | 126 | //! match records.next() { |
127 | 127 | //! Some(Ok(mut record)) => { |
128 | 128 | //! //println!("next record"); |
@@ -309,7 +309,7 @@ macro_rules! embl { |
309 | 309 | for rec in reader.records() { |
310 | 310 | match rec { |
311 | 311 | Ok(r) => { |
312 | | - println!("this is r {:?}", &r); |
| 312 | + //println!("this is r {:?}", &r); |
313 | 313 | vec.push(r); |
314 | 314 | } |
315 | 315 | Err(e) => panic!("Error reading record: {:?}", e), |
@@ -722,37 +722,36 @@ where |
722 | 722 | //collects the DNA sequence and translations on the correct strand |
723 | 723 | if stra == -1 { |
724 | 724 | if cod > 1 { |
725 | | - println!( |
726 | | - "reverse strand coding start more than one {:?}", |
727 | | - &iterablecount |
728 | | - ); |
| 725 | + //println!("reverse strand coding start more than one {:?}", |
| 726 | + // &iterablecount |
| 727 | + //); |
729 | 728 | if sto < record.sequence.len() { |
730 | 729 | sliced_sequence = &record.sequence[sta + cod..sto + 1]; |
731 | 730 | } else { |
732 | 731 | sliced_sequence = &record.sequence[sta + cod..sto]; |
733 | 732 | } |
734 | 733 | } else { |
735 | | - println!("record sta {:?} sto {:?} cod {:?} stra {:?} record.seq length {:?}", &sta, &sto, &cod, &stra, &record.sequence.len()); |
736 | | - println!( |
737 | | - "sliced sta {:?} sliced sto {:?} record.id {:?}", |
738 | | - sta, sto, &record.id |
739 | | - ); |
740 | | - println!( |
| 734 | + //println!("record sta {:?} sto {:?} cod {:?} stra {:?} record.seq length {:?}", &sta, &sto, &cod, &stra, &record.sequence.len()); |
| 735 | + //println!( |
| 736 | + // "sliced sta {:?} sliced sto {:?} record.id {:?}", |
| 737 | + // sta, sto, &record.id |
| 738 | + //); |
| 739 | + //println!( |
741 | 740 | "iterable count is {:?} reverse strand codon start one", |
742 | 741 | &iterablecount |
743 | | - ); |
744 | | - println!("this is the sequence len {:?}", &record.sequence.len()); |
| 742 | + //); |
| 743 | + //println!("this is the sequence len {:?}", &record.sequence.len()); |
745 | 744 | if sto < record.sequence.len() { |
746 | 745 | sliced_sequence = &record.sequence[sta..sto + 1]; |
747 | 746 | } else { |
748 | 747 | sliced_sequence = &record.sequence[sta..sto]; |
749 | 748 | } |
750 | | - println!("iterable count after is {:?}", &iterablecount); |
| 749 | + //println!("iterable count after is {:?}", &iterablecount); |
751 | 750 | } |
752 | 751 | let cds_char = sliced_sequence; |
753 | 752 | let prot_seq = translate(&revcomp(cds_char.as_bytes())); |
754 | 753 | let parts: Vec<&str> = prot_seq.split('*').collect(); |
755 | | - println!("this is the prot_seq {:?}", &prot_seq); |
| 754 | + //println!("this is the prot_seq {:?}", &prot_seq); |
756 | 755 | record |
757 | 756 | .seq_features |
758 | 757 | .set_counter(key.to_string()) |
|
0 commit comments