Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1250TA42.5% 39.7 / 27.8 40intergenic (‑330/+60)ybgD/gltAputative fimbrial protein YbgD/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (21/2);  new base A (0/17);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.93e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCA  >  NC_000913/752987‑753251
                                                                                                                                          |                                                                                                                               
ggTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCAtata                                                                                                                                >  1:347187/1‑140 (MQ=255)
    ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaaaa                                                                                                                            >  2:158283/1‑140 (MQ=255)
          gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                      >  2:724865/1‑140 (MQ=255)
          gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                      >  1:91215/1‑140 (MQ=255)
            gACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                                                    >  1:633482/1‑140 (MQ=255)
             attACCATGAATTTTTTAATTAAAAGTTTTGGTTTGCTGATAATTTGAGCTGTTTTTTTTTTTTTAAATATTTATATAGGTTTTTTAATGTATTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCACCCCGc                                                                                                                   <  2:633482/138‑1 (MQ=255)
                        ttttttATTTAAAATTTTTGGGTGGCTAATAATTTGAGCTTTTCTATTTTTTTAAAATATCTATATTGGTCTTTTAATGGATTTTATTTTTACAATTTTTTTTGTTTAGGCATATAAAAATCAACCCGCCATATGAAcgg                                                                                                        <  2:206359/139‑1 (MQ=255)
                                                 cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGc                                                                              >  1:15064/1‑140 (MQ=255)
                                                      aaTTTGAGCTTTTCTTTTTTTTTAAAATATCTATATAGGTCTGTTAATGGTTTTTTTTTTTAAAATTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCa                                                                         <  2:543709/140‑1 (MQ=255)
                                                           gAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAAccct                                                                    >  2:735107/1‑138 (MQ=255)
                                                                 ttCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGc                                                              >  2:700706/1‑140 (MQ=255)
                                                                      ttCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGat                                                         >  1:286958/1‑140 (MQ=255)
                                                                        ttttttAAATATTTATATAGGTCTTTTAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                                       <  2:15064/139‑1 (MQ=255)
                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTTGCAATCAACCATTAAAGCTTGATATTGc                                                    >  2:269153/1‑140 (MQ=255)
                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                    >  1:502187/1‑140 (MQ=255)
                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                    >  1:14410/1‑140 (MQ=255)
                                                                                atatCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTTGCAATTAACCATTAACGCTTGGTATCGCttttt                                               >  2:841580/1‑139 (MQ=255)
                                                                                   tCTATATAGGTCTTTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGAAATAAAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAg                                            <  1:626721/140‑1 (MQ=255)
                                                                                         taGGTTTGTTAATGTTTTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                      <  2:816965/140‑1 (MQ=255)
                                                                                         taGGTTTGTTAATGGTTTTTATTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                      <  2:91215/140‑1 (MQ=255)
                                                                                           ggTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGttt                                    >  2:734136/1‑140 (MQ=255)
                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                   >  1:62978/1‑140 (MQ=255)
                                                                                              tttttAATGTTTTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                 <  2:206320/137‑1 (MQ=255)
                                                                                              ttGTTAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                 <  2:62978/139‑1 (MQ=255)
                                                                                              ttGTTAATGGTTTTTATTTTTACATTTTTTTTTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                 <  2:167775/139‑1 (MQ=255)
                                                                                               tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                >  1:816965/1‑140 (MQ=255)
                                                                                               tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                >  1:80228/1‑140 (MQ=255)
                                                                                               tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAAACATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                >  1:167775/1‑140 (MQ=255)
                                                                                               tGTTAATGGATTTTA‑TTTTACAATTTTTTGGGTTTAAGGATATAAAAATCAACCCGGCATATTAACGGCGGGGTTAAAAAATTTCAAACTTACAAATAAACAATTAACCCTT‑ATATTGCCTTTTAAATCCCCttttttct                               >  2:7663/1‑137 (MQ=255)
                                                                                                gTTAATGATTTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                               <  2:347187/140‑1 (MQ=255)
                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGTTTTTCata                              >  2:556484/1‑140 (MQ=255)
                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTtata                              >  1:642950/1‑140 (MQ=255)
                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                              >  1:554937/1‑140 (MQ=255)
                                                                                                   aaTGTATTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                            <  1:713753/140‑1 (MQ=255)
                                                                                                   aaTGGTTTTTATTTTTACATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                            <  1:232321/140‑1 (MQ=255)
                                                                                                   aaTGGATTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                            <  1:466869/140‑1 (MQ=255)
                                                                                                   aaTGGATTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                            <  1:817803/140‑1 (MQ=255)
                                                                                                           ttttttttAAATTTTTTTTTTTTTAGGAAAAAAAAAACCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAc                    <  2:218722/137‑1 (MQ=255)
                                                                                                                   aaatttttttGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcg            <  2:554937/136‑1 (MQ=255)
                                                                                                                     atttttttGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcgcg          <  1:158283/138‑1 (MQ=255)
                                                                                                                       ttttttGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGa        <  2:286958/140‑1 (MQ=255)
                                                                                                                        tttttGTGTTTAGAAATAAAAAAATCACCCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGAc       <  2:642950/140‑1 (MQ=255)
                                                                                                                            tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc   <  1:700706/140‑1 (MQ=255)
                                                                                                                             gtgtTTAGGCATAAAAAAACCACCCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCa  <  1:269153/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                               
GGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGC‑GGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCA  >  NC_000913/752987‑753251

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