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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE607,8840GA98.3% 130.6 / ‑5.9 41R209* (CGA→TGA) dckDeoxyadenosine/deoxycytidine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/19);  minor base C (1/0);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTT  >  USA300TCH1516_ALE/607750‑608013
                                                                                                                                      |                                                                                                                                 
atGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAAt                                                                                                                             >  1:480955/1‑141 (MQ=255)
atGTGAATTAATGTACTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAAt                                                                                                                             <  1:326098/141‑1 (MQ=255)
    gAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTg                                                                                                                         <  1:239046/141‑1 (MQ=255)
           tGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCa                                                                                                                  >  1:608359/1‑141 (MQ=255)
           tGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCa                                                                                                                  >  1:402895/1‑141 (MQ=255)
            gTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTCTCGTGTATCTACTTTTCGATATGTTTGAATAATTCCAGCATTTTTGTTTATCAt                                                                                                                 >  1:364404/1‑141 (MQ=255)
            gTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCAt                                                                                                                 >  1:349740/1‑141 (MQ=255)
                  gATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTTTCTCCTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGAtt                                                                                                           >  1:282945/1‑141 (MQ=255)
                    tAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGGTTTa                                                                                                         >  1:498114/1‑141 (MQ=255)
                    tAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATCTTGTTTATCATAGGATTTa                                                                                                         <  1:477024/141‑1 (MQ=255)
                     aGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTaa                                                                                                        >  1:835250/1‑141 (MQ=255)
                     aGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTaa                                                                                                        <  1:838965/141‑1 (MQ=255)
                     aGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATCTTGTTTATCATAGGATTTaa                                                                                                        >  1:815861/1‑141 (MQ=255)
                        aTTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGa                                                                                                     >  1:806257/1‑141 (MQ=255)
                           gATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTc                                                                                                  >  1:655713/1‑141 (MQ=255)
                             aaaaaTAATCGATACTATCTCTTGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCgg                                                                                                <  1:942/140‑1 (MQ=255)
                                aataacCAATAATATCTAAGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGt                                                                                             <  1:95670/141‑1 (MQ=255)
                                    atCGATAATATCTATGGAATATTGATTTAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCtt                                                                                         >  1:208893/1‑141 (MQ=255)
                                        ataataTCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATgg                                                                                     >  1:309551/1‑141 (MQ=255)
                                               cTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCat                                                                              <  1:365034/141‑1 (MQ=255)
                                                tATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCata                                                                             >  1:696943/1‑141 (MQ=255)
                                                                      aaCATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTAcc                                                       >  1:95429/1‑140 (MQ=255)
                                                                      aaCATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACa                                                       >  1:507157/1‑141 (MQ=255)
                                                                         atatACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACt                                                    <  1:675106/141‑1 (MQ=255)
                                                                                cTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATg                                             <  1:468672/141‑1 (MQ=255)
                                                                                    gTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGAGTTAAAGATTCGGCGTCCTTATGGATATCATATTCATTAACATTGATACGTCCAACTGGACATGCAtt                                         >  1:197701/1‑141 (MQ=255)
                                                                                     tCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCCTATGGATATCATATTCCTTAATATTGATACGTACAACTGGACATGCATTa                                        <  1:576854/141‑1 (MQ=255)
                                                                                       tttAGTCTTTTATCGGGTATCTACTTGTCGATATCTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGCGGTCCATATGGATGTCATATTCATTAATATTGATACGTACAACTGGACATGCATTaaa                                      <  1:528518/141‑1 (MQ=255)
                                                                                             cTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTAtt                                >  1:572086/1‑141 (MQ=255)
                                                                                              ttttATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTa                               >  1:395105/1‑141 (MQ=255)
                                                                                                     gtgtATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAAt                        >  1:627584/1‑141 (MQ=255)
                                                                                                     gtgtATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAAt                        >  1:729907/1‑141 (MQ=255)
                                                                                                     gtgtATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAAt                        >  1:410246/1‑141 (MQ=255)
                                                                                                        tATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGt                     <  1:824702/141‑1 (MQ=255)
                                                                                                          tCTACTTGGCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCa                   <  1:470019/141‑1 (MQ=255)
                                                                                                              cTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTCTATCATAGGATTTAAAGCTTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAATTGGACATGCATTAAAGCTCTTTATCCAATCGTCATAGc               <  1:486226/141‑1 (MQ=255)
                                                                                                              cTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTATATCATAGGATTTAAAGTTTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAAATGGACATGCATTAAAGCTTTTCATCCAATCGTCATAGc               <  1:486343/141‑1 (MQ=255)
                                                                                                                   cGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTa          <  1:470174/141‑1 (MQ=255)
                                                                                                                      tatGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAAt       <  1:769155/141‑1 (MQ=255)
                                                                                                                      tatGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAAt       <  1:361709/141‑1 (MQ=255)
                                                                                                                           ttGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCtt  <  1:441997/141‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                 
ATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTT  >  USA300TCH1516_ALE/607750‑608013

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