Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 111,512 A→G 19.0% 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,5120AG19.0% 71.9 / 13.7 37S194S (TCA→TCGUSA300HOU_0110possible staphylococcal tandem lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (21/9);  new base G (4/3);  total (25/12)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

TGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAA  >  CP000730/111377‑111651
                                                                                                                                       |                                                                                                                                            
tGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCa                                                                                                                                         <  1:111946/141‑1 (MQ=255)
   ggAAGTTACAGAGTACAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAata                                                                                                                                      <  1:119935/141‑1 (MQ=255)
          aaaGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTa                                                                                                                               >  1:176528/1‑141 (MQ=255)
           aagaagAAAACAAAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                                                              <  1:185391/141‑1 (MQ=255)
            agaagaAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATGTCAATTAAc                                                                                                                             <  2:192596/141‑1 (MQ=255)
             gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACt                                                                                                                            <  2:121003/141‑1 (MQ=255)
                aaaaaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTGTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATAGATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                                                                         <  1:141815/140‑1 (MQ=17)
                   aaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                                                                      >  1:77934/1‑141 (MQ=255)
                   aaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                                                                      >  1:124611/1‑141 (MQ=255)
                   aaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgat                                                                                                                      >  1:130446/1‑141 (MQ=255)
                      aTAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATgatgat                                                                                                                   <  1:105364/141‑1 (MQ=255)
                                 aaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaaa                                                                                                        >  2:30002/1‑141 (MQ=255)
                                      aaaaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATGTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATACAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                                                   >  2:78455/1‑141 (MQ=255)
                                        aaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTa                                                                                                 >  1:257237/1‑141 (MQ=255)
                                          cTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACg                                                                                               >  2:192658/1‑141 (MQ=255)
                                             tAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa                                                                                            <  2:193443/141‑1 (MQ=255)
                                                             tataGCAATTTTAAAGATTTAATGGATTGTAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGAAAAAGATATGATATAccg                                                                            >  1:171737/1‑140 (MQ=14)
                                                                 gCAAATTTAAAGATTTAATGAATTATAAAGCTGGAGATATATCATATAATCCAGAGGTGCCGAGTTATACAGCTCAATATCAATTACCTAAAGATGATTATAATGTAAAACAATTACGTAAAAGATATGCCATACcaacaa                                                                       >  2:164571/1‑141 (MQ=17)
                                                                         aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcg                                                               >  2:247191/1‑141 (MQ=255)
                                                                           aGATTGAATGAATTATAAAGATGGATCTAGACCATAAAATCTAGAGGTCCCGAGTAATTCAGAAAAATATGAATTAAGTAATGATGCTTATAGTGAAAGTCAATTAAGAAAAAGATAAGACATACCAAAAACTaaacactc                                                             >  2:130938/1‑136 (MQ=255)
                                                                            gATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACACTTACGTAAAAGATATGACAGACGAACAAATAAAGCGcca                                                            >  1:32374/1‑140 (MQ=255)
                                                                            gATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCg                                                            >  2:102886/1‑141 (MQ=255)
                                                                                          taAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGaa                                              >  2:265664/1‑141 (MQ=255)
                                                                                            aaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGACGTTGTTGAAAg                                            >  2:27738/1‑141 (MQ=255)
                                                                                              gATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGt                                          <  2:170075/141‑1 (MQ=255)
                                                                                                  gagaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAg                                      <  2:7770/141‑1 (MQ=255)
                                                                                                       atatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAt                                 >  1:79668/1‑141 (MQ=255)
                                                                                                         atCATATAAGCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATTATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTTTGGTTGAAAGGTACAGGGAAttt                               <  2:58161/141‑1 (MQ=255)
                                                                                                               taATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAgg                         <  1:33190/141‑1 (MQ=255)
                                                                                                                 aTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGGACAGGGAATTTAAAAGGtt                       >  1:263931/1‑141 (MQ=255)
                                                                                                                         gTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGtt               >  1:238450/1‑141 (MQ=255)
                                                                                                                          tGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTg              >  1:247070/1‑141 (MQ=255)
                                                                                                                             cGAGTTAGTCAGCTCAATATCAATTAACTAATGATGATTATAATATAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTCGat           >  2:190423/1‑141 (MQ=255)
                                                                                                                                  tATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaa      >  2:205055/1‑141 (MQ=255)
                                                                                                                                  aaTTCAGCTCAATATCAATTTCTTCTTGTTGTTGATAATGTAAAGCAATTACGTAAAAGAAATGAAATACAAGACAATTAAACGCCAAAGAAGTTTTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaa      >  2:70287/2‑141 (MQ=9)
                                                                                                                                     tCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaaaa   >  1:178906/1‑141 (MQ=255)
                                                                                                                                      cAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCAGAAGCTGTGGTTGAAAGGAACAGGGAATTTAAAGGGTTCATCAGTTGGATATaaaaaaa  >  2:88258/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                            
TGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAA  >  CP000730/111377‑111651

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: