Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,929,8621.T53.7% 78.3 / 77.0 41intergenic (‑945/‑223)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (9/10);  new base T (15/7);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCT  >  CP000730/1929725‑1929994
                                                                                                                                          |                                                                                                                                      
aGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑Taa                                                                                                                                     <  1:206301/141‑1 (MQ=255)
aGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑Taa                                                                                                                                     <  1:49618/141‑1 (MQ=255)
 gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑g                                                                                                                                   >  1:74743/1‑141 (MQ=255)
  tgtAGACGATTCAATCTATCAGCTAGAATCATATAACTCATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑gg                                                                                                                                  >  1:119146/1‑141 (MQ=255)
                                             atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑AT‑TGCTAATCATTTTCGTTCACATCATCGCACCAGTCATCAGTGGc                                                                                     >  2:111074/15‑141 (MQ=255)
                                              tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                    >  2:100389/14‑141 (MQ=255)
                                                 gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCACTT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                 >  2:68253/11‑141 (MQ=255)
                                                 gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                 >  2:68252/11‑141 (MQ=255)
                                                 gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                 <  2:185548/131‑1 (MQ=255)
                                                  aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                <  2:185691/132‑1 (MQ=255)
                                                    gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                              >  1:40664/8‑139 (MQ=255)
                                                         tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCACTT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                                         >  2:96276/3‑141 (MQ=255)
                                                         tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                                         >  2:142208/3‑141 (MQ=255)
                                                          gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                                                                         <  1:83854/141‑1 (MQ=255)
                                                                   tattcacagatggtATAATAATAGAGTCACCGATCGCTCTGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAGTCATTTTCGTTCACATCATAGCATCAGTCATCAGTGGCTGTGATGTTGCGTATTATACtt                                                               <  2:165032/133‑1 (MQ=255)
                                                                       cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                           >  2:148565/5‑141 (MQ=255)
                                                                       cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                           >  1:102059/5‑141 (MQ=255)
                                                                         ccGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCGGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgg                                                          <  1:113353/141‑1 (MQ=255)
                                                                            atggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                      <  2:48029/137‑3 (MQ=255)
                                                                             tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                     >  2:97832/4‑138 (MQ=255)
                                                                                tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                  <  2:111736/141‑2 (MQ=255)
                                                                                 ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                 >  2:197019/1‑139 (MQ=255)
                                                                                   aataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCGCCATTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGAc                                                <  1:164047/141‑1 (MQ=255)
                                                                                     tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                              <  1:18467/141‑1 (MQ=255)
                                                                                         agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                         >  1:26210/1‑138 (MQ=255)
                                                                                           agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATTAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                       >  1:71741/1‑136 (MQ=255)
                                                                                           agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                       >  1:90606/1‑136 (MQ=255)
                                                                                           agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                       >  1:89878/1‑136 (MQ=255)
                                                                                               gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                   >  1:116797/1‑132 (MQ=255)
                                                                                                  tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                <  1:89683/141‑13 (MQ=255)
                                                                                                  tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                <  1:215934/141‑13 (MQ=255)
                                                                                                  tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                >  2:64242/1‑129 (MQ=255)
                                                                                                  tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                <  1:152639/141‑13 (MQ=255)
                                                                                                   aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                               >  1:146787/1‑128 (MQ=255)
                                                                                                      tctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGTGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGAc                             >  1:171619/1‑141 (MQ=255)
                                                                                                      tctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGAc                             <  1:49610/141‑1 (MQ=255)
                                                                                                         cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                          <  2:105924/141‑1 (MQ=255)
                                                                                                             cGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGc                      <  1:175054/141‑1 (MQ=255)
                                                                                                                              agagGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATcc     >  1:38663/1‑141 (MQ=255)
                                                                                                                                 ggaggTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCt  >  2:21130/1‑141 (MQ=255)
                                                                                                                                 ggaggTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCt  >  2:162118/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                      
AGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCT  >  CP000730/1929725‑1929994

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