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	96.0%	63.0 / ‑3.3	26	R209* (CGA→TGA)	dck	deoxynucleoside kinase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (7/17);  total (8/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

TCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTA‑ATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCG  >  CP000730/607766‑607993                                                                                                                       |                                                                                                               tCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGa                                                                                          <  2:383970/141‑1 (MQ=255)         aTTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGa                                                                                  <  2:156671/141‑1 (MQ=255)           ttaTAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGAtt                                                                                <  2:453648/139‑1 (MQ=255)            gATCAAATAATGCATAATATCTATGGAATATACATTGAATATACACATGTACGTTGTTGAAAGAATCTGGTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCTAGCAATTTTGTTTATCATAGGATTTAAAGATTc                                                                               >  2:226080/1‑141 (MQ=255)            gATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTGTTTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTc                                                                               <  2:27879/141‑1 (MQ=255)            gATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTc                                                                               <  2:2828/141‑1 (MQ=255)            gATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTc                                                                               <  2:130115/141‑1 (MQ=255)               aaaaTAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATAGGTTTGAATAATTCAAGCAATTTTATTTATCATAGGATTTAAAGATTCtgg                                                                            >  2:90369/1‑141 (MQ=255)                   taatCGATAATATCTATGGAATATACATTGAATATAAACATATAGGCTCTGTCTATAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTcc                                                                        >  2:116611/1‑141 (MQ=255)                      tCGATAATATCTATGGAATTTACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGGAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTa                                                                     <  2:145283/141‑1 (MQ=255)                      tCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGGCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTa                                                                     >  2:493446/1‑141 (MQ=255)                        gATAATATCTATGGAATATACATTGAATATAAACATATACGCTATTTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATg                                                                   <  2:71852/141‑1 (MQ=255)                          taataTCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGa                                                                 <  2:170377/141‑1 (MQ=255)                             tatCTATGGAATATACATTGAATATAAACAGATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGatat                                                              <  2:249111/141‑1 (MQ=255)                              atCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATc                                                             >  2:362558/1‑141 (MQ=255)                               aatATGTAATATACATTAAAGATAAACATATACGCTATGGCTTTCGGCTTTTATCGTGTATCTACTTGTGGATATGTTTTAATAATTCAAGCAATTTTGTTCCTCACAGGATTTAAACATTCGGGCTCCTTATGGATATCa                                                            <  2:147248/139‑1 (MQ=255)                                               ttaatagtaaCATTGACGCTATGGCTTTAGTCTTTGATCGTGTATCTAATTGTAGATATGTTTGAATAATTCAAGCAATTTTGTCTATCATAGCACTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTA‑ATATTGa                                           <  2:187650/133‑1 (MQ=255)                                                gAATATAAACATATACGCTATGGCTTTAGTCTTTTATCGTGTGTCCACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTA‑ATATTGAt                                          <  1:5210/141‑1 (MQ=255)                                                         cATATACGCTATGTCTTTAGTCTTTTATCGTGTATTTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAGATATT‑TTACGTACAAc                                 <  2:229963/141‑1 (MQ=255)                                                                     gTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTA‑ATATTGATACGTACAACTGGACATGCAtt                     <  2:196711/141‑1 (MQ=255)                                                                      tCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATGTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTG‑ATATTGATACGTACAACTGGACATGCATTa                    >  2:28491/1‑141 (MQ=255)                                                                       cTTTAGTCTTTTATCGTGTATCTGCTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTA‑ATATTGATACGTACAACTGGACCTGCATTaa                   >  2:65838/1‑141 (MQ=255)                                                                        tttAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTA‑ATATTGATACGTACAACTGGACATGCATTaaa                  >  2:15684/1‑141 (MQ=255)                                                                            gttttttATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTA‑ATATTGATACGTACAACTGGACATGCATTAAAGCTa              <  2:91651/138‑1 (MQ=255)                                                                              cTTTTATCGTGTATGTACTTGTCGATATTGTTGACCCATTCAAGTAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTA‑ATATTGATACGTACAACTGGACATGCATTAAAGCTAtt            <  2:102718/141‑1 (MQ=255)                                                                                        gtATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTA‑ATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCg  <  2:263210/141‑1 (MQ=255)                                                                                                                       |                                                                                                               TCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTA‑ATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCG  >  CP000730/607766‑607993

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