【bwa mem比对问题汇总】

简介: 【bwa mem比对问题汇总】

bwa mem比对出现:skip orientation RF as there are not enough pairs

This is normal. The Illumina reads orientation is FR. 
You got majority FR and skipped other three.


mapping生成sam文件时出现[mem_sam_pe] paired reads have different names错误

方法一(这种方法我在使用后依旧报错):

bbrename.sh in1=read1.fq in2=read2.fq out1=renamed1.fq out2=renamed2.fq
#To install this package with conda run one of the following: conda install -c bioconda bbmap
对于多个fq文件,可以用以下命令:
while read line
do
      nohup bbrename.sh in=${line}_R1.fq in2=${line}_R2.fq out=${line}_R1.sh.fq out2=${line}_R2.sh.fq &    
done < name.txt #####其中,name.txt指的是包含sample名字文件

方法二 使用repair.sh进行修复:(使用这种方法后fastq文件变小,可能是合并部分文件造成的。

此外, fastq.gz也可以,另外就是记得输出文件要和源文件命名不同,以防覆盖。

while read line
do
      nohup repair.sh in=${line}_R1.fastq in2=${line}_R2.fastq out=${line}_R1.sh.fastq out2=${line}_R2.sh.fastq &    
done < name.txt


相关文章
|
算法 索引
ChIP-seq 分析:数据比对(3)
ChIP-seq 分析:数据比对(3)
277 0
|
数据库 开发者 索引
mem 100%场景优化|学习笔记
快速学习 mem 100%场景优化
136 0
mem 100%场景优化|学习笔记
|
网络协议 Linux API
PF_RING实现分析(2)
4、mmap操作 用户态的接下来调用: ring->buffer = (char *)mmap(NULL, PAGE_SIZE, PROT_READ|PROT_WRITE, MAP_SHARED, ring->fd, 0); 复制代码 进行内存映射。 同样地,内核调用相应的ring_mmap进行处理。 Ring选项结构通过ring_sk宏与sk 建立关联 struct ring_opt *pfr = ring_sk(sk); 复制代码
455 0
|
存储 Linux API
PF_RING实现分析(1)
内核版本:Linux 2.6.30.9 PF_RING版本:4.1.0 最近看了一个PF_RING的实现,看了个大概,发上来大家讨论讨论,共同学习。 一、什么是PF_RING PF_RING是一个第三方的内核数据包捕获接口,类似于libpcap。 二、为什么需要PF_RING 一切为了效率,按照其官方网站上的测试数据,在Linux平台之上,其效率至少高于libpcap 50% - 60%,甚至是一倍。更好的是,PF_RING提供了一个修改版本的libpcap,使之建立在PF_RING接口之上。这样,原来使用libpcap的程序,就可以自然过渡了。
754 0
PF_RING实现分析(3)
刚刚接触PF_RING,学习了,O(∩_∩)O 有个问题请教一下九贱前辈:内核中的PACKET_MMAP跟PF_RING有什么不同呢? 感觉PACKET_MMAP跟transparent_mode=0时的PF_RING原理上是一样的( ⊙ o ⊙ )啊...?
175 0
|
机器学习/深度学习 人工智能 算法
|
SQL NoSQL Oracle
[20180305]手工模拟buffer busy wait.txt
[20180305]手工模拟buffer busy wait.txt --//一般出现buffer busy wait原因,主要是对热块,大量dml操作. --//一种提法:oracle读不会阻塞写,写不会阻塞读,实际上写一定程度会阻塞读,只不过时间很短罢了.
1145 0
|
存储 TensorFlow 算法框架/工具