Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,929,8120GA54.4% 3.7 / 73.4 57intergenic (‑895/‑273)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (17/9);  new base A (15/16);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGAC  >  CP000730/1929675‑1929948
                                                                                                                                         |                                                                                                                                           
ccTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagtag                                                                                                                                          >  2:243121/1‑141 (MQ=255)
      ttCTTAAATTAACAATATCTATTCTCATAGAATTTATCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGc                                                                                                                                    >  1:338193/1‑141 (MQ=255)
        cTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCt                                                                                                                                  <  1:98492/141‑1 (MQ=255)
          tAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTAt                                                                                                                                <  1:379868/141‑1 (MQ=255)
                        cTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCaa                                                                                                                  >  2:24127/1‑141 (MQ=255)
                                      tttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagaga                                                                                                    <  1:133087/141‑1 (MQ=255)
                                          tCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATATAACAAACGTTGATTTAGACGCatattagcg                                                                                                <  2:338282/141‑9 (MQ=255)
                                                      tAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTg                                                                                 >  1:252854/1‑141 (MQ=255)
                                                      tAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTg                                                                                 >  2:59624/1‑141 (MQ=255)
                                                      tAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTg                                                                                 >  2:85903/1‑141 (MQ=255)
                                                       aGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGATTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCATAGAGGAGGTGTA‑T‑AA‑GGTGa                                                                                <  2:172418/141‑1 (MQ=255)
                                                         aCGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTGATg                                                                              <  2:74269/141‑1 (MQ=255)
                                                                           tAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTGATGCTTATTTTCGTTCAcatc                                                            >  2:143265/1‑141 (MQ=255)
                                                                           tAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTGATGCTTATTTTCGTTCAcatc                                                            >  2:186334/1‑141 (MQ=255)
                                                                                  atatAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGGTGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTGATGCTTATTTTCGTTCACATCATAGCAc                                                     >  1:186959/1‑141 (MQ=255)
                                                                                               gAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTgg                                        >  1:62668/1‑141 (MQ=255)
                                                                                               ataggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCCTT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                       <  2:182513/127‑1 (MQ=255)
                                                                                                tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                      >  2:183733/14‑141 (MQ=255)
                                                                                                tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                      >  1:103069/14‑141 (MQ=255)
                                                                                                gaggaagtagatatgCCTGGATATTGACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGACAATT‑‑GACGAAGAGAGGGGGAGCATT‑AC‑CGTGCTAATCATTTTCGTTCACATCATAGCACCAGACATCAGTGGCt                                      >  2:138943/14‑141 (MQ=255)
                                                                                                 cggaagtagatatgCTTGTATATTCACAGCTGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                     >  2:217139/13‑141 (MQ=255)
                                                                                                 cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                     <  1:152687/129‑1 (MQ=255)
                                                                                                   gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                   <  1:57667/131‑1 (MQ=255)
                                                                                                    aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCTTCAGTGTCTGTGc                                  <  2:326714/132‑1 (MQ=255)
                                                                                                    aaggggatatgCGTGTATATTCACAGATGGTATAATAAGAGAGTCGCCTATCTCTCGGGCGTCAATT‑‑GGCGAAGAGAGGAGGTGCA‑TGAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                  <  2:272346/132‑1 (MQ=255)
                                                                                                     agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                 <  1:100559/133‑2 (MQ=255)
                                                                                                     agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                 <  1:59555/133‑2 (MQ=255)
                                                                                                     aGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                  >  2:278177/1‑141 (MQ=255)
                                                                                                      gTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGcc                                 <  2:277437/141‑1 (MQ=255)
                                                                                                      gTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGcc                                 >  2:295535/1‑141 (MQ=255)
                                                                                                       tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                               >  2:246800/7‑138 (MQ=255)
                                                                                                       tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACTTTGCTAATCATTTTCGTTCACATCATAGCACCTGTCATCAGTGGCTGTGctgt                               >  1:62661/7‑138 (MQ=255)
                                                                                                       tagatatgCTTGTATATGCGCAGAAGGTATAATAATAGAGGCGCCTATCTCGCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                               <  1:67919/135‑4 (MQ=255)
                                                                                                       tGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGGCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCa                                >  1:240362/1‑141 (MQ=255)
                                                                                                         gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TGTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                             <  1:299236/137‑6 (MQ=255)
                                                                                                         aCTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                              >  1:363617/1‑141 (MQ=255)
                                                                                                          cTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTg                             <  1:245386/141‑1 (MQ=255)
                                                                                                          cTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTg                             <  2:328488/141‑1 (MQ=255)
                                                                                                          atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                            >  2:292395/4‑135 (MQ=255)
                                                                                                           tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgcg                           >  1:330656/3‑134 (MQ=255)
                                                                                                           tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgcg                           >  1:97334/3‑134 (MQ=255)
                                                                                                           tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                           >  1:65867/3‑141 (MQ=255)
                                                                                                            atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                          >  2:82501/2‑141 (MQ=255)
                                                                                                             tgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGta                         <  1:47551/141‑2 (MQ=255)
                                                                                                               cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                       <  1:380227/141‑4 (MQ=255)
                                                                                                               cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                       >  1:225136/1‑141 (MQ=255)
                                                                                                                 tgtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                     >  1:369883/13‑139 (MQ=255)
                                                                                                                 tGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                     <  1:311268/141‑6 (MQ=255)
                                                                                                                    atatTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTC‑T‑AT‑GGTGCTGCTTATTTTATTTCCCATCATTCTAACAACAATCAGTAGCTGTGGCGTTTCGtttttttt                   >  2:94342/1‑140 (MQ=255)
                                                                                                                         tGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAtt              >  2:153752/1‑141 (MQ=255)
                                                                                                                          acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg            <  2:232418/138‑1 (MQ=255)
                                                                                                                          acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg            <  1:201304/138‑1 (MQ=255)
                                                                                                                           cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc           >  1:130518/3‑141 (MQ=255)
                                                                                                                           cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc           <  1:73559/139‑1 (MQ=255)
                                                                                                                              atggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta        >  2:31345/5‑139 (MQ=255)
                                                                                                                                gatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AC‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg       >  2:10248/1‑141 (MQ=255)
                                                                                                                                 gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa     <  1:120606/140‑1 (MQ=255)
                                                                                                                                   ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa   <  1:157982/141‑3 (MQ=255)
                                                                                                                                     aataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGTCTAAGTAggc  >  2:191871/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                           
CCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGAC  >  CP000730/1929675‑1929948

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