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% 110.6 / NA 33A124P (GCG→CCG) prkCSerine/threonine‑protein kinase PrkC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (17/16);  total (17/16)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAGATGACTGTTACTACTTAGTAATGGAATATATCGAAGGTCCGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATGCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCTTTAAGTGAGACGTCTTTAACTCAGACTAATCA  >  USA300TCH1516_ALE/1232331‑1232598
                                                                                                                                      |                                                                                                                                     
aagaTGACTGTTACAACATAATAATGAAATATATCGAAGGTCCGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATg                                                                                                                                 <  1:547782/141‑1 (MQ=255)
        tGTTACTACTTAGTAATGGAATATATCGAAGGTCCGACTTTGTCTGGGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGt                                                                                                                         >  1:663794/1‑141 (MQ=255)
               actTATTAATGGAATATATCGAAGGTCCGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTAc                                                                                                                  <  1:386258/141‑1 (MQ=255)
                        tGGAATATAGCGAAGGCACGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGata                                                                                                         <  1:159539/141‑1 (MQ=255)
                              atatCGAAGGTCCGACGTTGTCTGATTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGc                                                                                                   <  1:377195/141‑1 (MQ=255)
                              agatCGAAGGTCCGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGc                                                                                                   <  1:261898/139‑1 (MQ=255)
                                  cGAAGGTCCGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCcaca                                                                                               >  1:276697/1‑141 (MQ=255)
                                  cGAAGGTCCGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCcaca                                                                                               >  1:298063/1‑141 (MQ=255)
                                          cGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAgtatat                                                                                       >  1:206896/1‑141 (MQ=255)
                                                tGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTg                                                                                 >  1:609193/1‑141 (MQ=255)
                                                tGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTg                                                                                 >  1:609113/1‑141 (MQ=255)
                                                    tGAGTATATTGAAAGTCTTGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATGTTAAGCCACAAAATATATTAATTGACAg                                                                             >  1:56531/1‑141 (MQ=255)
                                                     gAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAAAATATTATTTGACAGc                                                                            >  1:508798/1‑141 (MQ=255)
                                                              gAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACg                                                                   <  1:384687/141‑1 (MQ=255)
                                                                 aGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGAAATTAAGCAAAAAAATATATTAATGGACAGCAATAAAACATTg                                                                >  1:534367/1‑141 (MQ=255)
                                                                     aTGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGCCAGCAATAAAACGTTGaaaa                                                            >  1:548602/1‑141 (MQ=255)
                                                                      tGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTCATTATGGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAat                                                           >  1:658681/1‑141 (MQ=255)
                                                                      tGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAAAGTTGaaaaa                                                           >  1:586697/1‑140 (MQ=255)
                                                                               aaGGGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCAAAAAATTTATTAATTGACAGCAATAAAACGTTGAAAATATTTGAttt                                                  <  1:159622/141‑1 (MQ=255)
                                                                                    ttGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGaa                                             <  1:563332/141‑1 (MQ=255)
                                                                                    ttGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGaa                                             <  1:563298/141‑1 (MQ=255)
                                                                                        cacaGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGc                                         <  1:331897/141‑1 (MQ=255)
                                                                                         acaGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCt                                        <  1:351731/141‑1 (MQ=255)
                                                                                                       tACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCTTTAAGTga                          >  1:78184/1‑141 (MQ=255)
                                                                                                       tACGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCATATTGTACATAGACATATTAAGCCACAAAATATATTAGTTGACAGCAATAAAACGTTGAAAGTATTTGAGTTTGGAATTGGTAAAGCTTTAAGTga                          >  1:408862/1‑141 (MQ=255)
                                                                                                        aCGAATCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCTTTAAGTgag                         >  1:506568/1‑141 (MQ=255)
                                                                                                             tCAAATATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCTTTAAGTGCGACGTc                    >  1:138706/1‑141 (MQ=255)
                                                                                                               aaaTATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCTTTAAGTGAGACGTCtt                  <  1:31690/141‑1 (MQ=255)
                                                                                                               aaaTATTGGATGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGCAATTGCTAAAGCTTCAAGTGAGACGTCtt                  <  1:359634/141‑1 (MQ=255)
                                                                                                                         tGGCATTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCTTTAAGTGAGACGTCTTTAACTCAGAc        <  1:364754/141‑1 (MQ=255)
                                                                                                                             aTTAAACATCCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCTTTAAGTGAGACGTCTTTACCTCAGACTAAt    >  1:52351/1‑141 (MQ=255)
                                                                                                                               tAAACATCCGCATGAGATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATAATTTATTTTGGAATTGCTAAAGCTTTAAGTGAGACGTCTTTAACTCAGACTAATCa  <  1:19612/141‑1 (MQ=255)
                                                                                                                               tAAACAGCCGCATGATATGAGAATTGTACATAGAAATATTAAGCCACAAAATATATGAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCGTTTAGTGATACGTCTTTAACTCAGACTAATCa  <  1:160262/141‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
AAGATGACTGTTACTACTTAGTAATGGAATATATCGAAGGTCCGACTTTGTCTGAGTATATTGAAAGTCATGGGCCATTAAGTGTTGACACAGCGATTAATTTTACGAATCAAATATTGGATGGCATTAAACATGCGCATGATATGCGTATTGTACATAGAGATATTAAGCCACAAAATATATTAATTGACAGCAATAAAACGTTGAAAATATTTGATTTTGGAATTGCTAAAGCTTTAAGTGAGACGTCTTTAACTCAGACTAATCA  >  USA300TCH1516_ALE/1232331‑1232598

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