Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1250TA25.5% 94.7 / 19.1 47intergenic (‑330/+60)ybgD/gltAputative fimbrial protein YbgD/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (27/8);  new base A (0/12);  total (27/20)
Fisher's exact test for biased strand distribution p-value = 2.41e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCA  >  NC_000913/752987‑753258
                                                                                                                                          |                                                                                                                                     
ggTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCAtata                                                                                                                                      >  1:1229294/1‑140 (MQ=255)
      aaaTGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATc                                                                                                                                >  1:46719/1‑140 (MQ=255)
        gtGCGACTACCATGATTTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTTTTTTTTTTTTTTTAAATATCTATATAGGTTTTTTAATGGATTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCaa                                                                                                                              <  2:221406/139‑1 (MQ=255)
          gCGATTACCATGAATTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGTTTTTTTATTTTTTTTAAATATTTATATAGGTTTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGAATATAAAAATCAAcc                                                                                                                            <  2:432154/140‑1 (MQ=255)
           cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAccc                                                                                                                           <  1:1223911/140‑1 (MQ=255)
           cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAA‑TTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                                                          >  1:409021/1‑140 (MQ=255)
           cGACTACCATGAAGTTTTTAATTGAAAGCATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAccc                                                                                                                           >  2:257695/1‑140 (MQ=255)
             aCTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGc                                                                                                                         >  1:8856/1‑140 (MQ=255)
                    tGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATgc                                                                                                                  >  1:221406/1‑139 (MQ=255)
                    tGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGa                                                                                                                  >  2:229592/1‑140 (MQ=255)
                         tttttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAAcggc                                                                                                             <  1:75143/140‑1 (MQ=255)
                           tttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAAcggcgg                                                                                                           >  2:74432/1‑140 (MQ=255)
                                 taaaaGTATTGGGTTGCTGATAATTTGAGCTTTTTTTTTCTTTTTAAATATTTATATAGGTTTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTaaa                                                                                                     <  2:1029151/138‑1 (MQ=255)
                                                 cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                                     >  1:749257/1‑140 (MQ=255)
                                                 cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                                     >  1:722669/1‑140 (MQ=255)
                                                     tAATTTGAGTTTTTCTATTTTTTTTAAATATCTATATAGGTCTGTTAATGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATc                                                                                 <  2:1415794/140‑1 (MQ=255)
                                                           gAGCTTTTCTATTTTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTAAAATTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                           <  1:1610784/140‑1 (MQ=255)
                                                             gCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTa                                                                         <  1:1355825/140‑1 (MQ=255)
                                                                        tttttaAAATATTTATATAGGTTTTTTAATGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                                              <  2:722669/139‑1 (MQ=255)
                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATCAACCATTAAAACTTGATATCGc                                                           >  2:1068201/1‑140 (MQ=255)
                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                           >  2:1045150/1‑140 (MQ=255)
                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                           >  2:1340673/1‑140 (MQ=255)
                                                                                         taGGTTTGTTAATGTATTTTTTTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                             <  1:797392/140‑1 (MQ=255)
                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                          >  1:1564997/1‑140 (MQ=255)
                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                          >  1:1477117/1‑140 (MQ=255)
                                                                                             tCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGttttt                                         >  1:1278802/1‑140 (MQ=255)
                                                                                             tCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGttttt                                         >  1:1398878/1‑140 (MQ=255)
                                                                                              tttttAATGTATTTTTTTTTTAAAATTTTTTGTTTTTAGAAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                        <  2:1564997/137‑1 (MQ=255)
                                                                                               tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                       >  1:613348/1‑140 (MQ=255)
                                                                                               tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                       >  1:1368826/1‑140 (MQ=255)
                                                                                               tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                       >  2:91669/1‑140 (MQ=255)
                                                                                                tttAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGAAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                      <  2:960921/139‑1 (MQ=255)
                                                                                                tttAATGGTTTTTATTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                      <  2:747742/139‑1 (MQ=255)
                                                                                                tttAATGGATTTTTTTTTTACAATTTTTTGTGTTTGGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                      <  1:463932/139‑1 (MQ=255)
                                                                                                tttAATGGATTTTTTTTTTAAATTTTTTTGTTTTTAGCAAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                      <  2:1368826/139‑1 (MQ=255)
                                                                                                tttAAGGGATTTTTTTTTTACATTTTTTTGTTTTTAGTCATAAAAAAACCATCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                      <  2:1229294/139‑1 (MQ=255)
                                                                                                gTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                      <  1:1045150/140‑1 (MQ=255)
                                                                                                gTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                      <  1:1340673/140‑1 (MQ=255)
                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                     >  1:1571159/1‑140 (MQ=255)
                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                     >  1:960921/1‑140 (MQ=255)
                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                     >  1:231749/1‑140 (MQ=255)
                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                     >  1:747742/1‑140 (MQ=255)
                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTACAACTTATCAATCAACCATTAAAGCTTGATATTGCTTTTTAAGTCGCGTTTTTCata                                     >  2:769744/1‑140 (MQ=255)
                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTACAACTTAGCAATCAACCATTAAAGCTTGATATTGCTTTTTAAGTCGCGTTTTTCata                                     >  2:1085574/1‑140 (MQ=255)
                                                                                                   aagattttttttttttACAATTTTTTGTGTTTAGGCATATAAAACACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                   <  1:710138/134‑1 (MQ=255)
                                                                                                    aTGGTTTTTTTTTTTAAAATTTTTTGTTTTTAGACATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATcc                                  <  2:688701/140‑1 (MQ=255)
                                                                                                           ttattttTACATTTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAATGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                             <  1:1026940/122‑1 (MQ=255)
                                                                                                           ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAATGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                             >  2:1026940/1‑122 (MQ=255)
                                                                                                           ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAc                           >  2:591725/1‑140 (MQ=255)
                                                                                                                   tacaTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcg                   <  2:409021/137‑1 (MQ=255)
                                                                                                                           acatgtTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACggg           <  1:1500096/137‑1 (MQ=255)
                                                                                                                            tgttTTTAGGAATATAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc          <  1:292939/140‑1 (MQ=255)
                                                                                                                                  tAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt    <  1:74432/140‑1 (MQ=255)
                                                                                                                                  tAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt    <  2:1278802/140‑1 (MQ=255)
                                                                                                                                    ggCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCa  <  1:229592/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                     
GGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCA  >  NC_000913/752987‑753258

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑