Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 111,560 T→C 18.6% 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,5600TC18.6% 58.0 / 6.9 32R210R (CGT→CGCUSA300HOU_0110possible staphylococcal tandem lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (14/12);  new base C (3/3);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

AGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATATTTACTTTAC  >  CP000730/111424‑111697
                                                                                                                                        |                                                                                                                                         
aGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATAATAAAGTTGGATATATATCATATAATCCAGATGTGCCGTGTTATTCGTCTTAATATCAATTAACTAATGATGATTATAATTGTAAACAATTAAGCaaaa                                                                                                                                       >  2:172315/1‑141 (MQ=17)
         cgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGata                                                                                                                              <  1:239978/141‑3 (MQ=255)
          tttttttaGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCTGAGGTTCCGAGTTATTCAGCTTAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGAAAt                                                                                                                             >  2:249710/7‑141 (MQ=18)
          ccAATACAGCAATTTTAAAGATTGAATGAATTATAAAGATGGAGATATAGCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACAt                                                                                                                             <  1:238420/140‑1 (MQ=255)
           tattttttcAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGAtata                                                                                                                            >  2:35195/9‑141 (MQ=14)
           gaagttaGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGAtata                                                                                                                            <  1:112888/136‑1 (MQ=14)
           cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATa                                                                                                                            >  1:48700/1‑141 (MQ=255)
           cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGATGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACACTTACGCAAAAGATATGAtata                                                                                                                            >  2:162939/1‑141 (MQ=17)
                   cAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAccgacgag                                                                                                                    <  1:272075/141‑7 (MQ=11)
                     aTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATa                                                                                                                  <  1:235006/141‑1 (MQ=255)
                          aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGAGATACCAACtaataaggtg                                                                                                             >  1:170815/1‑137 (MQ=255)
                          aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcg                                                                                                             >  1:52846/1‑141 (MQ=255)
                           aaGATTTAATGAATTATAAAGAAGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgagc                                                                                                            >  2:158638/1‑141 (MQ=255)
                               tttAATGAAATATAAAGATAGAGATACATCATATAATCGTAAGCGGCTGATTTATACAGCTCCATATCAATTAACTAGTGATGATTATAATGTAAAACAATTACGTAAAATATATGACATCCCAACGAATATAGCGCCgca                                                                                                        <  2:31356/141‑3 (MQ=18)
                                               gATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGt                                                                                        >  2:264193/1‑141 (MQ=255)
                                                      atatatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGa                                                                                 >  2:86735/1‑141 (MQ=255)
                                                                 aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATCAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGt                                                                      <  2:215687/141‑1 (MQ=255)
                                                                  aTCCAGAGGTGCCGAGTTATTCAGCGCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGtt                                                                     <  1:104066/141‑1 (MQ=255)
                                                                   accAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtc                                                                    <  1:187061/140‑1 (MQ=255)
                                                                    ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtca                                                                   >  1:219207/1‑141 (MQ=255)
                                                                    ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtca                                                                   <  1:58270/141‑1 (MQ=255)
                                                                    ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtca                                                                   <  2:280794/141‑1 (MQ=255)
                                                                          gTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGAAAAAAAGTTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGtt                                                             >  1:95522/1‑141 (MQ=255)
                                                                                ggTTATTCAGCTCAATATCAATTAACTAATGATGATTATATTGTAGTACGATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGatat                                                       <  2:201037/140‑1 (MQ=255)
                                                                                                 tCAATTAACTAATGATGATTATAATGTAAAACAATTACATAAAAGATATGACATACCAACAAATAAAGAGCCGAAGCTGTAGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTGGGATATAAAAAAATTGAATTTAc                                      >  2:36993/1‑141 (MQ=255)
                                                                                                  cAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACt                                     <  1:50210/141‑1 (MQ=255)
                                                                                                         cTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTag                              <  1:169175/141‑1 (MQ=255)
                                                                                                         cTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACAAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGCATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTag                              >  2:94145/1‑141 (MQ=255)
                                                                                                               atgatTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATa                        >  2:51499/1‑141 (MQ=255)
                                                                                                                  atTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAGAAAATTGAATTTACTTTTTTAGAGAATaaaa                     >  2:3280/1‑141 (MQ=255)
                                                                                                                   ttATAATGTAAAACAATTACGTAAAAGATATGACCTACCAACAAATAAAGCGCCGAAGCTGTTGTAGCAAGGTACAGGGAATTTAAAAGGTTCAACAGTTGGATAGAAAAAAATTAAATTAACGTTTTTCGAGCAtgcaaa                    >  2:150215/1‑136 (MQ=14)
                                                                                                                                     aCGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATAtttactttac  <  2:95522/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                         
AGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATATTTACTTTAC  >  CP000730/111424‑111697

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