Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,552,031:1 +G 6.2% coding (1211/1734 nt) USA300HOU_RS13150 ← putative ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,552,0311.G 6.2% 340.9 / 16.5 81coding (1211/1734 nt)USA300HOU_RS13150putative ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (31/45);  new base G (3/2);  total (34/47)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

CCTGGATTGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGACATAATGTAG  >  USA300TCH1516_ALE/2551941‑2552121
                                                                                           |                                                                                          
ccTGGATTGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTAc                                                                                             >  1:1450245/1‑91 (MQ=255)
ccTGGATTGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTAc                                                                                             >  1:969867/1‑91 (MQ=255)
ccTGGATTGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTAc                                                                                             >  1:4126273/1‑91 (MQ=255)
   ggATTGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTa                                                                                         >  1:3336604/1‑91 (MQ=255)
    gATTGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTaa                                                                                        >  1:3730247/1‑91 (MQ=255)
      ttGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATg                                                                                      <  1:3813505/91‑1 (MQ=255)
      ttGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATg                                                                                      <  1:2952711/91‑1 (MQ=255)
      ttGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATg                                                                                      <  1:3813600/91‑1 (MQ=255)
       tGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGt                                                                                     <  1:839263/91‑1 (MQ=255)
       tGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAAGTACGTTAATGt                                                                                     <  1:1835847/91‑1 (MQ=255)
        gCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTc                                                                                    >  1:79310/1‑91 (MQ=255)
        gCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAAGTACGTTAATGTc                                                                                    <  1:1140292/91‑1 (MQ=255)
         ccATACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCa                                                                                   <  1:3268152/91‑1 (MQ=255)
          cAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCat                                                                                  <  1:453235/91‑1 (MQ=255)
            aaCAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCatat                                                                                <  1:4059708/91‑1 (MQ=255)
             aCAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATc                                                                               <  1:3579964/91‑1 (MQ=255)
                atatGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCttt                                                                            <  1:2830783/91‑1 (MQ=255)
                 tatGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTa                                                                           <  1:548379/91‑1 (MQ=255)
                  atGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTaa                                                                          >  1:3162543/1‑91 (MQ=255)
                  atGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTaa                                                                          <  1:2957951/91‑1 (MQ=255)
                  atGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTCTCGACATTAACGTACTTAATGTCATATCTTTaa                                                                          >  1:3681097/1‑91 (MQ=255)
                   tGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAat                                                                         >  1:3930688/1‑91 (MQ=255)
                   tGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAat                                                                         >  1:3664234/1‑91 (MQ=255)
                    gTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAata                                                                        >  1:1338931/1‑91 (MQ=255)
                    gTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAata                                                                        >  1:3399885/1‑91 (MQ=255)
                      tATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATc                                                                      >  1:1395962/1‑91 (MQ=255)
                      tATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATc                                                                      >  1:1395813/1‑91 (MQ=255)
                         tttCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAac                                                                   <  1:1984648/91‑1 (MQ=255)
                          ttCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAaca                                                                  <  1:1713815/91‑1 (MQ=255)
                           tCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAacac                                                                 <  1:1137866/91‑1 (MQ=255)
                           tCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAacac                                                                 <  1:3563014/91‑1 (MQ=255)
                           tCAATCGTATCATGCAATAAATACTCCTTTTGAAATACTGCACTAATTTTCGACATTAATGTACTTAATGTCATAGCTTTAATATCAacac                                                                 <  1:3743470/91‑1 (MQ=255)
                            cAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACAcc                                                                >  1:2709994/1‑91 (MQ=255)
                             aaTCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCa                                                               <  1:3573638/91‑1 (MQ=255)
                              aTCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCAt                                                              <  1:1929785/91‑1 (MQ=255)
                                cGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATAGCTTTAATATCAATACCATCg                                                            >  1:2068747/1‑91 (MQ=255)
                                 gTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGa                                                           <  1:2791733/91‑1 (MQ=255)
                                  tATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGAt                                                          <  1:3132988/91‑1 (MQ=255)
                                   aTCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATg                                                         >  1:3142443/1‑91 (MQ=255)
                                   aTCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATg                                                         <  1:1076345/91‑1 (MQ=255)
                                      aTTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGa                                                      <  1:3142511/91‑1 (MQ=255)
                                       ttAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGaa                                                     >  1:1495778/1‑91 (MQ=255)
                                        taaataaaTACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAat                                                    >  1:3458573/1‑91 (MQ=255)
                                         aaataaatACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAata                                                   <  1:1744259/91‑1 (MQ=255)
                                           ataaataCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATAtt                                                 <  1:1699798/91‑1 (MQ=255)
                                             aaataCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTc                                               <  1:2496161/91‑1 (MQ=255)
                                             aaataCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTc                                               <  1:443285/91‑1 (MQ=255)
                                             aaataCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTc                                               >  1:1317878/1‑91 (MQ=255)
                                             aaataCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTc                                               <  1:3758573/91‑1 (MQ=255)
                                                 acacCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCcatc                                           <  1:3245528/91‑1 (MQ=255)
                                                 acacCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCcatc                                           <  1:3979175/91‑1 (MQ=255)
                                                   acCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCcatcat                                         <  1:3140965/91‑1 (MQ=255)
                                                    ccTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCatcatc                                        <  1:1629977/91‑1 (MQ=255)
                                                    ccTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCatcatc                                        >  1:3498700/1‑91 (MQ=255)
                                                    ccTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCatcatc                                        >  1:3498779/1‑91 (MQ=255)
                                                      ttttGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGa                                      >  1:115065/1‑91 (MQ=255)
                                                       tttGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGat                                     <  1:3436909/91‑1 (MQ=255)
                                                        ttGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGata                                    <  1:736357/91‑1 (MQ=255)
                                                         tGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGatat                                   >  1:1097932/1‑91 (MQ=255)
                                                             aTACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATa                               <  1:308280/91‑1 (MQ=255)
                                                              tACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATaa                              >  1:3078569/1‑91 (MQ=255)
                                                               aCTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATaaa                             <  1:2328961/91‑1 (MQ=255)
                                                                   cACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAgcg                         <  1:521764/91‑1 (MQ=255)
                                                                     cTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCa                       <  1:680058/91‑1 (MQ=255)
                                                                       aaTTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAag                     >  1:2990927/1‑91 (MQ=255)
                                                                       aaTTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAag                     >  1:2988257/1‑91 (MQ=255)
                                                                         ttttCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAagag                   <  1:2049338/91‑1 (MQ=255)
                                                                          atatgACATTAAGTACGTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGt                  >  1:1421800/5‑91 (MQ=38)
                                                                           ttCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTa                 <  1:859142/91‑1 (MQ=255)
                                                                           ttCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTa                 <  1:286913/91‑1 (MQ=255)
                                                                            tCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTaa                >  1:1379002/1‑91 (MQ=255)
                                                                            tCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTaa                <  1:3034211/91‑1 (MQ=255)
                                                                            atgACATTAAGTACGTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTaa                >  1:2990835/3‑91 (MQ=255)
                                                                             cGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAg               >  1:1412549/1‑91 (MQ=255)
                                                                             cGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAg               >  1:3406551/1‑91 (MQ=255)
                                                                              gACATTAAGTACGTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATGATCGATATCATAAAAGCGCAAGAGTAAGt              >  1:2305122/1‑91 (MQ=255)
                                                                               aCATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTg             <  1:1482074/91‑1 (MQ=255)
                                                                               aCATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGCAAGTg             >  1:2827923/1‑91 (MQ=255)
                                                                                 aTTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGAc           >  1:3920919/1‑91 (MQ=255)
                                                                                   tAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGACAt         <  1:183151/91‑1 (MQ=255)
                                                                                   tAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGACAt         >  1:1521082/1‑91 (MQ=255)
                                                                                       gTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGACATAATg     >  1:3190352/1‑91 (MQ=255)
                                                                                         aCTTAATGTCATATCTTTAATATCAACCCCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGACATAATGTa   <  1:1481576/91‑1 (MQ=255)
                                                                                         aCTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGACATAATGTa   <  1:728081/91‑1 (MQ=255)
                                                                                          cTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGACATAAAGTAg  >  1:780948/1‑91 (MQ=255)
                                                                                           |                                                                                          
CCTGGATTGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGACATAATGTAG  >  USA300TCH1516_ALE/2551941‑2552121

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