Codice Allineamento FASTA corretto

cd /home/genomica/DATA/GENOLICS/RUN/GENOLICS_R_Leccino_Confronti
cp index* /home/genomica/DATA/GENOLICS/Sonde_FARGA

#Librerie da importare
library(Rsubread)

##Individuazione file da utilizzare per l'allineamento a Leccino
fastqPath_R1 <- list.files("//home/genomica/DATA/GENOLICS/Sonde_FARGA", pattern="fasta$", full=TRUE)

##Creazione dei .BAM di allineamento
align(index="index", maxMismatches=1, readfile1=fastqPath_R1, input_format="FASTA", output_format="BAM")

samtools sort SondeFASTA.fasta.subread.BAM > SondeFASTA.sorted.BAM
samtools index SondeFASTA.sorted.BAM
bamCoverage -b SondeFASTA.sorted.BAM --normalizeUsing RPKM -o SondeFASTA.bw --verbose --numberOfProcessors max/2 --binSize 1