近期使用RNA-seq,在NCBI上获取原始数据后,想使用fastqc检查二代测序数据是否存在问题。
于官网下载fastqc的软件包,并解压至Biosofts文件夹。
运行
fastqc -h
报错,显示无此命令。
起初认为没有配置环境变量,但配置完成后依旧报错。
通过网络查其原因是解压后的FastQC文件夹里虽然存在fastqc文件,但是使用ls -lh 代码发现该文件并无执行权限。
解决办法如下:
在fastqc文件所在目录下运行下列代码
sudo chmod 777 fastqc
再次运行依旧报错
报错信息如下:
Approx 95% complete for output.fastq
Analysis complete for output.fastq
Failed to process file output.fastq
java.lang.IllegalArgumentException: No key called gc_sequence:ignore in the config data
at uk.ac.babraham.FastQC.Modules.ModuleConfig.getParam(ModuleConfig.java:148)
at uk.ac.babraham.FastQC.Modules.PerSequenceGCContent.ignoreInReport(PerSequenceGCContent.java:57)
at uk.ac.babraham.FastQC.Report.HTMLReportArchive.startDocument(HTMLReportArchive.java:331)
at uk.ac.babraham.FastQC.Report.HTMLReportArchive.<init>(HTMLReportArchive.java:84)
at uk.ac.babraham.FastQC.Analysis.OfflineRunner.analysisComplete(OfflineRunner.java:155)
at uk.ac.babraham.FastQC.Analysis.AnalysisRunner.run(AnalysisRunner.java:110)
at java.lang.Thread.run(Thread.java:745)
询问得知虽然我从官网上下载了fastqc软件包解压后,并且修改文件可执行权限,但是由于我曾经用过apt安装系统自带的fastqc,但是系统自带的fastqc是存在bug的,上述情况就是运行系统自带的fastqc是出现的bug
解决这样的问题,我们要先把系统自带的fastqc卸载掉
sudo apt remove fastqc
再将下载的fastqc配置环境变量
echo 'export PATH=~/Biosofts/FastQC:/$PATH'>>~/.bashrc
source ~/.bashrc
再次运行没有报错
问题解决了!
网友评论