Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000731 25,669 (A)7→6 coding (322/333 nt) USA300HOU_pUSA300HOUMR0031 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073125,6630A.97.8% 192.6 / ‑1.0 46coding (316/333 nt)USA300HOU_pUSA300HOUMR0031hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (24/20);  total (25/21)
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.04e-01

TGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATA  >  CP000731/25529‑25790
                                                                                                                                       |                                                                                                                                
tGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCATGACAGACAATGGGATCTTTGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                                                                  <  2:219417/135‑1 (MQ=255)
  ttATTTACACTAGGGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTATCTTTATAATTGTAGCACTATTTCTTATGCATTTTAATTACGCTACTCCAAATaaa                                                                                                                                >  2:68146/1‑135 (MQ=255)
   tATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaa                                                                                                                               <  2:60884/135‑1 (MQ=255)
    aTTTACACTAGTTCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGTACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaa                                                                                                                              <  1:200132/135‑1 (MQ=255)
         cacTAGTGCTTAAATTTACAGGTTTGTAACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatata                                                                                                                         >  1:3728/1‑131 (MQ=255)
           cTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTATTATGCTGTTTAATTACTATACTACAAATAAAAAATATACa                                                                                                                      >  2:179941/1‑135 (MQ=255)
               tGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATTCAGAAATTTTCTTTATAATTGTAGTACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTa                                                                                                                  <  1:88434/135‑1 (MQ=255)
               tGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTa                                                                                                                  >  1:63317/1‑135 (MQ=255)
               tGCTGATATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATACAGAAATTTTCTTTATACTTTTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTa                                                                                                                  <  2:181262/135‑1 (MQ=255)
                gCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTCTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTaa                                                                                                                 <  1:115911/135‑1 (MQ=255)
                     aaTTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATCCCTCTAAAGGtt                                                                                                            >  2:202386/1‑135 (MQ=255)
                      aTTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTa                                                                                                           >  2:108134/1‑135 (MQ=255)
                       gttACACGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTAAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTaa                                                                                                          <  1:80285/134‑1 (MQ=255)
                           cAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTATTATACTACAAATAAAAAATATACAACTAAAGGTTAAATTa                                                                                                      <  1:72765/135‑1 (MQ=255)
                           cAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACTAATATAAAATATACATCTAAAGGTTAAATTa                                                                                                      <  1:170313/135‑1 (MQ=255)
                           cAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTa                                                                                                      <  1:34481/135‑1 (MQ=255)
                               tttGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGa                                                                                                  >  1:147066/1‑135 (MQ=255)
                                  gTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGaaaa                                                                                               <  1:196361/135‑1 (MQ=255)
                                    cACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATATAAAGGTTAAATTAATGAAAAAt                                                                                             >  2:64546/1‑135 (MQ=255)
                                       cAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATTTTATAAATAAAAAATATATATCTAAAGGTTAAATTAATGAAAAATAGa                                                                                          <  1:10419/135‑1 (MQ=255)
                                       cAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTATAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGa                                                                                          <  1:27434/135‑1 (MQ=255)
                                        aaGTATTTAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCATTATTTCTTATGCTGTTTAATTACTATATTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATATaa                                                                                         >  2:74192/1‑135 (MQ=255)
                                        aaGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGATTAAATTAATGAAAAATAGaa                                                                                         >  2:6072/1‑135 (MQ=255)
                                         aGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAAt                                                                                        <  1:80342/135‑1 (MQ=255)
                                           tATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTa                                                                                      <  1:222340/135‑1 (MQ=255)
                                           tATTGAAGGATCTTGATCTTCCGCAAATCCATAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACCTCTAAAGGTTAAATTAATGAAAAATAGAATTa                                                                                      <  2:205720/135‑1 (MQ=255)
                                               gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATACAACAGAGATTGTTAAATTGATGAAAAATACAATAagagg                                                                                  >  2:218219/1‑134 (MQ=255)
                                                 aGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTAAAAATAAAAAAAATACATCTAAAGATTAAATTAATGGAAAAAAGAATTAGACAat                                                                                >  1:22971/1‑135 (MQ=255)
                                                    aTCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAtata                                                                             >  2:89507/1‑135 (MQ=255)
                                                                tcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAg                                                                 >  1:148962/1‑135 (MQ=255)
                                                                         aGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtctc                                                        <  1:32897/135‑1 (MQ=255)
                                                                              ttttCTTTATAATTGTAGCACTATTTCTTATGCGGTTTAATAACTATAGTAAAAAGAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGaa                                                   <  2:11576/135‑1 (MQ=255)
                                                                               tttCTTTATAATTGTAGCACTATTTCTTATGCTGATTAATTAAAATAATACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGaaa                                                  >  1:107605/1‑135 (MQ=255)
                                                                                ttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAAc                                                 <  2:63317/135‑1 (MQ=255)
                                                                                 tCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt                                                >  2:61153/1‑135 (MQ=255)
                                                                                 tCTTTATAATTGTAGCACTATTTCTTATACTATTAAAATACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt                                                >  1:56403/1‑135 (MQ=255)
                                                                                  cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACtt                                               >  2:104218/1‑135 (MQ=255)
                                                                                        aaTTGTAGCACTATTTATCATGTGCCAGATTAACTATACTACAAATAAAAAATATACATCTAAATGTTAAATTAATGAAAAATAGAATTATAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTa                                        >  2:133448/1‑135 (MQ=255)
                                                                                               gCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAAtt                                 >  2:106104/1‑135 (MQ=255)
                                                                                                   tATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTaaaa                             <  1:32359/135‑1 (MQ=255)
                                                                                                   tATTTCTTATACTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTaaaa                             <  1:212327/135‑1 (MQ=255)
                                                                                                       tCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTAAAATTAATGAAAAATAGAAATAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTATAAGTAAAAAATAaaaatttt                         >  1:37716/1‑131 (MQ=255)
                                                                                                          tATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTATCTa                      >  2:80645/1‑135 (MQ=255)
                                                                                                              cTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCa                  >  1:29516/1‑135 (MQ=255)
                                                                                                              cTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAACATATAGAAAACAAAAAGGAATTTCTCAAGAAACCATACACGAATGAAAAAAAGTGTCTAGGCa                  >  1:211690/1‑135 (MQ=255)
                                                                                                               tGTTTGATGTCTATACTACTAGTAAACAATAGACATCTAAAGGTTAGATTAATGAAAACTAGAATTAGCGAATATAGAAACCTATAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCaa                 <  2:109120/135‑1 (MQ=255)
                                                                                                                      aacactatactaCAAATAAAAAATATACATCTAAACGTTACATTAATGAAAAATAGAAATAGAGAATATAGAAAACAAAAAGGAATTTCTTAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTAtt           <  1:51129/132‑1 (MQ=255)
                                                                                                                            atactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCa     >  1:207911/1‑135 (MQ=255)
                                                                                                                             tactacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCaa    >  1:159948/1‑135 (MQ=255)
                                                                                                                              actacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAAt   >  1:74075/1‑135 (MQ=255)
                                                                                                                               ctacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATa  <  1:169339/135‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                
TGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATA  >  CP000731/25529‑25790

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