Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 110,651 C→T 11.6% D170D (GAC→GAT USA300HOU_0109 → staphylococcal tandem lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730110,6510CT11.6% 144.6 / 7.1 58D170D (GAC→GATUSA300HOU_0109staphylococcal tandem lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (23/28);  new base T (2/5);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

AAAAGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGG  >  CP000730/110511‑110788
                                                                                                                                            |                                                                                                                                         
aaaaGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTGAGTTCTTTGTTAAATATGGCGACTGTAAAAAATTGAAAAATAATAAAGAc                                                                                                                                           >  1:284767/1‑141 (MQ=18)
 aaaGATAACCGGTTAAAATGGTCGATAATAAAAACATTGCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACg                                                                                                                                          <  2:72112/141‑1 (MQ=255)
        tCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGatat                                                                                                                                   <  1:642055/141‑1 (MQ=255)
             taaaaaTGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCat                                                                                                                              <  2:411542/139‑1 (MQ=255)
                    ggTCGATAATAAAATCCTTCAAGGAGTAGAAGTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAAACGGAGATATTTCATATAAtgc                                                                                                                       >  2:143927/1‑139 (MQ=17)
                       cGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCaga                                                                                                                    <  2:100273/141‑1 (MQ=255)
                       cGATAATAAAATCATGCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCaga                                                                                                                    >  1:494552/1‑141 (MQ=255)
                         ataataAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGTGATATTTCATATAATCCAGAgg                                                                                                                  >  1:249549/1‑141 (MQ=255)
                             taaAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                                                              >  1:434445/1‑141 (MQ=255)
                             aaaaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAACTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                                                              <  1:261769/140‑1 (MQ=255)
                               aaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAATAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGa                                                                                                            <  1:471607/141‑1 (MQ=255)
                               aaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGa                                                                                                            >  1:547439/1‑141 (MQ=255)
                               aaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGa                                                                                                            <  1:266579/141‑1 (MQ=255)
                               aaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTACCGa                                                                                                            >  1:326498/1‑141 (MQ=255)
                               aaaCCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCGTATAATCCAGAGGTGCCGa                                                                                                            <  1:353909/141‑1 (MQ=255)
                                      tCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCATTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTCttg                                                                                                     >  1:29738/1‑140 (MQ=255)
                                       ccAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCg                                                                                                    >  2:160864/1‑141 (MQ=255)
                                         aacaaAAGAAATTAAAGAGGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTGTGTTCAATCTGGAGACGTTAAAAATTTGAAAAATTATAAAGACTGAGATATTTCATATAATCCATAGGTGCCGAGTTATTCGGc                                                                                                  >  2:370967/1‑141 (MQ=17)
                                         aacaaAAAAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGc                                                                                                  <  2:585252/141‑1 (MQ=255)
                                              aaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                                                             >  2:425814/1‑141 (MQ=255)
                                              aaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                                                             >  2:543355/1‑141 (MQ=255)
                                              aaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                                                             <  2:331842/141‑1 (MQ=255)
                                              aaGAAATTAAAGATGAAAAAAAAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTACAAGATGGGGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                                                             <  2:401760/141‑1 (MQ=255)
                                               aGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGCGATATTTCACATAATCCAGAGGTGCCGAGTTATTCGGCTAAatt                                                                                            >  1:355628/1‑140 (MQ=255)
                                               aGAAATTAAAGAATAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                                                            <  2:487829/141‑1 (MQ=255)
                                                 aaaTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatata                                                                                          >  2:360582/1‑140 (MQ=255)
                                                    ttAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAAACGGAGATATTTCATATAATCCAGAGGTGACGAGTTATTAGGCTAAATATCAAt                                                                                       >  2:304346/1‑141 (MQ=255)
                                                         gATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                  <  1:509236/141‑1 (MQ=255)
                                                               aaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                            <  2:237911/141‑1 (MQ=255)
                                                               aaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                            <  1:140110/141‑1 (MQ=255)
                                                                     aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTCGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCAGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                      >  2:108035/1‑141 (MQ=255)
                                                                     aaaaaaGAAATCGAAAACTTTAAGTGCTGTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATCATtat                                                                      <  2:577117/141‑1 (MQ=255)
                                                                      aaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAATTAATGATGATtata                                                                     >  1:549574/1‑141 (MQ=255)
                                                                      aaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCACAGGTGCCGAGTTATCAGGCTACACAGCGACTGACGAGTGATGATtata                                                                     >  1:12938/1‑141 (MQ=17)
                                                                        aaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAAt                                                                   <  2:114113/141‑1 (MQ=255)
                                                                        aaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAAt                                                                   <  1:405086/141‑1 (MQ=255)
                                                                            aaaTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                               <  2:502173/141‑1 (MQ=255)
                                                                                  aaaaCTTTAAGTCCTTTGTTCAATGTCGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                         <  2:92939/141‑1 (MQ=255)
                                                                                   aaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                                                        <  2:344812/141‑1 (MQ=255)
                                                                                       tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                    <  1:238334/141‑1 (MQ=255)
                                                                                       tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATACTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                    <  2:168001/141‑1 (MQ=255)
                                                                                                   gTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatgat                                        <  2:194548/141‑1 (MQ=255)
                                                                                                   gTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatgat                                        <  2:194297/141‑1 (MQ=255)
                                                                                                                 cTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                          >  2:229484/1‑141 (MQ=255)
                                                                                                                 cTTTAAAAATTTGAAAAAAAATCAGGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAGACAATTACGCAAAAGATATGATATACCGACGAGTaa                          >  1:322945/1‑141 (MQ=18)
                                                                                                                  tttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaaa                         >  2:359919/1‑141 (MQ=255)
                                                                                                                    taaagatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTATTGTAAGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGc                       >  1:654563/14‑141 (MQ=11)
                                                                                                                     aaagatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCt                      >  2:242991/13‑141 (MQ=11)
                                                                                                                       aaaTTTGACAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                    <  2:255132/141‑1 (MQ=255)
                                                                                                                             gAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGtt              >  2:50196/1‑141 (MQ=255)
                                                                                                                             gAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGtt              >  1:61476/1‑141 (MQ=255)
                                                                                                                               tgaaTTATAAAGATGGACATATACCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACCATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt            <  2:169282/139‑1 (MQ=17)
                                                                                                                                  aTTATCAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGt         <  2:568288/141‑1 (MQ=255)
                                                                                                                                  aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGt         <  1:415612/141‑1 (MQ=255)
                                                                                                                                  aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGt         <  1:522279/141‑1 (MQ=255)
                                                                                                                                    tataAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTa       <  2:310084/141‑1 (MQ=255)
                                                                                                                                     ataAAGCTGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTATGTCAATTACGCACCTGTCACTGTATACCGACGAGTAAAGCTCCAAAGTTATTGTTaa      <  1:657503/141‑1 (MQ=14)
                                                                                                                                       aaaGACGGAGCTATTTCATATAATCCAGAGGTGCGGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATCATGTACAACAATTACGCAAAAGCTATGCTATACCGACGAGGAACGCTCCACAGATATTGCTgaaa    >  2:580345/1‑141 (MQ=17)
                                                                                                                                         agagGGAGATATCTCAGAGAATCCAGAGGTTCCGAATTAATCGGCAAAAAATCAATTAACTAATGATGATTATAATGTAAAACAATTACGAAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAgg  <  2:233752/141‑1 (MQ=17)
                                                                                                                                         acaCGCAGATATTCCATATAAACCAGAGGTCCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAgg  <  1:210954/139‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                         
AAAAGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGG  >  CP000730/110511‑110788

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: