Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 110,660 T→A 9.7% I173I (ATT→ATA)  USA300HOU_0109 → staphylococcal tandem lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730110,6600TA 9.7% 101.4 / 3.4 40I173I (ATT→ATA) USA300HOU_0109staphylococcal tandem lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (13/22);  new base A (2/2);  total (15/25)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

CGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGG  >  CP000730/110534‑110795
                                                                                                                              |                                                                                                                                       
tgATTATAATAGCATTCCAACAAAAGAAAATAATGATGAAAAAATAAAAAAAGAAATTGAACACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGCAGGTATTTCATATAATCCaga                                                                                                                           <  2:155878/140‑1 (MQ=17)
          aTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGt                                                                                                                 <  1:164736/141‑1 (MQ=255)
             aTTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATACTCCAGAGGTGCCgcataac                                                                                                              <  2:135663/141‑7 (MQ=255)
             aTTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTAt                                                                                                              >  2:17485/1‑141 (MQ=255)
             aTTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTAt                                                                                                              >  2:17486/1‑141 (MQ=255)
               tCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTc                                                                                                            >  1:153809/1‑141 (MQ=255)
                  aacaaAAGACATTAAAGACGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAAAACTATCTTATAATCCAGAGGTGCCGAGTTATTCGGc                                                                                                         <  1:10941/141‑1 (MQ=17)
                  aacaaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCATAGGTGCCGaataactcggc                                                                                                         <  1:168071/141‑11 (MQ=255)
                   ctaaagGAAATTAAACCTGAAAAAATAAAAAAAGAAATCTAAATCTCTAAGTGCGTTGTTCAGTATGGCGTCTTTAAACATTTTAAAAATTATAAAGACCGAGATATTTCATTTAATCCAGAGGTGCCGAGTTATTCGGCt                                                                                                        <  2:264791/139‑1 (MQ=17)
                        aGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGATAAata                                                                                                   >  2:23789/1‑141 (MQ=255)
                        aGAAATTAAAGATGAAAAAATAAAAAAAAAAATCGAAAACTATAAGTACTATGTTAAATATGACAACTTTAAAAATATGAAAAAATATAAAGACGGAGAAATTTAATATAATCCAGAGGTGCGGAGTGATTCGACTAAata                                                                                                   >  2:180150/1‑141 (MQ=255)
                            aTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGATGTGCCGAGTGATTCGGCTAAATATCaa                                                                                               >  1:250810/1‑141 (MQ=255)
                            aTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGATAATTATCaa                                                                                               >  2:254904/1‑141 (MQ=255)
                               aaaGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                                                                            <  1:25479/141‑1 (MQ=255)
                                aaGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                           >  1:227201/1‑141 (MQ=255)
                                 aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                          <  2:76750/141‑1 (MQ=255)
                                   aTGAAAAAATAGAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAAGTTGAAAAATTATAAAGACGGAGATATGTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                        <  1:195669/141‑1 (MQ=255)
                                       aaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGCCTTTAACAATTTGAAACATTATAAAGACGGAGATAGTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                                    >  2:31019/1‑141 (MQ=255)
                                       aaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                                    <  1:269519/141‑1 (MQ=255)
                                              aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                             <  1:23519/141‑1 (MQ=255)
                                                aaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtatat                                                                           >  1:143540/1‑140 (MQ=255)
                                                    gAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                                                                       <  1:13004/141‑1 (MQ=255)
                                                     ataTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                                      <  1:72803/139‑1 (MQ=255)
                                                                tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGATGCAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                           <  2:153484/141‑1 (MQ=255)
                                                                tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTAGTCGGCTAAATATCAATTAACTAACGATGATTATAATGTAAAACAATTACGc                                                           <  2:199014/141‑1 (MQ=255)
                                                                                 atatGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCCAAAGATATGAGATAcg                                          >  2:131995/1‑140 (MQ=255)
                                                                                  tatGGTGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcg                                         <  1:10879/141‑1 (MQ=255)
                                                                                  tatGGCGACTTTAAAAATTTGAAAAATTATAAAGCCGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATTAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcg                                         >  2:162975/1‑141 (MQ=255)
                                                                                      gCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacga                                     <  2:227201/141‑1 (MQ=255)
                                                                                          cTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                                 <  2:118998/141‑1 (MQ=255)
                                                                                               aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTc                            <  1:273080/141‑1 (MQ=255)
                                                                                                    ttGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTAATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATAGACCGACGAGTAAAGCTCCaaaa                       >  2:154303/1‑140 (MQ=255)
                                                                                                        aaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt                   <  1:59259/141‑1 (MQ=255)
                                                                                                        aaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt                   <  1:91063/141‑1 (MQ=255)
                                                                                                           aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGTTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATCATCTACCGACGAGTAACGCTCCAAAGTTATTGt                <  2:185630/141‑1 (MQ=17)
                                                                                                           aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGt                >  1:31798/1‑141 (MQ=255)
                                                                                                             tataACAACGGAGATATTTCATATAAGCCAGGGTGGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGTGTAAAGCTCCAAAGTTATTGTTa              <  2:216399/141‑1 (MQ=18)
                                                                                                             tataAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTa              <  1:17486/141‑1 (MQ=255)
                                                                                                             tataAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTa              <  1:17485/141‑1 (MQ=255)
                                                                                                                       gagaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGATATTGTTAAAAGGTTCAg    >  1:245672/1‑141 (MQ=255)
                                                                                                                         aatatatCATATAAGCCAGAGGTGCCGAATTATTCGGCTAAATATCAATTAACTAATTATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAggg  <  1:196364/140‑1 (MQ=255)
                                                                                                                              |                                                                                                                                       
CGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGG  >  CP000730/110534‑110795

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