Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,386,433 A→G 8.2% I9I (ATT→ATC)  USA300HOU_2245 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,386,4330AG 8.2% 167.8 / 4.2 58I9I (ATT→ATC) USA300HOU_2245hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/22);  new base G (2/3);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTAC  >  CP000730/2386294‑2386566
                                                                                                                                           |                                                                                                                                     
cTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATaatgaa                                                                                                                                      <  2:573085/141‑1 (MQ=255)
    aTCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAatgat                                                                                                                                  >  2:395774/1‑141 (MQ=255)
     tCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGGatgatg                                                                                                                                 >  1:82195/1‑141 (MQ=255)
     tCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGGatgatg                                                                                                                                 >  1:197666/1‑141 (MQ=255)
          tagCTCAATATAATTTTCATTCTTATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATc                                                                                                                            <  1:403157/141‑1 (MQ=255)
           agCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCa                                                                                                                           <  1:763439/141‑1 (MQ=255)
               cAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAAtt                                                                                                                       <  1:759815/141‑1 (MQ=255)
                aaTATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAAttt                                                                                                                      <  1:444621/141‑1 (MQ=255)
                  tgtgaTTTTCAGTGACATAACCTAATCTATTGAATGTTCCATATTTAGGATGAGACCCACCTGAATGTAATTTTATAACATATCATGCAAAACTAAACCGATCTTGCTGATCTATATAATGGATGATGTGATAAAATTTCa                                                                                                                    <  2:378475/137‑1 (MQ=255)
                  tataATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCa                                                                                                                    <  2:483102/141‑1 (MQ=255)
                   ataATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCaa                                                                                                                   >  1:468387/1‑141 (MQ=255)
                      aTTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAatat                                                                                                                >  1:190602/1‑141 (MQ=255)
                       ttttCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATc                                                                                                               >  2:784697/1‑141 (MQ=255)
                          tCATTGATATCACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCAtt                                                                                                            <  2:376435/141‑1 (MQ=255)
                           cATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTa                                                                                                           >  1:490710/1‑141 (MQ=255)
                              tGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTa                                                                                                        >  2:385278/1‑141 (MQ=255)
                                      ccTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTa                                                                                                >  1:99442/1‑141 (MQ=255)
                                      ccTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTa                                                                                                >  1:432971/1‑141 (MQ=255)
                                      ccTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTa                                                                                                >  1:757933/1‑141 (MQ=255)
                                            ttATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACa                                                                                          <  2:84724/141‑1 (MQ=255)
                                                tGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTa                                                                                      >  1:346790/1‑141 (MQ=255)
                                                 gAATGTTCCATATTTATGATGATACCCACCTGAATGTAGTTTTATAACATCGCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTaa                                                                                     >  2:200713/1‑141 (MQ=255)
                                                   aTGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAAtt                                                                                   >  1:500637/1‑141 (MQ=255)
                                                     gTTCCATATTTATGATGATACCCGCCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCGATTTACATGTAATTAc                                                                                 <  1:556017/141‑1 (MQ=255)
                                                        ccATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCt                                                                              >  1:673665/1‑141 (MQ=255)
                                                        ccATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCACTAGTATCCCCCTATTTACATGTAATTACGCt                                                                              >  2:314345/1‑141 (MQ=255)
                                                        ccATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATCTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCt                                                                              >  1:286785/1‑141 (MQ=255)
                                                         cATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCtt                                                                             >  1:730059/1‑141 (MQ=255)
                                                             tttATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTAttt                                                                         <  1:523254/141‑1 (MQ=255)
                                                              ttATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTAtttt                                                                        <  1:569464/141‑1 (MQ=255)
                                                                 tgatgaTAGCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCGAACTGATCTATATAATGAATGATGCGATCAAATTTCAATATCATTAGTAGCCCCCTATTTACATGTAATTACGCTTATTTTAAc                                                                     <  2:223740/141‑1 (MQ=255)
                                                                 tgatgaTACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAAc                                                                     >  2:22688/1‑141 (MQ=255)
                                                                   atgatACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACaa                                                                   <  1:756735/141‑1 (MQ=255)
                                                                   atgatACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACaa                                                                   >  1:166440/1‑141 (MQ=255)
                                                                      atACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGt                                                                >  1:644158/1‑141 (MQ=255)
                                                                       tACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGGATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGta                                                               <  2:197666/141‑1 (MQ=255)
                                                                       tACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGGATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGta                                                               <  2:82195/141‑1 (MQ=255)
                                                                             ccTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAatta                                                         >  1:64643/1‑141 (MQ=255)
                                                                                    gTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGc                                                  <  1:133303/141‑1 (MQ=255)
                                                                                     tAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGct                                                 <  2:103702/141‑1 (MQ=255)
                                                                                     tAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGct                                                 <  1:219180/141‑1 (MQ=255)
                                                                                     tAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGct                                                 >  2:104480/1‑141 (MQ=255)
                                                                                          ttATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCTCTtaa                                            >  2:723611/1‑141 (MQ=255)
                                                                                               aCATCTCCTGGAAAACTTAACAAATCTAACTGTTCTATATTTTGAATGATGTGATCTAATTTTAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCGCTTGATAAtt                                       <  2:116715/141‑1 (MQ=255)
                                                                                                       tGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTA‑TTTTGCTCTTAATAATTATATATTGa                              <  1:201328/141‑1 (MQ=255)
                                                                                                        ggAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCTCTTAATAATTATATATTGa                              <  2:767890/141‑1 (MQ=255)
                                                                                                         gAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCTCTTAATAATTATATATTGAt                             >  2:487824/1‑141 (MQ=255)
                                                                                                             aCTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCt                         >  1:541092/1‑141 (MQ=255)
                                                                                                              cTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTa                        <  1:39250/141‑1 (MQ=255)
                                                                                                                aaaCCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAAt                      >  1:570805/1‑141 (MQ=255)
                                                                                                                aaaCCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATCAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAAt                      <  2:624848/141‑1 (MQ=255)
                                                                                                                  accactCTAACTGATCTATATAATGAATGATGTGATCTAATTTCAATATTATTAGTATCCTATTATTTACATGTAATTACGTGTATTTTAACAAAGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTg                    <  2:189695/136‑1 (MQ=255)
                                                                                                                  aCCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTg                    >  2:610328/1‑141 (MQ=255)
                                                                                                                    cGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCt                  >  2:672368/1‑141 (MQ=255)
                                                                                                                    cGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCt                  >  2:435775/1‑141 (MQ=255)
                                                                                                                       tCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTccc               >  2:650632/1‑141 (MQ=255)
                                                                                                                       tCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTccc               >  1:388317/1‑141 (MQ=255)
                                                                                                                          aaCTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGtt            <  2:500637/141‑1 (MQ=255)
                                                                                                                                    tataATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTAc  >  2:565805/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                     
CTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTACATGTAATTACGCTTATTTTAACAAAGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTAC  >  CP000730/2386294‑2386566

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