NGS: How to use Stacks to remove low quality reads from RNA-seq data

Stacks is a software pipeline for building loci from short-read sequences, such as those generated on the Illumina platform. 




bedtools bamtofastq is a conversion utility for extracting FASTQ records from sequence alignments in BAM format.

Comments

Post a Comment

Popular posts from this blog

gspread error:gspread.exceptions.SpreadsheetNotFound

Error message: Traceback (most recent call last):   File "read_google_sheet.py", line 20, in <module>     wks = gc.open_by_key("1jUJOlICOH6VnyuALdU2VAc_Ci_ZVwdW-xc5kr_qkJHY")   File "/usr/lib/python2.7/site-packages/gspread/client.py", line 113, in open_by_key     raise SpreadsheetNotFound gspread.exceptions.SpreadsheetNotFound Solved : In the json file, there is a "client_email". You need to share your file with this email. 
Read more

Miniconda installation problem: concurrent.futures.process.BrokenProcessPool: A process in the process pool was terminated abruptly while the future was running or pending.

We got error message shown as below: [/home/omicsacademy/miniconda3] >>> PREFIX=/home/omicsacademy/miniconda3 Unpacking payload ... concurrent.futures.process._RemoteTraceback: ''' Traceback (most recent call last):   File "concurrent/futures/process.py", line 368, in _queue_management_worker   File "multiprocessing/connection.py", line 251, in recv TypeError: __init__() missing 1 required positional argument: 'msg' ''' The above exception was the direct cause of the following exception: Traceback (most recent call last):   File "entry_point.py", line 69, in <module>   File "concurrent/futures/process.py", line 484, in _chain_from_iterable_of_lists   File "concurrent/futures/_base.py", line 611, in result_iterator   File "concurrent/futures/_base.py", line 439, in result   File "concurrent/futures/_base.py", line 388, in __get_result concurrent.futures.process.BrokenProc...
Read more

Extract note from google slide using google-slide-api in Python

    This article is one after the previous post ( https://omicsacademy.blogspot.com/2021/03/save-all-google-presentation-slides-as.html ).      In this post, you are going to see how we can extract the note information from Google slide using  google-slide-api in Python. The first part of the code is the same as the previous post .  import urllib.request import json import os.path from googleapiclient.discovery import build from google_auth_oauthlib.flow import InstalledAppFlow from google.auth.transport.requests import Request from google.oauth2.credentials import Credentials # If modifying these scopes, delete the file token.json. SCOPES = ['https://www.googleapis.com/auth/presentations.readonly'] # The ID of a sample presentation. PRESENTATION_ID = '1-aTBNXcSIqlMRzn-FHnRmRPbGlh5eY8MgZNaBwo15IM' creds = None # The file token.json stores the user's access and refresh tokens, and is # created automatically when the authorization flow completes...
Read more