diff --git a/wgs.sh b/wgs.sh index be9ba9c..245793f 100644 --- a/wgs.sh +++ b/wgs.sh @@ -28,4 +28,12 @@ soft/megahit --k-list 21 --mem-flag 0 -1 EchE_R1.fastq.gz -2 EchE_R2.fastq.gz -o # Question 4 -soft/prodigal -i contigs/final.contigs.fa -f fasta > genes.fna +soft/prodigal -i contigs/final.contigs.fa -d genes.fna + +# Question 5 + +sed "s:>:*\n>:g" genes.fna | sed -n "/partial=00/,/*/p"|grep -v "*" > genes_full.fna + +# Question 6 + +