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% 118.2 / NA 36L48L (CTA→CTTlyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (14/22);  total (14/22)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAATTATTC  >  USA300TCH1516_ALE/2454491‑2454762
                                                                                                                                       |                                                                                                                                        
aaaTAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCAcc                                                                                                                                     <  1:738484/140‑1 (MQ=255)
aaaTAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCAcc                                                                                                                                     <  1:628957/140‑1 (MQ=255)
  aTAAAGTTTCTGGTTCTCAATGCGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCAt                                                                                                                                   <  1:644802/141‑1 (MQ=255)
     aaaTTGCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTTGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTaa                                                                                                                                <  1:482924/141‑1 (MQ=255)
       aTTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAAtg                                                                                                                              >  1:371341/1‑140 (MQ=255)
              gTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCaa                                                                                                                       <  1:404627/141‑1 (MQ=255)
                    aaTGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCAttt                                                                                                                 <  1:653091/141‑1 (MQ=255)
                        ggCAATGACGATTTTCTTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTACTGCAATCATTTGtat                                                                                                             >  1:405262/1‑141 (MQ=255)
                             tGACGCTTTTCGTCTTCTTTGTCATTACAATGGCGTTAGCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTGGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCg                                                                                                        <  1:655599/141‑1 (MQ=255)
                              gACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGt                                                                                                       >  1:246458/1‑141 (MQ=255)
                                    tttCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGAAATCACTTGTATACTCGttttaca                                                                                                 >  1:684259/1‑138 (MQ=255)
                                      tCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCaaaa                                                                                               >  1:4218/1‑140 (MQ=255)
                                      tCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAat                                                                                               <  1:120160/141‑1 (MQ=255)
                                       cGTCTTCTTTGTCATCACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAata                                                                                              <  1:91323/141‑1 (MQ=255)
                                                   cATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCa                                                                                  >  1:125393/1‑141 (MQ=255)
                                                   cATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCa                                                                                  >  1:696848/1‑141 (MQ=255)
                                                                  aTCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTaaa                                                                   <  1:538466/141‑1 (MQ=255)
                                                                            cTCAGAGAGTGCCAGGCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATc                                                         <  1:132263/141‑1 (MQ=255)
                                                                                        cAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTa                                             <  1:189713/141‑1 (MQ=255)
                                                                                        cAGTCTATAATTGAGGTCAAACACTCTATAATTGAAGCTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTa                                             <  1:32574/141‑1 (MQ=255)
                                                                                          gTCTATAATTGGTGTCAAACCCTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTGTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAAt                                           <  1:154467/141‑1 (MQ=255)
                                                                                                        tCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGc                             <  1:562579/141‑1 (MQ=255)
                                                                                                           aaCACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAg                          >  1:56747/1‑141 (MQ=255)
                                                                                                            acacTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTaag                         >  1:331994/1‑139 (MQ=255)
                                                                                                             cacTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAACCGTCTCAGTCATCCTGAAAGTAATTGAACGTCTATTGCTAGCt                        >  1:757703/1‑141 (MQ=255)
                                                                                                               cTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCttt                      <  1:338220/141‑1 (MQ=255)
                                                                                                                         gAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAGATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCt            >  1:353169/1‑141 (MQ=255)
                                                                                                                         gAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCACTCAGCCTGAAAGTAATTGAACGTATATTGCTAGATTTAATTTTACCt            >  1:54885/1‑141 (MQ=255)
                                                                                                                         gAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACTTCTATTGCTAGCTTTAATTTTCCCt            >  1:231193/1‑141 (MQ=255)
                                                                                                                           aagtcctgATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTGAAATATAAAAAGGTCCAACTTGCAGGTTTCGAATGCTCAATCAGCCTGAAAGTAACTGAACGTCTATTGCTAGCTTTAATTTTACCttt          <  1:355193/134‑1 (MQ=255)
                                                                                                                           aGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCttt          <  1:212168/141‑1 (MQ=255)
                                                                                                                            gTTACAGATCTTGCAGCATTAATTGCCGCAATCACTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAACATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGTTAGCTTTAATTTTACCTTTa         <  1:494289/141‑1 (MQ=255)
                                                                                                                            gTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTa         <  1:428078/141‑1 (MQ=255)
                                                                                                                               aCAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAatt      >  1:670632/1‑141 (MQ=255)
                                                                                                                                  gATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAattatt   <  1:733065/141‑1 (MQ=255)
                                                                                                                                   aTCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAATTATTc  <  1:235478/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
CAATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAATTATTC  >  USA300TCH1516_ALE/2454491‑2454762

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: