Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,455,698 2 bp→AT 100% coding (1216‑1217/1260 nt) lyrA → Lysostaphin resistance protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,455,6980CA90.3% 81.6 / 3.4 31coding (1216/1260 nt)lyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (11/16);  minor base . (2/1);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.
*USA300TCH1516_ALE2,455,6990AT100.0% 94.7 / ‑6.5 32H406L (CAT→CTT) lyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (13/18);  total (13/19)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACC  >  USA300TCH1516_ALE/2455565‑2455831
                                                                                                                                     ||                                                                                                                                    
accAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaa                                                                                                                               <  2:166275/139‑1 (MQ=255)
            aGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaagaa                                                                                                                   <  2:268206/141‑1 (MQ=255)
            aGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaagaa                                                                                                                   >  1:91967/1‑141 (MQ=255)
              caccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGCAGAACa                                                                                                                 >  2:332674/1‑141 (MQ=255)
                       gTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                                                                        <  1:58064/141‑1 (MQ=255)
                        tGACTATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGAt                                                                                                       <  1:368645/141‑1 (MQ=255)
                         gACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATa                                                                                                      <  1:87077/141‑1 (MQ=255)
                                      aCGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAataata                                                                                         >  1:307108/1‑141 (MQ=255)
                                             aaaTGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACGATCTTCAGATAAAAACAAATAATAACAAAgc                                                                                  >  1:186263/1‑141 (MQ=255)
                                                tGAAGCTGAAATAAAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATAATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCt                                                                               <  2:150933/141‑1 (MQ=255)
                                                tGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgcgca                                                                               >  1:93457/1‑140 (MQ=255)
                                                  aaaaTGAAATAGAGCAATAATTAAAAGAACCGGCGAATTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTg                                                                             <  2:135418/138‑1 (MQ=255)
                                                          aTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAACGCGCTTGCTAGTAcc                                                                     <  1:155564/141‑1 (MQ=255)
                                                          aTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAcc                                                                     >  1:246652/1‑141 (MQ=255)
                                                              agCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTGGCTAGTACCTCtt                                                                 >  1:310940/1‑141 (MQ=255)
                                                                 aagcATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaa                                                              <  1:2847/138‑1 (MQ=255)
                                                                 aaTCATTAAAAGAACCTGCGACTTATAAAGAAGACATACGTTCATCCGTTATAATTTATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaa                                                              <  1:366188/141‑1 (MQ=255)
                                                                      tgaaaaGAACCCGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCATCAGACAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTACAAAgat                                                         <  1:284957/139‑1 (MQ=255)
                                                                         aaaGAACCTGCGACGGATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAATACGTCTTAAAAAgatgat                                                      >  2:32288/1‑141 (MQ=255)
                                                                                cTGCGACTTATATATAAGACAGCTCTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACACATAATAACAACGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCa                                               <  2:333580/141‑1 (MQ=255)
                                                                                    gACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcg                                           <  1:158621/141‑1 (MQ=255)
                                                                                     aCTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcgc                                          <  1:275586/141‑1 (MQ=255)
                                                                                      attATAAAGAAGAGAGACGTACATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCt                                         <  2:71328/140‑1 (MQ=255)
                                                                                            aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTcta                                   >  2:120959/1‑141 (MQ=255)
                                                                                                agacagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGATAGTCCCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatact                               >  2:141821/1‑141 (MQ=255)
                                                                                                  acagacGTTCATCAGTTGTAATTGATATAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctaaactaa                             >  2:80134/1‑140 (MQ=255)
                                                                                                      acGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTatatat                         >  1:229182/1‑141 (MQ=255)
                                                                                                      acGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTatatat                         <  2:94726/141‑1 (MQ=255)
                                                                                                         ttCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTAGAtattat                      >  2:232885/1‑141 (MQ=255)
                                                                                                            atcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTTCCACGTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTAtttt                   <  1:332674/141‑1 (MQ=255)
                                                                                                               aGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAACAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCt                <  2:237513/141‑1 (MQ=255)
                                                                                                                              aaaaaaaaTTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTATCAcc  <  2:258185/139‑1 (MQ=255)
                                                                                                                                     ||                                                                                                                                    
ATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACC  >  USA300TCH1516_ALE/2455565‑2455831

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: