Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073076,2950TC53.3% 20.9 / 82.1 60intergenic (+353/‑781)USA300HOU_0072/USA300HOU_0075universal stress protein/lysophospholipase
Reads supporting (aligned to +/- strand):  ref base T (16/12);  new base C (16/16);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATAAAAAGCTAATT  >  CP000730/76158‑76424
                                                                                                                                         |                                                                                                                                 
aGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTAt                                                                                                                                >  2:39931/1‑141 (MQ=1)
          gCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGatatat                                                                                                                      >  2:19722/1‑141 (MQ=35)
                    atttataatataaTTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATAtt                                                                                                            >  2:62295/4‑141 (MQ=0)
                       tataatataaTTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAg                                                                                                         >  2:23897/1‑141 (MQ=17)
                       tataatataaTTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAg                                                                                                         <  2:75326/141‑1 (MQ=17)
                       tataatataaATATAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAg                                                                                                         <  2:75340/141‑1 (MQ=255)
                        ataatataATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTCATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAga                                                                                                        <  1:28309/137‑2 (MQ=255)
                         taGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAggg                                                                                                       >  2:77493/1‑141 (MQ=17)
                              taATTTTAAAAAAAAGGAGTCTTCTGTCTGAACTATCTCCGACCTCAACCATACAACAACGATGTTAAGACCGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACgt                                                                                                  <  1:51955/141‑1 (MQ=11)
                                  tttAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggtg                                                                                              >  2:58495/1‑141 (MQ=17)
                                  tttAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACgtggtg                                                                                              <  2:49927/141‑1 (MQ=12)
                                   ttAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                                                                                             <  1:7325/141‑1 (MQ=17)
                                   ttAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                                                                                             <  1:7316/141‑1 (MQ=17)
                                    tAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTa                                                                                            >  2:11479/1‑141 (MQ=17)
                                     aaCAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTaa                                                                                           >  1:42017/1‑141 (MQ=32)
                                      aCAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTGCTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTaaa                                                                                          >  2:69722/1‑141 (MQ=12)
                                         aaaaGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGt                                                                                       >  2:77633/1‑141 (MQ=12)
                                         aaaaGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCNGAAATATTAAGANAACGTGGTGTAAACGt                                                                                       <  1:51662/141‑1 (MQ=12)
                                             ggAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTtcat                                                                                   >  2:56307/1‑141 (MQ=17)
                                                gTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTcatcat                                                                                >  1:20819/1‑141 (MQ=32)
                                                 tcttctGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCAtt                                                                               >  1:48721/1‑141 (MQ=12)
                                                 tcttctGTATGAACTATTTCAGAGATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCAtt                                                                               >  2:71814/1‑141 (MQ=12)
                                                             aCTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTAcc                                                                   <  2:7565/141‑1 (MQ=17)
                                                               tATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGt                                                                 <  1:39931/141‑1 (MQ=17)
                                                               tATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATACTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGt                                                                 <  2:73196/141‑1 (MQ=17)
                                                               tATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGt                                                                 <  1:14636/141‑1 (MQ=255)
                                                                aTTTCAGATATAAACAATTGAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAAGATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGtt                                                                <  1:35341/141‑1 (MQ=11)
                                                                 tttCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTg                                                               >  2:64597/1‑141 (MQ=255)
                                                                   tCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTggg                                                             <  1:56307/141‑1 (MQ=17)
                                                                      gATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTc                                                          >  2:2733/1‑141 (MQ=255)
                                                                          taAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGa                                                      >  2:71865/1‑141 (MQ=17)
                                                                                 tttAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGcc                                               <  1:77030/141‑1 (MQ=255)
                                                                                 tttAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGcc                                               <  1:23897/141‑1 (MQ=255)
                                                                                   tAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGcccc                                             >  2:84235/1‑141 (MQ=12)
                                                                                     aCAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCaa                                           >  1:79555/1‑141 (MQ=17)
                                                                                     aCAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCaa                                           >  2:8250/1‑141 (MQ=17)
                                                                                     aCAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGGTCAAGAATATGCCCCaa                                           >  2:23480/1‑141 (MQ=17)
                                                                                             gTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCATTTTtatgt                                   >  2:3249/1‑139 (MQ=21)
                                                                                             gTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTatat                                   <  1:30770/141‑1 (MQ=255)
                                                                                             gTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTatat                                   >  2:63918/1‑141 (MQ=255)
                                                                                                     tGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTg                           <  1:62295/141‑1 (MQ=255)
                                                                                                        aGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGgaa                        <  1:20310/141‑1 (MQ=255)
                                                                                                        aGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGgaa                        <  1:51154/141‑1 (MQ=255)
                                                                                                         gctgTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGgaag                       <  1:16030/138‑1 (MQ=25)
                                                                                                          ccGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGgaaga                      >  1:72924/1‑141 (MQ=255)
                                                                                                              tGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGaaaaa                  >  1:58130/1‑141 (MQ=255)
                                                                                                               ggCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAg                 >  2:32264/1‑141 (MQ=255)
                                                                                                                 ctactaTCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCa               <  1:77493/141‑1 (MQ=35)
                                                                                                                 ctactaTCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCa               <  1:77633/141‑1 (MQ=255)
                                                                                                                 ctactaTCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCa               >  2:82407/1‑141 (MQ=255)
                                                                                                                  tactaTCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCAt              >  1:12651/1‑141 (MQ=35)
                                                                                                                  tactaTCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCAt              >  1:64007/1‑141 (MQ=255)
                                                                                                                   actaTCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATa             >  1:38932/1‑141 (MQ=255)
                                                                                                                    ctaTCTAAGATATGCATTGAGTTATCGTGATATATCTGAAAAATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGaaaaaaaaaaaaaa            <  2:24750/141‑7 (MQ=11)
                                                                                                                     taTCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATaaa           <  1:71865/141‑1 (MQ=35)
                                                                                                                      aTCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATaaaa          <  1:18763/141‑1 (MQ=255)
                                                                                                                         tAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATaaaaaag       <  2:12485/141‑3 (MQ=18)
                                                                                                                           aGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATaaaaaagct     >  1:79029/1‑137 (MQ=17)
                                                                                                                           aGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATaaaaaagct     >  2:35599/1‑137 (MQ=11)
                                                                                                                              tatGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATAAAAAGCTAAtt  <  2:42017/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                 
AGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATAAAAAGCTAATT  >  CP000730/76158‑76424

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