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

GCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGT  >  NC_000913/752997‑753264
                                                                                                                                |                                                                                                                                           
gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                                  >  1:1097312/1‑140 (MQ=255)
ccgccTACCATGAAGTTTTTAATTGAAAGTTGTGGGTGGCGAATATTTTGATGTTTTTTTTTTTTTTTAAATATTTTTATAGGTCTTTTAATGAATTTTATTTTTCAAATTTTTTGTGTTTAGGCATAAAAAAATCAAcc                                                                                                                                  <  2:733308/136‑1 (MQ=255)
 cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAccc                                                                                                                                 >  2:454263/1‑140 (MQ=255)
            aaGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAAc                                                                                                                      >  1:149429/1‑140 (MQ=255)
              ttttttAATTAAAAGTATTGGGTTGCTAATATTTTGAGCTTTTCTTTTCTTTTTAAATATTTATATAGGTCTTTTAAGGGATTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAAcgg                                                                                                                    <  2:277922/139‑1 (MQ=255)
                       tGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTTTTTTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACATTTTTTTGTGTTTAGGCATATAAAAATCACCCCGCCATATGAACGGCGGGTTaaa                                                                                                           <  2:1283806/140‑1 (MQ=255)
                            gttttGGGTTGCTGATATTTTGAGTTTTTTTATTTTTTTTAAATATCTATATAGGTTTGTTAATGTTTTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATAtt                                                                                                      <  2:483211/137‑1 (MQ=255)
                                   tttGCTTATATTTTGAGTTTTTTTTTTTTTTTTAAATATCTATATAGGTTTTTTAATGGATTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACt                                                                                               <  2:1344704/139‑1 (MQ=255)
                                       cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGc                                                                                           >  2:552010/1‑140 (MQ=255)
                                       cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                                           >  1:382949/1‑140 (MQ=255)
                                         gATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCaa                                                                                         >  1:1344704/1‑140 (MQ=255)
                                           tAATTTGAGCTTTTCTTTTTTTTTTAAATATCTATATAGGTCTGTTAATGGTTTTTTTTTTTACATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATc                                                                                       <  1:354652/140‑1 (MQ=255)
                                            aaTTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCa                                                                                      >  1:910497/1‑140 (MQ=255)
                                                 gatcttttCTATTCTTTTAAAATATCTATATAGGTCTTTTAAGGTTTTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                                 <  1:339676/134‑1 (MQ=255)
                                                        ttttttttttttAAATATCTATATAGGTCTGTTAATGTTTTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCt                                                                          <  1:454263/133‑1 (MQ=255)
                                                               tttttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATc                                                                   <  1:219741/140‑1 (MQ=255)
                                                                 tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  1:185889/1‑140 (MQ=255)
                                                                 tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  2:530256/1‑140 (MQ=255)
                                                                 tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  1:822070/1‑140 (MQ=255)
                                                                 tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  1:456395/1‑140 (MQ=255)
                                                                             tataGGTTTTTTAATGTTTTTTTTTTTTACATTTTTTTGTTTTTAGAAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcgc                                                     <  2:822070/140‑1 (MQ=255)
                                                                               taGTTCTGTTACTGGTTTTTATTTTTACAATTTTTTGTTTTTAGGCATATAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                   <  1:230614/140‑1 (MQ=255)
                                                                               taGGTCTGTTAATGTATTTTATTTTTAAATTTTTTTGTTTTAAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                   <  2:901692/140‑1 (MQ=255)
                                                                                  gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAAACATTAAAGCTTGATATCCCTTTTAAAGTCGCGtttt                                                >  2:622190/1‑140 (MQ=255)
                                                                                     tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                             >  1:901692/1‑140 (MQ=255)
                                                                                     tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                             >  2:884201/1‑140 (MQ=255)
                                                                                     tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                             >  2:711789/1‑140 (MQ=255)
                                                                                      tttAATGTTTTTTTTTTTTACATTTTTTTTTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  2:1615637/139‑1 (MQ=255)
                                                                                      tttAATGGTTTTTTTTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  2:1097312/139‑1 (MQ=255)
                                                                                      tttAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  2:28212/139‑1 (MQ=255)
                                                                                      tttAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  2:973444/139‑1 (MQ=255)
                                                                                      gTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  1:884201/140‑1 (MQ=255)
                                                                                       ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  2:1368638/1‑140 (MQ=255)
                                                                                       ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  1:1496135/1‑140 (MQ=255)
                                                                                       ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  2:115620/1‑140 (MQ=255)
                                                                                       ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  2:821847/1‑140 (MQ=255)
                                                                                       ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  1:1615637/1‑140 (MQ=255)
                                                                                       ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAACAAGCATTAAAGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  2:1227801/1‑140 (MQ=255)
                                                                                         aaTGTATTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                         <  2:149429/140‑1 (MQ=255)
                                                                                         aaTGGTTTTTATTTTTATATTTTTTTGTTTTTAGACATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                         <  1:622190/140‑1 (MQ=255)
                                                                                            ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTg                                      >  1:607892/1‑140 (MQ=255)
                                                                                            ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTg                                      >  1:973444/1‑140 (MQ=255)
                                                                                             gATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAAAATTTTCAACTTAGCAATCAACCTTTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGt                                     >  2:517756/1‑140 (MQ=255)
                                                                                              tttttattttTAAATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGta                                    <  1:1227801/139‑1 (MQ=255)
                                                                                                 ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAc                                 >  1:108323/1‑140 (MQ=255)
                                                                                                       ttACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGa                           >  1:1467966/1‑140 (MQ=255)
                                                                                                        tACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAc                          >  1:1514275/1‑140 (MQ=255)
                                                                                                         aCAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcg                         >  2:935454/1‑140 (MQ=255)
                                                                                                             ttttttGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGa                     >  2:571144/1‑140 (MQ=255)
                                                                                                              tttttGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACAGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGAc                    <  1:955021/140‑1 (MQ=255)
                                                                                                                 ttGTTTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACggg                 <  2:662984/140‑1 (MQ=255)
                                                                                                                            aataaAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATa      <  2:1474008/139‑1 (MQ=255)
                                                                                                                                taAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGt  >  2:1090440/1‑140 (MQ=255)
                                                                                                                                taAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAACAACCATTAACGCTTGGTATCGCTTTTAAAGTCGCGTTTTTCATATTCCGTATGCAGCTAACGCGGGCGGGCCATCTTCATACCGt  >  2:1308834/1‑140 (MQ=255)
                                                                                                                                |                                                                                                                                           
GCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGT  >  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: ‑