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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	607,884	0	G	A	97.9%	167.1 / NA	47	R209* (CGA→TGA)	dck	Deoxyadenosine/deoxycytidine kinase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (23/23);  total (24/23)
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

CATAGTGTATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAGT  >  USA300TCH1516_ALE/607740‑608021                                                                                                                                                 |                                                                                                                                          cATAGTGTATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGc                                                                                                                                        <  1:623583/148‑1 (MQ=255)      tgtATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCACttt                                                                                                                                   >  1:695065/1‑148 (MQ=255)      tgtATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAAttt                                                                                                                                   >  1:966187/1‑148 (MQ=255)      tgtATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAAttt                                                                                                                                   >  1:103660/1‑148 (MQ=255)      tgtATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAAttt                                                                                                                                   >  1:205294/1‑148 (MQ=255)       gtATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAAtttt                                                                                                                                  >  1:772865/1‑148 (MQ=255)        tataTGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTg                                                                                                                                 <  1:114641/148‑1 (MQ=255)             gttaattaattTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTAt                                                                                                                            <  1:1261018/145‑1 (MQ=255)             gtgAATTAATTTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTAt                                                                                                                            <  1:107188/148‑1 (MQ=255)                     aTGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCACTTTTGTTTATCATAGGAt                                                                                                                    >  1:599747/1‑148 (MQ=255)                        tactccaTAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTa                                                                                                                 <  1:256792/142‑1 (MQ=255)                        tATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTa                                                                                                                 <  1:430181/148‑1 (MQ=255)                        tATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTa                                                                                                                 <  1:430154/148‑1 (MQ=255)                        tATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTa                                                                                                                 <  1:12824/148‑1 (MQ=255)                        caTTCGATAGAATTGAAAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTa                                                                                                                 <  1:37757/147‑1 (MQ=255)                             gATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGAt                                                                                                            <  1:432465/148‑1 (MQ=255)                                   agtGATAAACTAATCGATAATATCTATGGAATATACAGTGAATATAAACATATACGCTATGTCTGTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCgggg                                                                                                      <  1:56686/146‑1 (MQ=255)                                   aTTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCgggg                                                                                                      <  1:673705/148‑1 (MQ=255)                                   aTTGACAACATAATATATTATCTCTATGGAATATACATTGAATATAAATATATACGCTATGTTTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCgggg                                                                                                      <  1:432686/148‑1 (MQ=255)                                   aTTGACAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGTTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGCTATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCgggg                                                                                                      <  1:170793/148‑1 (MQ=255)                                              aatCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGat                                                                                           <  1:446811/148‑1 (MQ=255)                                                   ataataTCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTGTCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCat                                                                                      >  1:173334/1‑148 (MQ=255)                                                       tatCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTc                                                                                  >  1:139117/1‑148 (MQ=255)                                                           tATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATCGCGCGAGCGATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTa                                                                              >  1:1238253/1‑148 (MQ=255)                                                           tATGGAATATACATTGAATATAAACACATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTa                                                                              >  1:938413/1‑148 (MQ=255)                                                               gAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAatat                                                                          >  1:258469/1‑148 (MQ=255)                                                               gAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAatat                                                                          <  1:1112404/148‑1 (MQ=255)                                                                            atatAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAAc                                                             >  1:1217262/1‑148 (MQ=255)                                                                            atatAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGATTCTAATGGATATCATATTCATTAATATTGATACGTACAAc                                                             >  1:969154/1‑148 (MQ=255)                                                                                  aCATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACa                                                       <  1:1104702/148‑1 (MQ=255)                                                                                      ataCGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCa                                                   <  1:503745/148‑1 (MQ=255)                                                                                       taCGCTATGTCTTTAGTCTTTTAACGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCAt                                                  >  1:936614/1‑148 (MQ=255)                                                                                           cTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTaaa                                              >  1:1072254/1‑148 (MQ=255)                                                                                           cTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTaaa                                              >  1:691419/1‑148 (MQ=255)                                                                                           cTATGTCTTTAGTCTTTTATCATGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTaaa                                              >  1:693507/1‑148 (MQ=255)                                                                                               gTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTa                                          <  1:821523/148‑1 (MQ=255)                                                                                                tCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTAt                                         >  1:270646/1‑148 (MQ=255)                                                                                                 cTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTGGGGTAGCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTAtt                                        >  1:979613/1‑148 (MQ=255)                                                                                                   ttAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTaa                                      <  1:830727/148‑1 (MQ=255)                                                                                                   ttAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTaa                                      <  1:1189590/148‑1 (MQ=255)                                                                                                        tttttATCGTGTATCTACTTGCCGATATGTTTGAATAATTCAAGCACTTTTGTTTACCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCaa                                 <  1:645068/147‑1 (MQ=255)                                                                                                        cTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCaa                                 >  1:949529/1‑148 (MQ=255)                                                                                                         ttttATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTGCAACTGGACATGCATTAAAGCTATTAATCCAAt                                >  1:1128769/1‑148 (MQ=255)                                                                                                                     tCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGtt                    <  1:199546/148‑1 (MQ=255)                                                                                                                        aCTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTaa                 >  1:351538/1‑148 (MQ=255)                                                                                                                        aCTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCACTAAAGCTATTAATCCAATCGTCATAGCGTTTaa                 >  1:852923/1‑148 (MQ=255)                                                                                                                                       tGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAGt  >  1:432022/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                          CATAGTGTATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAGT  >  USA300TCH1516_ALE/607740‑608021

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: ‑