Skip to content

Commit

Permalink
fmt
Browse files Browse the repository at this point in the history
  • Loading branch information
alethomas committed Nov 14, 2023
1 parent a8616d5 commit fdb401b
Show file tree
Hide file tree
Showing 2 changed files with 1 addition and 5 deletions.
2 changes: 0 additions & 2 deletions workflow/scripts/masking.py
Original file line number Diff line number Diff line change
Expand Up @@ -48,7 +48,6 @@ def get_base_count(pileupcolumn):
for pileupread in pileupcolumn.pileups:
# TODO Check pileupread for missing bases
if not pileupread.is_del and not pileupread.is_refskip:

read_base = pileupread.alignment.query_sequence[pileupread.query_position]

bases.append(read_base)
Expand Down Expand Up @@ -102,7 +101,6 @@ def mask_sequence(sequence, coverages, base_counts):
covered_postions = coverages.keys()

for position, base in enumerate(sequence):

if position not in covered_postions:
# TODO Check why there are postions that are not covered by any reads and are not Ns
# sequence[position] = "N"
Expand Down
4 changes: 1 addition & 3 deletions workflow/scripts/quality-filter.py
Original file line number Diff line number Diff line change
Expand Up @@ -31,9 +31,7 @@ def get_identity(quast_report_paths: List[str]) -> dict:
sample = path.dirname(report_path).split("/")[-1]

# load report
report_df = pd.read_csv(
report_path, delimiter="\t", names=["name", "value"]
)
report_df = pd.read_csv(report_path, delimiter="\t", names=["name", "value"])
report_df.set_index("name", inplace=True)
# select genome fraction (%)
try:
Expand Down

0 comments on commit fdb401b

Please sign in to comment.