Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,622,8160GA52.2% 6.3 / 32.1 23intergenic (+75/+144)yneM/dgcZinner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (10/1);  new base A (0/12);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 9.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-04

GGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACT  >  NC_000913/1622681‑1622966
                                                                                                                                         |                                                                                                                                                       
gggAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGca                                                                                                                                           <  1:214621/149‑2 (MQ=255)
           tattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACcggcgg                                                                                                                              >  1:165885/1‑151 (MQ=255)
               ttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGa                                                                                                                           >  2:337419/1‑150 (MQ=255)
                             ccggattttccacccAAAAAGGGAATGAATAAGGAAGGGAGATACCCCCCCACCAACCCGCCCCTTTGTTAAT‑TTTGTGAAAAGGGGTCC‑TTTTTTTTAGGAGGACA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa                                                                                                                  <  1:534302/132‑1 (MQ=255)
                                  tATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt                                                                                                       >  2:86700/1‑151 (MQ=255)
                                            aaaagggggATGACTAATGAAGGGAGATACCCCCCCCCCTACCCGCCCCCTTTTTACA‑GTTGTGAAAAGGGGCCC‑TTTTTTTTATGCGGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTg                                                                                               <  1:381340/143‑1 (MQ=255)
                                             aaaaGGGGATGATAAATGAAGGGAAATAACCCCCCCCCTACCCGCCCCTTTTTTACG‑TTTTTGTACAGGGGGCT‑TGTTTTTTAGGCCGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtgtg                                                                                                    <  1:219053/142‑1 (MQ=255)
                                                        aCTAAGAAAGGGAAAAAACCCCCCACCAACCCGCCCCTTTTTTACA‑TTGGTGAAAAGGGGGTG‑TGTTTTTAAGGCGGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCt                                                                                     <  1:578656/147‑1 (MQ=255)
                                                                 cGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTc                                                                        >  2:242094/1‑151 (MQ=255)
                                                                    agaTAACCCCTCACCAACCCGCCCCTTTTTTACA‑GTTGTGTAAAGGGGGCG‑TGTTTTTTAGGAGGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTc                                                                        <  2:80149/148‑1 (MQ=255)
                                                                           cccTCCCCTACCCGCCCCTTTTTTACA‑TTTGGTAACAGGGGGCT‑TGTTTTTTAGGACGAAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                              <  2:203795/151‑1 (MQ=255)
                                                                                 cccaaCCCGCCCCTTTTTAAATGTTGTG‑AAAAGGGGCCTTTTTTTTTAGGAGGGAA‑AAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATa                                                             <  2:133067/144‑1 (MQ=255)
                                                                                    acccGCCCCTTTTTTACA‑TTTGTAAACAGGGGCCT‑TTTTTTTTATGACGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa                                                         <  1:220763/145‑1 (MQ=255)
                                                                                       cgcccctTTGTTACT‑GTGGTGAACAAGGGGCT‑GTTTTTTTATGAGGAAA‑AAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCg                                                     <  2:225149/145‑1 (MQ=255)
                                                                                        ggCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGtt                                                   >  2:375815/1‑149 (MQ=255)
                                                                                                    cA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc                                     >  2:191625/1‑151 (MQ=255)
                                                                                                         tgtgTACAGGGGCCT‑TGTTTTTTAGGACGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGc                                        <  1:470157/144‑1 (MQ=255)
                                                                                                                  ggggCT‑TGTTTTTTAGGACGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                           <  1:406842/148‑1 (MQ=255)
                                                                                                                       cttgttttttAGGAGGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaat                     <  1:604003/145‑1 (MQ=255)
                                                                                                                         tGATTTTTATGACGGCGAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGt                   >  2:369444/1‑151 (MQ=255)
                                                                                                                             ttttATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaaa               >  2:143351/1‑150 (MQ=255)
                                                                                                                                   gACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCt        >  1:557367/1‑151 (MQ=255)
                                                                                                                                        cG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAAc   >  2:626405/1‑151 (MQ=255)
                                                                                                                                         g‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCCTGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACt  >  2:74302/1‑151 (MQ=255)
                                                                                                                                         |                                                                                                                                                       
GGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACT  >  NC_000913/1622681‑1622966

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: ‑