basic version of depth has been tested and works!

flatgfa/src/cmds.rs

@@ -303,11 +303,11 @@ pub fn depth(gfa: &flatgfa::FlatGFA) {
             };
         }
     }
-    for seg in gfa.segs {
-        let seg_id: u32 = seg.name as u32;
-        match depths.get(&seg_id) {
-            Some(depth) => println!("{}\t{}", seg_id, depth),
-            None => println!("{}\t{}", seg_id, 0)
+    for (id, seg) in gfa.segs.iter().enumerate() {
+        let name: u32 = seg.name as u32;
+        match depths.get(&(id as u32)) {
+            Some(depth) => println!("{}\t{}", name, depth),
+            None => println!("{}\t{}", name, 0)
         };
     }
 }

flatgfa/src/main.rs

@@ -64,6 +64,7 @@ enum Command {
     Stats(cmds::Stats),
     Position(cmds::Position),
     Extract(cmds::Extract),
+    Depth(cmds::Depth),
 }
 
 fn main() -> Result<(), &'static str> {
@@ -126,6 +127,9 @@ fn main() -> Result<(), &'static str> {
             let store = cmds::extract(&gfa, sub_args)?;
             dump(&store.view(), &args.output);
         }
+        Some(Command::Depth(_)) => {
+            cmds::depth(&gfa);
+        }
         None => {
             // Just emit the GFA or FlatGFA file.
             dump(&gfa, &args.output);