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% 78.1 / NA 28R209* (CGA→TGA) USA300HOU_0550deoxynucleoside kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/16);  total (12/16)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCC  >  CP000730/607752‑608018
                                                                                                                                    |                                                                                                                                      
gtgAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTCTGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCACGCAAttt                                                                                                                                >  2:100371/1‑141 (MQ=255)
          ttATTCGATAGTATTGAGAAAAAAATCGCTAATTTCTATGGAACAAACATGGAATATAAGCATATACGCTACGTCTTTAGTCTTGTATCGTGTATCTACATGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCAt                                                                                                                      <  2:262138/140‑1 (MQ=255)
          ttATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTGTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCAt                                                                                                                      <  1:288117/140‑1 (MQ=255)
                    gAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTaaa                                                                                                            >  2:431173/1‑141 (MQ=255)
                                     gATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTCTCTACTTGTCGATATGTTTGAATAATTCAAGCAATTATGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATg                                                                                           <  2:344570/141‑1 (MQ=255)
                                             ttATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCat                                                                                   <  1:301719/140‑1 (MQ=255)
                                                  gAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTct                                                                              >  2:151442/1‑140 (MQ=255)
                                                   aaTATACATTGAATATAAACATATACGCTATGTCGTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTGCTTATGGTTATCATATTGAt                                                                             >  1:180010/1‑141 (MQ=255)
                                                    atatACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCAtt                                                                            <  1:266207/141‑1 (MQ=255)
                                                        aCATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTAGGGATATCATAttagttagta                                                                        >  1:76531/1‑133 (MQ=255)
                                                             gAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGAt                                                                   <  1:431173/141‑1 (MQ=255)
                                                               atatAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATAc                                                                 <  1:23901/141‑1 (MQ=255)
                                                                   aaaCATATACGCGATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGAAATTTTGTTTGTCATAGGATTTAAAGATTCGGGGTCCTTATGGATATAATATTCTTTAATATTGTTACGTAc                                                             >  1:155370/1‑141 (MQ=255)
                                                                      cATATACGCTATGTCTTTAGTCTTTTATCTTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAAc                                                          <  2:372281/141‑1 (MQ=255)
                                                                       atatACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTAAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATAAGTACAACt                                                         >  1:297345/1‑141 (MQ=255)
                                                                        tataCGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTg                                                        >  1:312471/1‑141 (MQ=255)
                                                                             gCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACAt                                                   <  2:190012/141‑1 (MQ=255)
                                                                                 tGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACCACTGGACATGCAt                                               >  2:268857/1‑141 (MQ=255)
                                                                                 tGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCAt                                               <  2:157832/141‑1 (MQ=255)
                                                                                     tttAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTaaa                                           >  1:329967/1‑141 (MQ=255)
                                                                                      ttAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAg                                          >  2:447495/1‑141 (MQ=255)
                                                                                           tttttATCTTGTATCTACGTGTCGATATGTTTTAATTAGTCAAGCTATTTTGTTTATCATAGGATTTAAAGATGCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCAGTAAAGCTAtt                                     <  2:263189/140‑1 (MQ=255)
                                                                                              ttATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAAt                                  >  1:453071/1‑141 (MQ=255)
                                                                                                            cTTGGCGATATGTTTGAATAATTCAATCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGc                    <  1:280292/141‑1 (MQ=255)
                                                                                                             atGTCGATATATTTGAATAATTCAAGCAATTTTGTTTAACATAGGAATTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATGAAAGCTATTAATCCAATCGCCATAGCg                   <  2:308345/140‑1 (MQ=255)
                                                                                                                tCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGttt                <  2:229700/141‑1 (MQ=255)
                                                                                                                             ataatTCAAGCAATTTCTTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTc   <  1:79433/141‑1 (MQ=255)
                                                                                                                              taatTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTcc  <  1:100371/141‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                      
GTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCC  >  CP000730/607752‑608018

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: