Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,581 C→T 34.4% 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
*CP00073036,5810CT34.4% 78.5 / 44.4 58intergenic (‑30/‑493)tpn/USA300HOU_0028IS431mec transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/11);  new base T (9/11);  total (36/22)
Fisher's exact test for biased strand distribution p-value = 8.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

TTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTAT  >  CP000730/36441‑36710
                                                                                                                                             |                                                                                                                                 
ttCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGCAGACTCCTTTTTGTTAAAATTATAc                                                                                                                                   <  2:353262/141‑1 (MQ=17)
     ttAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaa                                                                                                                              >  2:206968/1‑141 (MQ=33)
        atatTTCAGATATATCACAATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAtt                                                                                                                           <  2:264455/141‑1 (MQ=25)
                      tCACGATAGCTCAATGAATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAg                                                                                                             >  1:441216/1‑141 (MQ=12)
                         cgagagcTCAATGCATATCTTAGTTAGTAGCCAACGTCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAGACTGAAGACTCCTTTTGGTTAAAATTATATGATAAATTCAACTTGGCAACAGGAc                                                                                                          <  2:339653/135‑1 (MQ=11)
                           aTAGCTCAATGAATCTCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcct                                                                                                        <  2:87522/141‑2 (MQ=255)
                            tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                                                                                                       <  1:194475/138‑1 (MQ=9)
                            gaaCTCAATGCAGATCTTAGATAGGAGCCAACAGCTACAGTGATAACATCCTTGTTAAAGGGTTTATGTGTGAAATAGGTCATACAGAAGACTCCGTTTTGTTAAAATTAGACTATAAATTCAACTTGGCAACAGAACCgc                                                                                                       >  2:11779/2‑140 (MQ=11)
                              gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGTACCAGAACCGtat                                                                                                     >  1:355051/2‑141 (MQ=255)
                              gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat                                                                                                     <  1:55524/140‑1 (MQ=12)
                              gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat                                                                                                     >  2:80198/2‑141 (MQ=12)
                              gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat                                                                                                     <  2:260222/140‑1 (MQ=12)
                                 cAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtattat                                                                                                  >  1:344436/1‑141 (MQ=255)
                                   atgcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATGTTATAAATTCAACTTTGCAACAGAACCGtattttat                                                                                                >  1:431513/5‑137 (MQ=255)
                                    tgcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatg                                                                                               >  1:206287/4‑136 (MQ=255)
                                     gcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatg                                                                                               >  2:46179/3‑136 (MQ=11)
                                     gcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatga                                                                                              >  1:286473/3‑135 (MQ=255)
                                       atatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATa                                                                                            <  1:5317/141‑1 (MQ=255)
                                        tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTTAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagta                                                                                           <  1:210698/141‑10 (MQ=11)
                                        tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtatta                                                                                           >  1:366780/1‑127 (MQ=11)
                                         atCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtac                                                                                          >  1:155771/1‑131 (MQ=255)
                                           cTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacac                                                                                        <  2:170464/141‑13 (MQ=11)
                                           cTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAgaga                                                                                        >  1:437625/1‑141 (MQ=255)
                                            ttAGATAGTAGCGAACGGCGAGAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTGTGCAACAGAACCGtattttatgagtacact                                                                                       <  2:409748/141‑14 (MQ=11)
                                             tAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacactc                                                                                      <  1:308050/141‑15 (MQ=11)
                                             tAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacactc                                                                                      <  1:129192/141‑15 (MQ=11)
                                             tAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATAGCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacactc                                                                                      <  1:243131/141‑15 (MQ=11)
                                                aTAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTg                                                                                   >  1:428972/1‑141 (MQ=255)
                                                   gtagCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTa                                                                                >  2:146942/1‑141 (MQ=255)
                                                     agCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAAc                                                                              >  1:194300/1‑141 (MQ=255)
                                                     agCCAACGGCTACAGTGATAACAACCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAAc                                                                              >  2:398298/1‑141 (MQ=255)
                                                      gCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGTACCGTATTATGGAATAGAGATGTTGGTAACa                                                                             >  2:207505/1‑141 (MQ=255)
                                                        cgttggCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACAtt                                                                           >  2:251960/5‑141 (MQ=255)
                                                         aaCGGCTACAGTGATAACATCCTTGTTAAATCGTGTATAGCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACAttt                                                                          >  1:448849/1‑141 (MQ=255)
                                                           cGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTtat                                                                        >  2:155299/1‑141 (MQ=255)
                                                            ggCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTATTTTTTTACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTtata                                                                       >  2:79414/1‑141 (MQ=18)
                                                               tACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAg                                                                    >  1:245196/1‑141 (MQ=255)
                                                                aCAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAgg                                                                   >  1:187579/1‑141 (MQ=255)
                                                                     gATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCAt                                                              <  2:142773/141‑1 (MQ=255)
                                                                          cATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATAc                                                         <  2:86266/141‑1 (MQ=255)
                                                                           aTCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACt                                                        >  1:342932/1‑141 (MQ=255)
                                                                           aTCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATGATACt                                                        >  2:370264/1‑141 (MQ=255)
                                                                             ccTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTa                                                      >  2:143967/1‑141 (MQ=255)
                                                                              cTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTaa                                                     <  2:347233/141‑1 (MQ=255)
                                                                               ttGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAg                                                    <  1:398298/141‑1 (MQ=255)
                                                                                 gTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGtt                                                  >  2:332940/1‑141 (MQ=255)
                                                                                 gTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGtt                                                  <  2:447634/141‑1 (MQ=255)
                                                                                   tAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTa                                                >  2:181424/1‑141 (MQ=255)
                                                                                        tGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTTTTTTTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTc                                           >  2:50617/1‑141 (MQ=37)
                                                                                                 cTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACa                                  >  2:261253/1‑141 (MQ=255)
                                                                                                  tGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAg                                 <  2:274089/141‑1 (MQ=255)
                                                                                                   gAAATAGTTCATACAGAAGACTCCTTTTTGTTTAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGa                                >  2:146815/1‑141 (MQ=255)
                                                                                                   gAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGa                                <  2:236106/141‑1 (MQ=255)
                                                                                                        aGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACcac                           <  1:35125/141‑1 (MQ=255)
                                                                                                                   aagaCTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAAcc                >  2:95702/1‑141 (MQ=255)
                                                                                                                   aagaCTCCTTGATGTAAAAAGTATACTATAAATTCAACTTTGCAACAGAACCATATTATGGAATAGAAATATTGGTAACATGTATATAAGATCATTATACTTAAGTTTAATTTGGTTATTACACAACCACACATTACAAcc                >  2:95645/1‑141 (MQ=37)
                                                                                                                         ccTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCagaaga          <  1:340178/141‑1 (MQ=255)
                                                                                                                            ttttGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGaa       >  1:419430/1‑141 (MQ=255)
                                                                                                                                 ttAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTAt  >  1:362577/1‑141 (MQ=255)
                                                                                                                                             |                                                                                                                                 
TTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTAT  >  CP000730/36441‑36710

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