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

TTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGG  >  NC_000913/1622667‑1622956
                                                                                                                                                        |                                                                                                                                            
tttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCg                                                                                                                                              <  2:433622/150‑1 (MQ=255)
  tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaa                                                                                                                                            >  1:782817/1‑150 (MQ=255)
  tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaa                                                                                                                                            >  1:53762/1‑150 (MQ=255)
   aTGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaa                                                                                                                                          <  1:533969/151‑1 (MQ=255)
      gCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaagaaaa                                                                                                                                       <  2:756372/151‑1 (MQ=255)
         gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAAcc                                                                                                                                    >  1:590582/1‑151 (MQ=255)
          tACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCg                                                                                                                                   >  2:606717/1‑151 (MQ=255)
              gggAATAATTTTATTTTTTGGTTTTGTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGcc                                                                                                                                 >  2:88439/1‑149 (MQ=255)
                     attttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcg                                                                                                                        >  2:249265/1‑151 (MQ=255)
                        ttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGCTGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggcg                                                                                                                     >  2:795489/1‑151 (MQ=255)
                         ttttttttggtttcctggCCGCTTTTTCCACCCAAAAAGGGGATAAAAAAGGAACGGAAAAAACCCCCCACCTACCCGCCCCTTTTTAACAGTTGTGAAAAGGGGGCTTGATTTTTAAGAGGGCAAAAAAAAACCGCCAGTAAACcggcg                                                                                                                     <  2:461415/148‑1 (MQ=255)
                           ttttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGACTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTg                                                                                                                  >  2:548640/1‑151 (MQ=255)
                            tttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTg                                                                                                                  >  2:317989/1‑150 (MQ=255)
                            tttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTg                                                                                                                  >  2:478238/1‑150 (MQ=255)
                                   ttttctggCCGCTTATTACAGCAAAAAGGGGGATAACTAAAAAGAGAAAAAAACCCCCCCCCAACCCGCCCCTTTTTTAAATTTGTGAAAAGGGGGCTTGTTTTTTAGGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt                                                                                                           <  2:550201/150‑1 (MQ=255)
                                   ttttcgggCGGTGTATTCCAGCCACAAAGGGGATGATTAATGAACGGAAAAAACCCCCCCCCTACCCGCCCCTTTGTAAAAGTTGTATACAGGGGGCTTGTTTTTTAGGGGGGGAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt                                                                                                           <  1:225691/150‑1 (MQ=255)
                                     ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                         >  1:655633/1‑150 (MQ=255)
                                       ctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                                      >  1:725995/1‑151 (MQ=255)
                                         ggcggcGTATTTCACCCACAAAGGGGATGAAAAAGGAACGAAAAAACCCCCCCCCCTACCCGGCCCTTTGTTACAGTTGTGTACAGGGGGCTTGTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGa                                                                                                    <  2:13572/151‑1 (MQ=255)
                                         ggCCGCGTTTTACAGCCAAAAAGGGGATGACTAATGAAGGGAAATAACCCCTCCCCTACCCGGCCCTTTGTTACAGTTGTGTAAAGGGGGCTTGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                                      <  1:611192/149‑1 (MQ=255)
                                           ccGCTTTTTTCAGCCAAAAAGGGGATAACAAATAAAGGGAAATAACCCCCCCCCAACCCGCCCCTTTTTTACGTTTGTGTACAGGGGGCTTTTTTTTTAGGCCGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                   <  1:637521/150‑1 (MQ=255)
                                                  tttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACCAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAGAGATTtgt                                                                                           >  1:113236/1‑151 (MQ=255)
                                                  ttccAGCCAAAAGTGGGATGAAAAATGAAGGGAAATACCCCCTCCCCTACCCGCCCCTTTTTTACAGTTGTGAAAAGGGGGCTTGTTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt                                                                                             <  2:612800/146‑1 (MQ=255)
                                                            aaaGGGAATAACAAATGAAGGAAAAAAACCCCCCCCCTACCCGCCCCTTTTTTACATTTGTGTAAAGGGGGCTTTTTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTg                                                                                     <  2:793446/147‑1 (MQ=255)
                                                             aaTGGGATGATTAATGAAGGGAAATATTCCCTCCCCGAGCCGCCCCTTTTTTAGATTGTTGCAAGGGGGGCTTTTTTTTTAGGCCGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt                                                                                <  1:161521/151‑1 (MQ=255)
                                                                 ggATAACTAATGAAGGGAAATAACCCCTCCCCTACCCGCCCCTTTGTTACAGTTGTGTACAGGGGGCTTGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGc                                                                            <  1:422080/151‑1 (MQ=255)
                                                                       cTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                        >  1:184233/1‑149 (MQ=255)
                                                                        tAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAgg                                                                     >  1:49708/1‑151 (MQ=255)
                                                                                ggAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtttt                                                               >  1:286386/1‑149 (MQ=255)
                                                                                   gaTAATCCCTCACCTATCCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc                                                          >  1:808933/1‑151 (MQ=255)
                                                                                         ccccccccTACCCGCCCCTTTTTTACAGTTGTTAAAAGGGGCCTTTTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc                                                          <  2:95506/139‑1 (MQ=255)
                                                                                                  acccGCCCCTTTTTTACATTTGTGAAAAGGGGGCCTGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCg                                           <  2:263828/149‑1 (MQ=255)
                                                                                                   aCCCGCCCCTTTTTAACTGTTGTTAAAAGGGGGCTTTTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAAc                                                 <  2:768841/145‑1 (MQ=255)
                                                                                                           cTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGacac                                    >  2:582028/1‑149 (MQ=255)
                                                                                                             ttGTTACAGTTGTGAACAGGGGGCTTGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAg                                  <  1:548695/150‑1 (MQ=255)
                                                                                                              tGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGt                                 >  2:240992/1‑150 (MQ=255)
                                                                                                                ttACAGTTGTTTAAAGGGGGCCTGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAg                               <  1:699298/150‑1 (MQ=255)
                                                                                                                ttACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAg                               >  2:407824/1‑150 (MQ=255)
                                                                                                                        tgtgTAAAAGGGGCTTGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGt                        <  2:782319/150‑1 (MQ=255)
                                                                                                                             aCAAGGGGCGTGTTTTTTAGGACGGAAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc                  <  2:26546/151‑1 (MQ=255)
                                                                                                                              aaaggggcgttgttttttAGGACGGTAAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct                   <  2:725924/137‑1 (MQ=255)
                                                                                                                               aaGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                 >  2:332238/1‑150 (MQ=255)
                                                                                                                                agggggctttttttttAGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGttt                      <  1:36506/134‑1 (MQ=255)
                                                                                                                                   ggCCTGATTTTTATGACGGCGAGAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga             >  1:651378/1‑150 (MQ=255)
                                                                                                                                    gCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaa            >  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)
                                                                                                                                                        |                                                                                                                                            
TTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGG  >  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: