Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 147,927 G→T 9.4% intergenic (+85/‑59) 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,9270GT 9.4% 163.3 / 2.9 52intergenic (+85/‑59)budC/USA300HOU_0139short chain dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/20);  new base T (2/3);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

ACATTACTGGACAAACAATTATTGTAGATGGTGGTATGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAACACCACTATCTAATTTAT  >  CP000730/147793‑148066
                                                                                                                                      |                                                                                                                                           
acattacTGGACAAACAATTATTGTAGATGGTGGTATGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTAtt                                                                                                                                       >  2:175251/1‑141 (MQ=255)
          acaaacaaTTATTGTAGATGGTGGTATGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAAt                                                                                                                             >  2:126849/1‑141 (MQ=255)
            aaacaaTTATTGTAGATGGTGGTATGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGTATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATta                                                                                                                           <  2:424972/141‑1 (MQ=255)
                aaTTATTGTAGATGGTGGTATGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATtaataa                                                                                                                       >  1:253637/1‑141 (MQ=255)
                     ttGTAGATGGTGGTATGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGaa                                                                                                                  >  2:328688/1‑141 (MQ=255)
                       gTAGATGGTGGTATGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGt                                                                                                                >  1:326/1‑141 (MQ=255)
                        tAGATGGTGGTATGAGATTCCGTTAATAATCATGCATTAGCGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGATTGTATTGGTAGAAAATTAATAAGCGAAGtt                                                                                                               <  2:404423/141‑1 (MQ=255)
                         aGATGGTGGTATGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGTAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTc                                                                                                              >  1:329016/1‑141 (MQ=255)
                         aGATGGTGGTATGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGGAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTc                                                                                                              >  1:98781/1‑141 (MQ=255)
                              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:349116/130‑1 (MQ=255)
                                                      catgcattagcgATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGATTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAATGAGTTATGATGTCACAGTCTAAt                                                                                 >  2:404591/12‑141 (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                                                        >  2:101910/1‑141 (MQ=255)
                                                                               tATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTa                                                        >  1:372743/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)
                                                                                                                                      |                                                                                                                                           
ACATTACTGGACAAACAATTATTGTAGATGGTGGTATGAGATTCCGTTAATAATCATCCACTAATGATAAATAAATCCTTATTGTTAAGTTTAATCACTTAGCAGTAAGGATTTTTTAGTGCACTTAGAAGGGAGTGTATTGGTAGAAAATTAATAAGCGAAGTTCTTAAGTGAGTTATGATGTCACAGTCTAATGCATCAGTTGAAAGCATTATTAGTATTAACACACCCAAGATATTATAAAACATCACAAAAACACCACTATCTAATTTAT  >  CP000730/147793‑148066

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: ‑