Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,814 Δ1 bp 17.2% intergenic (+45/+63) USA300HOU_0025 → / ← USA300HOU_0026 hypothetical protein/IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073035,8140A.17.2% 142.7 / 30.2 58intergenic (+45/+63)USA300HOU_0025/USA300HOU_0026hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base A (23/24);  new base . (2/8);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

TGCAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTC  >  CP000730/35676‑35946
                                                                                                                                          |                                                                                                                                    
tGCAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGa                                                                                                                                    >  1:319277/1‑141 (MQ=255)
tGCAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGa                                                                                                                                    >  1:119152/1‑141 (MQ=255)
 gCAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGaa                                                                                                                                   <  2:102269/141‑1 (MQ=255)
   aaaTTATTCTGAATTTAAAGACGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAgg                                                                                                                                 <  2:351340/141‑1 (MQ=255)
    aaTTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGt                                                                                                                                >  1:283624/1‑141 (MQ=255)
     attattCTGAATTTAAAGATGTAGATGATCTCATTGACAACGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAGAAATAGAAGGtt                                                                                                                               <  1:264123/141‑1 (MQ=255)
      ttattCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTc                                                                                                                              <  1:300537/141‑1 (MQ=255)
           cTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTg                                                                                                                         >  1:283779/1‑141 (MQ=255)
                   aaaGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTaa                                                                                                                 <  1:85517/141‑1 (MQ=255)
                    aaGACGTAGACGACCCCATTGACAAAGGATTTAATGTAGAAGCTACTCAGTATAAAGTTTCTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTaaa                                                                                                                <  1:51857/141‑1 (MQ=255)
                     aGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTaaaa                                                                                                               >  2:281277/1‑141 (MQ=255)
                     aGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAAAAGTTTGTGTTGCAAGGTaaaa                                                                                                               >  1:156609/1‑141 (MQ=255)
                            gatgatATCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAACAAAtcta                                                                                                        >  2:76564/1‑141 (MQ=255)
                              tgatCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAga                                                                                                      >  2:346096/1‑140 (MQ=255)
                                    cATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCa                                                                                                <  1:123037/141‑1 (MQ=255)
                                         aCAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAAt                                                                                           <  2:283624/141‑1 (MQ=255)
                                             aGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTAt                                                                                       >  1:108939/1‑141 (MQ=255)
                                                  ttAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGt                                                                                  >  1:244234/1‑141 (MQ=255)
                                                     aTGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAg                                                                               >  2:288150/1‑141 (MQ=255)
                                                       gTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAgtg                                                                             >  1:95266/1‑141 (MQ=255)
                                                       gTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGAAAAGTAAAAACATATAGCTAACCACTAATTTATCATGTCAgtg                                                                             >  2:87334/1‑141 (MQ=255)
                                                            aGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGc                                                                        <  1:157007/141‑1 (MQ=255)
                                                                aTTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGCTCTGATAAAAATAGAAGGTTCTGTTGCAACGAAAGAGAAGAAAGCTAACCAATAATTTATCATGTCCGTGTTCGCTTaa                                                                    >  2:78070/1‑141 (MQ=255)
                                                                  tCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACt                                                                  >  1:10051/1‑141 (MQ=255)
                                                                        taAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAg                                                            <  1:13430/141‑1 (MQ=255)
                                                                           aGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCat                                                         <  2:231336/141‑1 (MQ=255)
                                                                                tatgaaaaaaaTGATTACTTGAACAAGAGCTAAAATAAAATTGTAATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGc                                                    <  2:104839/135‑1 (MQ=255)
                                                                                     gAAAAATGATTACTTGAACAAGAACTAAAAAAAAATTGTGAACTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAAtt                                               <  2:55627/141‑1 (MQ=255)
                                                                                            gATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGc                                        >  2:75650/1‑141 (MQ=255)
                                                                                               tACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATg                                     <  2:302780/141‑1 (MQ=255)
                                                                                                 cTTGAACAAGATTTAAAATAAAATTTTGATCTAATAAGAATAGAAGGTTCTGTTGCAAAGTAAAATAATATAGCTTACCTCTAATTTCTCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGc                                   <  2:158653/141‑1 (MQ=17)
                                                                                                 cTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGc                                   <  2:208903/141‑1 (MQ=255)
                                                                                                 cTTAAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCGGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGc                                   <  2:243140/141‑1 (MQ=255)
                                                                                                      acaacaGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGaaaa                              >  2:226083/1‑141 (MQ=255)
                                                                                                      aCAAGAGCTAAAATAAAATTGTCATCTAATAAAAATAGACAGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATCTATCATTGGAGGGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGaaaa                              <  1:311465/141‑1 (MQ=17)
                                                                                                          gagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATTCTAATTTCGTGGCATGGCGAAAATCCg                          >  2:310540/1‑141 (MQ=255)
                                                                                                            gCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTa                        >  2:213002/1‑141 (MQ=255)
                                                                                                             cTAAAATAAAATTGTGATCTAATAAAAATAAAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTCATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATTGCGAAAATCCGTAg                       >  2:178940/1‑141 (MQ=255)
                                                                                                              taaaataaaaTTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGa                      >  2:187007/1‑141 (MQ=255)
                                                                                                                  acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                >  1:261188/15‑141 (MQ=11)
                                                                                                                   caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag               <  1:205985/128‑1 (MQ=11)
                                                                                                                      aagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAACCTTTAGCTAATCACTAATATATCATGTTAGTGTTCGTTTAACTTGCTAGCATGATGATCATTTCGGGGCATGGCGAAAATCCGTAGATCTGAAgtga            <  2:115899/131‑1 (MQ=2)
                                                                                                                       agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc           <  2:338356/132‑1 (MQ=11)
                                                                                                                         tgtgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc           <  2:108939/141‑1 (MQ=255)
                                                                                                                           tgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt         <  2:283779/141‑1 (MQ=255)
                                                                                                                           tgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt         >  1:60335/1‑141 (MQ=255)
                                                                                                                           tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAGTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc       <  2:342107/136‑1 (MQ=11)
                                                                                                                           tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc       <  1:81528/136‑1 (MQ=11)
                                                                                                                           tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc       <  2:213429/136‑1 (MQ=11)
                                                                                                                            gATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTCGATCTGTAAAGACCTg        >  2:10712/1‑141 (MQ=255)
                                                                                                                            attttaataaAAATGATGGTTCTGTTGCAAAGAAAAAAATATAGCTAATCACTAATATATCATGTCAGTGGACGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAAACCGTAGATATGAAGAGACCTGCg      >  2:365966/5‑141 (MQ=11)
                                                                                                                             ttttaataaAAATGCTGGTTCTGTTGCAAAGTAAAAAATATAGCCAACCACTAATATATCCTGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTGCGTTGCATGGCGAAAATCCGGAGATCTGAAGAGACCTGCgg     <  1:321997/138‑1 (MQ=11)
                                                                                                                             aTCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc       <  2:186990/141‑1 (MQ=255)
                                                                                                                             aTCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc       <  2:10051/141‑1 (MQ=255)
                                                                                                                              tCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg      <  1:213002/141‑1 (MQ=255)
                                                                                                                               ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt   <  2:340933/140‑1 (MQ=11)
                                                                                                                                taagaaAAATTGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATCTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt    <  2:61688/141‑1 (MQ=255)
                                                                                                                                  ataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc  <  2:156609/141‑1 (MQ=38)
                                                                                                                                          |                                                                                                                                    
TGCAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTC  >  CP000730/35676‑35946

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: