Skip to content

Yingying-Buynak Submission#5

Open
mbuynak wants to merge 27 commits into
lyy005:masterfrom
yunluyingying:master
Open

Yingying-Buynak Submission#5
mbuynak wants to merge 27 commits into
lyy005:masterfrom
yunluyingying:master

Conversation

@mbuynak

@mbuynak mbuynak commented Nov 17, 2017

Copy link
Copy Markdown

No description provided.

Comment thread Exercise11.py
outfile1.close()
outfile2.close()
outfile3.close()
vcffile.close() No newline at end of file

Copy link
Copy Markdown
Owner

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Good job

Comment thread Search Codes
cat transporter.hits | grep -v "#" | wc -l

./hmmsearch --tblout sporecoat.hits sporecoat.hmm Verrucomicrobia.fasta
cat sporecoat.hits | grep -v "#" | wc -l

@lyy005 lyy005 Dec 11, 2017

Copy link
Copy Markdown
Owner

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

The goal of this exercise is to practice for loop. It would be good to use a loop in Linux command:
Also, make sure you only print out the proteome, e-value of the alignment
-0.1 pts

for file in *.fasta
do
outSigma=$file.Sigma.hits
outSpore=$file.Spore.hits
outTrans=$file.Trans.hits
hmmsearch --tblout $outSigma sigma.hmm $file
hmmsearch --tblout $outSpore sporecoat.hmm $file
hmmsearch --tblout $outTrans transporter.hmm $file
done

cat *.hits | grep -v "#" | awk '{print $1,$3,$5}' | sed -E 's/tr|[A-Z0-9]+|[A-Z0-9]+_9//g' > hmmOut.txt

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants