Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1250TA40.7% 30.3 / 17.4 27intergenic (‑330/+60)ybgD/gltAputative fimbrial‑like adhesin protein/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (14/2);  new base A (0/11);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 5.98e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAG  >  NC_000913/752956‑753339
                                                                                                                                                                           |                                                                                                                                                                                                                           
aGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGGTAGAAATATTTTCACATAGCAATCAAACATCAAAGCATg                                                                                                                                                     >  2:19881/1‑239 (MQ=255)
                 cGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTTTAAATTATAAATCAACCCTTAAAGCTTGCAATCGCCTTTTaaatc                                                                                                                                    >  2:190018/1‑237 (MQ=255)
                   tCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGGTAAAATATTTTAAACTTAGTAATAAACCATTAACGCTTGATATAGCTTTTTAAATcca                                                                                                                                 >  2:315730/1‑237 (MQ=255)
                                                        ttttaAATTGAAAGTTTTGGGTTGCTGATAATTTGAGCTTTTCTATTCTTTTTAAATTTCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                                                                                           <  1:315730/239‑1 (MQ=255)
                                                          tttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATTAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGTTTTTCTTATCCCTGTTAACGCTGAAGCGGACGGGCaa                                                                                         >  2:708276/1‑239 (MQ=255)
                                                              aTTGAAATTATTGGTTTGCTGATAATTGGAGTTTTTTTATTTTTTTAAATTTTCTAAATAGGTTTGTTAATGGATTTTTTTTTTAAATTTTTTTGTGTTTAGGCATAAAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt                                                                                     <  1:190018/239‑1 (MQ=255)
                                                                                                                   ttatataGTTTTGTAAATGTATTTTATTTTACAATTTTTTTGTTTTTAGGCAAAAAAAAATCA‑‑CCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                                                                                               <  2:248120/144‑1 (MQ=255)
                                                                                                                    tatataGCTATGTAAATGATTTTTATTTTTAAATTTTTTTGTTTTTAGGCAAAAAAAAATCA‑‑CCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                                                                                           <  2:435047/179‑1 (MQ=255)
                                                                                                                      tatgGGTTTGTAAATGTATTTTTTTTTAAAATTTTTTTTTGTTTAGGAAAAAAAAAACCA‑‑CCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTATTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCa                                                                                   <  2:663371/185‑1 (MQ=255)
                                                                                                                      tataGGTTTGTTAATGAATTTTATTTTTAAATTTTTTTGTGTTTAGGCATAAAAAAATCA‑‑CCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCa                                              <  1:558902/222‑1 (MQ=255)
                                                                                                                      tataGGTCTTTTAATGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAACCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                                                                                           <  2:199198/177‑1 (MQ=255)
                                                                                                                           gttttttAATTGATTTTTATTTTAAAATTTTTTTGGTTTAAGAAATAAAAAAATCA‑‑CCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCa                                              <  2:319214/213‑1 (MQ=255)
                                                                                                                           gtttGTAAAGGGTTTTTATTTTAAAATTTTTTTTTTTTTAGGAATAAAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGcc                                           <  1:618352/217‑1 (MQ=255)
                                                                                                                           gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCa                                                                           >  1:416479/1‑188 (MQ=255)
                                                                                                                             cTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                                                                                                    >  1:248120/1‑145 (MQ=255)
                                                                                                                              tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATa                                                                                 >  1:435047/1‑179 (MQ=255)
                                                                                                                                ttAATGGTTTTTATTTTTACATTTTTTTGTTTTTAGGAATAAAAAAACCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                                                                                                        <  2:168811/139‑1 (MQ=255)
                                                                                                                                ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACt                                                                         >  1:663371/1‑185 (MQ=255)
                                                                                                                                ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGGTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATa                                                                                 >  1:199198/1‑177 (MQ=255)
                                                                                                                                ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAATCGCGGTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAAACTTCATACCGTCACTTGGCATTTTGCTCCAATGGGCGATCCACCCAACGGTACGTGCCATTTCGGAAAt                   >  2:306089/1‑239 (MQ=255)
                                                                                                                                     ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGGTTTTCATATCCTGTATACAGCTGACGCGGACGGGCACTCTTCATAACGTCACTGTGCATTTTGCTCCAGTGGGCGGTCCAGCCAACCGTACGt                                 >  2:618352/1‑220 (MQ=255)
                                                                                                                                     ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGGGATCCAGCCAACggga                                    >  1:319214/1‑215 (MQ=255)
                                                                                                                                          ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtata                                                                                                             >  1:168811/1‑139 (MQ=255)
                                                                                                                                          ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAg                                                                                                          >  2:193542/1‑142 (MQ=255)
                                                                                                                                          ttattttTACAATTTTTTGTGTTTAGGCATAAAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAg                                                                                                          <  1:193542/142‑1 (MQ=255)
                                                                                                                                               tttAAATTTTTTTGTTTTTAGGCATAAAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGa    <  1:306089/239‑1 (MQ=255)
                                                                                                                                                  aCAATTTTTTGTGTTAAAGATATAAAAAATATCACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCCCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAg  <  2:558843/239‑1 (MQ=255)
                                                                                                                                                                           |                                                                                                                                                                                                                           
AGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAG  >  NC_000913/752956‑753339

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: