Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,232,465 G→C 100% A124P (GCG→CCG)  prkC → Serine/threonine‑protein kinase PrkC

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

AAGAAGATGACTGTTACTACTTAGTAATGGAATATATCGAAGGTCCGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATGCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCTTTAAGTGAGACGTCTTTAACTCAGACTAATCAT  >  USA300TCH1516_ALE/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)
                                   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  >  USA300TCH1516_ALE/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: