Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 401,629 T→G 8.8% T58T (ACT→ACG)  USA300HOU_0370 → possible transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730401,6290TG 8.8% 123.0 / 2.7 45T58T (ACT→ACG) USA300HOU_0370possible transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (16/25);  new base G (2/2);  total (18/27)
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.62e-01

ATGGCTTAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGTGGAATTTTAGGTAATTTTGATGGATTTGGTATTGAGAACTTGACGTTTGC  >  CP000730/401490‑401768
                                                                                                                                           |                                                                                                                                           
aTGGCTTAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTg                                                                                                                                            <  1:374419/141‑1 (MQ=255)
      tAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTtta                                                                                                                                      <  1:279326/141‑1 (MQ=255)
      tAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTAGAAACTGTTTtta                                                                                                                                      >  2:62302/1‑141 (MQ=255)
       aaaCCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTttat                                                                                                                                     <  2:103806/141‑1 (MQ=255)
         atcAAACGCAGCTTGCCAAACTAGCTGGTGTTTCAAGACAAACCATTTCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAGATTGCTCGCATTTTCAATGAAACGGTGGAAACGGTTTttatta                                                                                                                                   <  1:39702/139‑1 (MQ=255)
          ccAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTtattat                                                                                                                                  <  1:103861/141‑1 (MQ=255)
             aaCGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTga                                                                                                                               <  2:261695/141‑1 (MQ=255)
               cGCAACTTGCTACACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTgagg                                                                                                                             <  2:366581/141‑1 (MQ=255)
                   aCTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAggagga                                                                                                                         >  1:188210/1‑141 (MQ=255)
                       gCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTTATCGATACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAACCGGGGGAAACTGTTTTTATTATTGAGGAGGATGAg                                                                                                                     >  2:386808/1‑141 (MQ=255)
                       gCTAAACAAGCGGGCGTTCCAAGACAAACCATATCGCTAATTGATCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAg                                                                                                                     >  2:134054/1‑141 (MQ=255)
                       gCTAAACAAGCGGGCGTTACAAGACAAACCATATCGCTACTTGAGCGAAACAATTTTATGCCATCGGTATTTACGACAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAg                                                                                                                     <  1:386641/141‑1 (MQ=255)
                           aaCATGCGGGCCTACAAAAACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAGCGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCAt                                                                                                                 <  2:116309/141‑1 (MQ=255)
                             cAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGa                                                                                                               <  2:59261/141‑1 (MQ=255)
                                 cGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAAt                                                                                                           >  2:246349/1‑141 (MQ=255)
                                  gggCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATa                                                                                                          >  2:106301/1‑141 (MQ=255)
                                       tttCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGAATGAGGCATGAAAATACTAAg                                                                                                     >  1:241879/1‑141 (MQ=255)
                                              aCAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATc                                                                                              >  1:306128/1‑141 (MQ=255)
                                                 aaCCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGa                                                                                           <  1:292928/141‑1 (MQ=255)
                                                   ccATATCGCTAATTGAGCGCAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAACCTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGata                                                                                         <  1:7697/141‑1 (MQ=255)
                                                     atatCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATc                                                                                       <  2:185775/141‑1 (MQ=255)
                                                      tatCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCt                                                                                      <  1:188822/141‑1 (MQ=255)
                                                       atCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGatatccgatatatctt                                                                                     <  1:356409/141‑13 (MQ=255)
                                                       atCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCtt                                                                                     <  1:167919/141‑1 (MQ=255)
                                                          gctcaattaGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTa                                                                                  >  2:286487/9‑141 (MQ=255)
                                                          gCTAATTGAGAGAAACAATTTTATGCCATCAGTATAAAAAAAAAAAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTa                                                                                  <  2:260655/141‑1 (MQ=255)
                                                           cTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTAc                                                                                 <  1:167699/141‑1 (MQ=255)
                                                           cTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTAc                                                                                 <  1:408292/141‑1 (MQ=255)
                                                                    cGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGAc                                                                        <  2:371862/141‑1 (MQ=255)
                                                                        aCAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGCGGACTTGt                                                                    >  1:61711/1‑141 (MQ=255)
                                                                              ttATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAggggg                                                              >  2:80985/1‑141 (MQ=255)
                                                                               tATGCCATCAGTATTAACGGCAATAAAGATTGCTCGCATTTTCAATGAAACGGTGGAAACGGTTTTTATTATTGAGGAGGATGAAGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGAATTGTAGGGGGt                                                             >  2:337179/1‑141 (MQ=255)
                                                                                        aGTATTAACGGCAATAAAGATTGCTCGCATTTTCAATGAAACGGTGGAAACGGTTTTTATTATTGAGGAGGATGAAGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGAATTGTAGGGGGTATCATTGGt                                                    <  1:419243/141‑1 (MQ=255)
                                                                                        aGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGt                                                    <  1:25278/141‑1 (MQ=255)
                                                                                        aGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGt                                                    <  1:25801/141‑1 (MQ=255)
                                                                                          tATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGTgg                                                  <  2:256725/141‑1 (MQ=255)
                                                                                              aaCGGCAATAAAGATTGCTCGCATTTTCAATGAAACGGTGGAAACGGTTTTTATTATTGAGGAGGATGAAGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGAATTGTAGGGGGTATCATTGGTGGAAtt                                              >  1:132/1‑141 (MQ=255)
                                                                                                     aTAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGTGGAATTTTAGGTa                                       >  1:438089/1‑141 (MQ=255)
                                                                                                         aaaTTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGTGGAATTTTAGGTAAttt                                   >  1:318580/1‑141 (MQ=255)
                                                                                                           aTTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGTGGAATTTTAGGTAATTTTg                                 >  2:21190/1‑141 (MQ=255)
                                                                                                                  gCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGTGGAATTTTAGGTAATTTTGATGGAtt                          <  1:176582/141‑1 (MQ=255)
                                                                                                                                   gTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGTGGAATTTTAGGTAATTTTGATGGATTTGGTATTGAGAACTTGa         >  1:446694/1‑141 (MQ=255)
                                                                                                                                      gAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGTGGAATTTTAGGTAATTTTGATGGATTTGGTATTGAGAACTTGACGt      <  1:260481/141‑1 (MQ=255)
                                                                                                                                        aaCTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATCTCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGTGGAATTTTAGGTAATTTTGATGGATTTGGTATTGAGAACTTGACGttt    >  1:351243/1‑141 (MQ=255)
                                                                                                                                          cTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGTGGAATTTTAGGTAATTTTGATGGATTTGGTATTGAGAACTTGACGTTTGc  <  2:257843/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                           
ATGGCTTAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGTGGAATTTTAGGTAATTTTGATGGATTTGGTATTGAGAACTTGACGTTTGC  >  CP000730/401490‑401768

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