Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,622,816	0	G	A	45.8%	19.6 / 51.3	48	intergenic (+75/+144)	yneM/dgcZ	inner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (25/1);  new base A (0/22);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 8.40e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAATGGGATGACTAAT‑GAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAA  >  NC_000913/1622669‑1622964                                                                                                                                                        |                                                                                                                                                      tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAATGGGATGACTAAT‑GAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGaa                                                                                                                                                     >  1:486670/1‑150 (MQ=255)     gCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAATGGGATGACTAAT‑GAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGaaaaaa                                                                                                                                                 >  1:809273/1‑150 (MQ=255)             gggAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAATGGGATGACTAAT‑GAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCa                                                                                                                                        <  2:439452/150‑1 (MQ=255)                     tttttttttGT‑GTTTTCCGGCCGTCTTATTACACCCAAAAAGGGGATGAAAAAA‑GAAGGGAAATATCCCCTCCCCAACCCCCCCCTTTGTTAAAT‑TTGTGAACAAGGGGCT‑TTTTTTTTATGACGGAAAAAAAAAA‑CCGCCAGTAAACcg                                                                                                                                <  2:764156/150‑1 (MQ=255)                        tattttCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAATGGGATGACTAAT‑GAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACcggcgg                                                                                                                            >  2:367901/1‑151 (MQ=255)                           tttctggttttctggCCG‑CGTATTTCAGCCACAAATGGGATGACTAAT‑GAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACcggcgg                                                                                                                            >  1:442849/1‑148 (MQ=255)                                 tttttctggCCG‑TTTTTTCCACCAAAAAAGGGGATGA‑TAATGGAAGGGAAATACCCCCCCCCCAACCCGCCCCTTTTTTAAAT‑TTGTGTA‑AAGGGGCCTTTTTTTTTAGGGGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGc                                                                                                                     <  1:655927/148‑1 (MQ=255)                                    ttctggCCG‑CGTATTTCAGCCACAAATGGGATGACTAAT‑GAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTg                                                                                                                  >  1:422757/1‑149 (MQ=255)                                       gggcgg‑cGTATTGCACCCACAAATGGATTTAATAAA‑GAAGGGATAAAAACCCACACAAAGCCGGCCCTTTTTTACGG‑GTGTATAAAGGGGGCG‑TGATTTTTATTAAGGAAAAAAAAAC‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                                             <  2:544614/150‑1 (MQ=255)                                          ccG‑CTTTTTTCACCCAAAAAGGGGATGACAAAT‑AAACGGAAATAACCCCCCCCCTACCCGCCCCCTTTTTACTG‑TTGTGAACAGGGGGCT‑TTTTTTTTAGGACGAAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                           <  1:558387/150‑1 (MQ=255)                                          ccG‑CTTTTTCCACCCAAAAAGGGAATGATTAAT‑GAACGGAAATACCCCCCCCCCAACCCGCCCCTTTGTTCCAG‑TTGTGAAAAGGGGGCT‑TTTTTTTTAGGACGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                           <  1:415521/150‑1 (MQ=255)                                          ccG‑CGTATTTCAGCCACAAATGGGATGACTAAT‑GAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGTCGGCGAAAACACA‑CCGCCAGTAAACCGGCGGCGAAAGCTTGCCTGGAt                                                                                                           >  1:394489/1‑150 (MQ=255)                                           cggcttATTCCACCCAAAAAGGGGATGAAAAAT‑AAAGGGAAAT‑ACCCCCCCCCCAACCGCCCCCTTTTTACCGTTTGTGAAAAGGGGGCG‑TGTTTTTTATGAGGCCAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                                             <  2:822708/143‑1 (MQ=255)                                            g‑cgTATTTCAGCCACAAATGGGATGACTAAT‑GAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCCTgg                                                                                                             >  2:795321/1‑146 (MQ=255)                                               gTATTTCAGCCACAAATGGGATGACTAAT‑GAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt                                                                                                      >  1:806821/1‑151 (MQ=255)                                                          aaaaaGGGGATGAATAAT‑GACCGGAAAAAACCCCCCCCCGACCCGGCCCTTTGTTAAGG‑TTGTGTAAAGGGGGCT‑TTTTTTTTTGGAGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCt                                                                                           <  1:577419/149‑1 (MQ=255)                                                            aaaTGGGATGACTAAT‑GAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCttt                                                                                         >  1:237863/1‑151 (MQ=255)                                                                gggATGACTAAT‑GAAGGGAAATAACCCCTCCCCAACCCGCCCCTTTTTTACAG‑TTGTGAACAGGGGGCC‑TGTTTTTTATGGCGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt                                                                                        <  1:157876/148‑1 (MQ=255)                                                                   aTGACTAAT‑GAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCt                                                                                   >  1:712087/1‑150 (MQ=255)                                                                    tGACTAAT‑GAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTa                                                                                  >  1:125247/1‑150 (MQ=255)                                                                     gACTAAT‑GAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                                >  1:646762/1‑151 (MQ=255)                                                                         aaT‑GAACGGAGATAATCCCTCACCTAACCGGCCACATGTTACAG‑TTGGGTACAAGGGGCC‑TGATTTTGATGACGGCGAAAAAAAA‑ACGCCAGTAAACCGGCGGTGAACGACTGCCTGGAACTAATTGTGTTTTGTTTGTAcaccaa                                                                                 >  2:339304/1‑141 (MQ=255)                                                                                 ggaaaaaaccccTCCCCAACCCGCCCCCTTGCTACGG‑TTGTGTACAGGGGGCT‑TGTTTTTTAGGCCGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCa                                                                           <  1:752156/136‑1 (MQ=255)                                                                                   agaTAATCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGGTAGATTTGTGTTTTGCTTTTACGCTAAAAGGCATTTTCCt                                                                    >  2:113353/1‑150 (MQ=255)                                                                                   agaTAACCCCCCACCTACCCGCCCCTTTTTTACTG‑TTGTTTACAAGGGGTT‑TGTTTTTTATGGCGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCt                                                                    <  1:250267/150‑1 (MQ=255)                                                                                       aaTCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACt                                                               >  1:737291/1‑151 (MQ=255)                                                                                           ccctcccctcccGGCCCCTTTTTAAGG‑TTTTTTAAAAGGGGCC‑TTTTTTTTAAGGCGGAAAAAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                            <  1:368018/141‑1 (MQ=255)                                                                                                     cgggCCCTTTGTTAAGT‑TTGTGAAAAGGGGGCT‑TGTTTTTTATGACGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGtt                                                 <  2:203069/149‑1 (MQ=255)                                                                                                      cGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                                >  2:634899/1‑151 (MQ=255)                                                                                                         ccttgtGTTACTG‑TTGGGAAAAAGGGGCT‑TTTTTTTTATGACGAAAAAAAAACT‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGaca                                             <  2:672190/146‑1 (MQ=255)                                                                                                         ccctTTTTTACAT‑TTGTGAAAAGGGGGCT‑TTTTTTTTATGACGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                                <  1:733942/148‑1 (MQ=255)                                                                                                                ttACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGc                                      >  1:129552/1‑151 (MQ=255)                                                                                                                ttACAG‑TTGTGAACAGGGGGCT‑TTTTTTTTATGACGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTa                                        <  1:175376/149‑1 (MQ=255)                                                                                                                        ttgtgaaaaGGGGGC‑TTTTTTTTTTGAAGAAAAAAAAAAC‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAg                                 <  2:347529/142‑1 (MQ=255)                                                                                                                         gtgaaaagggggCG‑TGTTTTTTATAACGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtt                               <  2:340528/142‑1 (MQ=255)                                                                                                                                 ggggCT‑TGTTTTTTAGGAGGAAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc                          <  1:740382/147‑1 (MQ=255)                                                                                                                                 ggggCCTTGTTTTTTAGGAGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtt                               <  1:600238/143‑1 (MQ=255)                                                                                                                                  gggccttgttttttAGGACGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatg                      <  1:673634/142‑1 (MQ=255)                                                                                                                                   ggCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga                     >  1:10234/1‑150 (MQ=255)                                                                                                                                   ggCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga                     >  2:687215/1‑150 (MQ=255)                                                                                                                                        tttttttttGGACGGCAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGt                 <  2:486522/147‑1 (MQ=255)                                                                                                                                          ttttttAGGAGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt                <  1:635093/148‑1 (MQ=255)                                                                                                                                           tttttATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaa              >  2:390213/1‑150 (MQ=255)                                                                                                                                              ttATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACgg          >  2:224358/1‑151 (MQ=255)                                                                                                                                               tATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACgg          >  1:616623/1‑150 (MQ=255)                                                                                                                                                aTGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGATTTTTCCTGCAATGATAACGAAACGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTCAACg           >  2:238811/1‑148 (MQ=255)                                                                                                                                                 tGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAAATTTGTGTTTTGATTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAg        >  2:847286/1‑150 (MQ=255)                                                                                                                                                      gCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGCATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTaa   >  2:402747/1‑150 (MQ=255)                                                                                                                                                        gAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTaaa  >  2:744948/1‑149 (MQ=255)                                                                                                                                                        |                                                                                                                                                      TATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAATGGGATGACTAAT‑GAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAA  >  NC_000913/1622669‑1622964

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