Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	607,884	G→A	R209* (CGA→TGA)	dck ←	deoxynucleoside kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	607,884	0	G	A	100.0%	52.8 / NA	19	R209* (CGA→TGA)	dck	deoxynucleoside kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/8);  total (11/8)

TGTATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTT  >  CP000730/607745‑608015                                                                                                                                            |                                                                                                                                    tgtATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTTTCTACTTGTCGATATGTTTGAATAATTCAAGCAAttt                                                                                                                             >  1:121186/1‑148 (MQ=255) tgtATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTAGTCGATATGTTTGAATAATTCAAGCAAttt                                                                                                                             >  1:295747/1‑148 (MQ=255)  gtATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAAtttt                                                                                                                            <  2:161041/148‑1 (MQ=255)        gtgAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATCATTCAAGCAATTTTGTTTAt                                                                                                                      <  2:63616/148‑1 (MQ=255)          gAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCa                                                                                                                    >  1:44187/1‑148 (MQ=255)            aTTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATa                                                                                                                  >  1:102703/1‑148 (MQ=255)                        gATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCGCTGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCCATTTTGTTTATCATAGGATTTAAAGAt                                                                                                      >  2:192210/1‑148 (MQ=255)                                  aTAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCt                                                                                            <  1:145315/148‑1 (MQ=255)                                                aataTCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCCTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCatat                                                                              <  2:216184/148‑1 (MQ=255)                                                     cTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCAtt                                                                         >  1:180381/1‑148 (MQ=255)                                                        tGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAat                                                                      <  1:30906/148‑1 (MQ=255)                                                                aCATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATAc                                                              <  2:215312/148‑1 (MQ=255)                                                                                    cgctcggTCTTTAGTCTTTTATCGTGTATCTGCTTGTCGATCTGTTTGAATAATTCAAGCGATTTTGTTTATCATAGGATTTAAAGATACGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTAGAATAGCATTa                                          >  2:183825/7‑148 (MQ=255)                                                                                    cGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTa                                          >  2:291462/1‑148 (MQ=255)                                                                                              ttAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTaa                                <  2:201841/148‑1 (MQ=255)                                                                                                  tCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCa                            >  1:4999/1‑148 (MQ=255)                                                                                                      ttATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCg                        >  2:226852/1‑148 (MQ=255)                                                                                                                           atatGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCttt   <  2:44187/148‑1 (MQ=255)                                                                                                                            tatGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCtttt  >  2:74576/1‑148 (MQ=255)                                                                                                                                            |                                                                                                                                    TGTATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTT  >  CP000730/607745‑608015

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑