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

ATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACC  >  NC_000913/752991‑753262
                                                                                                                                      |                                                                                                                                         
ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaaaa                                                                                                                                      >  1:315558/1‑140 (MQ=255)
        gACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                                                              >  1:387927/1‑140 (MQ=255)
                    ttttttAATTGAAAGTTTTGGGTTGCTAATATTTTGAGTTTTTTTATTTTTTTTAAATTTTTATATAGGTTTGTTAATGTTTTTTTTTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAAcgg                                                                                                                  <  2:315558/139‑1 (MQ=255)
                                  ttATTGGGTTGCTAATAATTTGAGTTTTTTTTTTTTTTTTAAATATTTATATAGGTTTTTTAATGGTTTTTATTTTTACAATTTTTTGTTTTTAGGAATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATAtt                                                                                                    <  2:837263/139‑1 (MQ=255)
                                      tGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACa                                                                                                >  2:478351/1‑140 (MQ=255)
                                          ttGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACtt                                                                                            <  1:115651/140‑1 (MQ=255)
                                                        aGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAtt                                                                              >  1:107243/1‑140 (MQ=255)
                                                             ttCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATCAACCATTAACGc                                                                         >  2:709962/1‑140 (MQ=255)
                                                             ttCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGc                                                                         >  1:301165/1‑140 (MQ=255)
                                                                  ttCTTTTTAAATATCTATATAGGCCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGATTAAAATATTTACAACCTAGCAATCAACCATTAACGCTTgtt                                                                    >  1:231179/1‑138 (MQ=255)
                                                                       tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                               >  1:944767/1‑140 (MQ=255)
                                                                       tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                               >  2:233944/1‑140 (MQ=255)
                                                                       tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                               >  1:843767/1‑140 (MQ=255)
                                                                       tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                               >  1:748069/1‑140 (MQ=255)
                                                                                     taGTTCTGTTAATGTTTTTTATTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                 <  2:105867/140‑1 (MQ=255)
                                                                                        gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGTTTTAAAATATTTAAAAAATAGCAATCAACCAATAAAGCTTGATATTTGATTTAAAGCAGCGtttt                                              >  2:7936/1‑140 (MQ=255)
                                                                                        gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAAAGCTTGATATCGCTTTTTAAGTCGCGtttt                                              >  2:920412/1‑140 (MQ=255)
                                                                                           tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                           >  1:588200/1‑140 (MQ=255)
                                                                                           tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                           >  1:105867/1‑140 (MQ=255)
                                                                                           tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                           >  1:360285/1‑140 (MQ=255)
                                                                                            tttATTGGTTTTTTTTTTTAAATTTTTTTGTGTTTAGAAAAAAAAAAACAACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                          <  2:843767/139‑1 (MQ=255)
                                                                                            tttAATGGTTTTTATTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                          <  1:219067/139‑1 (MQ=255)
                                                                                            gTTAATGGATTTTATTTTTATATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                          <  1:404034/140‑1 (MQ=255)
                                                                                             ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                         >  2:442634/1‑140 (MQ=255)
                                                                                             ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAATCGCGTTTTTCatt                                         >  2:608034/1‑139 (MQ=255)
                                                                                               aatgttttttttttttAAATTTTTTTGTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                       <  2:301165/134‑1 (MQ=255)
                                                                                                  ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTGTAATATTTTCAACATAACAATCAACCAATAACGGTTGAACTCCGGTTTAAAGGCGTGTTTTTTATATCCTg                                    >  2:219067/1‑140 (MQ=255)
                                                                                                       ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAc                               >  1:945860/1‑140 (MQ=255)
                                                                                                                        tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTTCAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGTTTTTTATATCCTGTATTCAGCTGACGCGGACgggg              >  2:256778/1‑139 (MQ=255)
                                                                                                                        tgtgTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc              <  1:233944/140‑1 (MQ=255)
                                                                                                                                    tataAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATAcc  >  1:152690/1‑140 (MQ=255)
                                                                                                                                    tataAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATAcc  >  1:830848/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                         
ATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACC  >  NC_000913/752991‑753262

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