Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 736,763 A→G 9.1% I415M (ATA→ATG)  bceB_1 → Bacitracin export permease protein BceB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE736,7630AG 9.1% 111.1 / 4.8 44I415M (ATA→ATG) bceB_1Bacitracin export permease protein BceB
Reads supporting (aligned to +/- strand):  ref base A (25/15);  new base G (2/2);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

ATACAATGTAACAATTACTAGTGATAAATACATCCCTAATACTGATTTGAAACGTGGGCAAGCTGATTTATTTGTAGCGGAAGGTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGAAATAAGAAAGTATACAAAGGC  >  USA300TCH1516_ALE/736631‑736903
                                                                                                                                    |                                                                                                                                            
aTACAATGTAACAATTACTAGTGCTAAATACATCCCTAATACTGATTTGAAACGTGGGCAAGCTGAATTATTTGTAGCGGAAGGATCTATCAGAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGaa                                                                                                                                      >  1:264013/1‑141 (MQ=255)
aTACAATGTAACAATTACTAGTGATAAATACATCCCTAATACTGATTTGAAACGTGGGCAAGCTGATTTATTTATCGGGGAAGGTTCTATCAAAGATTTAGTGAAACAGAAGAAGCATGGTAAGGCAATTATAGGAACGaa                                                                                                                                      >  1:296895/1‑141 (MQ=255)
            aaTTACTAGTGATAAATACATCCCTAATACTGATTTGAAACGTGGGCAAGCTGATTTATTTGTAGCGGAAGGTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAGCATCATGt                                                                                                                          >  2:270285/1‑141 (MQ=255)
              ttACTAGTGATAAATACATCCCTAATACTGATTTGAAACGTGGGCAAGCTGATTTATTTGTAGCGGAAGGTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTa                                                                                                                        <  1:159325/141‑1 (MQ=255)
               tACTAGTGATAAATACATCCCTAATACTGATTTGAAACGTGGGCAGTCTGATTTATTTGTAGCGGAAGGTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTaa                                                                                                                       >  1:245988/1‑141 (MQ=255)
                  tAGTGATAAATACATCCCTAATACTGATTTGAAACGTGGGCAAGCTGATTTATTTGTAGCGGAAGGTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAatat                                                                                                                    >  1:314457/1‑141 (MQ=255)
                   aGTGATAAATACATCCCTAATACTGATTTGAAACGTGGGCAAGCTGATTTATTTGTAGCGGAAGGTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATAtt                                                                                                                   <  1:358938/141‑1 (MQ=255)
                      gaaaaaTACATCCCTAATACTGATTTGAAACGTGGGCAAGCTGATTTATTTGTAGCGGAAGGTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAg                                                                                                                <  2:266084/138‑1 (MQ=255)
                               aTCCCTAATACTGATTTGAAACGTGGGCAAGCTGATTTATTTGTAGCGGAAGGTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTaaa                                                                                                       <  1:327982/141‑1 (MQ=255)
                                tCCCTAATACTGATTTGAAACGTGGGCAAGCTGATTTATTTGTAGCGGAAGGTTCTATCAAAGATTTAGTGAAACATAAGCAGCATGTTAGGGCAATGATAAAAAAGAAAAAACATCATGTTAATATTAAGTTACGTAAAg                                                                                                      <  2:163383/141‑1 (MQ=255)
                                        aCTGATTTGAAACGTGGGCAAGCTGATTTATTTGTGGCGGAAGGTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATtaat                                                                                              >  2:256539/1‑141 (MQ=255)
                                             tttGAAACGTGGGCAAGCTGATTTATTTGTAGCGGAAGGTTGTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAAAGAAAAAACATCATGTTAATATTAAGTTACGTCAAGATATTAATAAAAt                                                                                         >  2:233023/1‑141 (MQ=255)
                                               tGAAACGGGGGCAAGCTGATTTATTTGTAGCGGAAGGTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCt                                                                                       <  1:346018/141‑1 (MQ=255)
                                                gAAACGTGGGCAAGCTGATTTATTTGTAGCGGAAGGTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTa                                                                                      >  1:64348/1‑141 (MQ=255)
                                                      tGGGCAAGCTGATTTATTTGTAGCGGAAGGTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTAt                                                                                >  1:128697/1‑141 (MQ=255)
                                                      tGGGCAAGCTGATTTATTTGCAGCGGAAGGTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTAt                                                                                >  1:216644/1‑141 (MQ=255)
                                                                           aGCGGAAGGTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTgg                                                           >  2:236401/1‑141 (MQ=255)
                                                                               gAAGGTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCa                                                       >  1:342099/1‑141 (MQ=255)
                                                                                   gTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATCACAGATGTTGATTTAGGTGGACCAACGt                                                   >  1:134946/1‑141 (MQ=255)
                                                                                      cTATCAAAGATTTAGTGAAACATAAGAAGCATAGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTg                                                >  1:256893/1‑141 (MQ=255)
                                                                                          cAAAGATTTAGGGAAACGGAAGAAGCAGGGGAAGGCAATTATAGGAACGAAAAAACATCATGGTGATATTAAGTGACGTAAAGATATTAATAAAATCTATTGTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCtt                                            <  1:232699/141‑1 (MQ=255)
                                                                                           aaaGATTTAGGGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTa                                           <  2:46916/141‑1 (MQ=255)
                                                                                              gATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAAt                                        >  2:174818/1‑141 (MQ=255)
                                                                                               aTTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATg                                       <  1:52196/141‑1 (MQ=255)
                                                                                               aTTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAAAGAAAAAACATTATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTCTTTTATTACAGATGTTGATTTAGGTTGACCAACGTTTGTCTTAAATg                                       <  1:295208/141‑1 (MQ=255)
                                                                                                  tAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACa                                    >  2:244479/1‑141 (MQ=255)
                                                                                                   aGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTGAAATGACaa                                   >  1:121682/1‑141 (MQ=255)
                                                                                                       aaaCATAAGAAGCATGGTAAGGCAATTAAAGAAAAGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGAc                               <  2:154131/141‑1 (MQ=255)
                                                                                                            tAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCa                          >  2:149824/1‑141 (MQ=255)
                                                                                                            tAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCa                          >  1:206479/1‑141 (MQ=255)
                                                                                                              agaagCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATaaaa                        >  2:292294/1‑140 (MQ=255)
                                                                                                               gaagCATGGTAAGGCAATTAAAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGa                       <  1:182050/141‑1 (MQ=255)
                                                                                                               aaagCATGGTAAGGCAATAAAAAAAAAAAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGa                       <  2:64348/140‑1 (MQ=255)
                                                                                                                   cATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGAAATa                   >  1:175260/1‑141 (MQ=255)
                                                                                                                      ggTAAGGCAATTATGGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGAAATAAGa                >  2:269615/1‑141 (MQ=255)
                                                                                                                      ggTAAGGCAATTATGGGAACGAAAAAACATCATGTTAATATGAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGAAATAAGa                >  2:269669/1‑141 (MQ=255)
                                                                                                                      ggTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGCCAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGAAATAAGa                >  2:210741/1‑141 (MQ=255)
                                                                                                                      ggTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGAAATAAGa                >  1:123778/1‑141 (MQ=255)
                                                                                                                       gTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGAAATAAGaa               >  2:44079/1‑141 (MQ=255)
                                                                                                                            gCAATTATAGGAACGAAAAAACATCATGTTAATATTACGTTGCGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGAAATAAGAAAGtat          >  1:246182/1‑141 (MQ=255)
                                                                                                                               aTTATGGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGAAATAAGAAAGTATACa       <  1:269669/141‑1 (MQ=255)
                                                                                                                               aTTATGGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGAAATAAGAAAGTATACa       <  1:269615/141‑1 (MQ=255)
                                                                                                                                ttATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGAAATAAGAAAGTATACaa      <  2:302594/141‑1 (MQ=255)
                                                                                                                                ttATAGGAACAAAAAAACATCATGTGAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGAAATAAGAAAGTATACaa      <  2:90379/141‑1 (MQ=255)
                                                                                                                                   taGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGAAATAAGAAAGTATACAAAgg   >  2:186175/1‑141 (MQ=255)
                                                                                                                                    aGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGAAATAAGAAAGTATACAAAGGc  >  2:118024/1‑141 (MQ=255)
                                                                                                                                    |                                                                                                                                            
ATACAATGTAACAATTACTAGTGATAAATACATCCCTAATACTGATTTGAAACGTGGGCAAGCTGATTTATTTGTAGCGGAAGGTTCTATCAAAGATTTAGTGAAACATAAGAAGCATGGTAAGGCAATTATAGGAACGAAAAAACATCATGTTAATATTAAGTTACGTAAAGATATTAATAAAATCTATTTTATGACAGATGTTGATTTAGGTGGACCAACGTTTGTCTTAAATGACAAAGACTATCAAGAAATAAGAAAGTATACAAAGGC  >  USA300TCH1516_ALE/736631‑736903

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