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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730607,8840GA100.0% 71.1 / NA 26R209* (CGA→TGA) USA300HOU_0550deoxynucleoside kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/16);  total (10/16)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAGT  >  CP000730/607747‑608021
                                                                                                                                         |                                                                                                                                         
tatttttgAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGc                                                                                                                                        <  1:29333/135‑1 (MQ=255)
               gTATTCGATAGAATTGATAAAATAATCGAGAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGGCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCAt                                                                                                                         <  2:361419/141‑1 (MQ=255)
                   tCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGa                                                                                                                     <  2:274340/141‑1 (MQ=255)
                                      aatCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCt                                                                                                  <  2:309897/141‑1 (MQ=255)
                                       atCGATAATATCTATGGAATATTCATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCtt                                                                                                 <  2:38425/141‑1 (MQ=255)
                                       atCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCtt                                                                                                 <  1:363898/141‑1 (MQ=255)
                                                  cTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCat                                                                                      >  2:219716/1‑141 (MQ=255)
                                                         atatACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGCTCTAATGATTCGGGGTCCTTATGGATATCCTATTCttt                                                                               >  1:99938/1‑141 (MQ=255)
                                                         atatACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCAtt                                                                               >  1:315719/1‑141 (MQ=255)
                                                              cATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATTGATATCCTATTCATTAatat                                                                          >  2:348728/1‑141 (MQ=255)
                                                                      ataAACATATACGCTATGTCTTTAGTCTTAGATCGTGTATCCACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGt                                                                  <  1:247599/141‑1 (MQ=255)
                                                                         aaCATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACa                                                               >  1:351154/1‑141 (MQ=255)
                                                                                 cGCTATGTCTTTAGTCTTTTATCGTTTATCTACTTGTCGATATGTTTGAATAATGCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGCATATCATATTCATTAATATTGCTACGTACAACTGGACa                                                       <  1:220850/141‑1 (MQ=255)
                                                                                  gCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACAt                                                      >  2:189751/1‑141 (MQ=255)
                                                                                   cTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACTACTGGACATg                                                     <  1:286912/141‑1 (MQ=255)
                                                                                                    tATCGTGTATCTACTTGTAGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGTACATGCATTAAAGCTATTAATc                                    <  2:289143/141‑1 (MQ=255)
                                                                                                     aTCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATcc                                   >  1:123226/1‑141 (MQ=255)
                                                                                                               tACTTGTCGTTATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATa                         <  1:116592/141‑1 (MQ=255)
                                                                                                                       gATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTaa                 <  2:226980/141‑1 (MQ=255)
                                                                                                                       gATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTaa                 <  2:227029/141‑1 (MQ=255)
                                                                                                                          atGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATa              >  1:351024/1‑141 (MQ=255)
                                                                                                                           tggctGAATAATTCAAGCAATTTTGTTTATCAGAGGATGTAAAGATTCGGGGTCCTCATGCATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAg             <  1:348728/137‑1 (MQ=255)
                                                                                                                           tGTTGGAATAATTCAAGTAATTTTGTTTATCATATGATTTAAAGATTCGGCGTCCTTGTGGATATCATATTCATTAATATTGTTACGTTCAACTGGACATGCATTAAAACTATTAATCAAATCGTCATAGCGTTTAAATAg             <  2:9021/141‑1 (MQ=255)
                                                                                                                              ttGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCtt          >  1:331066/1‑141 (MQ=255)
                                                                                                                                 aataatTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGTATATCATATTCATTAATATTGATACGTACAACTGGCCATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAACAGCttttt       <  2:291667/141‑1 (MQ=255)
                                                                                                                                      ttCAAGCAATTTTGTTTATAATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGAGGTTAAACAGCTTATTCTAga  >  2:284854/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                         
TATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAGT  >  CP000730/607747‑608021

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