| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | CP000730 | 36,581 | C→T | 36.6% | intergenic (‑30/‑493) | tpn ← / → USA300HOU_0028 | IS431mec transposase/hypothetical protein | 
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 36,581 | 0 | C | T | 36.6% | 71.6 / 55.8 | 63 | intergenic (‑30/‑493) | tpn/USA300HOU_0028 | IS431mec transposase/hypothetical protein | 
| Reads supporting (aligned to +/- strand): ref base C (19/21); new base T (15/8); total (34/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.99e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01 | |||||||||||
TCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTATGTCTATT  >  CP000730/36442‑36717                                                                                                                                           |                                                                                                                                        tCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATAc                                                                                                                                          <  1:132764/140‑1 (MQ=33)       atatTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAt                                                                                                                                   <  2:75347/140‑1 (MQ=35)       atatTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAt                                                                                                                                   <  1:12640/140‑1 (MQ=34)                      cACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAg                                                                                                                    <  1:207486/140‑1 (MQ=35)                       aCGATAGCTCAATGCATTTCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGa                                                                                                                   <  1:9725/140‑1 (MQ=255)                           tAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACATGAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCg                                                                                                               <  1:61128/140‑1 (MQ=255)                             gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGGTTGTATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta                                                                                                             >  1:140894/2‑140 (MQ=11)                             gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta                                                                                                             >  2:218620/2‑140 (MQ=9)                              cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAAGAGCCGGTGTGAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGTACCGtat                                                                                                            >  1:168956/1‑140 (MQ=11)                              cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat                                                                                                            >  1:77866/1‑140 (MQ=12)                               tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGGTCATACAGAAGACGCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt                                                                                                           >  1:299583/1‑140 (MQ=11)                               tCAATGCATATCGTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAGAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt                                                                                                           >  1:226513/1‑140 (MQ=11)                                cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt                                                                                                          >  2:115346/1‑139 (MQ=255)                                  atgcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttta                                                                                                        >  2:154754/5‑137 (MQ=255)                                  aTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttt                                                                                                        >  2:66915/1‑137 (MQ=255)                                   tgcatttCTTAAATAGTAGCCAACGTCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttat                                                                                                       <  1:3994/134‑5 (MQ=11)                                    gcatttCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGCTAATATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatg                                                                                                      <  2:291464/135‑6 (MQ=11)                                    gcatatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttacta                                                                                                      <  2:113179/138‑11 (MQ=11)                                     catatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTTTTTTTTTACTATAAATTCAACTTTGCAACAGAAcctcatttactat                                                                                                     >  2:272284/2‑129 (MQ=11)                                      atatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAATTTGTTTAGATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtat                                                                                                    <  2:94594/140‑13 (MQ=11)                                      atatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgag                                                                                                    >  2:113924/1‑133 (MQ=255)                                       tatCTTAGATAGTAGCCAACGGCTGCAGGGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtatt                                                                                                   >  1:203328/1‑127 (MQ=11)                                       tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtatt                                                                                                   >  2:22239/1‑127 (MQ=11)                                       tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtatt                                                                                                   >  1:298721/1‑127 (MQ=11)                                       tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatgttatgagt                                                                                                   >  1:5573/1‑131 (MQ=255)                                       tatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTGAATTTCGAAACAGAActcaatttactatat                                                                                                   >  2:290199/1‑126 (MQ=11)                                        atCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATag                                                                                                  >  1:119098/1‑140 (MQ=255)                                         tCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtac                                                                                                 >  2:237468/1‑130 (MQ=11)                                         tCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtac                                                                                                 <  2:230519/140‑11 (MQ=11)                                          cTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATagag                                                                                                >  2:116202/1‑140 (MQ=255)                                            tAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacact                                                                                              <  2:299826/140‑14 (MQ=11)                                             aGATAGTAGCCAACGGCTACAGTTATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacactc                                                                                             <  1:20775/140‑15 (MQ=11)                                             aGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacactc                                                                                             <  1:145448/140‑15 (MQ=11)                                                tagtagCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTg                                                                                          >  1:231866/1‑140 (MQ=255)                                                   aagCCAACGGCTACAGTGATAACTTCCTTTTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTCTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTa                                                                                       <  1:320513/139‑1 (MQ=255)                                                           ggCTACAGTGATCAACTCTTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACGCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTCTGGAATAGAGATGTTGGGAACCTTtat                                                                               >  1:63689/1‑140 (MQ=17)                                                            gCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTtata                                                                              <  1:274268/140‑1 (MQ=255)                                                             cTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATAc                                                                             >  1:91964/1‑140 (MQ=255)                                                              tACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTACAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACa                                                                            <  1:10331/140‑1 (MQ=255)                                                               aCAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAg                                                                           <  2:60004/140‑1 (MQ=255)                                                                    gATAACATCCTTGTTAAATCGTTTAGATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCa                                                                      >  1:175485/1‑140 (MQ=255)                                                                    gATAACATCCTTGATAAATTGTGTATATCTGGAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGAATCa                                                                      >  1:248189/1‑140 (MQ=255)                                                                           cccTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTTCTTTTTTTTAAAATTATACTATAAATTCAACTTCGCAACACAGCCGTATGATGGAATAGAGATGTTGGTAACATTTATACACGATCATTATACt                                                               <  2:4373/139‑1 (MQ=21)                                                                               tGTTTAATTGTTTATATCTGGAAAAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAGTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAg                                                           <  1:52059/140‑1 (MQ=255)                                                                               tGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAg                                                           <  1:137111/140‑1 (MQ=255)                                                                               tGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAACGTACGAAGGAATAGAGATGTTGGTAACATTAATACAGGATCATtagaattcag                                                           >  1:299443/1‑132 (MQ=18)                                                                                gTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGt                                                          <  2:91964/140‑1 (MQ=255)                                                                                  tAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAACATTAGACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGttt                                                        <  1:42945/140‑1 (MQ=255)                                                                                   aaaTTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGAAACATTTATACAGGATCATTATACTTAAGTTTa                                                       >  2:83501/1‑140 (MQ=255)                                                                                          ttATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGt                                                >  1:78843/1‑140 (MQ=255)                                                                                             tatCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATCATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATCATTAAGTTTAATTTCGttat                                             <  2:3088/140‑1 (MQ=255)                                                                                                    aaTAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGaa                                      <  1:66217/140‑1 (MQ=255)                                                                                                    aaTAGGTCCTAATGAAGGCTCCTTTTTGTTAAAATTATACTATAAAATCAACGTTGCAACGGACGCGGTTTTTGGGATGGAGGTGTTGGTAACATTTATACAGTATCATTATACTTAAGTTTAATTTCGTTATTACAGaa                                      <  2:39651/140‑1 (MQ=25)                                                                                                     aTAGTTCATACAGAAGACTCCTTTTTGTTAAAATTGTACTATAAAATCAAGGGTGCAGCAGAACAGTATTATGGAATAGAGATGTTGGTACCATTTATACAGGATCATTATACTTAAGTTTAGTTTCGTTATTCCGGaaa                                     >  1:140548/1‑139 (MQ=25)                                                                                                        gTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACcac                                  <  1:275793/140‑1 (MQ=255)                                                                                                         ttCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACcaca                                 <  1:205027/140‑1 (MQ=255)                                                                                                              aCAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTc                            >  2:78198/1‑140 (MQ=255)                                                                                                                   agaCTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAAAACCACACATTCCAAcc                       >  1:83662/1‑140 (MQ=255)                                                                                                                          tttttGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACAAATTACAACCAGAagag                >  2:135450/1‑140 (MQ=255)                                                                                                                             ttGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACAAGAAGAGaaa             >  2:297395/1‑140 (MQ=255)                                                                                                                                  aaaaTTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAATTagc        >  1:103044/1‑138 (MQ=255)                                                                                                                                        ataCTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTATGTCTAtt  <  2:245500/140‑1 (MQ=255)                                                                                                                                        ataCTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTATGTCTAtt  >  2:270156/1‑140 (MQ=255)                                                                                                                                        ataCTAAAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTATGTCTAtt  <  2:292517/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                        TCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTATGTCTATT  >  CP000730/36442‑36717 | 
| 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: ‑ |