Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,622,8150CA46.0% 39.3 / 21.2 38intergenic (+74/+145)yneM/dgcZinner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base C (19/1);  new base A (1/16);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 2.14e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGG  >  NC_000913/1622675‑1622956
                                                                                                                                                 |                                                                                                                                             
cGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAAc                                                                                                                                     >  1:1079555/1‑151 (MQ=255)
 gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAAcc                                                                                                                                    >  2:687624/1‑151 (MQ=255)
    cTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGcc                                                                                                                                 >  2:1098071/1‑151 (MQ=255)
      gggAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCa                                                                                                                                <  1:258755/150‑1 (MQ=255)
                     ttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGa                                                                                                                 >  2:299659/1‑150 (MQ=255)
                             ttctggCCGTTTTTTCCACCCAAAAAGGGAATAATAAAAGACGGGAAATAACCCCCCCCCAACCCGCCCCTTTT‑TT‑ACTT‑TTTTGTACAGGGGG‑CCTTT‑TTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                         <  2:1137844/150‑1 (MQ=255)
                              tctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCAt                                                                                                       >  2:917781/1‑151 (MQ=255)
                                tggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                                     >  2:759798/1‑151 (MQ=255)
                                 ggCCGCTTTTTCCAGCCAAAAGGGGGAGAACAAATGAGGGAAAATACCCCCCCCCCAACCCGCCCCTTTG‑TT‑ACAG‑TTGTGTAAAGGGGG‑CCTTT‑TTTTTATGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCAt                                                                                                       <  2:946541/148‑1 (MQ=255)
                                 ggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                                     >  2:412953/1‑150 (MQ=255)
                                    cgcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGa                                                                                                    >  2:750155/1‑148 (MQ=255)
                                     gcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                   >  2:703908/1‑148 (MQ=255)
                                        ttttcaacccAAAAAGGGGATAAAAAATAAACGGAAATACCCCCCCCCCAACCCGCCCCTTTT‑TT‑ACGG‑TTGTAAACAGGGGG‑CCTGT‑TTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt                                                                                              <  2:569536/142‑1 (MQ=255)
                                        tATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑CTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt                                                                                             >  1:1023607/1‑151 (MQ=255)
                                         aTTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGGGAATGCTTGCATGGATAGAttt                                                                                             >  2:699437/1‑150 (MQ=255)
                                                        gggATGAAAAATGAACGGAGATACCCCCTCCCCTACCCGCCCCTTTT‑TT‑ACAT‑TTGTGTAAAGGGGG‑CTTTT‑TTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt                                                                                <  2:772306/148‑1 (MQ=255)
                                                        gggATAATTAATGAAGGGAAAT‑ACCCCCCCCCCAACCCGCCCCTTT‑TT‑TCAGTTTGTGAAAAGGGGG‑CTTTT‑TTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt                                                                                <  2:785273/148‑1 (MQ=255)
                                                         ggtggaaaagtGAAGGGATATACCCCCCCCCCTACCCGCCCCCTTT‑TT‑ACTG‑TTG‑GCAAAAGGGGTTTTGT‑TTTTTATGGGGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGc                                                                            <  1:703494/141‑1 (MQ=255)
                                                              aCTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                        >  1:968943/1‑150 (MQ=255)
                                                                       cGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtttt                                                               >  2:989082/1‑150 (MQ=255)
                                                                           gaTAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc                                                          >  2:704536/1‑151 (MQ=255)
                                                                                 ccccccccTACCCGCCCCTTTG‑TT‑ACAG‑TTGTTTACAGGGGG‑CTTTT‑TTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGa                                                     <  2:436960/144‑1 (MQ=255)
                                                                                  cccccaccaaccctccccTTGTTT‑AAAG‑TTGGTTACAGGGGT‑CTTTA‑TTTTTAGGGGGAAAAAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTg                                                      <  1:418999/135‑1 (MQ=255)
                                                                                       ccTAACCGCCCCTTTG‑TT‑ACAT‑TTGTGAACAGGGGG‑CTTGT‑TTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGaa                                              <  2:817180/151‑1 (MQ=255)
                                                                                        cTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGaa                                              >  2:550006/1‑150 (MQ=255)
                                                                                               gCCCCTTT‑TTAACGG‑TTGTGAAAAGGGGG‑CTTTT‑TTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                        <  2:894433/150‑1 (MQ=255)
                                                                                               gCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGac                                      >  2:838453/1‑151 (MQ=255)
                                                                                                  ccTTG‑TT‑ACAG‑TTTTGTACAAGGGG‑CCTGA‑TTTTTATGACGGGGCTAAAAAAACGCCAGTCAACCGGCGGTGAATGCTAGTACGGAAAGATTTGTGTTTTGCTTTTACGCCAACAGGGATTTTCCTCCACTGATAACGAATCGTTTactc                                    >  1:626010/1‑148 (MQ=255)
                                                                                                       ttttaaaatttttgtAAAGGGGG‑CTTGT‑TTTTTAGGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGAcaca                                   <  2:542073/137‑1 (MQ=255)
                                                                                                        tt‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGc                              >  1:111429/1‑151 (MQ=255)
                                                                                                               gttttgaagaaGGGG‑CTTGT‑TTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt                            <  1:1007457/139‑1 (MQ=255)
                                                                                                                ttgtaaaaaaGGGC‑CTTGT‑TTTTTAAGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt                            <  2:26329/139‑1 (MQ=255)
                                                                                                                 ttgggaaaaGGGG‑CCTTTTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc                           <  2:563835/141‑1 (MQ=255)
                                                                                                                        aggggg‑CCTTT‑TTTTTATGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                 <  2:3807/148‑1 (MQ=255)
                                                                                                                          gggg‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGGATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCCGTTTTCTCa                 >  1:858683/1‑148 (MQ=255)
                                                                                                                           ggg‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatg              >  2:687958/1‑150 (MQ=255)
                                                                                                                                      ttttAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaaa     <  2:750953/150‑1 (MQ=255)
                                                                                                                                        ttATGACGGAGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGAGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACgg  >  1:409342/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGG  >  NC_000913/1622675‑1622956

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