Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,232,465 G→C 100% A124P (GCG→CCG)  USA300HOU_1157 → non‑specific serine/threonine protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,232,4650GC100.0% 67.0 / NA 23A124P (GCG→CCG) USA300HOU_1157non‑specific serine/threonine protein kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (12/11);  total (12/11)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAGAAGATGACTGTTACTACTTAGTAATGGAATATATCGAAGGTCCGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATGCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCTTTAAGTGAGACGTCTTTAACTCAGACTAATCAT  >  CP000730/1232328‑1232599
                                                                                                                                         |                                                                                                                                      
aagaagATTACTATTACTACTTAGTAATGGAATATATCGAAGGTCCGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATccgc                                                                                                                                     <  1:91069/141‑1 (MQ=255)
  gaagaTGACTGTTACTACTTAGTAATGGAATATATCGAAGGTCCGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCAt                                                                                                                                   <  2:12870/141‑1 (MQ=255)
         aCTGTTACTACTTAGTAATGGAATATATCGAAGGTCCGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAACGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGc                                                                                                                            >  2:52200/1‑141 (MQ=255)
            gTTACTACTTAGTAATGGAATATATCGAAGGTCCGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTa                                                                                                                         >  1:85766/1‑141 (MQ=255)
                                tatatCGAAGGTCCGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAg                                                                                                     >  1:241477/1‑141 (MQ=255)
                                  aatcgaattgcacaCTTTGTGTGAATATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGcc                                                                                                   <  1:239498/128‑1 (MQ=255)
                                   atCGAAGGTCCGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCTGCATGGTATGCGTCTTGTTCATAGAGATATTaacaca                                                                                                  >  2:187557/1‑137 (MQ=255)
                                        aGGGCCGAATTTGTATAAGTAAAATGAAAGACATGGGCCATTAAGTGTTCACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACaaaa                                                                                             <  1:156368/141‑1 (MQ=255)
                                                          gTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTATACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGAAAACaa                                                                           >  2:208168/1‑141 (MQ=255)
                                                          gTATATTGAAAGACATGGGACATTAAGTGTTGAAACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCaa                                                                           <  1:125536/141‑1 (MQ=255)
                                                                  aaaGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATCTTGCGGGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGt                                                                   <  1:275075/141‑1 (MQ=255)
                                                                      tCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGaa                                                               <  2:24713/141‑1 (MQ=255)
                                                                       cATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGaaa                                                              >  2:23204/1‑141 (MQ=255)
                                                                       cATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATCAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACTGCAATAAAACGtgcaaa                                                              >  2:187801/1‑136 (MQ=255)
                                                                                 tAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGCACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGAtt                                                    <  1:994/141‑1 (MQ=255)
                                                                                       ttGACACAGCGATTAATTTTACGAACAAAATATTGGATGCCAATAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGaa                                              <  1:258952/141‑1 (MQ=255)
                                                                                        tGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAAt                                             <  2:50957/141‑1 (MQ=255)
                                                                                                cGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATCAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAg                                     >  2:126384/1‑141 (MQ=255)
                                                                                                 gATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGc                                    >  1:160167/1‑141 (MQ=255)
                                                                                                     aaTTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCTTTa                                >  1:179792/1‑141 (MQ=255)
                                                                                                                   aaTATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCTTTAAGTGAGACGTCttt                  >  1:85325/1‑141 (MQ=255)
                                                                                                                          gATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCTTTAAGTGACACGTCTTTAACTCAg           <  1:208109/141‑1 (MQ=255)
                                                                                                                              gCATTAAACATCCGGATGATATGCGGATTGTAGATAGAGATAGTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCTTTAAGTGAGACGTCTTTAGCTCAGACTa       <  2:127053/141‑1 (MQ=255)
                                                                                                                                   aaaCATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCTTTAAGTGAGACGTCTTTAACTCAGACTAATCAt  >  1:209476/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                      
AAGAAGATGACTGTTACTACTTAGTAATGGAATATATCGAAGGTCCGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATGCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCTTTAAGTGAGACGTCTTTAACTCAGACTAATCAT  >  CP000730/1232328‑1232599

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: