Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,243,499 T→A 10.2% intergenic (+77/+32) USA300HOU_2114 → / ← pyrG hypothetical protein/CTP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,243,4990TA10.2% 181.1 / 12.8 59intergenic (+77/+32)USA300HOU_2114/pyrGhypothetical protein/CTP synthase
Reads supporting (aligned to +/- strand):  ref base T (20/33);  new base A (3/3);  total (23/36)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

ATGTATACTATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAGCAATAAAG  >  CP000730/2243364‑2243639
                                                                                                                                       |                                                                                                                                            
aTGTATACTATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACaa                                                                                                                                         <  1:669455/141‑1 (MQ=255)
 tGTATACTATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACatt                                                                                                                                        <  2:254513/141‑3 (MQ=255)
   tataCTATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACATTGa                                                                                                                                      <  1:2041/141‑1 (MQ=255)
     taCTATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGaaa                                                                                                                                    >  2:222490/1‑141 (MQ=255)
       cTATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGt                                                                                                                                  >  2:482049/1‑141 (MQ=255)
        tATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTa                                                                                                                                 >  2:763960/1‑141 (MQ=255)
           gACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCg                                                                                                                              <  1:161249/141‑1 (MQ=255)
             cATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGt                                                                                                                            >  2:366576/1‑141 (MQ=255)
                 tataCACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTtta                                                                                                                        >  1:572303/1‑141 (MQ=255)
                   taCACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTttatt                                                                                                                      >  1:241085/1‑141 (MQ=255)
                    acacAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACat                                                                                                                                         <  1:758140/121‑2 (MQ=255)
                    acacAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACat                                                                                                                                         >  2:758140/1‑120 (MQ=255)
                      acaAATGATTGATATGACAATAGACTTAAAAATATAAAAGTCACTATTGTTATTTCATCCCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTATATTAGc                                                                                                                   <  2:545852/141‑1 (MQ=255)
                      acaAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACAGGTTTTATTAGc                                                                                                                   >  1:768023/1‑141 (MQ=255)
                           taaTTGAAACTACAAGAAACTTATAATGATAAAATTCACTATTATTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatt                                                                                                              <  2:467993/139‑1 (MQ=255)
                            gATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAattt                                                                                                             <  2:696420/141‑1 (MQ=255)
                              ttGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttat                                                                                                           <  1:52148/141‑1 (MQ=255)
                              ttGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttat                                                                                                           <  2:547380/141‑1 (MQ=255)
                               tGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatt                                                                                                          <  1:192936/141‑1 (MQ=255)
                                  aaTGATAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatttat                                                                                                       >  2:165823/1‑141 (MQ=255)
                                  aaTGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGTAAatttatttat                                                                                                       >  2:224730/1‑141 (MQ=255)
                                        aaGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgtt                                                                                                 <  1:630452/141‑1 (MQ=255)
                                           agaCTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGat                                                                                              <  1:567530/141‑1 (MQ=255)
                                             aCTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGatat                                                                                            <  1:361025/141‑1 (MQ=255)
                                                 aGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATCTTAAAATTATACGAAAAACCGGAACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTACTTTGTTGATATTTTa                                                                                        >  2:755566/1‑141 (MQ=255)
                                                 aGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTa                                                                                        <  2:249556/141‑1 (MQ=255)
                                                    aTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATg                                                                                     <  2:242754/141‑1 (MQ=255)
                                                        taAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGc                                                                                 <  2:718837/141‑1 (MQ=255)
                                                         aaaaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCt                                                                                >  1:45606/1‑141 (MQ=255)
                                                                 cTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGCACTTTCACTGTTTAACAATGAAAGTACAGGCATTATTAACAAGTATATTTATTTTGTTGATAATTTAATGAAGCTTCAATaaa                                                                        >  1:508521/1‑141 (MQ=255)
                                                                         aTTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaa                                                                >  2:204548/1‑141 (MQ=255)
                                                                              aTACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATc                                                           <  2:500513/141‑1 (MQ=255)
                                                                                aCCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCgg                                                         >  1:152015/1‑141 (MQ=255)
                                                                                      tCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCgg                                                   <  1:96442/141‑1 (MQ=255)
                                                                                       cACTAACATAATATTAAAATTATACGAACAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATAAATTAAAAATCGGGTGCGGa                                                  <  2:393081/141‑1 (MQ=255)
                                                                                          tAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGa                                               <  1:595806/141‑1 (MQ=255)
                                                                                          tAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGa                                               <  2:241085/141‑1 (MQ=255)
                                                                                            aCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGAtt                                             <  2:301437/141‑1 (MQ=255)
                                                                                               taataTTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTgg                                          <  2:145176/141‑1 (MQ=255)
                                                                                                   atTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGtata                                      >  1:166588/1‑141 (MQ=255)
                                                                                                         aTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAg                                <  1:251071/141‑1 (MQ=255)
                                                                                                          ttATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGa                               <  2:302018/141‑1 (MQ=255)
                                                                                                               cGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCt                          >  2:98747/1‑141 (MQ=255)
                                                                                                                gAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTg                         >  1:682917/1‑141 (MQ=255)
                                                                                                                  aaaaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGTAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTggg                       <  1:224730/141‑1 (MQ=255)
                                                                                                                      ccGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCa                                                                             <  2:260916/83‑1 (MQ=255)
                                                                                                                      ccGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCa                                                                             >  1:260916/1‑83 (MQ=255)
                                                                                                                      ccGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGa                   <  2:681851/141‑1 (MQ=255)
                                                                                                                         gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAAtt                >  1:396670/1‑141 (MQ=255)
                                                                                                                            cTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGAc             <  2:56888/141‑1 (MQ=255)
                                                                                                                              ttCATTGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACaa           >  2:288026/1‑141 (MQ=255)
                                                                                                                               tCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAg          <  2:103745/141‑1 (MQ=255)
                                                                                                                               tCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAg          <  2:532565/141‑1 (MQ=255)
                                                                                                                                cAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAGc         <  2:545849/141‑1 (MQ=255)
                                                                                                                                cAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAGc         <  2:768023/141‑1 (MQ=255)
                                                                                                                                 attGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAGCa        >  2:147466/3‑141 (MQ=255)
                                                                                                                                 actgtgTAACCATGAAAGTTACGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGGGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAGCa        <  2:166588/139‑1 (MQ=255)
                                                                                                                                 aaTGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAGCa        >  1:619178/1‑141 (MQ=255)
                                                                                                                                    gTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAGCAATa     <  2:360805/141‑1 (MQ=255)
                                                                                                                                     tttAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAGCAATaa    >  1:226842/1‑141 (MQ=255)
                                                                                                                                      ttAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTATGGGTGGAATTGACAAGCAATaaa   >  2:407897/1‑141 (MQ=255)
                                                                                                                                       tAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAGCAATAAAg  <  1:551126/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                            
ATGTATACTATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAGCAATAAAG  >  CP000730/2243364‑2243639

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