Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 39,378 C→A 100% D586Y (GAT→TAT)  pbp ← Beta‑lactam‑inducible penicillin‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE39,3780CA100.0% 50.3 / NA 16D586Y (GAT→TAT) pbpBeta‑lactam‑inducible penicillin‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (12/4);  total (12/4)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCATTA  >  USA300TCH1516_ALE/39256‑39514
                                                                                                                          |                                                                                                                                        
gTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTGATTTTAAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTtttatt                                                                                                                        <  1:387499/141‑1 (MQ=255)
   tGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTAc                                                                                                                     <  1:92137/141‑1 (MQ=255)
               aTCATATGATATAAACCACGCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGGGTTCTGCAGTGCCGGATTTGCCGATTAAGTTTGCATAAGATCTGTAAGTATATTCTTTATGTGTTTTATTTACGACTTGTTGcgt                                                                                                         >  1:430203/1‑139 (MQ=255)
                 cATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATAc                                                                                                       >  1:268408/1‑141 (MQ=255)
                        tataAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGt                                                                                                >  1:178250/1‑141 (MQ=255)
                        tataAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTGGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTGTATTTACGACTTGTTGCATACCATCAGt                                                                                                >  1:325335/1‑141 (MQ=255)
                                                           gTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATAtttt                                                             <  1:311819/141‑1 (MQ=255)
                                                                          cTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGt                                              >  1:268446/1‑141 (MQ=255)
                                                                          cTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGt                                              >  1:435742/1‑141 (MQ=255)
                                                                                 aCCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCgg                                       >  1:213394/1‑140 (MQ=255)
                                                                                         tGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTta                               >  1:375993/1‑141 (MQ=255)
                                                                                                                  ataAATATATTCTTTTTGTGTTTTATTTTCGTCTTTTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGcc      <  1:242458/141‑1 (MQ=255)
                                                                                                                  ataAATATATTCTTTATGTGTTTTGTTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGATGTGCGTTAATATTGca      >  1:537856/1‑140 (MQ=255)
                                                                                                                   taAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCa     >  1:342233/1‑141 (MQ=255)
                                                                                                                   taAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATAACATCAGTTCATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAGCTTGTTTTTTGGTGTATTTTAAGCAGTGAGGAGCGGTAATATTGCCa     >  1:333255/1‑141 (MQ=255)
                                                                                                                     aaTATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAACAAGTGAGGTGCGGTAATATTGCCatt   >  1:159671/1‑141 (MQ=255)
                                                                                                                      atatatTGTTTATGTGTTTTATTTACGACTTGTTGCATACCTTCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTCCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCatta  <  1:251706/141‑1 (MQ=255)
                                                                                                                          |                                                                                                                                        
GTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCATTA  >  USA300TCH1516_ALE/39256‑39514

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