Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 623,744 A→G 8.8% I242V (ATT→GTT)  USA300HOU_0558 → glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730623,7440AG 8.8% 128.3 / 4.0 46I242V (ATT→GTT) USA300HOU_0558glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base G (2/2);  total (22/24)
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.09e-01

TCAACAATTTTTTGATACAGAAGATGAATTAGTTCAATATTTTCTCCATCAATTATGTAAAAATAATGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTAACACGTTATAAAGCGATTGTTGTATCAACAGAAAAGCAATGCCAAGATATTTCACAA  >  CP000730/623605‑623881
                                                                                                                                           |                                                                                                                                          
tCAACAATTTTTTGATACAGAAGATGAATTAGTTCAATATTTTCTCCATCAATTATGTAAAAATAATGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTAt                                                                                                                                           <  2:303910/141‑1 (MQ=255)
   acaaTTTTTTGATACAGAAGATGAATTAGTTCAATATTTTCTCCATCAATTATGTAAAAATAATGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTAttgt                                                                                                                                        >  1:88099/1‑141 (MQ=255)
   acaaTTTTTTGATACAGAAGATGAATTAGTTCAATATTTTCTCCATCAATTATGTAAAAATAATGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTAttgt                                                                                                                                        <  1:472611/141‑1 (MQ=255)
            tGATACAGAAGATGAATTAGTTCAATATTTTCTCCATCAATTATGTAAAAATAATGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCa                                                                                                                               >  2:238916/1‑141 (MQ=255)
                        tGAATTAGTTCAATATTTTCTCCATCAATTATGTAAAAATAATGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTAATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACAttt                                                                                                                   <  1:75445/141‑1 (MQ=255)
                                 tCAATATTTTCTCCATCAATTATGTAAAAATAATGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGc                                                                                                          <  2:398289/141‑1 (MQ=255)
                                             ccATCAATTATGTAAAAATAATGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTCTTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGtat                                                                                              <  1:585603/141‑1 (MQ=255)
                                             ccATCAATTATGTAAAAATAATGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGtat                                                                                              >  1:312716/1‑141 (MQ=255)
                                             ccATCAATTATGTAAAAATAATGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGtat                                                                                              <  2:552060/141‑1 (MQ=255)
                                               aTCAATTATGTAAAAATAATGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATaa                                                                                            <  1:647888/141‑1 (MQ=255)
                                                     tATGTAAAAATAATGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGtt                                                                                     <  2:306951/141‑1 (MQ=255)
                                                         tAAAAATAATGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTta                                                                                 >  2:143314/1‑141 (MQ=255)
                                                         tAAAAATAATGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTta                                                                                 <  1:666582/141‑1 (MQ=255)
                                                         tAAAAATAATGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATACCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTta                                                                                 >  2:143497/1‑141 (MQ=255)
                                                                  tGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAgt                                                                        >  2:433017/1‑141 (MQ=255)
                                                                   gATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAgtg                                                                       >  1:649463/1‑141 (MQ=255)
                                                                   gATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAgtg                                                                       <  1:7771/141‑1 (MQ=255)
                                                                    aTCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAgtgt                                                                      <  2:376738/141‑1 (MQ=255)
                                                                     tCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGtt                                                                     >  1:393736/1‑141 (MQ=255)
                                                                          tCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTtaata                                                                <  1:563704/141‑1 (MQ=255)
                                                                           cATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTGATGGTATAAAAAGTTTTTATAAAACAGTGTTtaataa                                                               <  2:502066/141‑1 (MQ=255)
                                                                            atatTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACGTTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTaataat                                                              >  1:398296/1‑141 (MQ=255)
                                                                             tatTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAAtt                                                             >  2:179767/1‑141 (MQ=255)
                                                                                tAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTaa                                                          <  1:141737/141‑1 (MQ=255)
                                                                                tAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTaa                                                          <  1:349567/141‑1 (MQ=255)
                                                                                           cATGAATTAGGAAATGTTATAGCGGGATTAAATCAACGTATTCCAGTTACTGTTGCACCCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTAACACGTTATaaa                                               <  1:709380/141‑1 (MQ=255)
                                                                                           cATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGtttt                                                                                   >  2:583513/1‑105 (MQ=255)
                                                                                           cATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGtttt                                                                                   <  1:583513/105‑1 (MQ=255)
                                                                                           cATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTAACACGTTATaaa                                               <  2:531403/141‑1 (MQ=255)
                                                                                                aTTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTAACACGTTATAAAGCGAt                                          >  2:577780/1‑141 (MQ=255)
                                                                                                   aGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTAACACGTTATAAAGCGAttgt                                       >  2:213425/1‑141 (MQ=255)
                                                                                                    ggAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTAACACGTTATAAAGCGAttgtt                                      <  1:389773/141‑1 (MQ=255)
                                                                                                        aTGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTAACACGTTATAAAGCGATTGTTGTAt                                  >  2:409119/1‑141 (MQ=255)
                                                                                                          gTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAATAGTGTTTAATAATTTAACACGTTATAAAGCGATTGTTGTATCa                                >  2:660026/1‑141 (MQ=255)
                                                                                                              taGCGGGATTAAATCAAAGTATTCCAGTTGTTGTTGTGCTCCATAGTACACATTTATCCGGTACCGGTaa                                                                                                    <  2:241532/70‑1 (MQ=255)
                                                                                                              taGCGGGATTAAATCAAAGTATTCCAGTTGTTGTTGTGCTCCATAGTACACATTTATCCGGTACCGGTaa                                                                                                    <  2:241400/70‑1 (MQ=255)
                                                                                                              taGCGGGATTAAATCAAAGTATTCCAGTTGTTGTTGTGCTCCATAGTACACATTTATCCGGTACCGGTaa                                                                                                    >  1:241532/1‑70 (MQ=255)
                                                                                                              taGCGGGATTAAATCAAAGTATTCCAGTTGTTGTTGTGCTCCATAGTACACATTTATCCGGTACCGGTaa                                                                                                    >  1:241400/1‑70 (MQ=255)
                                                                                                                    gATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTAACACGTTATAAAGCGATTGTTGTATCAACAGAAAAGc                      <  1:442843/141‑1 (MQ=255)
                                                                                                                    gATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTAACACGTTATAAAGCGATTGTTGTATCAACAGAAAAGc                      >  2:149170/1‑141 (MQ=255)
                                                                                                                          aTCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTAACACGTTATAAAGCGATTGTTGTATCAACAGAAAAGCAATGcc                <  2:623585/141‑1 (MQ=255)
                                                                                                                           tCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTTTTTTTTTTTTTTTTAATAATTTAACACGTTATAAAGCGATTGTTGTATCAGCAGAAAAGCAATGcc                >  2:488802/1‑141 (MQ=255)
                                                                                                                                 tATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTAACACGTTATAAAGCGATTGTTGTATCAACAGAAAAGCAATGCCAAGatat         >  1:338377/1‑141 (MQ=255)
                                                                                                                                   ttCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTAACACGTTATAAAGCGATTGTTGTATCAACAGAAAAGCAATGCCAAGATAttt       >  1:525399/1‑141 (MQ=255)
                                                                                                                                     ccAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTAACACGTTATAAAGCGATTGTTGTATCAACAGAAAAGCAATGCCAAGATATTTca     >  1:15897/1‑141 (MQ=255)
                                                                                                                                        gTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTAACACGTTATAAAGCGATTGTTGTATCAACAGAAAAGCAATGCCAAGATAttgtgaaa  <  2:51104/141‑7 (MQ=255)
                                                                                                                                           |                                                                                                                                          
TCAACAATTTTTTGATACAGAAGATGAATTAGTTCAATATTTTCTCCATCAATTATGTAAAAATAATGATCAAATCATATTAGATCGTCCTCATGAATTAGGAAATGTTATAGCGGGATTAAATCAAAGTATTCCAGTTATTGTTGTACTCCATAGTACACATTTATCCGGTGCCGGTAATGGTATAAAAAGTTTTTATAAAACAGTGTTTAATAATTTAACACGTTATAAAGCGATTGTTGTATCAACAGAAAAGCAATGCCAAGATATTTCACAA  >  CP000730/623605‑623881

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