G***G
发帖数: 16778 | 1 what kind of data are they?
is it possible to extract mitochondra sequence from this kind of data?
@HWI-1KL138:2:1101:1092:1984#GCCAATGCCAAT/1
NGGGCACAAAGAATTACAGCAAACTAAAATCTCCCTAGTTCAAGTATGGAAGACAATGACCATTTATCTTTGGCAG
GATTCACCCTTTTGCAAGTNGNNGA
+HWI-1KL138:2:1101:1092:1984#GCCAATGCCAAT/1
BOWZ^^^
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBB
@HWI-1KL138:2:1101:1138:1988#GCCAATGCCAAT/1
NTAAATTTAATGTGACCGTTTATCGCAATCTGCCGACCACTCGCGATTCAATCATGACTTCGTGATAAAAGATTGA
GTGTGAGGTTATAACGCCGAAGCGG
+HWI-1KL138:2:1101:1138:1988#GCCAATGCCAAT/1
BP\a`ccceceaY`deehcJ`_K`bfe`P^[IP^H^e^cZIH[FFZZZV\`d]d`Y_ZZHZTU^_`ca`_`_
aYYTYTY`]``^GJY`cd_WETO[_aaR] | l**********1
发帖数: 5204 | 2 Samtool file
flag inside sam file
SEQanswers - Bioinformatics — I got a read aligned as below: ~ HWI-1KL138:2
:2105:12847:125331#GCCAATGCCAAT 161 chr1 12036 1 74M = 12645 1188
CTGTGCCAGGGTGCAAGCTGAGCACTGGAGTGGAGTTTTCCTGTGGAGAGGAGCCATGCCTAGAGTGGGATGGG
hhhhhhhhhhhghhhhhhhhhhhhhhhhhheffffdffffhhhghhehefhhhhhghhfhfffWdbfffbffff
NM:i:0 NH:i:3 CC:Z:chr15 CP:i:102519061 HI:i:0 ~ The flag here is 161=128+32
+1. 128 means 2nd pair; 32 means mated reverse strand; 1 means paired read.
I am wondering why there is no information about the alignment in this flag.
The alignment info can be either both aligned, read unmapped, mate unmapped
indicated by 0x2, 0x4, 0x8 in sam files. Thanks advance for any answer or
comment. ~ fanping
//www.buzzbox.com/ws?u=seqanswers.com&id=85531 |
|
