|
|
@@ -585,7 +585,7 @@ impl Fasta {
|
|
|
Fasta { fa: IndexedFasta::from_file(path).expect("Error opening fa") }
|
|
|
}
|
|
|
fn get_sequence (&self, contig: &str, start: usize, end: usize) -> String { // end is included
|
|
|
- let start = start - 1;
|
|
|
+ // let start = start - 1;
|
|
|
let end = end - 1;
|
|
|
let chr_index = self.fa.fai().tid(contig).expect("Cannot find chr in index");
|
|
|
if (end as i64 - start as i64) < (0 as i64) {
|
