Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,680,137 T→C 16.5% intergenic (+189/+343) cbpD → / ← USA300HOU_2537 possible secretory antigen/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,680,1370TC16.5% 153.1 / 21.7 66intergenic (+189/+343)cbpD/USA300HOU_2537possible secretory antigen/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/31);  new base C (5/6);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTCAAAACC  >  CP000730/2679999‑2680273
                                                                                                                                          |                                                                                                                                            
aaGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCg                                                                                                                                            >  2:71232/1‑141 (MQ=255)
  gTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCt                                                                                                                                          <  2:335580/141‑1 (MQ=255)
     tttattttatGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCtata                                                                                                                                       >  1:268760/1‑141 (MQ=255)
         ttttatGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACtt                                                                                                                                   >  1:396059/1‑141 (MQ=255)
         ttttatGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACtt                                                                                                                                   >  1:90242/1‑141 (MQ=255)
            tatGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAAc                                                                                                                                <  1:19903/141‑1 (MQ=255)
            tatGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAAc                                                                                                                                >  1:194153/1‑141 (MQ=255)
            tatGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAAc                                                                                                                                <  1:232305/141‑1 (MQ=255)
             atGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACa                                                                                                                               <  1:304368/141‑1 (MQ=255)
                 aaTTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACg                                                                                                                           >  2:114913/1‑141 (MQ=255)
                         tattaCTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGat                                                                                                                   >  1:81191/1‑141 (MQ=255)
                            taCTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGAtata                                                                                                                <  1:681520/141‑1 (MQ=255)
                             tcTTCTACAAATTTCAAATTGCCCTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATaa                                                                                                               <  2:643580/140‑1 (MQ=255)
                              tttCTACAATTTTCAGATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATACTATTATTTAAATAATTCGCTATAACTTAACATGCGTTTTCGATATaaa                                                                                                              <  2:317542/140‑1 (MQ=255)
                                 cTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCt                                                                                                           >  2:386494/1‑141 (MQ=255)
                                 cTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCt                                                                                                           <  1:165389/141‑1 (MQ=255)
                                   aCAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTg                                                                                                         <  2:665491/141‑1 (MQ=255)
                                     aaaTTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGtt                                                                                                       <  1:604382/141‑1 (MQ=255)
                                      aaTTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAAAAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTc                                                                                                      >  1:180995/1‑141 (MQ=255)
                                            aaaTTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAAtc                                                                                                <  1:679995/141‑1 (MQ=255)
                                                         tGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAAtttttt                                                                                <  2:194153/141‑1 (MQ=255)
                                                               tataTTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGtt                                                                         <  2:191955/141‑1 (MQ=255)
                                                                   tttCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCa                                                                     <  2:51644/141‑1 (MQ=255)
                                                                             aaCTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTAGATTTATACCTTGTTCTATCTCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAg                                                           >  1:491093/1‑141 (MQ=255)
                                                                               cTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGtt                                                         <  1:71232/141‑1 (MQ=255)
                                                                                caaattaccATCTTCGGCATCGTCTAAGTTTAAAATTTTATGCTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACTAACCTTTTTCTAAATCACT‑‑TACAATTTCTTGATTTTTTTCGTCGTCATCAGt                                                          <  2:185798/133‑1 (MQ=255)
                                                                                  ttattTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTaaa                                                      >  2:425276/1‑141 (MQ=255)
                                                                                  ttattTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTaaa                                                      >  1:243904/1‑141 (MQ=255)
                                                                                     tttCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTaaaaaa                                                   <  1:662660/141‑1 (MQ=255)
                                                                                     tttCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTaaaaaa                                                   <  1:662637/141‑1 (MQ=255)
                                                                                     tttCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTaaaaaa                                                   <  1:662634/141‑1 (MQ=255)
                                                                                         aTCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAaataa                                               >  1:474850/1‑141 (MQ=255)
                                                                                              tAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTaa                                          >  1:662603/1‑141 (MQ=255)
                                                                                               aGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAAc                                         >  2:427991/1‑141 (MQ=255)
                                                                                               aGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAAc                                         >  2:295574/1‑141 (MQ=255)
                                                                                                gCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACt                                        <  1:401871/141‑1 (MQ=255)
                                                                                                   tAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAATTGAg                                     >  1:394452/1‑141 (MQ=255)
                                                                                                   tAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAg                                     <  1:295574/141‑1 (MQ=255)
                                                                                                   tAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAg                                     <  1:583869/141‑1 (MQ=255)
                                                                                                   tAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAg                                     >  1:525950/1‑141 (MQ=255)
                                                                                                     aTCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGtt                                   >  1:87450/1‑141 (MQ=255)
                                                                                                       ctaagtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAAcc                                                                                                            >  2:6147/8‑70 (MQ=255)
                                                                                                       ctaagtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAAcc                                                                                                            <  1:6147/63‑1 (MQ=255)
                                                                                                        taagtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGata                                                                                                                  >  1:401215/7‑63 (MQ=255)
                                                                                                        taagtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGata                                                                                                                  <  2:401215/57‑1 (MQ=255)
                                                                                                         aagtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGatac                                                                                                                 <  2:408277/58‑2 (MQ=255)
                                                                                                         aagtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGatac                                                                                                                 >  1:408277/6‑62 (MQ=255)
                                                                                                          agtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAAtcactt                                                                                            <  2:523125/79‑5 (MQ=255)
                                                                                                          agtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAAtcactt                                                                                            >  1:523125/5‑79 (MQ=255)
                                                                                                           gtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAAtcactta                                                                                           <  1:669334/80‑6 (MQ=255)
                                                                                                           gtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAAtcactta                                                                                           >  2:669334/4‑78 (MQ=255)
                                                                                                             ataAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGAc                           <  1:126136/141‑1 (MQ=255)
                                                                                                                aaaTTTTATGCTATTATTTAAATAATTAGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTT‑‑‑AATAATTTTTTTCTGTTTTCATAGTAATAAGATAAAAAAATAATTTAACTGAGTTTTTGAGACCTg                        <  2:383740/141‑1 (MQ=255)
                                                                                                                 aaTTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGAcc                          <  1:154507/141‑1 (MQ=255)
                                                                                                                    tttATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTa                    <  2:708587/141‑1 (MQ=255)
                                                                                                                     ttATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGAAAAAAAAAAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTaa                   <  2:512610/141‑1 (MQ=255)
                                                                                                                       aTGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAAt                 >  2:381549/1‑141 (MQ=255)
                                                                                                                       aTGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAANTTAACTGAGTTTNTGAGACCTGTTTAAAt                 <  1:394475/141‑1 (MQ=255)
                                                                                                                        tGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAAtt                >  1:610512/1‑141 (MQ=255)
                                                                                                                            attattTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATtgatg            <  1:640241/141‑1 (MQ=255)
                                                                                                                            attattTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATtgatg            <  1:599351/141‑1 (MQ=255)
                                                                                                                              tattTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGAc          >  1:103019/1‑141 (MQ=255)
                                                                                                                               attTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACt         >  1:658330/1‑141 (MQ=255)
                                                                                                                                tttAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTc        <  2:369981/141‑1 (MQ=255)
                                                                                                                                   aaaTAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTTAATCTA‑‑‑AAAACATGTTTGCTTGTTTTATTTACATTATGTAAAATAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTCaaa     >  1:137004/1‑141 (MQ=255)
                                                                                                                                      taatTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGAAAAAAAATAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTCAAAAcc  <  2:118726/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                            
AAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTCAAAACC  >  CP000730/2679999‑2680273

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: