Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000731 25,663 Δ1 bp 100% coding (316/333 nt) USA300HOU_pUSA300HOUMR0031 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073125,6630A.94.3% 215.9 / 1.6 52coding (316/333 nt)USA300HOU_pUSA300HOUMR0031hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base . (28/21);  minor base G (1/1);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACAAT  >  CP000731/25529‑25799
                                                                                                                                        |                                                                                                                                        
tGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaat                                                                                                                                    >  1:49900/1‑140 (MQ=255)
   tATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatata                                                                                                                                 <  1:90345/141‑5 (MQ=255)
    aTTTACACTAGTGCTTAAATTTACAGGTTTGTCACCCCGTA‑‑AACGGGTCAGGGGGGTCCAAGTCACCATATATATGGTCTTTATAATTGTAGAACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATAc                                                                                                                               >  1:55830/1‑141 (MQ=255)
     tttACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACa                                                                                                                              >  2:53869/1‑141 (MQ=255)
      ttACACTAGTGCTTAAATTTACAGGTTTGTCATCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACAt                                                                                                                             <  2:167036/141‑1 (MQ=255)
             aatGCTTAAATTTAAAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATGGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAgg                                                                                                                      <  2:152881/139‑1 (MQ=255)
                gCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTa                                                                                                                   >  1:126061/1‑141 (MQ=255)
                 cTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTaa                                                                                                                  <  1:184278/141‑1 (MQ=255)
                        ttACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAACAAATATACATCTAAAGGTTAAATTAATg                                                                                                           <  1:118815/141‑1 (MQ=255)
                          aCAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGaa                                                                                                         <  2:54787/141‑1 (MQ=255)
                             agTTTGTCACCAAGTATTAAAGAATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TATTTATAATTGTAGCAATATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGaaaaa                                                                                                      <  2:171760/140‑1 (MQ=255)
                                ttGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTCTTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAg                                                                                                   >  2:76581/1‑141 (MQ=255)
                                        aaGTATTGAAGGATATTTGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTGTGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTagag                                                                                           <  2:208786/141‑1 (MQ=255)
                                          gTATTGAAGGATCTTTGTTGTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATGTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                                                                         <  1:116993/141‑1 (MQ=255)
                                               gAAGGATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATCAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGCAtataa                                                                                    >  2:31091/1‑140 (MQ=255)
                                               gAAGGATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGATAAtataa                                                                                    >  2:37267/1‑140 (MQ=255)
                                               gAAGGATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                                                                    >  1:106954/1‑141 (MQ=255)
                                               aaaGGAACTTGGTTTTTTTAATATACAGAAATTT‑‑TATTTATAATTGTAGCAATATTTCTTATGATGTTTAATTACTAAAATAAAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                                                                    <  2:89008/140‑1 (MQ=255)
                                                aatgATCTTTGTTTTCTCAATATACAGAAATTT‑‑TATTTATGAATGTAGGAATAATTATTATGCTGTTTAATTACTATAATAAAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGa                                                                                   <  2:167158/138‑1 (MQ=255)
                                                aaGGATCTTAGATTACTCAATATCAAGAAATTT‑‑TCTTTATAATTATAGAACTATTTCTTATGCTTTATAATTACTATAATAAAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGa                                                                                   <  2:38528/141‑1 (MQ=255)
                                                 aGGATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaa                                                                                  >  1:10909/1‑141 (MQ=255)
                                                    aTCGTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAAc                                                                               <  1:202069/141‑1 (MQ=255)
                                                        tGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaa                                                                           >  1:131615/1‑141 (MQ=255)
                                                        tGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATAAGAGAATATAGAAAACaaaa                                                                           >  2:53139/1‑141 (MQ=255)
                                                        tGGTTTTCTCAATATCCAGAAATTT‑‑TATTTAAAATTGTAGCACTATGTCTTATGCTGTATAATTACTATACTACAAATAAAAAATAAACATCTAAAGGTTAAATTAATGAAAACTAGAATTAGAGCATATAAAACACaaaa                                                                           >  2:64426/1‑141 (MQ=255)
                                                         ggTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaaa                                                                          >  1:22040/1‑141 (MQ=255)
                                                         ggTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaaa                                                                          >  1:200640/1‑141 (MQ=255)
                                                         ggTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaaa                                                                          >  2:147140/1‑141 (MQ=255)
                                                           ttttCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTCGCGAATATAGAAAACAAAAAgg                                                                        >  1:201269/1‑141 (MQ=255)
                                                             ttCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGaa                                                                      <  1:124285/141‑1 (MQ=255)
                                                                tcAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAAttt                                                                   >  2:62080/1‑141 (MQ=255)
                                                                 cAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaaa                                                                          >  2:143472/1‑133 (MQ=255)
                                                                 cAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaaa                                                                          <  1:143472/133‑1 (MQ=255)
                                                                        cAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTCCAAAGAAAAAATATACCTCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGaa                                                           >  2:96688/1‑141 (MQ=255)
                                                                            aaTTT‑‑TCTTTATAATTGAAGCACTATTTCTTATACTGTTTAATTACTATACTAAAAATAAAAAATATACATCTAAATGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACtt                                                       <  2:61011/141‑1 (MQ=255)
                                                                              ttt‑‑tATTTATAATTGTAGCACTATTTATTATGCTGTTTAATTACTCGACGACAAAAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTa                                                     <  2:49900/141‑1 (MQ=255)
                                                                                         taATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGaa                                            >  1:47241/1‑141 (MQ=255)
                                                                                           atagTAGAACTATTTATTATGCTGTTAAATTACTATACAAAAAATAAAAAATATAAATATAAAAGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAAtt                                          <  2:132874/138‑1 (MQ=255)
                                                                                            ttGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTa                                         <  2:22040/141‑1 (MQ=255)
                                                                                            ttGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACATTAAGTAAAGAATTa                                         >  2:194921/1‑141 (MQ=255)
                                                                                              gTAGCACTATTTTTTATGCTATAAAAATAAGATAATACAAAGAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTaaa                                       <  1:42977/141‑1 (MQ=255)
                                                                                              gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAACGGTTAAATTAATGAAAAATAGATTTAGATAATAGAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTaaa                                       >  1:101803/1‑141 (MQ=255)
                                                                                              gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTaaa                                       >  2:9933/1‑141 (MQ=255)
                                                                                                   aCTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAGTAAAGAATTAAAAGTGTc                                 >  1:24336/1‑141 (MQ=255)
                                                                                                       tttCTCATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAg                              >  2:233414/1‑141 (MQ=255)
                                                                                                         tCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTGAAAGTGTCTAGGc                            <  1:113890/141‑1 (MQ=255)
                                                                                                         tCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGc                            <  1:147140/141‑1 (MQ=255)
                                                                                                            tATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGACAAAATATATGAATATCTCCAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCaaa                         >  1:233065/1‑141 (MQ=255)
                                                                                                               gCTGTTTAATGACTATAATACAAATAAAAAATATACATCTAAAGGTTAAACTAATGACAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGTAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTa                      <  1:197547/141‑1 (MQ=255)
                                                                                                                 tGTTTAATTACTATACTACAAATAAAAAATATTTATCTGCGGCTGAAAATAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTAtt                    <  1:161189/141‑1 (MQ=255)
                                                                                                                   tttAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGGATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTaa                  >  1:169218/1‑141 (MQ=255)
                                                                                                                       aTTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCa              >  2:224781/1‑141 (MQ=255)
                                                                                                                            tatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGa         >  1:191725/1‑141 (MQ=255)
                                                                                                                                   cAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACaat  >  2:100732/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                        
TGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTT‑‑TCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACAAT  >  CP000731/25529‑25799

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