Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,143,110 G→A 77.2% V163V (GTG→GTA USA300HOU_2026 → hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,143,1100GA77.2% 118.6 / 25.8 63V163V (GTG→GTAUSA300HOU_2026hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base G (10/4);  new base A (23/24);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCATATGTTTGGACAAAACGGCTTGCTGTT  >  CP000730/2142974‑2143245
                                                                                                                                        |                                                                                                                                       
aGCTAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCAt                                                                                                                                     <  2:773920/141‑1 (MQ=255)
aGCTAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCAt                                                                                                                                     <  2:20535/141‑1 (MQ=255)
  cTAAGTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCAttt                                                                                                                                   >  2:132006/1‑141 (MQ=255)
    aaTTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATttat                                                                                                                                 >  1:241747/1‑141 (MQ=255)
    aaTTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTATCATttat                                                                                                                                 >  2:418523/1‑141 (MQ=255)
      ttattattGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATttatta                                                                                                                               >  1:752465/1‑141 (MQ=255)
              gCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGTTAATCGTATCGTAGCATTTATTATTGCTaca                                                                                                                       >  2:422566/1‑141 (MQ=255)
                         aCCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTc                                                                                                            >  1:664027/1‑141 (MQ=255)
                          ccTATATTTATCTTTTTTGTATCTAATATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCa                                                                                                           <  2:490756/141‑1 (MQ=255)
                          ccTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTACGttga                                                                                                           >  1:249792/1‑139 (MQ=255)
                           cTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTCTTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGCCATCTATTATTGCTGCATCAGTAAGTTCAt                                                                                                          >  1:738892/1‑141 (MQ=255)
                           cTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCAt                                                                                                          >  1:131769/1‑141 (MQ=255)
                            tataTTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCAtt                                                                                                         >  1:28745/1‑141 (MQ=255)
                             atatTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTa                                                                                                        <  1:233226/141‑1 (MQ=255)
                                tttATCTTCTTTGTATCGATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATAGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAata                                                                                                     <  2:509148/141‑1 (MQ=255)
                                 ttATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAatat                                                                                                    <  2:778884/141‑1 (MQ=255)
                                    tcttctTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTg                                                                                                 >  1:637230/1‑141 (MQ=255)
                                         tttGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAAAAAAATAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATAGCTAAAAAAGTAAGTTCAGTAATATTTCCATTa                                                                                            >  1:749153/1‑141 (MQ=255)
                                            gTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGca                                                                                         >  2:234358/1‑141 (MQ=255)
                                                  attattGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAAt                                                                                   <  2:628896/141‑1 (MQ=255)
                                                          tCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCaa                                                                           >  2:644672/1‑141 (MQ=255)
                                                           ccTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCaaa                                                                          >  1:199205/1‑141 (MQ=255)
                                                              ttattaGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAAttt                                                                       <  2:769790/141‑1 (MQ=255)
                                                              ttattaGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAAttt                                                                       <  1:675876/141‑1 (MQ=255)
                                                                  taGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTc                                                                   >  2:396051/1‑141 (MQ=255)
                                                                  taGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACGACAGTAAGTTCATTAATATTTGCATTATCACATAATGCTTTCACATTTAATc                                                                   >  2:307546/1‑141 (MQ=255)
                                                                   agaagaATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTcc                                                                  >  1:310797/1‑141 (MQ=255)
                                                                     aagaaTATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTCAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAg                                                                <  2:507659/141‑1 (MQ=255)
                                                                     aagaaTATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAg                                                                <  2:231588/141‑1 (MQ=255)
                                                                      agaaTATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGt                                                               >  1:96553/1‑141 (MQ=255)
                                                                         atatGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTa                                                            <  1:111799/141‑1 (MQ=255)
                                                                           atGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTAtt                                                          <  1:768638/141‑1 (MQ=255)
                                                                             gTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTAtttt                                                        <  1:555070/141‑1 (MQ=255)
                                                                             gTATTCAGAAAAGTAATCGTTGGATAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTAtttt                                                        <  2:537996/141‑1 (MQ=255)
                                                                             gTATTCAGAAAAGGAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATAGTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTAtttt                                                        <  2:344020/141‑1 (MQ=255)
                                                                                ttCAGAAAAGTAATCTTTGGAGACTTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTGGCATTAGCACATAATGATTTCACATTTATTCCATTTTATTTTGGt                                                     >  1:636511/1‑141 (MQ=255)
                                                                                ttCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGt                                                     <  1:646567/141‑1 (MQ=255)
                                                                                       aaGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGgtg                                              <  1:666361/141‑1 (MQ=255)
                                                                                       aaGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGgtg                                              <  1:226881/141‑1 (MQ=255)
                                                                                        aGTGATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTAGTGCTATAACAGTAAGTTCATTAGTATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGgtgt                                             <  2:534144/141‑1 (MQ=255)
                                                                                        aGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGgtgt                                             >  1:467576/1‑141 (MQ=255)
                                                                                            aTCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGTTATGGTTTTCAtt                                         >  2:57710/1‑141 (MQ=255)
                                                                                             tCTTTGGAGAATTATTTAATGAGATTGTATTTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCAAGTTTATTTTGGTGTGGGTGTAAttt                                        >  2:749879/1‑141 (MQ=255)
                                                                                              cTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCAtttt                                       >  1:3220/1‑141 (MQ=255)
                                                                                               tttGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCAttttt                                      <  1:503469/141‑1 (MQ=255)
                                                                                                 tGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTc                                    <  1:228117/141‑1 (MQ=255)
                                                                                                  ggAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCa                                   >  1:337396/1‑141 (MQ=255)
                                                                                                      aaTTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAg                               <  2:724345/141‑1 (MQ=255)
                                                                                                          attTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCtta                           >  2:731026/1‑141 (MQ=255)
                                                                                                          attTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAATTTATTTTGGTGTGGGTGTCATTTTTTCATTataata                           >  2:481512/1‑136 (MQ=255)
                                                                                                           tttAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCatat                          >  2:384825/1‑141 (MQ=255)
                                                                                                                ttagactAAAGGTATTCTTATCTTTGTATTTATTATTGATATAACAGTAAGTTCATTAATATTTGCATTAGCACATAATTATTTCTAATTTATTCCAGTTTATTTTGGTATGGGTGTAATTTTTTCATTAGCATATGTTTg                     <  2:428822/135‑1 (MQ=255)
                                                                                                                    aTTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTTATTAGCATATGTTTGGACa                 >  2:588563/1‑141 (MQ=255)
                                                                                                                    aTTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCATATGTTTGGACa                 <  1:165730/141‑1 (MQ=255)
                                                                                                                      tAAAGGTAATCGTATCGTCACATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCATATGTTTGGACaaa               >  1:222640/1‑141 (MQ=255)
                                                                                                                      tAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCATATGTTTGGACaaa               >  1:320660/1‑141 (MQ=255)
                                                                                                                        aaGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCCTTTTTTCATTGGCATATGTTTGGACaaaag             >  1:65648/1‑140 (MQ=255)
                                                                                                                           gTAATCGTATCGTAGCTTTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCATATGTTTGGACAAAACGGc          >  1:385876/1‑141 (MQ=255)
                                                                                                                            tAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCATATGTTTGGACAAAACGGCt         <  2:67933/141‑1 (MQ=255)
                                                                                                                                cgtatcgtAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCATATGTTTGGACAAAACGGCTTGCt     <  1:366568/141‑1 (MQ=255)
                                                                                                                                   atcgtGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCATATGTTTGGACGAAACGGCTTGCTGtt  <  1:358920/141‑1 (MQ=255)
                                                                                                                                   atcgtAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCATATGTTTGGACAAAACGGCTTGCTGtt  <  2:772315/141‑1 (MQ=255)
                                                                                                                                   atcgtAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTAATTTGGTATGGGTGTCATTTTTTCATTAGCATATGTTTGGACAAAACGGCTTGCTGtt  >  2:107920/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                       
AGCTAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCATATGTTTGGACAAAACGGCTTGCTGTT  >  CP000730/2142974‑2143245

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: