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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,4790TC30.4% 75.8 / 30.2 55T25A (ACA→GCA) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (13/25);  new base C (9/8);  total (22/33)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATG  >  CP000730/36340‑36615
                                                                                                                                           |                                                                                                                                        
taataaGCTTTTTTATGCCTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTATACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGta                                                                                                                                         <  1:92299/141‑1 (MQ=255)
   taaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTatat                                                                                                                                      <  2:336082/141‑1 (MQ=40)
      gCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCtt                                                                                                                                   >  1:60057/1‑141 (MQ=38)
          ttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAt                                                                                                                               <  1:337162/141‑1 (MQ=35)
           ttttTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGCACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAta                                                                                                                              <  1:152762/141‑1 (MQ=25)
           tttATGTTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAta                                                                                                                              <  1:152820/141‑1 (MQ=25)
           tttATGCGTTTTCTTCCAAGTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTATTCCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAta                                                                                                                              <  1:152768/141‑1 (MQ=255)
              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                           >  1:237534/1‑141 (MQ=35)
                        ttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCt                                                                                                                 <  1:54437/141‑1 (MQ=35)
                        ttcCAAATTGGAGACAAAATTGGGGCCTATTCTTGAACCCAACGGGAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCt                                                                                                                 <  2:168775/141‑1 (MQ=255)
                             aaTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGt                                                                                                            <  2:155771/141‑1 (MQ=34)
                             aaTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGt                                                                                                            <  2:155579/141‑1 (MQ=35)
                                     aCAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTg                                                                                                           >  1:301000/1‑134 (MQ=33)
                                     aCAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTg                                                                                                           <  2:301000/134‑1 (MQ=34)
                                      cAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATc                                                                                                   >  2:149589/1‑141 (MQ=33)
                                          aTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                                                               <  2:337420/141‑1 (MQ=32)
                                            tGGGGCATATTCTTGAACCCACCGGTAGACCGTTGAATGATGACCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGtt                                                                                             <  1:174269/141‑1 (MQ=12)
                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaaa                                                                                          >  2:274496/1‑141 (MQ=32)
                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAAATCCTTGTTaaa                                                                                          >  2:314020/1‑141 (MQ=17)
                                                    atTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGtt                                                                                     >  2:224855/1‑141 (MQ=32)
                                                       cTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTtat                                                                                  <  2:161697/141‑1 (MQ=32)
                                                           aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                              <  2:60057/141‑1 (MQ=32)
                                                                        aCCGTTGAATGATGAAAGCTTAGACCAGGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATa                                                                 >  1:393003/1‑141 (MQ=255)
                                                                            ttcaaTGAGGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACaga                                                             <  2:376735/138‑1 (MQ=17)
                                                                                   atgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTcc                                                      <  1:314492/141‑1 (MQ=32)
                                                                                        cGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttttt                                                 <  2:237534/141‑1 (MQ=17)
                                                                                                ccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTACATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAt                                         >  2:269971/1‑141 (MQ=17)
                                                                                                ccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAt                                         >  2:60889/1‑141 (MQ=32)
                                                                                                 cacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATCCAGAAGACTCCTTTTTGTTAAAAtt                                        >  2:364218/1‑141 (MQ=17)
                                                                                                            aaTATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAAt                             <  2:120899/141‑1 (MQ=21)
                                                                                                               atTTCAGTTATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                          <  2:341181/141‑1 (MQ=25)
                                                                                                               aattCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                          <  2:313371/139‑1 (MQ=25)
                                                                                                                         tatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTACAAc                >  2:390579/1‑141 (MQ=25)
                                                                                                                         aatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAAc                <  1:334007/140‑1 (MQ=25)
                                                                                                                           ttacaatatCTGAATGCATATCTTAGATAATAGCCAACTGCTACAGTTATAACATCCTTGTTAAATTTTTTATATCTGAAATAGTTCATACATAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAg              <  1:132207/136‑1 (MQ=1)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACATTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt        >  2:231374/4‑141 (MQ=9)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt        <  1:114297/138‑1 (MQ=9)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt        <  2:266808/138‑1 (MQ=9)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt        <  1:341964/138‑1 (MQ=12)
                                                                                                                                 tAACTCAATGTATATCTTAGATAGTAGCCAACGTCTACAGTCATAACATACTATTTAAATTGTTTACATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt        <  2:295754/141‑1 (MQ=11)
                                                                                                                                 tAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt        <  2:312671/141‑1 (MQ=37)
                                                                                                                                 tAACTCAATGTATATCTTAGATAGTAGACAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAGACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt        <  2:391034/141‑1 (MQ=255)
                                                                                                                                 aaaCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt        <  2:245004/140‑1 (MQ=37)
                                                                                                                                    gtCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt     >  2:102875/2‑141 (MQ=16)
                                                                                                                                    cgcAATGCGTATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt     <  2:187865/139‑1 (MQ=255)
                                                                                                                                    cTCAATGCATATCTTAGATAGTCGCCAACGGCTACAGTGATAACATCCTTGTTAAGTTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACAGtatt     >  2:288262/1‑141 (MQ=11)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGCAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt     >  1:150085/1‑141 (MQ=12)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAACTTATAGTATAAATTCAACTTTGCAACATAACCGtgat     >  1:138063/1‑138 (MQ=11)
                                                                                                                                     tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt   >  1:26978/1‑140 (MQ=12)
                                                                                                                                      cAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtattat   >  1:25919/1‑141 (MQ=40)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt   >  1:32293/1‑139 (MQ=12)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTCAAAATTAGATTATAAATTCAACTTTGCAACAGAACCGtatttt   <  1:21100/141‑3 (MQ=11)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTAAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTTCAACAGAACCGtatttt   >  1:261939/1‑139 (MQ=11)
                                                                                                                                       aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATTCTTGTTAAATTGTTTATAGCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttta  <  1:380664/141‑4 (MQ=11)
                                                                                                                                       aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatttta  >  2:321788/1‑138 (MQ=21)
                                                                                                                                       aaTGCATATCGTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTACATTGTTTGTAACTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGATCCGtatttta  <  1:91112/141‑4 (MQ=11)
                                                                                                                                           |                                                                                                                                        
TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATG  >  CP000730/36340‑36615

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