Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 547,929 A→G 11.6% K191K (AAA→AAG ftsH → ATP‑dependent zinc metalloprotease FtsH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE547,9290AG11.6% 105.3 / 4.7 42K191K (AAA→AAGftsHATP‑dependent zinc metalloprotease FtsH
Reads supporting (aligned to +/- strand):  ref base A (20/17);  new base G (2/3);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TAAATCTAAAGCAAAAATGTACGATAATAATAAACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACTTTGTAGAGATGT  >  USA300TCH1516_ALE/547800‑548068
                                                                                                                                 |                                                                                                                                           
aaaaTCTAAAGCAAAAATGTACGATAATAATAAACGTCGTGTTCGTCTCTCTGATGTAGCGGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATc                                                                                                                                  <  2:106484/140‑1 (MQ=255)
        aaGCAAAAATGTACGATAATAATAGACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTc                                                                                                                          >  2:179122/1‑141 (MQ=255)
        aaGCAAAAATGTACGATAATAATAAACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTCATTTCTTGAAAGATAATAAAAAATTCGAAGAGATGGGATCTAGGATTc                                                                                                                          >  1:49017/1‑141 (MQ=255)
          gCAAAAATGTACGATAATAATAAACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCt                                                                                                                        <  1:134403/141‑1 (MQ=255)
           cAAAAATGTACGATAATAATAAACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTa                                                                                                                       >  1:23293/1‑141 (MQ=255)
              aaaTGTACGATAATAATAAACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGCTGATTTCATGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAg                                                                                                                    >  1:282240/1‑141 (MQ=255)
                aTGTACGATAATAATAAACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGgt                                                                                                                  >  1:317149/1‑141 (MQ=255)
                     cGATAATAATAAACGGCGGGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCtt                                                                                                             <  1:100245/141‑1 (MQ=255)
                       ataataataAACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTat                                                                                                           >  1:326344/1‑140 (MQ=255)
                           taataaACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACttgt                                                                                                       <  2:55915/141‑1 (MQ=255)
                                 aCGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGGCTTACTTGTTGGAcc                                                                                                 >  1:39868/1‑141 (MQ=255)
                                        gtTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTAGGACCTCCAggg                                                                                          >  2:55879/1‑140 (MQ=255)
                                           cGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAGGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACt                                                                                       >  1:316277/1‑141 (MQ=255)
                                           cGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAGGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACt                                                                                       >  1:143045/1‑141 (MQ=255)
                                                 tctGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCAAGGTACTGGTaaa                                                                                 >  1:144692/1‑141 (MQ=255)
                                                   tGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAAc                                                                               <  2:62359/141‑1 (MQ=255)
                                                          gCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACtt                                                                        >  1:169961/1‑141 (MQ=255)
                                                          gCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACtt                                                                        >  1:158912/1‑141 (MQ=255)
                                                          gCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTATTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACtt                                                                        >  1:170053/1‑141 (MQ=255)
                                                             ggggCAGATGCAGAAAAACAAGAATGAATTGAAATTGTTGATTGCGTGAAAGATAACAAAAAATGCAAAGAGATGGGATCTAGGATTCCTAAAGGTGGCTTACTTGTTGGACCTCCAGGTACTGGTAAAAAGTTACTTGCt                                                                     >  2:90693/1‑141 (MQ=255)
                                                               ggCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAGTGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTag                                                                   <  1:10587/141‑1 (MQ=255)
                                                               ggCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTag                                                                   <  1:10585/141‑1 (MQ=255)
                                                                  aGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGc                                                                <  1:301898/141‑1 (MQ=255)
                                                                           aaaaCAAGAATTAATTGAAATTGTTACTGGTTTTGAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGATGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATGACTTGCTAGAGCGGTTGGAgg                                                       >  2:175489/1‑141 (MQ=255)
                                                                             aaCAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTg                                                     >  2:233847/1‑141 (MQ=255)
                                                                                           gAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCa                                       <  1:253596/141‑1 (MQ=255)
                                                                                           gAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATGACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCa                                       <  2:157486/141‑1 (MQ=255)
                                                                                                            gAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTgg                      >  1:383277/1‑141 (MQ=255)
                                                                                                              aaGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGtt                    >  2:369853/1‑141 (MQ=255)
                                                                                                                   aataaAACATTTAAGGAAAGGGGATCTAGGAGTCCTAAAGTTGTCCTACTAGTGGCTCCTCCAGGTACTTTTAACATTTTGCTTGCTAGAGCGGTTGCAGGTGAAGCTCGCGCACCATTCTTCTCTATTAGTGGTTCATAc               <  2:185787/141‑1 (MQ=255)
                                                                                                                   aataaAAAATTCAAGGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGAc               <  2:143045/141‑1 (MQ=255)
                                                                                                                   aataaAAAATTCAAGGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGAc               <  2:316277/141‑1 (MQ=255)
                                                                                                                    ataaAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACt              <  2:169961/141‑1 (MQ=255)
                                                                                                                    ataaAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACt              <  2:170053/141‑1 (MQ=255)
                                                                                                                    ataaAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCACTCCATTCTTCTCTATTAGTGGTTCAGACt              <  2:158912/141‑1 (MQ=255)
                                                                                                                     taaAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCGGGTGACGCTGGCGCACCATTCTTCTCTATTGGTGGTTCAggagt             >  2:111892/1‑137 (MQ=255)
                                                                                                                     taaAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACtt             <  2:139241/141‑1 (MQ=255)
                                                                                                                     taaAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACtt             <  2:139266/141‑1 (MQ=255)
                                                                                                                     taaAAAATTCAAAGAAATGGGATCGAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACtt             <  2:49017/141‑1 (MQ=255)
                                                                                                                      aaaaaaTTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTGTTCTCTATTAGTGGTTCAGACttt            >  2:35277/1‑141 (MQ=255)
                                                                                                                       aaaaaTTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACTTTg           <  1:20681/141‑1 (MQ=255)
                                                                                                                       aaaaaTTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAAGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACttgg           >  2:96800/1‑139 (MQ=255)
                                                                                                                                aaGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACTTTGTAGAGAtgt  >  1:67746/1‑141 (MQ=255)
                                                                                                                                 |                                                                                                                                           
TAAATCTAAAGCAAAAATGTACGATAATAATAAACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACTTTGTAGAGATGT  >  USA300TCH1516_ALE/547800‑548068

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: