Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,479 T→C 21.7% T25A (ACA→GCA)  USA300HOU_0026 ← IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,4790TC21.7% 75.9 / 14.6 43T25A (ACA→GCA) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (15/17);  new base C (4/5);  total (20/23)
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.89e-01

TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGA  >  CP000730/36340‑36617
                                                                                                                                           |                                                                                                                                          
taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGta                                                                                                                                           <  2:333585/141‑1 (MQ=40)
  ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtata                                                                                                                                         <  1:118692/141‑1 (MQ=40)
           tttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAta                                                                                                                                >  1:66598/1‑141 (MQ=35)
             tATGCTTTTTCCTCCAAATTTGATACAAAATTGGGGCATATTCTTGACCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                              <  1:258564/141‑1 (MQ=25)
             tATGCTTGTGGTTCTACATTTGCTACAAGAGGGGGGCATATTCTTGAACCCAACGGTTGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                              <  1:240746/141‑1 (MQ=11)
              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                             >  2:236258/1‑141 (MQ=35)
                      tcttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACgg                                                                                                                     <  2:207615/141‑1 (MQ=35)
                       cttcCATATTTGATACAAAATTGGGGCATATTCTTGAACCCACCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTCATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                    <  2:108500/141‑1 (MQ=255)
                        ttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAATGGCt                                                                                                                   >  1:291155/1‑141 (MQ=25)
                         tcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTa                                                                                                                  <  2:113635/141‑1 (MQ=35)
                             aaTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGt                                                                                                              >  1:118378/1‑141 (MQ=35)
                                 tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATa                                                                                                          >  1:323081/1‑141 (MQ=35)
                                             ggggCATATTCTTGAACCCAACGGTAGACCGTTGACTGATGAACGATTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAAGGAAGCCAACGGCTACAGTGATAACATCCTCGTTa                                                                                              >  1:172309/1‑141 (MQ=11)
                                             ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACCGATCATTTATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTa                                                                                              >  1:108224/1‑141 (MQ=11)
                                                   tatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGt                                                                                        <  1:125598/141‑1 (MQ=32)
                                                    atTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTCCCCTTAATATGTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGtt                                                                                       <  2:217931/141‑1 (MQ=17)
                                                      tCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTta                                                                                     <  1:13723/141‑1 (MQ=32)
                                                         tGAACCCAACGGTAGACCGTTGTATTATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTatat                                                                                  <  1:116132/141‑1 (MQ=17)
                                                           aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                                <  2:237456/141‑1 (MQ=32)
                                                           aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                                <  2:102456/141‑1 (MQ=32)
                                                            aCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTg                                                                               <  2:4894/141‑1 (MQ=32)
                                                                aaCGGTAGACCGTTGAATGATGAACGTTTACACCCCGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAAt                                                                           <  2:172309/141‑1 (MQ=17)
                                                                            ttGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTAATAAaga                                                               >  2:5408/1‑141 (MQ=17)
                                                                              gAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaag                                                             >  1:7897/1‑141 (MQ=32)
                                                                                atgatgAACGTTTACACCTCGTTCCATTAATATTTCACTTATATCACGATAACTCAATGTATAGCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTCAGATCATAAATAGTTCATACAGAAGAc                                                           <  1:161215/141‑1 (MQ=11)
                                                                                         gTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTg                                                  <  2:57195/141‑1 (MQ=32)
                                                                                            tACGCCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTa                                               <  2:282617/141‑1 (MQ=17)
                                                                                                 cacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAtt                                          >  1:162333/1‑141 (MQ=32)
                                                                                                       cccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACta                                    <  2:7897/141‑1 (MQ=32)
                                                                                                              tacttCAGATATATCACGATAACTCAATGCATATCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTc                             >  1:304285/4‑141 (MQ=2)
                                                                                                                  tctgATAGATCACGATAGAACAACGCATATCTTAGATAGGATCCAACGGCTACAGTTATAACAACCTTTTTAAATTGTGTATATCTGAAATAGTTCATACAGACGACTCCTTTTTGTTAAAAATATATTATAAATTCAACt                         <  2:271922/138‑1 (MQ=1)
                                                                                                                    aGATATATCACGATAGCTCAATGGATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAACTTCAACttt                       <  2:214127/141‑1 (MQ=255)
                                                                                                                        atatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGGTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCaa                   >  1:210175/1‑141 (MQ=25)
                                                                                                                         tatCACGATAGCTCAATGAATATATTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTAATACAGAAGCCTCCTTTTTATTCAAATTATGTTATAAATTCAACTTTGCaag                  >  1:130499/1‑140 (MQ=11)
                                                                                                                          atCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACa                 >  1:97397/1‑141 (MQ=35)
                                                                                                                          atCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACa                 >  2:272158/1‑141 (MQ=35)
                                                                                                                           acacggcaaCTCAATGCATATCTTAGATAGTAGCCAACAGCTACAGTGACAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAg                <  1:183352/134‑1 (MQ=1)
                                                                                                                               ggtagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTGAAAATTATATTATAAATTCAACTTTGCAACAGAAcc            <  1:47459/136‑1 (MQ=11)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt          >  2:154003/4‑141 (MQ=2)
                                                                                                                                 acgcttaaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGAAAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt          <  2:186749/135‑1 (MQ=11)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat        >  1:320322/2‑141 (MQ=12)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat        >  2:95185/2‑141 (MQ=12)
                                                                                                                                      cAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtattat     >  1:316680/1‑141 (MQ=255)
                                                                                                                                      cAATGCATATCGTAGATAGTAGCCAACGGCTATAGTGATAACATACTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGCCTCTGTTTTGTTAAAATTATATTATTTATGCAACTTTGCAACAGAACCGtatttt     <  2:56836/141‑3 (MQ=11)
                                                                                                                                         tGTATATCTTAGATAGTAGCCACCGGCTACAGTGATAACATCCTTGTTAAAGTGTTTATATCTGAAATAGTTCATACAGACGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCCACAGAACGGTCTTATGGa  >  2:51190/1‑141 (MQ=14)
                                                                                                                                           |                                                                                                                                          
TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGA  >  CP000730/36340‑36617

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: