Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,620,970 A→G 7.3% K42E (AAA→GAA)  USA300HOU_2484 → hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,620,9700AG 7.3% 154.8 / 3.7 55K42E (AAA→GAA) USA300HOU_2484hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base A (32/19);  new base G (2/2);  total (34/21)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

ATCAGGAAGGTGCAACAATGACCGTCAACATCGTATTTAGTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATGATTTTATTGCATTAGCATTCTTATTTTTAACGATAGCATTAT  >  CP000730/2620830‑2621102
                                                                                                                                            |                                                                                                                                    
aTCAGGAAGGTGCAACAATGACCGTCAACATCGTATTTAGTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATa                                                                                                                                      >  1:116542/1‑141 (MQ=255)
    ggAAGGTGCAACAATGACCGTCAACATCGTATTTAGTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAaca                                                                                                                                  >  2:267786/1‑141 (MQ=255)
     gAAGGTGCAACAATGCCCGTCAACATCGTATTTAGTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAacac                                                                                                                                 <  2:505273/141‑1 (MQ=255)
      ggggTGCAACAAGGGCCGTCAAGATCGTATTGGGTATCGTATTCTGCATAGGGAGGGTGATATTAGGCGTTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAAtaca                                                                                                                                <  2:607990/139‑1 (MQ=255)
        ggTGCAACAATGACCGTCAACATCGTATTTAGTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAg                                                                                                                              >  1:616559/1‑141 (MQ=255)
        ggTGCAACAATGACCGTCAACATCGTATTTAGTATCGTATTCTGCATAAATATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACAATAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAg                                                                                                                              >  2:568485/1‑141 (MQ=255)
         gTGCAACAATGACCGTCAACATCGTATTTAGTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGa                                                                                                                             >  2:642191/1‑141 (MQ=255)
         gTGCAACAATGACCGTCAACATCGTATTTAGTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGa                                                                                                                             >  2:663992/1‑141 (MQ=255)
         gTGCAACAATGACCGTCAACATCGTATTTAGTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGa                                                                                                                             >  1:88448/1‑141 (MQ=255)
           gCAACAATGACCGTCAACATCGTATTTAGTATCGTATTCTGCATAAGTATGGTCATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAAc                                                                                                                           <  2:276114/141‑1 (MQ=255)
                      cGTCAACATCGTATTTAGTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCata                                                                                                                <  2:615643/141‑1 (MQ=255)
                        tCAACATCGTATTTAGTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCatata                                                                                                              <  1:538321/141‑1 (MQ=255)
                           acgtCGTATTTAGTATCGTATTCTGCATGAGTATGGTTAGATTAGGCAGGTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAAATAAACACAAGAACCAAATTGCagatatat                                                                                                           >  2:600991/4‑141 (MQ=255)
                               cGTATTTAGTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTAc                                                                                                       >  2:627248/1‑141 (MQ=255)
                                gTATTTAGTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACg                                                                                                      >  1:588538/1‑141 (MQ=255)
                                   tttAGTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTa                                                                                                   >  2:297911/1‑141 (MQ=255)
                                    ttAGTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTAACACTCATATCTTTTATAAATCAAACGGACATTGCAGATAAACACAAGAACCACATTGCATATATATTTACGatat                                                                                                  >  1:708379/1‑141 (MQ=255)
                                     tAGTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATg                                                                                                 <  2:440288/141‑1 (MQ=255)
                                       gTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGta                                                                                               <  2:241344/141‑1 (MQ=255)
                                       gTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGAAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGta                                                                                               <  2:126440/141‑1 (MQ=255)
                                        tATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGtat                                                                                              >  2:658850/1‑141 (MQ=255)
                                        tATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGtat                                                                                              >  2:658860/1‑141 (MQ=255)
                                          tCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATTTCGTAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATaa                                                                                            >  2:37998/1‑141 (MQ=255)
                                                       aaGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGt                                                                               <  2:682951/141‑1 (MQ=255)
                                                            tGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTAtttt                                                                          >  1:452752/1‑141 (MQ=255)
                                                             ggTTATATTAGGCATATATGTCGCAATAACTAAAGATTTCACACTAACCTCTTTTATAAATCAAACGGCCATTGCAGCTACACACAAGAACCAAATTGCATATATATTTACGCTATGCACACGTTTGAGTGCAGTAttttt                                                                         >  2:276472/1‑141 (MQ=255)
                                                               ttATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTaa                                                                       <  1:642191/141‑1 (MQ=255)
                                                                   atTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTc                                                                   >  1:442614/1‑141 (MQ=255)
                                                                     tAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCtt                                                                 <  1:534696/141‑1 (MQ=255)
                                                                     tAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCtt                                                                 >  2:87474/1‑141 (MQ=255)
                                                                            ttATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCt                                                          <  1:627248/141‑1 (MQ=255)
                                                                               tGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAg                                                       >  1:259838/1‑141 (MQ=255)
                                                                                gTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATGAACACAAg                                                                                                                              <  2:687771/69‑1 (MQ=255)
                                                                                gTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATGAACACAAg                                                                                                                              >  1:687788/1‑69 (MQ=255)
                                                                                gTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATGAACACAAg                                                                                                                              <  2:687788/69‑1 (MQ=255)
                                                                                gTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATGAACACAAg                                                                                                                              >  1:687771/1‑69 (MQ=255)
                                                                                    cAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTg                                                  <  2:667988/141‑1 (MQ=255)
                                                                                         aCTAAAGATTTCACACTAATCTCATTTATAAATCAAACGGTCATTGCAGATAAACATAAGAATTAATTTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATGTTTAATGTCTTCAATGCTAAGCTTTGAatag                                             <  2:453258/141‑2 (MQ=255)
                                                                                          cTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATg                                            >  2:247452/1‑141 (MQ=255)
                                                                                          cTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATg                                            >  2:247436/1‑141 (MQ=255)
                                                                                            aaaGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATGAt                                          >  1:667434/1‑141 (MQ=255)
                                                                                            aaaGATTTCACACTAATCTCTTTTATAAATC‑AACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATGAtt                                         >  1:248108/1‑141 (MQ=255)
                                                                                             aaGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATGAtt                                         >  2:80661/1‑141 (MQ=255)
                                                                                               gATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATGAtttt                                       >  2:403983/1‑141 (MQ=255)
                                                                                                    cacaCTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATGATTTTATTGc                                  >  1:491622/1‑141 (MQ=255)
                                                                                                     acacTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATGATTTTATTGCa                                 <  2:616559/141‑1 (MQ=255)
                                                                                                        cTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATGATTTTATTGCATTa                              >  2:504419/1‑141 (MQ=255)
                                                                                                              tctTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATGATTTTATTGCATTAGCATTc                        >  1:306080/1‑141 (MQ=255)
                                                                                                                ttttATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGATTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATGATTTTATTGCATTAGCATTCtt                      >  2:609172/1‑141 (MQ=255)
                                                                                                                   tataAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATGATTTTATTGCATTAGCATTCTTAtt                   <  1:37998/141‑1 (MQ=255)
                                                                                                                     taAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATGATTTTATTGCATTAGCATTCTTAtttt                 >  1:545646/1‑141 (MQ=255)
                                                                                                                      aaaTCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCatatatat                                                                                                           >  1:23546/1‑50 (MQ=255)
                                                                                                                      aaaTCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCatatatat                                                                                                           <  2:23546/50‑1 (MQ=255)
                                                                                                                       aaTCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATGATTTTATTGCATTAGCATTCTTATTTTTa               <  2:554768/141‑1 (MQ=255)
                                                                                                                                 ccATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCGATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATGATTTTATTGCATTAGCATTCTTATTTTTAACGATAGCat     >  1:404584/1‑141 (MQ=255)
                                                                                                                                    ttGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATGATTTTATTGCATTAGCATTCTTATTTTTAACGATAGCattat  <  2:562938/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                    
ATCAGGAAGGTGCAACAATGACCGTCAACATCGTATTTAGTATCGTATTCTGCATAAGTATGGTTATATTAGGCATTTATGTCGCAATAACTAAAGATTTCACACTAATCTCTTTTATAAATCAAACGGCCATTGCAGATAAACACAAGAACCAAATTGCATATATATTTACGCTATGTATAAGTTTGAGTGCAGTATTTTTAATGTCTTCAATTCTAAGCTTTGAATATGATTTTATTGCATTAGCATTCTTATTTTTAACGATAGCATTAT  >  CP000730/2620830‑2621102

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