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

TTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGT  >  NC_000913/752989‑753264
                                                                                                                                        |                                                                                                                                           
ttATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaa                                                                                                                                          >  2:201439/1‑140 (MQ=255)
        gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                                  >  1:1244667/1‑140 (MQ=255)
        gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                                  >  1:1419034/1‑140 (MQ=255)
        gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                                  >  2:500133/1‑140 (MQ=255)
        gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                                  >  1:600995/1‑140 (MQ=255)
         cGACTACAATGAAGTTTTTAATTAAAAGTTTTGGGTTGCTGATATTTTGAGCTGTTTTTTTTTTTTTAAATATCTATATAGGTTTGTTAATGGATTTTTTTTTTACATTTTTTTGTGTTTAGGAATATAAAAATCAAccc                                                                                                                                 <  2:920391/140‑1 (MQ=255)
          gACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                                                                >  2:1129418/1‑140 (MQ=255)
                       tttttAATTGAAAGTTTTGGTTTGCTGATAATTTGAGCTTTTTTTTTTTTTTTAAATATTTATATAGGTCTGTTAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAAcggc                                                                                                                   <  2:32988/140‑1 (MQ=255)
                                      tttGGTTTGCTAATATTTTGAGTTTTTTTTTTCTTTTTAAATATTTATATAGGTTTGTTAATGGTTTTTTTTTTTACATTTTTTTTTTTTTAGGCATAAAAAAAACACCCCGCCATATGAACGGCGGGTTAAAATATTTa                                                                                                    <  2:131814/139‑1 (MQ=255)
                                      tttGGGTTGTTGATATTTGGTTTTTTTTTTTTTTTTTTAAATTTTTATATGGGTTTTTTAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGTATATAAAAAAAAACCCGCCATATGAACGGCGGGTTAAAATATTTa                                                                                                    <  2:564623/139‑1 (MQ=255)
                                             tGCTGATAATTTGAGTTTTTTTTTTTTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACATTTTTTTTTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTa                                                                                             <  2:421136/140‑1 (MQ=255)
                                             tGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACATTTTTTTGTTTTTAGGCATATAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTa                                                                                             <  1:1079176/140‑1 (MQ=255)
                                               cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                                           >  1:59383/1‑140 (MQ=255)
                                               cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                                           >  2:729015/1‑140 (MQ=255)
                                                 gATAATTTGAGCTGTTCTTTTTTTTTTAAATATCTATATAGGTCTGTTAATGGTTTTTTTTTTTAAAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCaa                                                                                         <  1:1129418/140‑1 (MQ=255)
                                                              tttCTTTTTTTTTTAAATATCTATATAGTTTTGTTAATGTTTTTTTTTTTTACATTTTTTTTTTTTTGGGAAAAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACg                                                                            <  2:704959/139‑1 (MQ=255)
                                                                         tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAAAGCTTGATATCGc                                                                 >  2:71492/1‑140 (MQ=255)
                                                                                           tCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGttttt                                               >  2:407303/1‑140 (MQ=255)
                                                                                           tCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGgtttt                                               >  1:768028/1‑140 (MQ=255)
                                                                                            cTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCAGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcgcg                                                    >  1:361595/1‑134 (MQ=255)
                                                                                            cTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCAGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcgcg                                                    <  2:361595/134‑1 (MQ=255)
                                                                                              tttAATGTTTTTTATTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  2:1382960/139‑1 (MQ=255)
                                                                                              tttAATGTTTTTTATTTTTACATTTTTTTGGTTTTAGACATAAAAAACTCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  2:1419034/139‑1 (MQ=255)
                                                                                              tttAATGGTTTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  1:274515/139‑1 (MQ=255)
                                                                                              gTTAGTGGTTTTTATTTTTACAATTTTTTGTTTTTGGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  1:854329/140‑1 (MQ=255)
                                                                                               ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                   >  2:847389/1‑132 (MQ=255)
                                                                                               ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                   <  1:847389/132‑1 (MQ=255)
                                                                                               ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCatt                                           >  1:1382960/1‑139 (MQ=255)
                                                                                               ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACCTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  1:55355/1‑140 (MQ=255)
                                                                                                tAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATCAAACATTAACGCTTGATATCGCTTTTAAAATCGCGTTTTTCatat                                          >  2:1007760/1‑140 (MQ=255)
                                                                                                 aatgttttttttttttAATTTTTTTTTTGTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                         <  2:1072640/134‑1 (MQ=255)
                                                                                                 aatgatttttttttttAAATTTTTTTGTTTTTAGGCAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                         <  2:413718/134‑1 (MQ=255)
                                                                                                 aaTGGTTTTTTTTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                         <  1:639684/140‑1 (MQ=255)
                                                                                                 aaTGGTTTTTTTTTTAACATTTTTTTGTTTTTAGACAAAAAAAATCCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                         <  2:1416538/140‑1 (MQ=255)
                                                                                                   tGGATTTTTTTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCt                                       <  1:142127/140‑1 (MQ=255)
                                                                                                      attttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGta                                    >  1:808444/1‑140 (MQ=255)
                                                                                                                    tttttttttttttAGGAAAAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCgg                      <  2:790526/139‑1 (MQ=255)
                                                                                                                          tgttTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                <  1:973316/140‑1 (MQ=255)
                                                                                                                                tAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt          <  2:718491/140‑1 (MQ=255)
                                                                                                                                 aGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTc         >  1:368579/1‑140 (MQ=255)
                                                                                                                                  ggCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCa        <  2:808444/140‑1 (MQ=255)
                                                                                                                                  ggCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCa        <  1:729015/140‑1 (MQ=255)
                                                                                                                                   gCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCAt       >  1:924002/1‑140 (MQ=255)
                                                                                                                                      tataAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAAAGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATAcc    >  2:1329449/1‑140 (MQ=255)
                                                                                                                                      tataAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATAcc    >  1:781619/1‑140 (MQ=255)
                                                                                                                                      tataAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATAcc    >  1:1110119/1‑140 (MQ=255)
                                                                                                                                        taAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGt  >  1:1290044/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                           
TTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGT  >  NC_000913/752989‑753264

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