Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,929,8610TC48.0% ‑1.5 / 67.4 50intergenic (‑944/‑224)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/12);  new base C (12/12);  total (26/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 = 9.03e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACC  >  CP000730/1929728‑1930001
                                                                                                                                     |                                                                                                                                              
gtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑t                                                                                                                                       >  2:129423/1‑141 (MQ=255)
     cGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGc                                                                                                                                  <  1:203316/141‑1 (MQ=255)
       aTTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCtt                                                                                                                                <  1:31467/141‑1 (MQ=255)
             tCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTNTTTTc                                                                                                                          >  1:32057/1‑141 (MQ=255)
             tCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTc                                                                                                                          <  2:33343/141‑1 (MQ=255)
                aTCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGtt                                                                                                                       >  2:11119/1‑141 (MQ=255)
                                taACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAg                                                                                                       <  2:306075/141‑1 (MQ=255)
                                     tATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGtcatc                                                                                                  <  2:190677/141‑1 (MQ=255)
                                          atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                            <  1:300624/127‑1 (MQ=255)
                                           tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                           >  2:331506/14‑141 (MQ=255)
                                            aGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                           >  1:104041/1‑141 (MQ=255)
                                              gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                        >  2:256877/11‑141 (MQ=255)
                                               aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                       <  1:138632/132‑1 (MQ=255)
                                                agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtcat                                                                                      >  1:346066/9‑138 (MQ=255)
                                                 gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTAAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTAATCAGTGGCTGTGctg                                                                                     >  1:291339/8‑139 (MQ=255)
                                                         gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                             <  1:150593/141‑3 (MQ=255)
                                                          cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                            <  1:5346/141‑4 (MQ=255)
                                                          cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                            <  2:226543/141‑4 (MQ=255)
                                                           ttGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttt                                                                            <  1:309807/141‑1 (MQ=255)
                                                           ttGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttt                                                                            >  1:247717/1‑141 (MQ=255)
                                                           ttGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTACGACATCTATTTAGCCAGAGACCGTGGCTCCA‑TGTGT‑TCTGGCTTATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttt                                                                            <  1:377989/141‑1 (MQ=255)
                                                            tGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                                          <  2:246795/141‑6 (MQ=255)
                                                               atatTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCCTTGCGTTTTTTTc                                                                        >  2:188694/1‑141 (MQ=255)
                                                                tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatactt                                                                      >  2:233152/9‑135 (MQ=255)
                                                                 atTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCCTTGCGTTTTTTTCtt                                                                      <  2:18674/141‑1 (MQ=255)
                                                                   tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                   <  1:2912/136‑1 (MQ=255)
                                                                      cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                <  2:24079/139‑1 (MQ=255)
                                                                          tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                            >  1:66995/4‑141 (MQ=255)
                                                                            gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                          >  2:66986/2‑141 (MQ=255)
                                                                              ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                        >  1:162890/1‑139 (MQ=255)
                                                                              ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                        >  2:349056/1‑139 (MQ=255)
                                                                              ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                        <  2:375338/141‑3 (MQ=255)
                                                                                 ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                     <  2:167903/137‑1 (MQ=255)
                                                                                 ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                     >  1:45614/5‑141 (MQ=255)
                                                                                     tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                 <  1:364281/141‑3 (MQ=255)
                                                                                     tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                 >  1:350120/1‑139 (MQ=255)
                                                                                       gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                >  1:176270/1‑141 (MQ=255)
                                                                                           cGCCTATCTCTCAGGCGTCAATTTCGTCGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaaaa                                            <  2:114941/141‑1 (MQ=255)
                                                                                            gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                          >  1:302587/1‑132 (MQ=255)
                                                                                                      cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                                 >  2:12750/1‑141 (MQ=255)
                                                                                                      cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                                 >  2:249047/1‑141 (MQ=255)
                                                                                                       aGGCGTCCATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACatat                                >  1:351445/1‑141 (MQ=255)
                                                                                                       aGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACTAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACatat                                >  1:163387/1‑141 (MQ=255)
                                                                                                       aGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACatat                                <  2:124022/141‑1 (MQ=255)
                                                                                                        ggCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACAtata                               <  2:376128/141‑1 (MQ=255)
                                                                                                          cGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGc                             >  1:142037/1‑141 (MQ=255)
                                                                                                                 ttAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCa                      >  2:195193/1‑141 (MQ=255)
                                                                                                                   aGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt                    <  1:195193/141‑1 (MQ=255)
                                                                                                                     aCGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaa                  <  2:298629/141‑1 (MQ=255)
                                                                                                                                aggTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCa       >  2:54310/1‑141 (MQ=255)
                                                                                                                                    gtA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTAc   >  1:308256/1‑141 (MQ=255)
                                                                                                                                     ta‑taAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTAcc  <  2:176270/141‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                              
GTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACC  >  CP000730/1929728‑1930001

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: ‑