Hi,
My HiSeq data for small RNA looks like this
@HWI-ST705:254:C0G8HACXX:5:1101:1674:1995 1:N:0:AGTCAA
TGAGATGAAGCACTGTAGCTCTGGAATTCTCGGGT
+
CCCFFFFFHHHHHJJJJJJJIJJJJJJJJJJJJJH
@HWI-ST705:254:C0G8HACXX:5:1101:1765:1986 1:N:0:AGTCAA
TGAGAACTGAATTCCATAGGCTGTTGGAATTCTCG
+
BCCFFFFDHFHHHIJJJIJJJJJHHJJEHIJIJJJ
@HWI-ST705:254:C0G8HACXX:5:1101:1785:1990 1:N:0:AGTCAA
GCTCTGTGATGAACCCTGGAATTCTCGGGTGCCAA
+
=?@DFFFA=AADFHIJJJJGAHHGIIIJI?CFGHC
@HWI-ST705:254:C0G8HACXX:5:1101:1825:1999 1:N:0:AGTCAA
TTTGGCAATGGTAGAACTCACACCTGGAATTCTCG
+
CCCFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJJJ
@HWI-ST705:254:C0G8HACXX:5:1101:2182:1985 1:N:0:AGTCAA
CAACNGAATCCCAAAAGCAGCTGTGGAATTCTCGG
+
@@@D#2=BDDHFHBGHHIIIGHHHIGGGHH<?DHI
@HWI-ST705:254:C0G8HACXX:5:1101:2106:1988 1:N:0:AGTCAA
TAGCTTATCAGACTGATGTTGACTTGGAATTCTCG
+
??@FFFD+=CFFFHGIJJGIHHHHHJCFHEHHHDH
@HWI-ST705:254:C0G8HACXX:5:1101:2543:1995 1:N:0:AGTCAA
TTCACAGTGGCTAAGTTCTGCTGGAATTCTCGGGT
+
CCCFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJJH
When I use the Adapter_trim.pl script from miRTools with format option "3" (for illumina format 1.3+) or evn "2" for the older formats ..I get a empty output file..
The previous illumina datasets had the complete IDs repeated befire the quality value lines and the script used to work good for me...
Any suggestions?
regards,
Nandan
The read grouping script from miRAnalyser also gives me issues for HiSeq dataset (again this worked well for GAII illumina data)
My HiSeq data for small RNA looks like this
@HWI-ST705:254:C0G8HACXX:5:1101:1674:1995 1:N:0:AGTCAA
TGAGATGAAGCACTGTAGCTCTGGAATTCTCGGGT
+
CCCFFFFFHHHHHJJJJJJJIJJJJJJJJJJJJJH
@HWI-ST705:254:C0G8HACXX:5:1101:1765:1986 1:N:0:AGTCAA
TGAGAACTGAATTCCATAGGCTGTTGGAATTCTCG
+
BCCFFFFDHFHHHIJJJIJJJJJHHJJEHIJIJJJ
@HWI-ST705:254:C0G8HACXX:5:1101:1785:1990 1:N:0:AGTCAA
GCTCTGTGATGAACCCTGGAATTCTCGGGTGCCAA
+
=?@DFFFA=AADFHIJJJJGAHHGIIIJI?CFGHC
@HWI-ST705:254:C0G8HACXX:5:1101:1825:1999 1:N:0:AGTCAA
TTTGGCAATGGTAGAACTCACACCTGGAATTCTCG
+
CCCFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJJJ
@HWI-ST705:254:C0G8HACXX:5:1101:2182:1985 1:N:0:AGTCAA
CAACNGAATCCCAAAAGCAGCTGTGGAATTCTCGG
+
@@@D#2=BDDHFHBGHHIIIGHHHIGGGHH<?DHI
@HWI-ST705:254:C0G8HACXX:5:1101:2106:1988 1:N:0:AGTCAA
TAGCTTATCAGACTGATGTTGACTTGGAATTCTCG
+
??@FFFD+=CFFFHGIJJGIHHHHHJCFHEHHHDH
@HWI-ST705:254:C0G8HACXX:5:1101:2543:1995 1:N:0:AGTCAA
TTCACAGTGGCTAAGTTCTGCTGGAATTCTCGGGT
+
CCCFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJJH
When I use the Adapter_trim.pl script from miRTools with format option "3" (for illumina format 1.3+) or evn "2" for the older formats ..I get a empty output file..
The previous illumina datasets had the complete IDs repeated befire the quality value lines and the script used to work good for me...
Any suggestions?
regards,
Nandan
The read grouping script from miRAnalyser also gives me issues for HiSeq dataset (again this worked well for GAII illumina data)
Comment