bowtie2, and python script
Hello, I am trying to do bowtie2 alignment for 100 fastq files. To run bowtie2 alignment I am trying to write a python script but I am stuck at writing ''for loop''. Please help me. Thank you
• 2,867 views
•
link
updated
by
Ram
•
written
by
qudrat.nii •
•
written
by
qudrat.nii •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
bowtie2, and python script
written by qudrat.nii •
-
For loop to go over all files in a folder to execute after qc
written by ravigenomeinfo •
-
how to run seqkit in a loop with multiple files
written by Kumar
-
bowtie2; using for multiple fastq files, linux loop code
written by Farzaneh •
-
How to calculate bowtie2 run time using python
written by John
-
Python script to find error rate
written by zaidnab •
-
Using trimmomatic on multiple single-end read files
written by armel.donkpegan •
-
Bowtie2 throws error no. 139
written by anu014
-
Linux for loop for 2 inputs and 4 outputs for Trimmomatic fastq quality trimming
written by Rashedul Islam
-
Bowtie2 allowing mismatches syntax issue
written by aggiegurl2009 •
A simple bash script would probably serve you better, or a Makefile. Alternatively a workflow manager such as Snakemake which is python-ish. Please show what you tried.
import subprocess folder = "/Desktop/RNA-Seq" fastq_file = [] try: for files in folder: if files.endswith('fastq'): fastq_files.append(files) subprocess.run("bowtie2 -x folder/bt2_idx -1 f1.q -2 f2.q -S f.sam") except: print('some error')
What is error you're getting here? What are the
f1.qandf2.q. Also, please format your code blocks.f1.q and f2.q is fastq file. Actually I am not able to write the script which can take pairs of files as input for bowtie2
Here is a snippet that can help you do that:
https://stackoverflow.com/questions/312443/how-do-i-split-a-list-into-equally-sized-chunks
However, you should listen to other's advices here and look into shell or workflow managers.
using bash for multiple files: bowtie2; using for multiple fastq files, linux loop code