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.

AGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAG  >  NC_000913/752956‑753339
                                                                                                                                                                           |                                                                                                                                                                                                                           
aGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGGTAGAAATATTTTC‑ACATAGCAATCAAACATCAAAGCATg                                                                                                                                                     >  2:19881/1‑239 (MQ=255)
                 cGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTTTAAATTATAAATCAACCCTTAAAGCTTGCAATCGCCTTTTaaatc                                                                                                                                    >  2:190018/1‑237 (MQ=255)
                   tCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGGTA‑AAATATTTTAAACTTAGTAATAAACCATTAACGCTTGATATAG‑CTTTTTAAATcca                                                                                                                                 >  2:315730/1‑237 (MQ=255)
                                                        ttttaAATTGAAAGTTTTGGGTTGCTGATAATTTGAGCTTTTCTATTCTTTTTAAATTTCTATATAGGTCTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGc                                                                                           <  1:315730/239‑1 (MQ=255)
                                                          tttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATTAACCATTAACGCTTGATATCG‑CTTTTTAAGTCGCGTTTTTCTTATCCCTGT‑TAACGCTGAAGCGGACGGGCaa                                                                                         >  2:708276/1‑239 (MQ=255)
                                                              aTTGAAATTATTGGTTTGCTGATAATTGGAGTTTTTTTATTTTTTTAAATTTTCTAAATAGGTTTGTTAATGGATTTTTTTTTTAAAT‑TTTTTTGTGTTT‑AGGCATAAAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCtt                                                                                     <  1:190018/239‑1 (MQ=255)
                                                                                                                   ttatataGTTTTGTAAATGTATTTTA‑TTTTACAATTTTTTTGTTTTT‑AGGCAAAAAAAAATCA‑‑CCCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGt                                                                                                                               <  2:248120/144‑1 (MQ=255)
                                                                                                                    tatataGCTATGTAAATGATTTTTATTTTTAAAT‑TTTTTTGTTTTT‑AGGCAAAAAAAAATCA‑‑CCCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGc                                                                                           <  2:435047/179‑1 (MQ=255)
                                                                                                                      tatgGGTTTGTAAATGTATTTTTTTTTAAAAT‑TTTTTTTTGTTT‑AGGAAAAAAAAAACCA‑‑CCCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTATTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCa                                                                                   <  2:663371/185‑1 (MQ=255)
                                                                                                                      tataGGTTTGTTAATGAATTTTATTTTTAAAT‑TTTTTTGTGTTT‑AGGCATAAAAAAATCA‑‑CCCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCa                                              <  1:558902/222‑1 (MQ=255)
                                                                                                                      tataGGTCTTTTAATGGTTTTTTTTTTTACAT‑TTTTTTGTTTTT‑AGGAATAAAAAAACCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGc                                                                                           <  2:199198/177‑1 (MQ=255)
                                                                                                                           gttttttAATTGATTTTTATTTTAAAA‑TTTTTTTGGTTTAAGAAATAAAAAAATCA‑‑CCCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCa                                              <  2:319214/213‑1 (MQ=255)
                                                                                                                           gtttGTAAAGGGTTTTTATTTTAAAAT‑TTTTTTTTTTTT‑AGGAATAAAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGcc                                           <  1:618352/217‑1 (MQ=255)
                                                                                                                           gTCTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCa                                                                           >  1:416479/1‑188 (MQ=255)
                                                                                                                             cTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑c                                                                                                                    >  1:248120/1‑145 (MQ=255)
                                                                                                                              tGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATa                                                                                 >  1:435047/1‑179 (MQ=255)
                                                                                                                                ttAATGGTTTTTATTTTTACAT‑TTTTTTGTTTTT‑AGGAATAAAAAAACCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCat                                                                                                                        <  2:168811/139‑1 (MQ=255)
                                                                                                                                ttAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACt                                                                         >  1:663371/1‑185 (MQ=255)
                                                                                                                                ttAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGGTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATa                                                                                 >  1:199198/1‑177 (MQ=255)
                                                                                                                                ttAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAATCGCGGTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAAACTTCATACCGTCACTTGGCATTTTGCTCCAATGGGCGATCCACCCAACGGTACGTGCCATTTCGGAAAt                   >  2:306089/1‑239 (MQ=255)
                                                                                                                                     ggATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTT‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTTAAGTCGCGGTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCACTCTTCATAACGTCACTGTGCATTTTGCTCCAGTGGGCGGTCCAGCCAACCGTACGt                                 >  2:618352/1‑220 (MQ=255)
                                                                                                                                     ggATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGGGATCCAGCCAACggga                                    >  1:319214/1‑215 (MQ=255)
                                                                                                                                          ttattttTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGtata                                                                                                             >  1:168811/1‑139 (MQ=255)
                                                                                                                                          ttattttTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAg                                                                                                          >  2:193542/1‑142 (MQ=255)
                                                                                                                                          ttattttTACAA‑TTTTTTGTGTTT‑AGGCATAAAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAg                                                                                                          <  1:193542/142‑1 (MQ=255)
                                                                                                                                               tttAAAT‑TTTTTTGTTTTT‑AGGCATAAAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGa    <  1:306089/239‑1 (MQ=255)
                                                                                                                                                  aCAA‑TTTTTTGTGTTA‑AAG‑ATATAAAAAATATCACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCCCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAg  <  2:558843/239‑1 (MQ=255)
                                                                                                                                                                           |                                                                                                                                                                                                                           
AGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAG  >  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: ‑