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

GCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGT  >  NC_000913/752997‑753264
                                                                                                                                |                                                                                                                                            
gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                                   >  2:464756/1‑140 (MQ=255)
  gACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                                                                 >  2:68532/1‑140 (MQ=255)
  gACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                                                                 >  1:938621/1‑140 (MQ=255)
   aCTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGc                                                                                                                                >  1:32345/1‑140 (MQ=255)
   aCTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGc                                                                                                                                >  1:87402/1‑140 (MQ=255)
              ttttttAATTGAAAGTATTGGGTTGTTGATAATTTGAGCTTTTCTTTTTTTTTAAAATATCTATATAGGTTTGTTAATGGATTTTATTTTTACATTTTTTTGTTTTTAGGAATATAAAAATCAACCCGCCATATGAAcgg                                                                                                                     <  2:930182/139‑1 (MQ=255)
                              tttGGTTTGCTGATATTTTGAGCTTTTCTTTTTTTTTAAAATATTTATATAGGTCTATTAATGGATTTTATTTTTCGATTTTTTTGTGTTTAGGAATATAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTa                                                                                                     <  2:158694/139‑1 (MQ=255)
                                   tttGCTAATAATTTGATCTTTTCTTTTTTTTTAAAATATCTATATAGGTCTTTTAATGGATTTTATTTTTACAATTTTTTTTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACt                                                                                                <  2:525642/139‑1 (MQ=255)
                                   gTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACt                                                                                                <  1:429626/140‑1 (MQ=255)
                                       cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTATAACTTAga                                                                                            >  2:717786/1‑139 (MQ=255)
                                       cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                                            >  2:188969/1‑140 (MQ=255)
                                           tattttGAGTTTTTTTTTTTTTTTTAAATTTTTATATAGGTCTGTTAATGGATTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATc                                                                                        <  2:861642/137‑1 (MQ=255)
                                             aTTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCaa                                                                                      >  1:782883/1‑140 (MQ=255)
                                                    cTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTaa                                                                               >  1:418558/1‑140 (MQ=255)
                                                              tttttaAAATTTTTATATAGGTCTGTTAATGGATTTTTTTTTTACAATTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                                                     <  1:355090/139‑1 (MQ=255)
                                                               tttttAAATTTCTATATAGGTCTGTTAATGGTTTTTTTTTTTACAATTTTTTGTTTTTAGGAATATAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATc                                                                    <  1:343330/140‑1 (MQ=255)
                                                                 tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATcac                                                                  >  2:153006/1‑138 (MQ=255)
                                                                 tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                  >  1:754665/1‑140 (MQ=255)
                                                                 tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                  >  1:268238/1‑140 (MQ=255)
                                                                      atatCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTa                                                             <  1:4481/140‑1 (MQ=255)
                                                                             tataGGTTTTTTAATGGTTTTTATTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTcgc                                                     <  2:87402/140‑1 (MQ=255)
                                                                               taGGTTTTTTAATGGTTTTTATTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGt                                                   <  2:418558/140‑1 (MQ=255)
                                                                               taGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGt                                                   <  1:493838/140‑1 (MQ=255)
                                                                               taGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGt                                                   >  1:224105/1‑140 (MQ=255)
                                                                                  gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGtttt                                                >  2:428357/1‑140 (MQ=255)
                                                                                  gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGtttt                                                >  1:233479/1‑140 (MQ=255)
                                                                                    tttttAATGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTc                                              <  1:68532/137‑1 (MQ=255)
                                                                                    ttttaaaTGTATTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTc                                              <  2:628301/135‑1 (MQ=255)
                                                                                     tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTCa                                             >  2:802915/1‑140 (MQ=255)
                                                                                     tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAAAGCTTGATATCGCTTTTTAA‑GTCGCGGTTTTCa                                             >  2:645584/1‑140 (MQ=255)
                                                                                      tttAATGTTTTTTATTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTCat                                            <  2:503699/139‑1 (MQ=255)
                                                                                      tttAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTCat                                            <  2:268238/139‑1 (MQ=255)
                                                                                       ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCTGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTCata                                           >  1:808559/1‑140 (MQ=255)
                                                                                       ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTATAACTTAACAATCAAACATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTtata                                           >  2:204878/1‑140 (MQ=255)
                                                                                       ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTCata                                           >  1:256555/1‑140 (MQ=255)
                                                                                       ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTCata                                           >  2:870592/1‑140 (MQ=255)
                                                                                         aatgttttttttttttAAATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTCATATc                                         <  2:32345/134‑1 (MQ=255)
                                                                                         aaTGTATTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTCATATc                                         <  2:782883/140‑1 (MQ=255)
                                                                                             gATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACCCTTGATATCGCTTTTAAA‑GTCGCGGTTTTCATATCCTGt                                     >  1:503699/1‑140 (MQ=255)
                                                                                                 ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTATCACTTAGCCAACAACCAATAACGCTTGAAATTGCTTTTAAAGGTTCCGTTTTTC‑TATCCttttcac                                 >  2:170658/1‑134 (MQ=255)
                                                                                                      tttACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTCATATCCTGTATACAGCTg                            <  1:513990/140‑1 (MQ=255)
                                                                                                             ttttttGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGa                     >  1:920994/1‑140 (MQ=255)
                                                                                                                  tgttTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                <  2:754665/140‑1 (MQ=255)
                                                                                                                  tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                >  2:333878/1‑140 (MQ=255)
                                                                                                                              tataAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATAcc    >  1:855936/1‑140 (MQ=255)
                                                                                                                                taAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGt  >  1:443321/1‑140 (MQ=255)
                                                                                                                                |                                                                                                                                            
GCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAA‑GTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGT  >  NC_000913/752997‑753264

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