Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,454,626 A→T 100% L48L (CTA→CTT lyrA → Lysostaphin resistance protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,454,6260AT100.0% 60.8 / NA 19L48L (CTA→CTTlyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (7/12);  total (7/12)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAA  >  USA300TCH1516_ALE/2454493‑2454756
                                                                                                                                     |                                                                                                                                  
aTAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACGGATCTTGCACCAt                                                                                                                             <  1:272342/141‑1 (MQ=255)
   aaaTTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTaa                                                                                                                          <  1:209803/141‑1 (MQ=255)
    aaTTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAAt                                                                                                                         >  1:357491/1‑141 (MQ=255)
                     gggCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATCATAATGTTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGta                                                                                                        >  1:220556/1‑141 (MQ=255)
                     gggCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAGTCATTTGta                                                                                                        >  1:220533/1‑141 (MQ=255)
                     gggCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGta                                                                                                        <  1:397915/141‑1 (MQ=255)
                                                           gCGTTATCCATTATGCTCAGAGTTACTAAGTCGATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGt                                                                  >  1:460882/1‑141 (MQ=255)
                                                               tATCCATTAGGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTaa                                                              <  1:180557/141‑1 (MQ=255)
                                                                                gaTTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTg                                             <  1:306155/141‑1 (MQ=255)
                                                                                         tCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAGTATAAAAAGGTCCAACTTGCAGGTTGAAAATTCGCAATCAGCCTGAAAGTAAtt                                    <  1:301211/141‑1 (MQ=255)
                                                                                              aaTTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACTATGAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACg                               <  1:85861/141‑1 (MQ=255)
                                                                                               aTTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGt                              <  1:241831/141‑1 (MQ=255)
                                                                                                    tgtCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTAtt                         <  1:179865/141‑1 (MQ=255)
                                                                                                      tCAAACAGTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATGCTGGTACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGGTTAAAACTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGc                       <  1:455970/141‑1 (MQ=255)
                                                                                                           cacTTTTTATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCt                  <  1:160927/141‑1 (MQ=255)
                                                                                                                  tatTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAACGTCATTGAACGTATATTGCTAGCTTTAAttt           <  1:93199/141‑1 (MQ=255)
                                                                                                                       gAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATGAGGATGAAGGTAAGTGAACGTCTATTGCTAGCTTTAATTTTACCt      >  1:433190/1‑141 (MQ=255)
                                                                                                                       gAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCt      >  1:155916/1‑141 (MQ=255)
                                                                                                                           ttACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAACGTAATTGAACGTCTATTGCTAGCTTGAATTTTACCTTTat  >  1:375269/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                                  
ATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAA  >  USA300TCH1516_ALE/2454493‑2454756

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