Fix stats (#20)

* fixed bam stats bug

* updated version

bamdash/__init__.py

@@ -1,3 +1,3 @@
 """interactively visualize coverage and tracks"""
 _program = "bamdash"
-__version__ = "0.3"
+__version__ = "0.3.1"

bamdash/command.py

@@ -300,7 +300,7 @@ def main(sysargs=sys.argv[1:]):
     # dump track data
     vcf_track_count, bed_track_count, gb_track_count = 0, 0, 0
     if args.dump:
-        pd.DataFrame.from_dict(stat_dict, orient="index").to_csv("bam_stats.tabular", sep="\t", header=False, index=True)
+        pd.DataFrame.from_dict(stat_dict, orient="index").to_csv(f"{args.reference}_bam_stats.tabular", sep="\t", header=False, index=True)
         if track_data:
             for track in track_data:
                 if track[1] == "vcf":

bamdash/scripts/data.py

@@ -51,6 +51,9 @@ def make_title_string(parsed_bam, coverage_df, reference, min_cov):
     stat_dict = {}
     # get bam stats for correct chrom
     bam_stats = parsed_bam.get_index_statistics()[0]
+    for bam_stats in parsed_bam.get_index_statistics():
+        if bam_stats[0] == reference:
+            break
     # pop dict
     stat_dict["reference"] = bam_stats[0]
     stat_dict["reference length (bp)"] = parsed_bam.get_reference_length(reference)