Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,455,692 Δ1 bp 69.7% coding (1210/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,6920G.69.7% 92.7 / 53.3 46coding (1210/1260 nt)lyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base G (5/5);  new base . (18/14);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATCAGCAGAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATT  >  USA300TCH1516_ALE/2455554‑2455823
                                                                                                                                          |                                                                                                                                   
gAATCAGCAGAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCaaaa                                                                                                                                   >  1:99372/1‑141 (MQ=255)
      acagAATCAGTTAAACAATCAGCACAAATTGACGATATAAAAAACGATTAAAATGAAGATAAAATAAAGCAATGATAAAAAGAGCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaaaaa                                                                                                                             <  2:522324/140‑4 (MQ=255)
                 tAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACt                                                                                                                 >  1:68235/1‑141 (MQ=255)
                 tAAACAAGCACCACAAAGTGAAGATTTAACAAACGATTCAAATGACGATGGAATAGAGCAATCATTAAAAGAACCTGCGAGTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATTCAAAAAAATTATCGAAGAAGCt                                                                                                                 >  1:583146/1‑141 (MQ=255)
                   aaCAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACAGATCGAAAAAACTg                                                                                                                >  2:402410/1‑141 (MQ=255)
                       aGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaagaa                                                                                                           <  1:28334/141‑1 (MQ=255)
                          accacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACa                                                                                                         <  2:27467/141‑1 (MQ=255)
                            cacaAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATc                                                                                                      >  1:542971/1‑141 (MQ=255)
                              caAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCtt                                                                                                    >  1:529689/1‑141 (MQ=255)
                              caAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCtt                                                                                                    >  2:467703/1‑141 (MQ=255)
                                   tGACGATTGAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGAt                                                                                               <  1:110813/141‑1 (MQ=255)
                                          ttAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAAc                                                                                         <  1:374998/141‑1 (MQ=255)
                                                accGATTCAAATGAAGACAAAATAGAGGAATCATTAAACGATCCTGCGACTTATAAAGAAGACAGATGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaataa                                                                                   <  2:455860/139‑1 (MQ=255)
                                                  agATTCAAATGAAGATGAAACAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATAATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAAtaataa                                                                                <  2:537021/140‑1 (MQ=255)
                                                   gATTCAAATGAAGATGAAACAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAAc                                                                               <  1:18539/141‑1 (MQ=255)
                                                     ttCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACa                                                                              >  1:145292/1‑141 (MQ=255)
                                                      tCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGGAGACAAACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaaa                                                                            <  2:218742/141‑1 (MQ=255)
                                                         aaTGAAGATGAAATAGCGCAATCATTAAAAGAAGCTGCGACTTATAAAGAAGATAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgc                                                                          <  2:500573/141‑1 (MQ=255)
                                                         aaTGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgc                                                                          >  2:427304/1‑141 (MQ=255)
                                                         aaTAAAGAAGAAAAAGAGCAATCATTAAAAGAACCTGCAACTTATAAAGAAGACAGACGTTCATCAGTTGAAATAAAAAAAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgcg                                                                         <  2:542832/141‑1 (MQ=255)
                                                           tGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCt                                                                       <  1:166431/141‑1 (MQ=255)
                                                               gATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCt                                                                   <  1:587174/141‑1 (MQ=255)
                                                                  gAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGGCAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGt                                                                <  1:322128/141‑1 (MQ=255)
                                                                        gagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctct                                                          >  1:421427/1‑141 (MQ=255)
                                                                            aaTCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTGGCTAGTACCTCTTaa                                                       >  1:443433/1‑141 (MQ=255)
                                                                              tCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaa                                                     <  1:520430/141‑1 (MQ=255)
                                                                              tCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaaa                                                    >  1:128690/1‑141 (MQ=255)
                                                                                 ttAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgaa                                                 >  2:215703/1‑140 (MQ=255)
                                                                                    aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgatga                                               >  1:251373/1‑141 (MQ=255)
                                                                                      aGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGc                                            >  1:176662/1‑141 (MQ=255)
                                                                                         aCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAGAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTAGCTAGTACCTCTTAAAAAGATGATGCTaa                                         >  2:411365/1‑140 (MQ=255)
                                                                                         aCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAg                                         <  2:68235/141‑1 (MQ=255)
                                                                                         aCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAg                                         >  1:284641/1‑141 (MQ=255)
                                                                                          ccTGCGACTTATAACGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGc                                        >  2:413992/1‑141 (MQ=255)
                                                                                          ccTGCGACTTATAAAGAAGACAGACGTACATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAGCAAATAATAACAAAGGGCTCGCTAGTACCTCTTAAAAAAATGATGCTaaa                                        >  2:416720/1‑139 (MQ=255)
                                                                                           ctacGACTCAAAAAGAAGAAAGACGTTCATCAGTTGAAAAAAAAAAAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCa                                       <  2:230934/138‑1 (MQ=255)
                                                                                              cGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACGCCTTAAAAAAATGATGCTAGCAAgg                                    >  1:107247/1‑140 (MQ=255)
                                                                                                aCTTATAAAGAAGACCGACGTTCATCAGTTGTAATTGATGCAATAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATATCAAAGCGCTTGCTAGTACCTCTTAAAAAGATGAGGCTAGCATgcgc                                  <  1:430522/141‑1 (MQ=255)
                                                                                                aCTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcgc                                  <  2:251828/141‑1 (MQ=255)
                                                                                                   tataAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcgtaat                               <  2:6040/141‑5 (MQ=255)
                                                                                                   tataAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCttt                               <  1:144505/141‑1 (MQ=255)
                                                                                                         gaagaCAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctata                         >  1:107657/1‑141 (MQ=255)
                                                                                                            gacagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGGTAGCAACCGATTTTctatacta                      >  2:393948/1‑141 (MQ=255)
                                                                                                                 acGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTatatat                 >  1:213092/1‑141 (MQ=255)
                                                                                                                         caGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTc         <  2:24142/141‑1 (MQ=255)
                                                                                                                             tGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTAAATATTATTTTTCTaaa     >  1:4996/1‑141 (MQ=255)
                                                                                                                              gTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTAGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTaaaa    >  2:188147/1‑141 (MQ=255)
                                                                                                                                aaTTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAAtt  >  1:376172/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                   
GAATCAGCAGAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATT  >  USA300TCH1516_ALE/2455554‑2455823

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: