Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 111,512 A→G 35.2% S194S (TCA→TCG)  USA300HOU_0110 → possible staphylococcal tandem lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730111,5120AG35.2% 57.6 / 37.7 48S194S (TCA→TCG) USA300HOU_0110possible staphylococcal tandem lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (14/17);  new base G (6/11);  total (20/28)
Fisher's exact test for biased strand distribution p-value = 5.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAAT  >  CP000730/111374‑111652
                                                                                                                                          |                                                                                                                                            
cacggAGGAAGTTAAAGCAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGAGGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGc                                                                                                                                            <  1:182209/137‑1 (MQ=255)
    gAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCAAGAGGTACGGAGGTATTCAggcaaa                                                                                                                                        >  2:394711/1‑136 (MQ=25)
     aGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAat                                                                                                                                       <  1:7031/141‑1 (MQ=255)
         aGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCa                                                                                                                                   >  1:62135/1‑141 (MQ=255)
          gTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCaa                                                                                                                                  >  2:409515/1‑141 (MQ=255)
          attAAATAAAGAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                                                                                                  <  2:75004/140‑1 (MQ=17)
               agaagaAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAAc                                                                                                                             <  2:426653/141‑1 (MQ=255)
                aaagaaAACATAAAAAAAGAAATCGAAAACTTGAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                                                            <  2:419000/140‑1 (MQ=255)
                 aagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATGTGGCGATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTa                                                                                                                           <  1:358231/141‑1 (MQ=255)
                 aagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTa                                                                                                                           <  1:325505/141‑1 (MQ=255)
                 aagaaAACATAAAAAAAGAAATCGAAAACGTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTa                                                                                                                           <  1:119614/141‑1 (MQ=255)
                        cATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTTTTTTTTTTTTCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTTCCGAGTTATTCGGCTAAATATAAAATAAGAAatgaaaa                                                                                                                    >  2:323264/1‑138 (MQ=14)
                        cATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAtgatga                                                                                                                    <  1:152171/141‑1 (MQ=255)
                         aaaaaaaaaaaaACCGAAAACTTTAAGTTTTTTGCCCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                                                                                                   <  1:234079/139‑1 (MQ=17)
                                  gAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATAAAATCCAGAGGTGCCGAATTATTAGGCTAAATATAAATTAACTAATGATGATTATAATGTa                                                                                                          >  2:216838/1‑141 (MQ=18)
                                    aaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATCTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaaa                                                                                                        <  2:434718/141‑1 (MQ=255)
                                     aTCGAAAACTTTAAGTTTTTTGCGCACTCAATCAGTTGTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAAGGATGATTATAATGTaaaa                                                                                                       >  1:110049/1‑141 (MQ=21)
                                     aTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaaaa                                                                                                       >  2:16338/1‑141 (MQ=255)
                                     aTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATGTAATGAATTATAAAGATGGAGATATATAATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                                                                                                       >  2:199268/1‑141 (MQ=255)
                                         aaaaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                                                   <  2:91731/141‑1 (MQ=255)
                                         aaaaCTTTAAGTGTTTGGCGCAATATAGCAATTGTAAAGAGTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                                                   <  2:49884/141‑1 (MQ=17)
                                          aaaCTTTAAGTTTTTTGCGCAATATAGCAATTTGAAAGAGTTAATGAATTATAAAGATTGAGATATATCATATAATGCAGAGGTGCCGAGTTATTAGGCTAAATATAAATTAAATAATGATGATTATAATGTAAAAGAAtt                                                                                                  >  2:160634/1‑141 (MQ=18)
                                            aCTTGAAGTTTTTTGCGAAATATAGAAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc                                                                                                <  1:36104/141‑1 (MQ=255)
                                             cTTTACGCTTTTTTCTCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGACGGATATATATTATATAATCCAGAGGTGCCGAGTTATTCAGATCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACg                                                                                               <  1:31238/141‑1 (MQ=17)
                                              tttAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGATATATATCATATAATCCAGAGGGGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACgc                                                                                              <  2:115438/141‑2 (MQ=255)
                                              tttAAGTTGTTTCAGCAATATATCAATTTTAAAGATTTAATGAATTAGAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACgc                                                                                              <  2:453678/141‑2 (MQ=17)
                                                 aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa                                                                                           >  1:3356/1‑141 (MQ=255)
                                                  aGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTaaaa                                                                                          >  1:160544/1‑141 (MQ=255)
                                                    ttttttGAGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGa                                                                                        <  2:381989/141‑1 (MQ=255)
                                                     tttttGGACAATATAGTAATTTTAAATATTCAAGTAATTTTAAAGATGTATTTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAACATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAgtt                                                                                       <  2:86902/141‑3 (MQ=17)
                                                     tttgtGAGCAATATAGCAATTTTAAATATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGat                                                                                       <  2:376437/141‑1 (MQ=14)
                                                                  taGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaac                                                                          >  2:194503/1‑141 (MQ=255)
                                                                   aGCAATTTTAAAGAGTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacg                                                                         <  2:257539/141‑2 (MQ=11)
                                                                    gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAccgacga                                                                        <  1:66206/141‑6 (MQ=11)
                                                                    gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAAAAATTACGAAAACGATTTGATATACcgacga                                                                        >  2:69291/1‑139 (MQ=14)
                                                                            aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAACAGATATGACCTACCAACAAATAAAgcc                                                                >  2:90681/1‑140 (MQ=255)
                                                                              acaTTTAATGAATTATAAAGATGGACATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGcc                                                              <  2:133931/139‑1 (MQ=255)
                                                                                   tAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGc                                                         <  2:454637/141‑1 (MQ=255)
                                                                                            ataAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGa                                                <  2:329536/141‑1 (MQ=255)
                                                                                                        atatatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGa                                    <  2:160544/141‑1 (MQ=255)
                                                                                                                     accAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtc                       <  2:19554/140‑1 (MQ=255)
                                                                                                                      ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTTTAATGTAAGACAATTACGTAAAAGATATGGCATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtca                      <  1:335465/141‑1 (MQ=255)
                                                                                                                       cAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAGATAAAGCGCCGAAGCTGTTGTTGGAAGGTACAGGGAATTTAAAAGGTtcat                     >  1:118523/1‑141 (MQ=255)
                                                                                                                         gagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatca                   >  1:388212/1‑141 (MQ=255)
                                                                                                                            gTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACTAATAAAGCGCCGCGGCTGTTGTTGAGAGGTACCGGGAATTTAAAAGGTTCGTCAGtt                >  1:132687/1‑141 (MQ=17)
                                                                                                                            gTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAGAGTTCATCAggt                >  1:123755/1‑139 (MQ=255)
                                                                                                                              gCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTgg              >  2:46410/1‑141 (MQ=255)
                                                                                                                                       ttCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGAACAGGGAAGTTAAAGGGTTCATCAGTTGGATATaaaaa     >  1:322440/1‑141 (MQ=255)
                                                                                                                                          aGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAAt  >  1:319851/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                            
TACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAAT  >  CP000730/111374‑111652

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: ‑