Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,455,698 2 bp→AT coding (1216‑1217/1260 nt) USA300HOU_2315 → hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,455,6980CA91.7% 183.8 / 11.1 60coding (1216/1260 nt)USA300HOU_2315hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (24/31);  minor base . (2/3);  total (26/34)
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.97e-01
*CP0007302,455,6990AT98.3% 191.6 / NA 60H406L (CAT→CTT) USA300HOU_2315hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (25/34);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTT  >  CP000730/2455567‑2455824
                                                                                                                                   ||                                                                                                                             
cAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAATACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaa                                                                                                                      <  1:1364591/141‑1 (MQ=255)
  gTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaaa                                                                                                                    >  1:94459/1‑141 (MQ=255)
           gCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATGGAAAAAACTGAAGAAc                                                                                                           <  1:818036/141‑1 (MQ=255)
           gCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAAc                                                                                                           <  1:1097068/141‑1 (MQ=255)
           gCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAAc                                                                                                           <  1:603341/141‑1 (MQ=255)
           gCACCACAAAGTCACGATTTAACAAACGATTCAAATAAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAAc                                                                                                           <  1:1230195/141‑1 (MQ=255)
           acaccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGACCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAAc                                                                                                           <  1:793492/140‑1 (MQ=255)
            caccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAGGAACCTGCGACTTATAAAGAAGACTGACGTTCATCAGTTGTAATTGATGCTAAAAAATGATCGAAAAAACTGAAGAACa                                                                                                          <  1:890633/141‑1 (MQ=255)
            caccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACa                                                                                                          <  1:444311/141‑1 (MQ=255)
            caccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCAGTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACa                                                                                                          <  1:862574/141‑1 (MQ=255)
                       gACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGCAAAAACTGAAGAACAATCTTCAGATa                                                                                               >  1:331410/1‑141 (MQ=255)
                        aCGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAACAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAGCTGAAGAACAATCTTCAGATaa                                                                                              <  1:1293767/141‑1 (MQ=255)
                        aCGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaa                                                                                              <  1:683569/141‑1 (MQ=255)
                        aCGATTTAACAAACGATTCAAATGAAGAGGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaa                                                                                              <  1:180822/141‑1 (MQ=255)
                          gATTTAACAAACGATTCAAATGAAGATGAACTAGAGCAATCATTCAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCCGATaaaa                                                                                            >  1:814360/1‑141 (MQ=255)
                          gATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaa                                                                                            <  1:201554/141‑1 (MQ=255)
                           aTTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaaa                                                                                           <  1:1190288/141‑1 (MQ=255)
                            tttAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAAc                                                                                          <  1:1004576/141‑1 (MQ=255)
                            tataataaACTATTAAACTGCAGAAGCAATAAATCATGCTTTAAAAGTACCATTGAATTATAAAGAAGACTGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAAc                                                                                          <  1:1320971/139‑1 (MQ=255)
                                aCAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaat                                                                                      <  1:413632/141‑1 (MQ=255)
                                aCAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaat                                                                                      <  1:499586/141‑1 (MQ=255)
                                 cAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaata                                                                                     >  1:1177740/1‑141 (MQ=255)
                                    aCGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAataata                                                                                  >  1:1256976/1‑141 (MQ=255)
                                          cAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAGAATGAATACGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaaaa                                                                            >  1:1069583/1‑140 (MQ=255)
                                          cAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAg                                                                            >  1:1262396/1‑141 (MQ=255)
                                                 aGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGc                                                                     <  1:862621/141‑1 (MQ=255)
                                                 aGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGGTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGc                                                                     <  1:862413/141‑1 (MQ=255)
                                                 aGATGAAATAGAGCAATCATGAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACACATAATAACAATGCGCTTGc                                                                     <  1:677483/141‑1 (MQ=255)
                                                   aTGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACATACGTTCATCAGTTGTAATTGATGCAGAAAAATTATCGAAAAAACTGAAGAATAATCTTCACATAAAAACATATAATAACAAAGCGCCTGCTa                                                                   <  1:1014930/141‑1 (MQ=255)
                                                     gAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTGGCTAGt                                                                 >  1:152947/1‑141 (MQ=255)
                                                           gagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCCGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctct                                                           <  1:258373/141‑1 (MQ=255)
                                                                    ttAAAAGAACCTGCGACTTATAAAGAAGACAGACGTACATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCATGCTAGTACCTCTTAAAAAgat                                                  >  1:954113/1‑141 (MQ=255)
                                                                      aaaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAcatga                                                >  1:1308501/1‑141 (MQ=255)
                                                                      aaaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTGGCTAGTACCTCTTAAAAAgatga                                                >  1:1284921/1‑141 (MQ=255)
                                                                         aGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGc                                             <  1:812025/141‑1 (MQ=255)
                                                                          gAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaaaaacgagacc                                            >  1:162422/1‑132 (MQ=255)
                                                                             ccTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACCAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGc                                         >  1:781849/1‑141 (MQ=255)
                                                                              cTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCa                                        <  1:614321/141‑1 (MQ=255)
                                                                               tGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCaa                                       >  1:1217678/1‑141 (MQ=255)
                                                                               tGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCaa                                       >  1:447364/1‑141 (MQ=255)
                                                                               tGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACCGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTACCaa                                       >  1:673007/1‑141 (MQ=255)
                                                                                 cGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgc                                     >  1:748490/1‑141 (MQ=255)
                                                                                  gACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAACAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcg                                    <  1:121736/141‑1 (MQ=255)
                                                                                   aCTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTCCCACTTAAAACGATGATGCTAGCAAgcgc                                   >  1:1197425/1‑141 (MQ=255)
                                                                                    cTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCt                                  >  1:847748/1‑141 (MQ=255)
                                                                                     ttATAAACAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCtt                                 <  1:23972/141‑1 (MQ=255)
                                                                                       ataAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACCACGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCtttt                               <  1:666729/141‑1 (MQ=255)
                                                                                        taAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTc                              <  1:1175306/141‑1 (MQ=255)
                                                                                              agacagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgggactttttatcct                        >  1:353521/1‑127 (MQ=255)
                                                                                              agacagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatact                        >  1:1117105/1‑141 (MQ=255)
                                                                                                 cagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGCTGCTAGCAAGCGCTTTTctatactata                     >  1:498587/1‑141 (MQ=255)
                                                                                                 cagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatactata                     >  1:1226773/1‑141 (MQ=255)
                                                                                                 cagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatactata                     >  1:943163/1‑141 (MQ=255)
                                                                                                 cagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGATTTTctatactata                     >  1:217996/1‑141 (MQ=255)
                                                                                                    acGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTatatat                  <  1:720646/141‑1 (MQ=255)
                                                                                                      gTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTTCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATAtatta                <  1:1006254/141‑1 (MQ=255)
                                                                                                      gTTCATCAGCTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATAtatta                <  1:797643/141‑1 (MQ=255)
                                                                                                           tcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTAttttt           <  1:19285/141‑1 (MQ=255)
                                                                                                              gTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTa        <  1:195026/141‑1 (MQ=255)
                                                                                                                    aTTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTATTCTATACTATATATTATTTTTCTAAACttt  >  1:1266620/1‑141 (MQ=255)
                                                                                                                                   ||                                                                                                                             
CAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTT  >  CP000730/2455567‑2455824

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: