Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,622,8160GA42.5% 40.9 / 36.4 47intergenic (+75/+144)yneM/dgcZinner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (25/2);  new base A (0/20);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 2.37e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGG  >  NC_000913/1622667‑1622956
                                                                                                                                                        |                                                                                                                                            
tttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCg                                                                                                                                              <  2:433622/150‑1 (MQ=255)
  tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGaa                                                                                                                                            >  1:782817/1‑150 (MQ=255)
  tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGaa                                                                                                                                            >  1:53762/1‑150 (MQ=255)
   aTGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGaaaa                                                                                                                                          <  1:533969/151‑1 (MQ=255)
      gCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGaagaaaa                                                                                                                                       <  2:756372/151‑1 (MQ=255)
         gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAAcc                                                                                                                                    >  1:590582/1‑151 (MQ=255)
          tACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCg                                                                                                                                   >  2:606717/1‑151 (MQ=255)
              gggAATAATTTTATTTTTTGGTTTTGTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGcc                                                                                                                                 >  2:88439/1‑149 (MQ=255)
                     attttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcg                                                                                                                        >  2:249265/1‑151 (MQ=255)
                        ttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGCTGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggcg                                                                                                                     >  2:795489/1‑151 (MQ=255)
                         ttttttttggtttcctggCCGCTTTTTCCACCCAAAAAGGGGATAAAAAAGGAACGGAAAAAACCCCCCACCTACCCGCCCCT‑TTTTA‑ACAGTTGTGAAAAGGGGGCT‑TGATTTTTAAGAGGGCAAAAAAAAACCGCCAGTAAACcggcg                                                                                                                     <  2:461415/148‑1 (MQ=255)
                           ttttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGACTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTg                                                                                                                  >  2:548640/1‑151 (MQ=255)
                            tttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTg                                                                                                                  >  2:317989/1‑150 (MQ=255)
                            tttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTg                                                                                                                  >  2:478238/1‑150 (MQ=255)
                                   ttttctggCCGCTTATTACAGCAAAAAGGGGGATAACTAA‑AAAGAGAAAAAAACCCCCCCCCAACCCGCCCCTTTTTT‑AAATTTGTGAAAAGGGGGCT‑TGTTTTTTAGGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt                                                                                                           <  2:550201/150‑1 (MQ=255)
                                   ttttcgggCGGTGTATTCCAGCCACAAAGGGGATGATTAATGAACGGAAAAAACCCCCCCCCTACCCGCCCCT‑TTGTA‑AAAGTTGTATACAGGGGGCT‑TGTTTTTTAGGGGGGGAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt                                                                                                           <  1:225691/150‑1 (MQ=255)
                                     ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                         >  1:655633/1‑150 (MQ=255)
                                       ctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                                      >  1:725995/1‑151 (MQ=255)
                                         ggcggcGTATTTCACCCACAAAGGGGATGAAAAAGGAACGAAAAAACCCCCCCCCCTACCCGGCCCT‑TTGTT‑ACAGTTGTGTACAGGGGGCT‑TGTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGa                                                                                                    <  2:13572/151‑1 (MQ=255)
                                         ggCCGCGTTTTACAGCCAAAAAGGGGATGACTAATGAAGGGAAATAACCCCTCCCCTACCCGGCCCT‑TTGTT‑ACAGTTGTGTAAAGGGGGCT‑TGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                                      <  1:611192/149‑1 (MQ=255)
                                           ccGCTTTTTTCAGCCAAAAAGGGGATAACAAATAAAGGGAAATAACCCCCCCCCAACCCGCCCCT‑TTTTT‑ACGTTTGTGTACAGGGGGCT‑TTTTTTTTAGGCCGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                   <  1:637521/150‑1 (MQ=255)
                                                  tttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACCAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAGAGATTtgt                                                                                           >  1:113236/1‑151 (MQ=255)
                                                  ttccAGCCAAAAGTGGGATGAAAAATGAAGGGAAATACCCCCTCCCCTACCCGCCCCT‑TTTTT‑ACAGTTGTGAAAAGGGGGCT‑TGTTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt                                                                                             <  2:612800/146‑1 (MQ=255)
                                                            aaaGGGAATAACAAATGAAGGAAAAAAACCCCCCCCCTACCCGCCCCT‑TTTTT‑ACATTTGTGTAAAGGGGGCT‑TTTTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTg                                                                                     <  2:793446/147‑1 (MQ=255)
                                                             aaTGGGATGATTAATGAAGGGAAATATTCCCTCCCCGAGCCGCCCCT‑TTTTT‑AGATTGTTGCAAGGGGGGCT‑TTTTTTTTAGGCCGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt                                                                                <  1:161521/151‑1 (MQ=255)
                                                                 ggATAACTAATGAAGGGAAATAACCCCTCCCCTACCCGCCCCT‑TTGTT‑ACAGTTGTGTACAGGGGGCT‑TGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGc                                                                            <  1:422080/151‑1 (MQ=255)
                                                                       cTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                        >  1:184233/1‑149 (MQ=255)
                                                                        tAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAgg                                                                     >  1:49708/1‑151 (MQ=255)
                                                                                ggAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtttt                                                               >  1:286386/1‑149 (MQ=255)
                                                                                   gaTAATCCCTCACCTATCCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc                                                          >  1:808933/1‑151 (MQ=255)
                                                                                         ccccccccTACCCGCCCCT‑TTTTT‑ACAGTTGT‑TAAAAGGGGCCTTTTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc                                                          <  2:95506/139‑1 (MQ=255)
                                                                                                  acccGCCCCT‑TTTTT‑ACATTTGTGAAAAGGGGGCC‑TGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCg                                           <  2:263828/149‑1 (MQ=255)
                                                                                                   aCCCGCCCC‑TTTTTAACTGTTGTTAAAAGGGGGCT‑TTTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAAc                                                 <  2:768841/145‑1 (MQ=255)
                                                                                                           c‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGacac                                    >  2:582028/1‑149 (MQ=255)
                                                                                                             ttGTT‑ACAGTTGTGAACAGGGGGCT‑TGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAg                                  <  1:548695/150‑1 (MQ=255)
                                                                                                              tGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGt                                 >  2:240992/1‑150 (MQ=255)
                                                                                                                tt‑ACAGTTGTTTAAAGGGGGCC‑TGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAg                               <  1:699298/150‑1 (MQ=255)
                                                                                                                tt‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAg                               >  2:407824/1‑150 (MQ=255)
                                                                                                                        tgtgTAAAAGGGGCT‑TGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGt                        <  2:782319/150‑1 (MQ=255)
                                                                                                                             aCAAGGGGCG‑TGTTTTTTAGGACGGAAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc                  <  2:26546/151‑1 (MQ=255)
                                                                                                                              aaaggggcgttgttttttAGGACGGTAAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct                   <  2:725924/137‑1 (MQ=255)
                                                                                                                               aaGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                 >  2:332238/1‑150 (MQ=255)
                                                                                                                                agggggctttttttttAGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGttt                      <  1:36506/134‑1 (MQ=255)
                                                                                                                                   ggCC‑TGATTTTTATGACGGCGAGAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga             >  1:651378/1‑150 (MQ=255)
                                                                                                                                    gCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaa            >  1:786770/1‑150 (MQ=255)
                                                                                                                                      cttgttttttATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaa            <  2:61244/144‑1 (MQ=255)
                                                                                                                                        tGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt        >  2:507548/1‑151 (MQ=255)
                                                                                                                                              ttATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACgg  >  1:884306/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                            
TTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGG  >  NC_000913/1622667‑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: ‑