For loop rna sequencing alignment
0
0
Entering edit mode
9 weeks ago
kcarery • 0

Hey guys... I am trying to do a for loop to align my samples to a reference file in terminal.....can someone help me...I am getting this error:

EXITING because of fatal input ERROR: could not open readFilesIn=Read1

This is the code I put in

for f1 in `cat MatchedFastqList.txt` ;do SAMPLE=$(echo ${f1}| sed "s/_1\.fastq\//"); echo ${SAMPLE}.bam; echo 
STAR --runThreadN 16 --genomeDir /folder/to/project/ncbi-genomes-2022-09-19/
GCF_000001405.40_GRCh38.p14_genomic.fna --sjdbGTFfile /folder/to/projects/GCF_000001405.40_GRCh38.p14_genomic.gtf 
--readFilesIn ${f1} --outReadsUnmapped Fastx ; done
Star • 325 views
ADD COMMENT
0
Entering edit mode

not an answer but : you'd better use a workflow manager like snakemake or nextflow.

..I am getting this error

what is "this error" ?

ADD REPLY
0
Entering edit mode

Thank you for your feedback. I am self-taught learning coding. What is a workflow manager? Is it super beneficial for longer pipelines?

ADD REPLY
0
Entering edit mode

Is it super beneficial for longer pipelines?

definitively

https://pubmed.ncbi.nlm.nih.gov/28398311/

Nextflow enables reproducible computational workflows

https://pubmed.ncbi.nlm.nih.gov/35037206/

Pipeline Automation via Snakemake

ADD REPLY
0
Entering edit mode

404: Error Not Found

(Now seriously, you forgot to copy the error ;) )

So far the only thing that I can see is that the done in the end of the loop is missing.

ADD REPLY
0
Entering edit mode

Sorry, I was rushing I got this error

EXITING because of fatal input ERROR: could not open readFilesIn=Read1

It also is not going down the list to each file...

ADD REPLY
0
Entering edit mode

Can you copy-paste the content of MatchedFastqList.txt?

ADD REPLY
0
Entering edit mode

It contacts samples that are set up like this "Sample345.fastq" and its all listed

ADD REPLY

Login before adding your answer.

Traffic: 1363 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6