Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 111,516 C→A 39.2% Q196K (CAA→AAA)  USA300HOU_0110 → possible staphylococcal tandem lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730111,5160CA39.2% 41.3 / 36.7 36Q196K (CAA→AAA) USA300HOU_0110possible staphylococcal tandem lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (12/10);  new base A (6/8);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

CATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTA  >  CP000730/111368‑111660
                                                                                                                                                    |                                                                                                                                                
cATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTc                                                                                                                                                  <  2:192790/149‑1 (MQ=255)
  tttCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCaa                                                                                                                                                >  1:147373/1‑149 (MQ=255)
            ggAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                                                                      <  1:24979/149‑1 (MQ=255)
                    aagaagAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                                                                              <  1:106873/149‑1 (MQ=255)
                     agaagaAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgatg                                                                                                                             >  2:204622/1‑149 (MQ=255)
                      gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAAtgatga                                                                                                                            >  2:32495/1‑149 (MQ=255)
                       aagaaAACATAAAAAGAGAAATCGAAAACTTTAAGTTCTTTGCGCAATATAGCAATCTTAAAGATTTAATGAATTATAAAGATGGAGAAATATCATATAATCCAGAGTTGCCGAGTTATTCAGCTCAATATCAATTAATTAATgatgat                                                                                                                           >  2:128299/1‑149 (MQ=18)
                         gaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATta                                                                                                                         >  2:79189/1‑149 (MQ=255)
                         gaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATta                                                                                                                         <  1:67639/149‑1 (MQ=255)
                           aaaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTTTTTTTTATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtata                                                                                                                       >  2:188146/1‑149 (MQ=17)
                           aaaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtata                                                                                                                       <  1:47852/149‑1 (MQ=255)
                            aaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTTTTTTTTTTATCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTTCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATaa                                                                                                                      >  2:71151/1‑149 (MQ=14)
                            aaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATaa                                                                                                                      <  2:84857/149‑1 (MQ=255)
                                   aaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                                                                                                               <  2:115540/149‑1 (MQ=255)
                                     aaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACa                                                                                                             <  1:238802/149‑1 (MQ=255)
                                      aaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACaa                                                                                                            <  2:113532/149‑1 (MQ=255)
                                        gAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                                                                                                          >  1:52164/1‑149 (MQ=255)
                                                aaaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTaaaa                                                                                                  <  1:151863/149‑1 (MQ=255)
                                                     ttAAGTTTTTTGAGCAATATAGACATTTTAGTGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCCAAATAACCATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGatat                                                                                             <  2:105430/149‑1 (MQ=18)
                                                         gTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACa                                                                                         <  2:288310/149‑1 (MQ=255)
                                                          ttttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACAt                                                                                        <  1:232603/149‑1 (MQ=255)
                                                                gcgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaac                                                                                  <  1:289342/149‑1 (MQ=255)
                                                                 cgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAccgacg                                                                                 <  2:49524/149‑5 (MQ=14)
                                                                 cgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaaca                                                                                 <  2:147373/149‑1 (MQ=255)
                                                                      tataGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCCAAAGATATGATATACCGACGAGTaa                                                                            >  2:208535/1‑149 (MQ=11)
                                                                         aGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAgc                                                                         >  2:106975/1‑149 (MQ=11)
                                                                          gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGTTATATCATATAATCCTGAGTTGCCGAGATATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACATTTACATAAAAGATATGACATACCAATAAATAAAgcg                                                                        >  2:105489/1‑149 (MQ=18)
                                                                           cAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcgc                                                                       <  1:291678/149‑1 (MQ=255)
                                                                           cAATTCTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCAAGAGGTGCCGAGATATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATAACAACAAATAAAgcga                                                                       >  2:34284/1‑148 (MQ=255)
                                                                                                   taAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAg                                               >  1:42354/1‑149 (MQ=255)
                                                                                                      aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGa                                            >  2:31024/1‑149 (MQ=255)
                                                                                                        aTGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAt                                          >  1:163878/1‑149 (MQ=255)
                                                                                                        aTGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAt                                          >  1:143920/1‑149 (MQ=255)
                                                                                                                         aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGt                         >  2:296106/1‑149 (MQ=255)
                                                                                                                            ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTgg                      <  2:60294/149‑1 (MQ=255)
                                                                                                                                       gAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaaaaa           >  1:103049/1‑149 (MQ=255)
                                                                                                                                                aGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTa  <  1:137635/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTA  >  CP000730/111368‑111660

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