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	82.8%	66.9 / 4.7	29	R209* (CGA→TGA)	dck	deoxynucleoside kinase
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (16/8);  total (19/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCG‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAG  >  CP000730/607747‑608020                                                                                                                                          |                                                                                                                                           tataTGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGc                                                                                                                                       <  2:288608/141‑1 (MQ=255)   tatGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCaa                                                                                                                                     >  2:244536/1‑141 (MQ=255)    atGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAAt                                                                                                                                    >  1:294793/1‑141 (MQ=255)               tGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCa                                                                                                                         >  1:149661/1‑141 (MQ=255)                         aGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTaa                                                                                                               <  1:122271/141‑1 (MQ=255)                                aTAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCG‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCg                                                                                                        >  2:385334/1‑141 (MQ=255)                                 tAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCgg                                                                                                       >  1:61674/1‑141 (MQ=255)                                     ataatCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTc                                                                                                   >  2:283277/1‑141 (MQ=255)                                       aatCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCt                                                                                                 <  2:243637/141‑1 (MQ=255)                                          cGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTAt                                                                                              >  2:376322/1‑141 (MQ=255)                                                  tCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCa                                                                                      >  2:154570/1‑141 (MQ=255)                                                    tATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCata                                                                                    <  1:413374/141‑1 (MQ=255)                                                             taCATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCG‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAat                                                                           <  1:276316/141‑1 (MQ=255)                                                              aCATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCG‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAata                                                                          >  1:44909/1‑141 (MQ=255)                                                               cATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAatat                                                                         <  1:417913/141‑1 (MQ=255)                                                                     atatAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATAc                                                                   >  1:42052/1‑141 (MQ=255)                                                                         aaaCATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTAc                                                               >  1:52297/1‑141 (MQ=255)                                                                          aaCATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCG‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACa                                                              >  2:148654/1‑141 (MQ=255)                                                                                           tttAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTaaa                                             >  1:174364/1‑141 (MQ=255)                                                                                                    ttATCGTGTATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAAt                                    >  1:300366/1‑141 (MQ=255)                                                                                                         gtgtATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAAt                               >  2:411575/1‑141 (MQ=255)                                                                                                         gtgtATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAAt                               >  1:232422/1‑141 (MQ=255)                                                                                                            tATCTACTTGTCGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGt                            >  2:384641/1‑141 (MQ=255)                                                                                                               cTACTTGTCGATATGTCT‑‑ACAATTGGCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCAt                         <  1:107261/141‑1 (MQ=255)                                                                                                                       cGATATGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTa                 <  2:372679/141‑1 (MQ=255)                                                                                                                           atGTTTGAATAATTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATa             >  1:204832/1‑141 (MQ=255)                                                                                                                             gTTTGAATAATTCA‑‑AGCAATTTTGTTTATGATAAGTTATAAATATTCTGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGc           <  1:184209/141‑1 (MQ=255)                                                                                                                                     aatTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCa   <  2:61674/141‑1 (MQ=255)                                                                                                                                      atTCA‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAg  >  2:19960/1‑141 (MQ=255)                                                                                                                                          |                                                                                                                                           TATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCG‑‑AGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAG  >  CP000730/607747‑608020

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