Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 111,512 A→G 20.1% 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,5120AG20.1% 45.9 / 9.1 25S194S (TCA→TCGUSA300HOU_0110possible staphylococcal tandem lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (13/7);  new base G (2/3);  total (15/10)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

CTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAA  >  CP000730/111376‑111650
                                                                                                                                        |                                                                                                                                          
cTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTCTAAGTTTTTTGCGCAATATAGCAATTTTAAAGAGTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTc                                                                                                                                        >  1:116770/1‑141 (MQ=255)
    ggAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAata                                                                                                                                    >  2:41897/1‑141 (MQ=255)
         ttAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                                                                                               >  1:16253/1‑141 (MQ=255)
                        tAAAAAAAGAAATCGAAAACTTTAAGTTTTTTTCTCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGAtt                                                                                                                >  2:8731/1‑141 (MQ=255)
                          aaaaaaGAAATCGAAATGTTTAAGGTTCTTGCGCAATATAGCAATTTTAAAATTTTAATGAATTATAACGCTGGACATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCGATTAACTATTGATGATtat                                                                                                              <  2:109264/141‑1 (MQ=18)
                               aGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGt                                                                                                         <  1:48138/141‑1 (MQ=255)
                                 aaaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                       <  1:84485/141‑1 (MQ=255)
                                        aaaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                                                                                                <  2:67761/141‑1 (MQ=255)
                                                   tttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGat                                                                                     <  1:30258/141‑1 (MQ=255)
                                                        gcgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGAc                                                                                >  1:12731/1‑141 (MQ=255)
                                                         cgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGata                                                                               <  2:99871/141‑3 (MQ=255)
                                                                    aaTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAAt                                                                    >  2:130098/1‑141 (MQ=255)
                                                                    aaTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATAAAATCCAGAGGTGCCGAGATATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAAt                                                                    >  2:23106/1‑141 (MQ=255)
                                                                           aaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcgc                                                             >  1:45105/1‑141 (MQ=255)
                                                                                               gATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGTTGTTGTTGAAAGGt                                         <  1:86640/141‑1 (MQ=255)
                                                                                                    agaTATATCATAGAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAATATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAgg                                    <  2:116770/141‑1 (MQ=255)
                                                                                                      atatatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGa                                  >  2:122960/1‑141 (MQ=255)
                                                                                                        atatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAt                                <  1:18403/141‑1 (MQ=255)
                                                                                                               ataATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAg                         >  1:46072/1‑141 (MQ=255)
                                                                                                                 aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGCAAGGTACAGGGAATTTAAAAgga                       >  2:58664/1‑140 (MQ=255)
                                                                                                                 aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGt                       >  1:16546/1‑141 (MQ=255)
                                                                                                                 aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGt                       <  2:108224/141‑1 (MQ=255)
                                                                                                                 aaTCCAGAGATGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGt                       >  2:33079/1‑141 (MQ=255)
                                                                                                                      agagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatc                  >  2:128339/1‑141 (MQ=255)
                                                                                                                           tGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAGAAGGTATAAATAACCAGAAAAGATAGCGCCGAAGCTGTAGATGAAAGGTACAGGGAATTAAAAAGGTTCATCAGTTg             >  1:104184/1‑141 (MQ=17)
                                                                                                                                      taaGCTCAATATAAATTAACTAATGATGATTATAATGTACAACAATTACGTAAAAGATATGAAATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTACAAGGATAATCAGTTGCATATaaaaaa  >  2:20220/3‑141 (MQ=17)
                                                                                                                                        |                                                                                                                                          
CTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAA  >  CP000730/111376‑111650

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: