Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 111,560 T→C 9.8% R210R (CGT→CGC USA300HOU_0110 → possible staphylococcal tandem lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730111,5600TC 9.8% 122.4 / 3.6 52R210R (CGT→CGCUSA300HOU_0110possible staphylococcal tandem lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (30/15);  new base C (3/2);  total (34/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATATTTACTTTACTGA  >  CP000730/111423‑111700
                                                                                                                                          |                                                                                                                                             
aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa                                                                                                                                            <  1:219038/141‑1 (MQ=255)
aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTaaa                                                                                                                                            >  2:307410/1‑141 (MQ=255)
   ttttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGa                                                                                                                                         >  1:457985/1‑141 (MQ=255)
      ttttttttttttttcAATTTTAAAGAGGAAATGAATTTTACAGATGGAGATATATCATATAATCCAGAGGTGCCGCGTTATTCAGCTTATTTTCTTTTATTTATTGATGTTACTTCGGTAAAACAATTACGATAAAGacat                                                                                                                                      >  2:581081/15‑141 (MQ=9)
       ttGCGCAATATAGCACTTTTAAAGATCTAATGAACTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCTGCTAACTATCAATAAACTAATGATGATTATAATGTAATACAATTACGCAAAAGATATg                                                                                                                                     <  1:372149/141‑1 (MQ=17)
         gcgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGaa                                                                                                                                   >  1:482008/1‑140 (MQ=255)
            ttttttttcAATATTAAAGAGTTAATGAATTTTAAAGATGGAGATAGATCATATAATCCAGAGGTGCCGAGTTATTTTTCTTTTTATCTATTAACTAATGATGATTATAATATAAAACAATTACGTAAAAGGTATGACATa                                                                                                                                >  2:374068/9‑141 (MQ=17)
            cattataGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACATTTCCGTAACAGATATTACAtt                                                                                                                                >  2:180524/4‑140 (MQ=18)
            cAATATAGCAATTTTAAAGCTTTAATGCATTATAAAGATAGAGATATATCATATAATCCGGATGTGCCTATTTATTCGACTAACTATCAATTAACTAATGATGATTGTCATCTAAAAAAATTAAGCAACatatattatatc                                                                                                                                >  2:174081/1‑135 (MQ=9)
            cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATa                                                                                                                                >  2:58867/1‑141 (MQ=255)
              atatAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAAATATGATATAcc                                                                                                                             >  2:409294/1‑141 (MQ=14)
                   gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATGCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATTATGATGATAATGTAAAACAATTACGTAAAAGATATGACATACcaacaa                                                                                                                        >  1:318573/1‑141 (MQ=255)
                          tAAAGATTTAATGAATTATAAAGATGGATATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgc                                                                                                                 <  1:348577/141‑1 (MQ=255)
                            aaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcgc                                                                                                               >  1:40986/1‑141 (MQ=255)
                                   aaTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCt                                                                                                        <  2:406040/141‑1 (MQ=255)
                                    aTGAATCATCAAGATGGAGATATATCATATAATCCACAGGTGCCAAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtg                                                                                                       <  1:14333/141‑1 (MQ=255)
                                          tataAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTCTGCACTCGAATAGCAAGTAACTAATGATGATTATAATGTCAAACAATTACGTCAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgttgtgg                                                                                                 >  2:68837/1‑139 (MQ=17)
                                               aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAgg                                                                                            <  1:233580/141‑1 (MQ=255)
                                                    gagaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAg                                                                                       >  2:577454/1‑141 (MQ=255)
                                                     agaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAgg                                                                                      >  1:502760/1‑141 (MQ=255)
                                                          tatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAtt                                                                                 <  1:273184/141‑1 (MQ=255)
                                                           atCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTACATCATGCATATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAttt                                                                                <  1:416476/141‑1 (MQ=17)
                                                               tataATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTaaaa                                                                            >  2:375369/1‑141 (MQ=255)
                                                                  aaTTCAGAGGGGGCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACACTTCTGGCGAAGAGATGAATGCTGAACGGAGAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGt                                                                         <  1:293812/141‑1 (MQ=12)
                                                                    tCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtc                                                                       <  1:527655/141‑1 (MQ=255)
                                                                     ccAGAGGAGCCGAGTTATTCAGCTCAACATCAAATAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtca                                                                      <  2:399062/141‑1 (MQ=255)
                                                                           gTGCCGAGTTATTCAGCTCAATATCAATCAACTAATGATGATTATACTGAAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGtt                                                                >  2:128778/1‑141 (MQ=255)
                                                                            tGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTGTAATGTAAAACAATTACGTAAAAGATATGACATACCGACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTg                                                               >  1:517505/1‑141 (MQ=255)
                                                                              ccGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAACGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGa                                                             >  1:224273/1‑141 (MQ=255)
                                                                                     aTTCAGCTCAATATCAAGTAACTAATGATAATTATAATGTAAAACAATTACGTAGTAGATATGACATACCAACAAGTAAAGCGCCGAAGCTGTTGGGGGGCGATCCAGGGAACGTAAGATGTTGATCCGTGGGATAGaaaa                                                      >  2:578794/1‑141 (MQ=12)
                                                                                       tCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCAACAGTTGGATATaaaaaa                                                    >  1:530392/1‑141 (MQ=255)
                                                                                               aaTATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAAt                                             >  2:560185/1‑141 (MQ=255)
                                                                                                atatCAATTAACTAATGATGATTATACTGTAAAACTATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCGGTTGATGCAAGGTACAGGGAAGTTAAAAGGTACATCCGTTGGATATAAAAAAAGTgcagt                                            >  1:296040/1‑137 (MQ=14)
                                                                                                atatCAAGGAGCGAATGATGATTATAATGTAAAACAAGGACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTGGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAAtt                                            <  1:580929/141‑1 (MQ=17)
                                                                                                 tatCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAAttt                                           >  2:566461/1‑141 (MQ=255)
                                                                                                   acAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTAc                                         <  1:357677/140‑1 (MQ=255)
                                                                                                    aaaTTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTatt                                        <  2:98960/140‑3 (MQ=255)
                                                                                                          aCTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATATCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTa                                  <  1:568601/141‑1 (MQ=255)
                                                                                                          aCTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTa                                  <  1:192930/141‑1 (MQ=255)
                                                                                                              atgatgATTGTAATGTAAAACAATTACGTAAAAGATATGACATACCGACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAgaga                              <  2:517505/141‑1 (MQ=255)
                                                                                                               tgatgaTTATAATGTAAAACAATTACGTAACAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGGTCATCAGTTGGATATAAAAAAATTGAATTTAATTTTTTagcaaa                             >  2:461084/1‑137 (MQ=255)
                                                                                                               tgatgaTTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGaa                             >  1:27345/1‑141 (MQ=255)
                                                                                                               tgatgaTTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGaa                             >  2:105028/1‑141 (MQ=255)
                                                                                                               tgatgaTTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGaa                             >  1:27339/1‑141 (MQ=255)
                                                                                                               tgatgaTTATAATGTAAAAAAATTACGTAAAAGATATGACAAACAAACAAATAACGCTCAGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGATGGATCTCAACACATTGAATGTACTTTTTTTgagag                             >  1:27344/1‑140 (MQ=14)
                                                                                                                gatgatTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAAt                            >  1:35790/1‑141 (MQ=255)
                                                                                                                       ataATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGGATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAGTTTACTGTGTTAGAGGATAAAGATg                     >  1:354689/1‑141 (MQ=17)
                                                                                                                         aaTGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAATTAAAAGGGTCATCAGTTGGAAATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGaa                   >  2:168727/1‑141 (MQ=255)
                                                                                                                               aaaCAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATAtt             >  1:577612/1‑141 (MQ=255)
                                                                                                                                  cAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAACAtttac          >  2:91469/1‑141 (MQ=255)
                                                                                                                                        cGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATATttactttact    >  1:26721/1‑141 (MQ=255)
                                                                                                                                        cGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATATttactttact    <  2:353839/141‑1 (MQ=255)
                                                                                                                                          tAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTGGGGAATTAAAATGGAAATATTTACTTTACTgg  >  1:220168/1‑140 (MQ=255)
                                                                                                                                          tAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATATTTACTTTACTGa  >  1:220638/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                             
AAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATATTTACTTTACTGA  >  CP000730/111423‑111700

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: