Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,622,8020AT40.6% 46.3 / 13.2 32intergenic (+61/+158)yneM/dgcZinner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base A (17/2);  new base T (0/13);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 3.02e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAATATGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAAT  >  NC_000913/1622653‑1622947
                                                                                                                                                          |                                                                                                                                                    
gTAATATGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGAt                                                                                                                                                     >  2:463502/1‑151 (MQ=255)
    tatGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCCGAttttt                                                                                                                                                 >  1:398673/1‑151 (MQ=255)
             ttttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCg                                                                                                                                      >  2:263432/1‑151 (MQ=255)
              tttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCg                                                                                                                                      >  2:397319/1‑150 (MQ=255)
                            gggAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGcc                                                                                                                        <  2:589289/149‑1 (MQ=255)
                                 taatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAAc                                                                                                                 >  2:292290/1‑151 (MQ=255)
                                    ttttttttttGGTTTTTCGGGCCGCTTTTTCAACCCCAAAATGGGTTGAATAAGGAAGGAAAAAAACCCCTCCCAAACCCGCCCCTTTTATAAATTGGTGAACAGGGCCCTGGATTTTTATTACGGCAAAAAAAAACCGCCAGTAAACcg                                                                                                               <  1:404381/150‑1 (MQ=255)
                                         ttttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGt                                                                                                          >  1:276957/1‑150 (MQ=255)
                                           ttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGt                                                                                                          >  1:279271/1‑148 (MQ=255)
                                           ttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGa                                                                                                        >  2:327650/1‑150 (MQ=255)
                                                tttttcgggCCGCGTTTTTCACCCAAAAAGGGGATGAATAATAAAGGGAAAAACCCCCTCCCCAACCCGCCCCTTTGTTACGTTTGTGAAAAGGGGGCTTTTTTTTTAGGGGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATg                                                                                                     <  2:573147/147‑1 (MQ=255)
                                                     ctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTATAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                             >  2:459816/1‑151 (MQ=255)
                                                      tggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                            >  2:293199/1‑151 (MQ=255)
                                                            tgTATTCCACCCAAAAAGGGAATAATTAATGAAGGAAAATATTCCCCCCCCTACCCGCCCCTTTGTTTTATTTGTGTACGAGGGCCTTTTTTTAAATGGGGCAAAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATa                                                                                         <  2:590612/147‑1 (MQ=255)
                                                                ttccacccAAAAAGGGAATGACAAATAAAGGGAAAAACCCCCTCCCCAACCCGCCCCTTTGTTACGTTTGTTTAAAGGGGCTTTTTTTTTTAGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa                                                                                               <  2:170812/132‑1 (MQ=255)
                                                                   cAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGtt                                                                               >  2:162548/1‑151 (MQ=255)
                                                                         caAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt                                                                         >  1:28035/1‑151 (MQ=255)
                                                                           aaagggaatcataaaTGATGGGAGATATTCCCCCACCAACCCGGCCTTTTGTTAGGTTTGTGAACAGGGGGCTTTTTTTTAATGAAGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGtt                                                                               <  2:561741/131‑1 (MQ=255)
                                                                           aaaGGGGATAACAAATGAAGGGAAATAACCCCTCCCCAACCCGCCCCTTTTTTACATTTGTGAACAGGGGGCTTTTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGtt                                                                               <  2:498001/144‑1 (MQ=255)
                                                                                                gagaTAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTcc                                                    >  1:577717/1‑150 (MQ=255)
                                                                                                   aTAACCCCCCCCCAACCCGCCCCTTTGTTACGGTTTTGTAAAGGGGCCTTGTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAt                                                         <  2:388949/142‑1 (MQ=255)
                                                                                                                     ggCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGtt                                >  1:107166/1‑149 (MQ=255)
                                                                                                                       ccctTTTTTACTTTTGTGTAAAGGGGCCTTTTTTTTGTGGAGGGAAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGtt                                <  1:376442/147‑1 (MQ=255)
                                                                                                                          ttttttACTGTTTTGTAAAGGGGCCTTTTTTTTTAGGCGGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGaca                            <  2:2966/144‑1 (MQ=255)
                                                                                                                             gTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGt                        >  1:400042/1‑149 (MQ=255)
                                                                                                                                    gTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc                  >  1:632947/1‑149 (MQ=255)
                                                                                                                                     tggtgaACAAGGGGCTTTTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt                   <  1:463255/145‑1 (MQ=255)
                                                                                                                                         gtaaagggggCTTGTTTTTTTTGACGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtt              <  1:459168/142‑1 (MQ=255)
                                                                                                                                          aacggggggCGTGTTTTTTAGGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGCGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtc           <  2:519815/145‑1 (MQ=255)
                                                                                                                                          aaCAGGGGGCTTGTTTTTTAGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct          <  2:675/150‑1 (MQ=255)
                                                                                                                                           aCAGGGGGCCTTTTTTTTATGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCACGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATCTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct          <  2:58228/150‑1 (MQ=255)
                                                                                                                                                  gCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaat  >  2:282210/1‑151 (MQ=255)
                                                                                                                                                          |                                                                                                                                                    
GTAATATGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAAT  >  NC_000913/1622653‑1622947

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