備忘錄
投影片放映
大綱
1
Genome Resequencing analysis
2
outline
  • Download & import data
  • Mapping reads to reference genome
  • SNP detect
  • DIP (InDel) detect


3
Rsequence sample data
  • Download data from http://163.25.92.61/course/454.zip



  • Extract the file


4
3 files are extracted from 454.zip
    • Ecoli.FLX.fna  (Reads sequence in fasta format)
    • Ecoli.FLX.qual  (Reads quelity in fasta format)
    • NC_010473.gbk (E. coli str. K-12 substr. DH10B, complete genome   sequence in Genbank format)

5
How many reads in a fasta file?
  • Extract lines with “>” character



  • And count it
6
Import 454 reads
7
Select 2 files via  Ctrl  key
8
 
9
Import Reference genome
10
 
11
Mapping reads to Reference
12
Select reads
13
Select NC_010473 as reference
14
 
15
Mapping result
16
Mapping report
17
SNP detection
18
 
19
filter of quality
20
Annotate SNP on reference/consensus sequence
21
show SNP table
22
show genome/reads mapping
23
 
24
 
25
DIP detection
26
 
27
 
28