Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1050AT22.0% 125.1 / 10.9 50intergenic (‑310/+80)ybgD/gltAputative fimbrial protein YbgD/citrate synthase
Reads supporting (aligned to +/- strand):  ref base A (28/11);  new base T (0/11);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 1.89e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTG  >  NC_000913/752972‑753228
                                                                                                                                     |                                                                                                                           
aCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttttt                                                                                                                       >  1:722039/1‑140 (MQ=255)
aCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttt                                                                                                                         <  2:793707/138‑1 (MQ=255)
aCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttt                                                                                                                         >  1:793707/1‑138 (MQ=255)
 cGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTtg                                                                                                                      >  2:241508/1‑140 (MQ=255)
 cGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTtg                                                                                                                      >  1:879241/1‑140 (MQ=255)
 cGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTtg                                                                                                                      >  2:463137/1‑140 (MQ=255)
     aTTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGtt                                                                                                                  >  1:873134/1‑140 (MQ=255)
         tCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAgg                                                                                                              <  1:865205/140‑1 (MQ=255)
         tCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAgg                                                                                                              <  1:463137/140‑1 (MQ=255)
         gagacaggTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAgg                                                                                                              <  2:260145/134‑1 (MQ=255)
              gggTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCatat                                                                                                         <  2:740610/140‑1 (MQ=255)
                   ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaaaa                                                                                                    >  1:260013/1‑140 (MQ=255)
                         gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                              >  1:523972/1‑140 (MQ=255)
                           gACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                            >  1:263768/1‑140 (MQ=255)
                                       ttttttAATTAAAATTATTGGTTTGTTGATATTTTGAGCTTTTCTTTTTTTTTTAAATATCTATATAGGTTTGTTAATGGATTTTATTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAAcgg                                                                                <  2:873134/139‑1 (MQ=255)
                                        tttttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAAcggc                                                                               <  1:241508/140‑1 (MQ=255)
                                        tttttAATTAAAAGTTTTGGTTTGCTAATAATTTGAGCTGTTCTATTTTTTTAAAATATTTATATGGTTCTTTTAATGGATTTTATTTTAACAATTTTTTGTTTTTAGGAAAAAAAAAATCACCCCGCCATATGAAcggc                                                                               <  2:774967/140‑1 (MQ=255)
                                          tttATTTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAGATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAAcggcgg                                                                             >  2:225286/1‑140 (MQ=255)
                                                taaaaGTTTTGGGTTGCTGATATTTTGAGTTTTTCTTTTTTTTTTAAATATTTATATAGGTTTGTTAATGGTTTTTATTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTaaa                                                                       <  2:514020/138‑1 (MQ=255)
                                                         tGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACa                                                              >  1:843884/1‑140 (MQ=255)
                                                                cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                       >  1:583174/1‑140 (MQ=255)
                                                                cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                       >  2:779094/1‑140 (MQ=255)
                                                                cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAAAATTTACAAATTAGc                                                       >  1:794631/1‑140 (MQ=255)
                                                                          gagtttttCTTTTTTTTTTAAATATCTATATAGGTTTTTAAATGGTTTTTATTTTTAAAATTTTTTGTTTTTAGATAAAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                             <  2:794631/134‑1 (MQ=255)
                                                                          gAGCTTTTCTATTTTTTTTAAATATCTATATGGGTTTGTTAATGGATTTTATTTTTACATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                             <  1:273145/140‑1 (MQ=255)
                                                                                        tttttAAATATCTATATAGGTTTGTTAATGGTTTTTTTTTTTAAATTTTTTTGTGTTTAGGCATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATc                               <  2:260013/140‑1 (MQ=255)
                                                                                          tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                             >  1:799174/1‑140 (MQ=255)
                                                                                          tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                             >  1:731717/1‑140 (MQ=255)
                                                                                                    tatataGGTCTTTTAATGGATTTTTTTTTTACAATTTTTTGTGTTTAGGAAAAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTc                   <  2:106937/140‑1 (MQ=255)
                                                                                                          gttCTTTTAATGGTTTTTATTTTTACAATTTTTTGTGTTTAGGAATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGttt             <  1:495709/138‑1 (MQ=255)
                                                                                                           gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt            >  2:98515/1‑140 (MQ=255)
                                                                                                           gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt            >  1:934771/1‑140 (MQ=255)
                                                                                                           gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGgttt            >  1:846735/1‑140 (MQ=255)
                                                                                                            tCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGttttt           >  1:144303/1‑140 (MQ=255)
                                                                                                             cTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc          >  2:803689/1‑140 (MQ=255)
                                                                                                              tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa         >  1:174497/1‑140 (MQ=255)
                                                                                                              tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa         >  1:231881/1‑140 (MQ=255)
                                                                                                               tttATTGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGACATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat        <  2:430312/139‑1 (MQ=255)
                                                                                                               tttAATGTTTTTTTTTTTTAAATTTTTTTTTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat        <  2:231881/139‑1 (MQ=255)
                                                                                                               tttAATGTATTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat        <  2:902199/139‑1 (MQ=255)
                                                                                                               tttAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat        <  2:846735/139‑1 (MQ=255)
                                                                                                               tttAATGGTTTTTATTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat        <  2:343934/139‑1 (MQ=255)
                                                                                                               ttaaaTGGTTTTTATTTTTAAAATTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat        <  2:583174/137‑1 (MQ=255)
                                                                                                                ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCTGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata       >  2:754869/1‑140 (MQ=255)
                                                                                                                ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata       >  1:8363/1‑140 (MQ=255)
                                                                                                                ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata       >  1:85195/1‑140 (MQ=255)
                                                                                                                 tAATGTTTTTTTTTTTTACATTTTTTTTTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCatat      <  2:548384/140‑1 (MQ=255)
                                                                                                                  aaTGGTTTTTATTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc     <  1:407304/140‑1 (MQ=255)
                                                                                                                    tGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCt   >  1:690877/1‑140 (MQ=255)
                                                                                                                    agcttttttttctTAAAATTTTTTGTTTTCAGCTTTATAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCt   <  1:225286/136‑1 (MQ=255)
                                                                                                                     ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTg  >  1:343934/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                           
ACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTG  >  NC_000913/752972‑753228

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: ‑