Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913635,0920AT22.2% 120.3 / 11.6 45D95E (GAT→GAA) ybdMSpo0J family protein, ParB‑like nuclease domain
Reads supporting (aligned to +/- strand):  ref base A (6/29);  new base T (10/0);  total (16/29)
Fisher's exact test for biased strand distribution p-value = 2.51e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCCAACTGGCTAAGTTCACGGACGATTTCCGACATGGCGGTGATTTGATGGCGTCCACGAGCGCGATTATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACT‑GGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGA  >  NC_000913/634877‑635305
                                                                                                                                                                                                                       |                                                                                                                                                                                                                       
cccaACTGGCTAAGTTCACGGACGATTTCCGACATGGCGGTGATTTGATGGCGTCCACGAGCGCGATTATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAAtggg                                                                                                                                                                                                  <  1:249928/239‑1 (MQ=255)
cccaACTGGCTAAGTTCACGGACGATTTCCGACATGGCGGTGATTTGATGGCGTCCACGAGCGCGATTATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAAtggg                                                                                                                                                                                                  <  1:323165/239‑1 (MQ=255)
cccaACTGGCTAAGTTCACGGACGATTTCCGACATGGCGGTGATTTGATGGCGTCCACGAGCGCGATTATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAAtggg                                                                                                                                                                                                  <  1:771671/239‑1 (MQ=255)
                         tttCCGACATGGCGGTGATTTGATGGCGTCCACGAGCGCGATTATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGa                                                                                                                                                                                                            <  1:826717/204‑1 (MQ=255)
                                                                  ttATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGc                                                                                                                                        <  1:440291/230‑1 (MQ=255)
                                                                  ttATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACgttgtt                                                                                                                               <  1:433165/239‑1 (MQ=255)
                                                                   tATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGttt                                                                                                                              <  1:14457/239‑1 (MQ=255)
                                                                           aTAGTCGTGGCGATGCTTTGAGTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTa                                                                                                                      <  2:724089/239‑1 (MQ=255)
                                                                                              taatACGTGTAACTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTgg                                                                                                   <  2:703108/236‑1 (MQ=255)
                                                                                               aTTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCAt                                                                                                                 <  1:268901/224‑1 (MQ=255)
                                                                                                         ccTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAAc                                                                                        <  1:918012/239‑1 (MQ=255)
                                                                                                           ttCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGTCTTTACCAATTTCATGGCGGTGAAAAACATTACAGATTTACCAGGCATTTTTTTTCGTCTGTGGGACCACAATAGGCTGGGTTAAACACTTTTTTTTCAA‑GGATTTTTTGAGAAGGTTTTTTTCGTGGGGTGCCACGACGCTTGGGTTGTATACATTATGGCTAACTTGGGTGGTTTTaatt                                                                                       >  2:909318/1‑236 (MQ=255)
                                                                                                           ttCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTa                                                                                                                                                                                                                   <  1:855998/115‑1 (MQ=255)
                                                                                                           ttCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGtttt                                                                                                                                                      <  1:636258/175‑1 (MQ=255)
                                                                                                             ccAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGa                                                                                                                                                                                                            >  1:618604/1‑120 (MQ=255)
                                                                                                             ccAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGa                                                                                                                                                                                                            <  2:618604/120‑1 (MQ=255)
                                                                                                                    ggTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGa                                                                                                                                                                                                            <  1:105631/113‑1 (MQ=255)
                                                                                                                    ggTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGa                                                                                                                                                                                                            >  2:105631/1‑113 (MQ=255)
                                                                                                                     gTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGGAAAACATCCACGATTTTCATTGCCTTTTTTTTCCTATTTGTGAccc                                                                                                                                                                                                         >  2:855998/1‑114 (MQ=255)
                                                                                                                     gTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTTCCTTGCATTTTTTTTCCTATGTGTGACCACAAAGGGGTGGGTGAAAC‑CAACTATTTTTCTGGGGTTTTTCAGCAGgtttttttttgggg                                                                                                                                            >  2:636258/1‑169 (MQ=255)
                                                                                                                     gTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGa                                                                                                                                                                                                            <  1:384039/112‑1 (MQ=255)
                                                                                                                     gTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCa                                                                                                                                      <  1:680187/181‑1 (MQ=255)
                                                                                                                               aaTAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCAACCACGATTTTCATGGCATTTTTTTCCCTATGTGTGACCACAATGGGGTGGGGAAAAA‑CATCTATTTTAATGGATTTTTTCAGCAGTTTTTTTTTTGGTGGTGCCACGTTGTTTgg                                                                                                                            >  2:680187/1‑181 (MQ=255)
                                                                                                                               aaTAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAACATCCACGATTTTCCTGGCATTTTTTTTCGTATGTGTGACCACAAtggg                                                                                                                                                                                                  >  2:384039/1‑112 (MQ=255)
                                                                                                                                     cTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATc                                                                                                                                                                                  <  1:80695/121‑1 (MQ=255)
                                                                                                                                      ttttAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGGGAAAACCATCCCAGATTTTCATGGCATTTTTTTCCGTTTGGGGGACCCAAATGGGTTGGGGAAAAC‑CATTTTTTTCAATGGGTTTTTGCAGCAGTTTTTTTTCTGGTGGTGCCACGTTGTTTGGGTTGTTATATTTTGGCATAAGTTGGCTGGTTTTTAACCACAGAAAGCAAACCACT‑TGGCCTTCAc                                                          >  2:169171/1‑239 (MQ=255)
                                                                                                                                               tAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACCATTTCCCATGCATTTTTTTCCGTATTGTTGACCACAACGGGGTGGGTAAAAC‑CATCTATTTTAATg                                                                                                                                                                        >  2:80695/1‑121 (MQ=255)
                                                                                                                                               tAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAACATCCCCGATTTTCCTTGCATTTTTTTTCGTATGTGTGAACACAATGGGGTGGGGAAAAA‑CATCTATTTTCATTGGTTTTTGTAGCAGGTTTTTTTTGGGGGGGGGCACGGTGTTTTGGTTGTAATTATTAAGGTTAAGTTGGCTGTTTTTTAACCCCCGAACGGCATCCAATAGGCTTTTTCGGAAAgggg                                                  >  2:499395/1‑239 (MQ=255)
                                                                                                                                                aaTTTTAATGATGATGGGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACT‑GGCTCTTCACGAAAGGGGCt                                                <  1:169171/239‑1 (MQ=255)
                                                                                                                                                aaTTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACT‑GGCTCTTCACGAAAGGGGCt                                                <  1:439916/239‑1 (MQ=255)
                                                                                                                                                  ttttAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACcaca                                                                                                                                                                                                       <  1:777465/88‑1 (MQ=255)
                                                                                                                                                            gatgaGCCTTTACCAATTTCATGGCGGTGAAAAACAACCACGATTTTCATGGCATTTTTTTCCCTATGTGTGACCACAAtgggttggg                                                                                                                                                                                             >  2:777465/1‑88 (MQ=255)
                                                                                                                                                                       aCCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACT‑GGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGa                         <  1:499395/239‑1 (MQ=255)
                                                                                                                                                                            tttCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCt                                                                                                                                                                                 <  1:302796/83‑1 (MQ=255)
                                                                                                                                                                               cATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATgg                                                                                                                                                                       >  1:773782/1‑90 (MQ=255)
                                                                                                                                                                               cATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATgg                                                                                                                                                                       <  2:773782/90‑1 (MQ=255)
                                                                                                                                                                                  ggCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAAtggg                                                                                                                                                                                                  >  1:726173/1‑61 (MQ=255)
                                                                                                                                                                                  ggCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAAtggg                                                                                                                                                                                                  <  2:726173/61‑1 (MQ=255)
                                                                                                                                                                                      gTGAAAACCATCCCCGATTTCCATGGCATTTTTTTCCGTATGTGTGACCACAAAGGGTTGGGGAAAAC‑CATCTATTTCAATgg                                                                                                                                                                       >  2:302796/1‑83 (MQ=255)
                                                                                                                                                                                       tGAAAACCATCCACGATTTCCATGGCATTTTTATc                                                                                                                                                                                                                       <  2:970128/35‑1 (MQ=255)
                                                                                                                                                                                        taaaaCCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACT‑GGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTg        <  1:909318/238‑1 (MQ=255)
                                                                                                                                                                                              cATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACT‑GGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGa  <  1:528522/239‑1 (MQ=255)
                                                                                                                                                                                                 ccACGATTTCCATGGCATTTTTATCCGTAtgtgtg                                                                                                                                                                                                             >  1:970128/1‑35 (MQ=255)
                                                                                                                                                                                                        ttCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGttt                                                                                                                              >  1:984249/1‑106 (MQ=255)
                                                                                                                                                                                                        ttCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGttt                                                                                                                              <  2:984249/106‑1 (MQ=255)
                                                                                                                                                                                                                       |                                                                                                                                                                                                                       
CCCAACTGGCTAAGTTCACGGACGATTTCCGACATGGCGGTGATTTGATGGCGTCCACGAGCGCGATTATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAAC‑CATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACT‑GGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGA  >  NC_000913/634877‑635305

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