Entering edit mode
20 months ago
hellokwmin
•
0
Hello,
After calculation of sequencing depth using samtools, I am trying to print out the result of sequencing depth using python script. But, I would like to make a limitation for y-axis to be "200".
The python script I am using as follosw: Probably, I need to revise one that I make a bold. Please, help me.
--script ----
import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
import sys
input_data = sys.argv[1]
output_data = sys.argv[2]
column_names
data = pd.read_csv(input_data, sep='\t', header=None, names=column_names)
x = data['基因组位置']
y = data['覆盖深度']
# ************
average_coverage_depth = y.mean()
maximal_depth = y.max()
minimal_depth = y.min()
total_genome_length = x.iloc[-1]
# ************
plt.figure(figsize=(12, 6))
plt.plot(x, y, color='blue', linewidth=1, linestyle='-')
plt.title('Sequencing Depth and Coverage Map', fontsize=18, pad=20)
plt.xlabel('Genomic Position (bp)', fontsize=14)
plt.ylabel('Sequencing Depth (×)', fontsize=14)
plt.xticks(fontsize=12)
plt.yticks(fontsize=12)
plt.grid(True, linestyle='--', alpha=0.5)