Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 631,628 Δ1 bp 7.8% intergenic (+254/‑248) gph_1 → / → proP Phosphoglycolate phosphatase/Proline/betaine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE631,6280A. 7.8% 170.9 / 8.8 51intergenic (+254/‑248)gph_1/proPPhosphoglycolate phosphatase/Proline/betaine transporter
Reads supporting (aligned to +/- strand):  ref base A (26/20);  new base . (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 = 9.86e-01

TATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATAAAATCATTGAATTTAATTCGATAAAGCGAAATTT  >  USA300TCH1516_ALE/631488‑631765
                                                                                                                                            |                                                                                                                                          
tATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTa                                                                                                                                            <  2:330018/141‑1 (MQ=255)
  tGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACa                                                                                                                                          <  2:265079/141‑1 (MQ=255)
    aaTTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATg                                                                                                                                        >  2:161114/1‑141 (MQ=255)
      tttttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGaa                                                                                                                                     <  2:228213/141‑1 (MQ=255)
       ttttAAAGAGTAATGTCTGACTCGTTGATAATTTATTATTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAg                                                                                                                                    >  1:571212/1‑141 (MQ=255)
        tttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGt                                                                                                                                   >  1:664443/1‑141 (MQ=255)
        tttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGt                                                                                                                                   <  1:646754/141‑1 (MQ=255)
            aaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTAtg                                                                                                                               >  1:1192/1‑141 (MQ=255)
               agTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGc                                                                                                                            <  1:617995/141‑1 (MQ=255)
                 tAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGCATAAACAATTAAAAACGCTGTAGATAATGGCGATAGCATAACGCTTACATAAAGAAATGTGCTa                                                                                                                          >  1:357613/1‑141 (MQ=255)
                  aaTGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTAt                                                                                                                         <  2:275371/141‑1 (MQ=255)
                   aTGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATAGCGATAGCACATAACTTCATGTAAGTGATATGCTATc                                                                                                                        <  2:84043/141‑1 (MQ=37)
                    tGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATAGCGATAGCACATAACTTCATGTAAGTGATATGCTATCg                                                                                                                       <  1:597899/141‑1 (MQ=37)
                     gTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGc                                                                                                                      <  1:223331/141‑1 (MQ=255)
                     gTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGc                                                                                                                      <  1:236795/141‑1 (MQ=255)
                       cTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTa                                                                                                                    >  2:154019/1‑141 (MQ=255)
                            tCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAttttt                                                                                                               <  2:93962/141‑1 (MQ=255)
                            tCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAttttt                                                                                                               <  1:319188/141‑1 (MQ=255)
                                tGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGtt                                                                                                           >  1:649754/1‑141 (MQ=255)
                                 gATAATTTATTTTTGTAAAAATAAATTAAAGTAATGAAAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAATTGGTGCTAGCATATCACTACCATGAAGTTATGTGCTATCGCGATTTTTAGttt                                                                                                          >  1:567767/1‑140 (MQ=255)
                                    aaTTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATaa                                                                                                       <  1:743352/141‑1 (MQ=255)
                                          tttttacaaaaaTAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCaa                                                                                                 >  2:117894/8‑141 (MQ=255)
                                             ttacaaaaaTAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCaaaaa                                                                                              >  1:446441/5‑141 (MQ=255)
                                               gTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGt                                                                                            <  1:341645/141‑1 (MQ=255)
                                               acaaaaaTAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGt                                                                                            >  1:24792/3‑141 (MQ=255)
                                                tAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGtt                                                                                           >  1:294331/1‑141 (MQ=255)
                                                      taaatTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGt                                                                                     >  2:236893/1‑141 (MQ=255)
                                                                 aaTGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGTGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTa                                                                          >  2:131153/1‑141 (MQ=255)
                                                                    gACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGa                                                                       <  2:254616/141‑1 (MQ=255)
                                                                          gTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattatt                                                                 <  1:272112/141‑1 (MQ=255)
                                                                                gAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGtt                                                           >  1:109389/1‑141 (MQ=255)
                                                                                      aaTTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTc                                                     >  1:527081/1‑141 (MQ=255)
                                                                                         tGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCaac                                                  <  1:746412/141‑2 (MQ=255)
                                                                                         tGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCaaa                                                  <  1:145062/141‑1 (MQ=255)
                                                                                                      ttAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGtat                                     >  2:65158/1‑141 (MQ=255)
                                                                                                      ttAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGtat                                     <  2:58103/141‑1 (MQ=255)
                                                                                                        aaaTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATaa                                   >  1:709729/1‑141 (MQ=255)
                                                                                                         aaTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATaaa                                  >  2:80608/1‑141 (MQ=255)
                                                                                                                    aaaaaTAGCGATAGCACATAACTTCATGTAAGTGATATGCTATCGCCATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATAAAATCATTGAAtt                       >  2:164727/2‑141 (MQ=21)
                                                                                                                       aaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATAAAATCATTGAATTTaa                    <  1:302460/141‑1 (MQ=255)
                                                                                                                       aaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATAAAATCATTGAATTTaa                    >  2:208574/1‑141 (MQ=255)
                                                                                                                       aaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATAAAATCATTGAATTTaa                    >  1:312547/1‑141 (MQ=255)
                                                                                                                        aTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATAAAATCATTGAATTTAAt                   <  1:276034/141‑1 (MQ=255)
                                                                                                                            cGATAGCACATAACTTCATGTAAGTGATATGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATAAAATCATTGAATTTAATTCGa               >  1:261818/1‑141 (MQ=37)
                                                                                                                             gATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATAAAATCATTGAATTTAATTCGAt              >  2:320153/1‑141 (MQ=255)
                                                                                                                             gATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATAAAATCATTGAATTTAATTCGAt              >  2:89035/1‑141 (MQ=255)
                                                                                                                              aTAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATAAAATCATTGAATTTAATTCGATa             >  2:310625/1‑141 (MQ=255)
                                                                                                                               tAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATAAAATCATTGAATTTAATTCGATaa            <  1:520871/141‑1 (MQ=255)
                                                                                                                                  cATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATAAAATCATTGAATTTAATTCGATAAAGc         <  1:204595/141‑1 (MQ=255)
                                                                                                                                    tatCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATAAAATCATTGAATTTAATTCGATAAAGCGa       >  2:10362/1‑141 (MQ=255)
                                                                                                                                     atCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATAAAATCATTGAATTTAATTCGATAAAGCGaa      >  1:644689/1‑141 (MQ=255)
                                                                                                                                         cTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATAAAATCATTGAATTTAATTCGATAAAGCGAAAttt  >  1:69764/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                          
TATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTATGAGGGTATAAAATCATTGAATTTAATTCGATAAAGCGAAATTT  >  USA300TCH1516_ALE/631488‑631765

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