|
@@ -167,7 +167,7 @@ pub fn get_start_end_qual(
|
|
|
let rend = record.reference_end() as i32;
|
|
let rend = record.reference_end() as i32;
|
|
|
|
|
|
|
|
if rstart >= start && rstart < stop {
|
|
if rstart >= start && rstart < stop {
|
|
|
- println!("start {rstart}");
|
|
|
|
|
|
|
+ // println!("start {rstart}");
|
|
|
if record.mapq() >= mapq {
|
|
if record.mapq() >= mapq {
|
|
|
let index = rstart - start;
|
|
let index = rstart - start;
|
|
|
*results.get_mut(index as usize).unwrap() += 1;
|
|
*results.get_mut(index as usize).unwrap() += 1;
|
|
@@ -175,7 +175,7 @@ pub fn get_start_end_qual(
|
|
|
}
|
|
}
|
|
|
|
|
|
|
|
if rend >= start && rend < stop {
|
|
if rend >= start && rend < stop {
|
|
|
- println!("{rend}");
|
|
|
|
|
|
|
+ // println!("{rend}");
|
|
|
if record.mapq() >= mapq {
|
|
if record.mapq() >= mapq {
|
|
|
let index = rend - start;
|
|
let index = rend - start;
|
|
|
*results.get_mut(index as usize).unwrap() += 1;
|
|
*results.get_mut(index as usize).unwrap() += 1;
|