Showing posts with label Research. Show all posts
Showing posts with label Research. Show all posts

Wednesday, 30 July 2014

How to parse MiRProf Results.

 MiRProf Results contains may unused information which might affect the effective analysis.

Here shows how to parse the MiRProf exported .csv files.

First get all the names of micro RNA from the file:

grep - All-Expression-hg19-MirBase.csv > mir.txt

 The get the corresponding Total number of reads:

grep Total All-Expression-hg19-MirBase.csv > Allnumber.txt

Then combine the two files on excel, then you are done.

Convert all fastq files in folder to fasta

for fq in *; do fn=$(echo $fq|sed s'/.fq//'); echo $fn;awk 'NR % 4 == 1 {print ">" $0 } NR % 4 == 2 {print $0}' $fq > ../MergedFasta/$fn.fa; done;

Merge two Replicates into one

for fq in 890/*; do fn=$( echo $fq | sed 's/_R1_001.fastq//'| sed 's/_R2_001.fastq//'); echo $fn; cat $fq >>$fn.fq; done;