Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 111,512 A→G 26.7% 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,5120AG26.7% 50.6 / 22.7 38S194S (TCA→TCGUSA300HOU_0110possible staphylococcal tandem lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (16/11);  new base G (5/5);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

CTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAA  >  CP000730/111373‑111650
                                                                                                                                           |                                                                                                                                          
ttACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAg                                                                                                                                           <  1:14360/140‑1 (MQ=255)
   ttgaagcagatAATGAAGAAAACAAAAAAAAATAAACCTAAAGCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAACGATGGAGATATATCATATAATCCAGAGGCGCTGAGTTATTCAGCTc                                                                                                                                        <  2:16829/131‑1 (MQ=25)
   cTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATGATCCACATGTGCCGAGTTATTCGGCta                                                                                                                                        <  1:206210/141‑2 (MQ=255)
             tAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                                                              >  1:31141/1‑141 (MQ=255)
                 gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                                                          >  1:96975/1‑141 (MQ=255)
                  aagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                                                         <  1:200391/141‑1 (MQ=255)
                    gaaAACATAAAAAAAGAAAGCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                                                                       >  2:25648/1‑141 (MQ=255)
                       aaCATAAAAAAAGAAATCGAAAACTTTAAGTACTTTGCACAATATCGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgat                                                                                                                    <  1:21698/141‑1 (MQ=255)
                        aCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAAGATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgatg                                                                                                                   >  1:24452/1‑141 (MQ=255)
                           tAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGCTGGAGTTAGATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGAtt                                                                                                                >  1:275589/1‑141 (MQ=18)
                           tAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAAGGATGAtt                                                                                                                >  1:26423/1‑141 (MQ=255)
                                   gAAATCGAAAATATTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                                                                                                        <  2:89578/141‑1 (MQ=255)
                                    aaaTCGAATACTTTAGTGTTTTTTCGCAATAGAGCAATTTTAAATATGTAATGAATTAGAAAGATGTATATATATCATATAAGCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                       <  1:57327/141‑1 (MQ=18)
                                    aaaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                       <  2:88039/141‑1 (MQ=255)
                                                  aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa                                                                                         >  1:224656/1‑141 (MQ=255)
                                                   aGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaaa                                                                                        >  1:267415/1‑141 (MQ=255)
                                                        tttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGAGGGATATAGATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGatat                                                                                   >  1:147768/1‑141 (MQ=255)
                                                        tttGCGCAATATAGCAATTTCAAAGATTTAATGAATTATAAAGATGGAAATATATCATATCATCCAGAGGTGCCGAGTTATTACGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGatat                                                                                   >  2:225665/1‑141 (MQ=18)
                                                          tGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGa                                                                                 <  1:160579/141‑1 (MQ=255)
                                                                          tttAAAGATTTAATCAATTATCAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCATTTAACTAATGATTATTATAATGTAAAACAATTACGTACAAGATATGACATAACAACCAATaaa                                                                 >  2:192282/1‑141 (MQ=18)
                                                                            tAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCACGAAAAAAAgc                                                               >  2:277134/1‑141 (MQ=255)
                                                                              aaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcgc                                                             >  1:40823/1‑141 (MQ=255)
                                                                               aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGcc                                                            >  1:78662/1‑141 (MQ=255)
                                                                                   ttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGAGCCGAAg                                                        >  1:187898/1‑141 (MQ=255)
                                                                                      aTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtg                                                     <  2:27933/141‑1 (MQ=255)
                                                                                      aTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtg                                                     <  1:36264/141‑1 (MQ=255)
                                                                                       tGAACTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgt                                                    <  1:107947/141‑1 (MQ=255)
                                                                                        gAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCCACCAGTAAAGCGCCGGAGCtgtt                                                   >  2:253474/1‑141 (MQ=255)
                                                                                              taAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGaa                                             <  1:108405/141‑1 (MQ=255)
                                                                                                  gATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGt                                         >  2:98161/1‑141 (MQ=255)
                                                                                                    tGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTAc                                       <  1:78284/141‑1 (MQ=255)
                                                                                                       agaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAgg                                    >  1:226208/1‑141 (MQ=255)
                                                                                                        gaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATATCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAggg                                   <  2:185032/141‑1 (MQ=255)
                                                                                                         atatatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGa                                  >  1:265090/1‑141 (MQ=255)
                                                                                                             atCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAttt                              <  1:225899/141‑1 (MQ=255)
                                                                                                                    aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAATATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGCAAGGTCCAGGGAATTTAAAAGGt                       >  2:202902/1‑141 (MQ=255)
                                                                                                                            ggTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGGATTTAAAAGGTTCATCAGt               >  1:27811/1‑141 (MQ=255)
                                                                                                                                         tCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaaaa  >  2:277479/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                          
CTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAA  >  CP000730/111373‑111650

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: