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.
Showing posts with label Research. Show all posts
Showing posts with label Research. Show all posts
Wednesday, 30 July 2014
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;
Subscribe to:
Posts (Atom)