【Linux blast比對(duì)】fastq轉(zhuǎn)fasta比對(duì)到參考序列

nohup less -S ../sample-1_L3_2.fq.gz |awk '{if(NR%4 == 1){print ">" substr($0, 2)}}{if(NR%4 == 2){print}}' > ./sample-1_L3_2.fasta &

$vi run_blast1.sh 

blastn \
 -query sample-1_L3_2.fasta \
 -db artificial.blastdb \
 -out sample_results2.xls \
 -outfmt 6 \
 -num_threads 8 

$nohup sh run_blast1.sh &

• 人面猴

  序言：七十年代末，一起剝皮案震驚了整個(gè)濱河市，隨后出現(xiàn)的幾起案子逞泄，更是在濱河造成了極大的恐慌涯竟，老刑警劉巖晃危，帶你破解...

  沈念sama閱讀 206,839評(píng)論 6贊 482
• 死咒

  序言：濱河連續(xù)發(fā)生了三起死亡事件，死亡現(xiàn)場(chǎng)離奇詭異裙品，居然都是意外死亡帝火，警方通過查閱死者的電腦和手機(jī)，發(fā)現(xiàn)死者居然都...

  沈念sama閱讀 88,543評(píng)論 2贊 382
• 救了他兩次的神仙讓他今天三更去死

  文/潘曉璐 我一進(jìn)店門风罩，熙熙樓的掌柜王于貴愁眉苦臉地迎上來糠排，“玉大人，你說我怎么就攤上這事超升∪牖拢” “怎么了？”我有些...

  開封第一講書人閱讀 153,116評(píng)論 0贊 344
• 道士緝兇錄：失蹤的賣姜人

  文/不壞的土叔 我叫張陵室琢，是天一觀的道長乾闰。 經(jīng)常有香客問我，道長盈滴，這世上最難降的妖魔是什么涯肩？ 我笑而不...

  開封第一講書人閱讀 55,371評(píng)論 1贊 279
• ?港島之戀（遺憾婚禮）

  正文 為了忘掉前任，我火速辦了婚禮雹熬，結(jié)果婚禮上宽菜，老公的妹妹穿的比我還像新娘谣膳。我一直安慰自己竿报，他們只是感情好，可當(dāng)我...

  茶點(diǎn)故事閱讀 64,384評(píng)論 5贊 374
• 惡毒庶女頂嫁案：這布局不是一般人想出來的

  文/花漫 我一把揭開白布继谚。 她就那樣靜靜地躺著烈菌，像睡著了一般。 火紅的嫁衣襯著肌膚如雪花履。 梳的紋絲不亂的頭發(fā)上芽世，一...

  開封第一講書人閱讀 49,111評(píng)論 1贊 285
• 城市分裂傳說

  那天，我揣著相機(jī)與錄音诡壁，去河邊找鬼济瓢。 笑死，一個(gè)胖子當(dāng)著我的面吹牛妹卿，可吹牛的內(nèi)容都是我干的旺矾。 我是一名探鬼主播，決...

  沈念sama閱讀 38,416評(píng)論 3贊 400
• 雙鴛鴦連環(huán)套：你想象不到人心有多黑

  文/蒼蘭香墨 我猛地睜開眼夺克，長吁一口氣：“原來是場(chǎng)噩夢(mèng)啊……” “哼箕宙！你這毒婦竟也來了？” 一聲冷哼從身側(cè)響起铺纽，我...

  開封第一講書人閱讀 37,053評(píng)論 0贊 259
• 萬榮殺人案實(shí)錄

  序言：老撾萬榮一對(duì)情侶失蹤柬帕，失蹤者是張志新（化名）和其女友劉穎，沒想到半個(gè)月后，有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體陷寝，經(jīng)...

  沈念sama閱讀 43,558評(píng)論 1贊 300
• ?護(hù)林員之死

  正文 獨(dú)居荒郊野嶺守林人離奇死亡锅很，尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...

  茶點(diǎn)故事閱讀 36,007評(píng)論 2贊 325
• ?白月光啟示錄

  正文 我和宋清朗相戀三年，在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了盼铁。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片粗蔚。...

  茶點(diǎn)故事閱讀 38,117評(píng)論 1贊 334
• 活死人

  序言：一個(gè)原本活蹦亂跳的男人離奇死亡，死狀恐怖饶火，靈堂內(nèi)的尸體忽然破棺而出鹏控，到底是詐尸還是另有隱情，我是刑警寧澤肤寝，帶...

  沈念sama閱讀 33,756評(píng)論 4贊 324
• ?日本核電站爆炸內(nèi)幕

  正文 年R本政府宣布当辐，位于F島的核電站，受9級(jí)特大地震影響鲤看，放射性物質(zhì)發(fā)生泄漏缘揪。R本人自食惡果不足惜，卻給世界環(huán)境...

  茶點(diǎn)故事閱讀 39,324評(píng)論 3贊 307
• 男人毒藥：我在死后第九天來索命

  文/蒙蒙 一义桂、第九天 我趴在偏房一處隱蔽的房頂上張望找筝。 院中可真熱鬧，春花似錦慷吊、人聲如沸袖裕。這莊子的主人今日做“春日...

  開封第一講書人閱讀 30,315評(píng)論 0贊 19
• 一樁弒父案溉瓶，背后竟有這般陰謀

  文/蒼蘭香墨 我抬頭看了看天上的太陽急鳄。三九已至，卻和暖如春堰酿，著一層夾襖步出監(jiān)牢的瞬間疾宏，已是汗流浹背。 一陣腳步聲響...

  開封第一講書人閱讀 31,539評(píng)論 1贊 262
• 情欲美人皮

  我被黑心中介騙來泰國打工触创， 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留坎藐，地道東北人。 一個(gè)月前我還...

  沈念sama閱讀 45,578評(píng)論 2贊 355
• 代替公主和親

  正文 我出身青樓哼绑，卻偏偏與公主長得像岩馍，于是被迫代替她去往敵國和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子凌那，可洞房花燭夜當(dāng)晚...

  茶點(diǎn)故事閱讀 42,877評(píng)論 2贊 345

推薦閱讀更多精彩內(nèi)容

• 2?? 雙序列比對(duì)(2)：BLAST詳細(xì)操作:web版和linux版

  序列比對(duì)和序列特征分析總目錄 網(wǎng)址：https://blast.ncbi.nlm.nih.gov/Blast.cg...

  Y大寬閱讀 24,290評(píng)論 1贊 36
• 本地化BLAST+進(jìn)行目標(biāo)序列比對(duì)

  舉個(gè)栗子??： 最近對(duì)擬南芥中的一個(gè)叫做CBF的基因非常感興趣兼雄，現(xiàn)在我想看看這個(gè)CBF基因在一個(gè)新測(cè)序的物種Am基因...

  Ruta閱讀 13,694評(píng)論 1贊 19
• linux服務(wù)器上的blast比對(duì)（my version：ncbi-blast-2.8.1+）

  一：安裝blast 1. 下載blast: $wget ftp://ftp.ncbi.nlm.nih.gov/bl...

  多啦A夢(mèng)的時(shí)光機(jī)_648d閱讀 8,004評(píng)論 1贊 23
• linux BLAST序列比對(duì) （nt/nr庫）

  安裝本地blast序列比對(duì)軟件，我們可以搜索一個(gè)查詢序列定制數(shù)據(jù)庫帽蝶，例如想研究一個(gè)新測(cè)序的基因組赦肋，或者感興趣的一組...

  像鳥一樣飛過你的高山閱讀 12,313評(píng)論 3贊 18
• 2021-03-04 blast+ 比對(duì)兩個(gè)denovo RNAseq序列

  問題描述：遇到一位同學(xué)他的SEQ數(shù)據(jù)是13年的块攒，illumina雙端測(cè)序，而且是無參物種佃乘，分別是病毒和昆蟲拼接好的...

  hopi_bai閱讀 938評(píng)論 0贊 1