Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 147,963 G→A 7.9% intergenic (+121/‑23) budC → / → USA300HOU_0139 short chain dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730147,9630GA 7.9% 145.4 / 2.7 51intergenic (+121/‑23)budC/USA300HOU_0139short chain dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/19);  new base A (2/2);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

GTGGTATGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAACACCACTATCTAATTTATCTCAATAAAAATTCACAAAGTTATCTCATTTTATTT  >  CP000730/147823‑148102
                                                                                                                                            |                                                                                                                                           
gtggtATGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAg                                                                                                                                             <  2:156796/141‑1 (MQ=255)
gtggtATGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAg                                                                                                                                             <  2:238026/141‑1 (MQ=255)
 tggtATGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGt                                                                                                                                            >  1:28893/1‑141 (MQ=255)
  ggtATGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTg                                                                                                                                           >  2:441246/1‑141 (MQ=255)
      tGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGtt                                                                                                                                       >  1:35867/1‑141 (MQ=255)
            tCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTatgatg                                                                                                                                 >  1:103044/1‑141 (MQ=255)
               gTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTca                                                                                                                              <  1:343151/141‑1 (MQ=255)
                ttAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTcac                                                                                                                             >  2:288058/1‑141 (MQ=255)
                ttAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTcac                                                                                                                             <  1:211808/141‑1 (MQ=255)
                    taatCATCCACTAATGATAAATAAATCCTTATTGTTANGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTc                                                                                                                         >  1:359338/1‑141 (MQ=255)
                     aatCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCt                                                                                                                        <  2:441208/141‑1 (MQ=255)
                     aatCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCt                                                                                                                        <  1:150870/141‑1 (MQ=255)
                     aatCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCt                                                                                                                        <  1:73570/141‑1 (MQ=255)
                        catgcattagcgATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGATTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAATGAGTTATGATGTCACAGTCTAAt                                                                                                                     >  2:404591/12‑141 (MQ=255)
                        catgcattagcgATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGATTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAATGAGTTATGATGTCACAGTCTAAt                                                                                                                     <  2:349116/130‑1 (MQ=255)
                        catgcattagcgATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGATTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAATGAGTTATGATGTCACAGTCTAAt                                                                                                                     >  2:255991/12‑141 (MQ=255)
                        catcCACTAATGATAAATAAATATTTATTGTTAAGTTTTATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTACTGGTAGAAAATTAATAAGCGTAGTTCTTAAGTGAGTTATGATGTCACAGTCTAAt                                                                                                                     <  2:319787/141‑1 (MQ=255)
                                 aTGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGtt                                                                                                            <  1:209936/141‑1 (MQ=255)
                                       aataaatCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGc                                                                                                      <  1:126849/141‑1 (MQ=255)
                                            atCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCattat                                                                                                 >  2:201604/1‑141 (MQ=255)
                                            atCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCattat                                                                                                 <  1:286136/141‑1 (MQ=255)
                                                 tATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTa                                                                                            >  1:372743/1‑141 (MQ=255)
                                                 tATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTa                                                                                            >  2:101910/1‑141 (MQ=255)
                                                     gTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTaa                                                                                        >  1:260327/1‑141 (MQ=255)
                                                      ttAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGATTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAATGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTACTAGTATTAac                                                                                       <  1:335053/141‑1 (MQ=255)
                                                      ttAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGACCATTAATAAGCGAAGTGCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAac                                                                                       <  1:200711/141‑1 (MQ=255)
                                                      ttAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAac                                                                                       <  1:175251/141‑1 (MQ=255)
                                                               aTCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAg                                                                              >  2:216365/1‑141 (MQ=255)
                                                                tCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGa                                                                             <  1:198525/141‑1 (MQ=255)
                                                                    ttAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGAtatt                                                                         >  1:160238/1‑141 (MQ=255)
                                                                     tAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGAtatta                                                                        <  2:189962/141‑1 (MQ=255)
                                                                      aGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGAtattat                                                                       >  1:97279/1‑141 (MQ=255)
                                                                       gCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATtata                                                                      >  2:17904/1‑141 (MQ=255)
                                                                         aGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATaaa                                                                    >  1:306945/1‑141 (MQ=255)
                                                                               gATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATca                                                              <  2:360227/141‑1 (MQ=255)
                                                                                  tttttAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACaa                                                           >  2:222324/1‑141 (MQ=255)
                                                                                    tttAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACaaaa                                                         <  2:256374/141‑1 (MQ=255)
                                                                                         tGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAAcacc                                                    >  1:434235/1‑141 (MQ=255)
                                                                                            aCTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAACACCACt                                                 >  2:356872/1‑141 (MQ=255)
                                                                                                aGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAATATAGTTATGATGTCACAGTCAAATGCATCAGTTTAAAGCAATTTTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAACACCACTATCt                                             <  2:21031/141‑1 (MQ=255)
                                                                                                  aaGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATGTTATAAAAAAACAAAAAAACACCACTATCTaa                                           <  2:326/141‑1 (MQ=255)
                                                                                                  aaGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAACACCACTATCTaa                                           >  2:173664/1‑141 (MQ=255)
                                                                                                       aGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAACACCACTATCTAATTTAt                                      >  1:416537/1‑141 (MQ=255)
                                                                                                              tGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAACACCACTATCTAATTTATCTCAATa                               >  2:40404/1‑141 (MQ=255)
                                                                                                                      aaTTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAACACCACTATCTAATTTATCTCAATAAAAATTca                       >  1:31075/1‑141 (MQ=255)
                                                                                                                          aataaGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCAAAAAAACACCACTAACTAATTTACATCAATAAAAATTCAaaaa                   >  2:415459/1‑141 (MQ=255)
                                                                                                                            taaGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAACACCACTATCTAATTTATCTCAATAAAAATTCACAAAGt                 >  2:28747/1‑141 (MQ=255)
                                                                                                                               gCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAACACCACTATCTAATTTATCTCAATAAAAATTCACAAAGTTAt              <  2:372743/141‑1 (MQ=255)
                                                                                                                                    gTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAACACCACTATCTAATTTATCTCAATAAAAATTCACAAAGTTATCTCat         <  1:101910/141‑1 (MQ=255)
                                                                                                                                     ttCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAACACCACTATCTAATTTATCTCAATAAAAATTCACAAAGTTATCTCatt        >  2:66616/1‑141 (MQ=255)
                                                                                                                                          aaGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAACACCACTATCTAATTTATCTCAATAAAAATTCACAAAGTTATCTCattttatt   >  1:89382/1‑141 (MQ=255)
                                                                                                                                          aaGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAACACCACTATCTAATTTATCTCAATAAAAATTCACAAAGTTATCTCattttatt   >  1:217217/1‑141 (MQ=255)
                                                                                                                                           aGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCGCAAAAACACCACTATCTAATTTATCTCAATAAAAATTCACAAAGTTATCTCattttattt  <  1:9863/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
GTGGTATGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAACACCACTATCTAATTTATCTCAATAAAAATTCACAAAGTTATCTCATTTTATTT  >  CP000730/147823‑148102

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