Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	523,387	0	A	G	43.8%	57.4 / 13.5	48	M43V (ATG→GTG)	rhsD	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (8/19);  new base G (21/0);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 3.77e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTT  >  NC_000913/523242‑523535                                                                                                                                                   |                                                                                                                                                       aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAAATTGGGCGCCCACCGGCGTGGCGGGCTCGGTGTGT‑CCGGGGGGGGTgg                                                                                                                                                     >  2:270189/1‑148 (MQ=255)  aCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGA‑Ctt                                                                                                                                                 <  1:268556/151‑1 (MQ=255)          ggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCa                                                                                                                                           <  1:320144/149‑1 (MQ=255)              cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCCGCGTGGCGGGCCCGGGGTGT‑CCGGGCGGGGTGA‑CTTCGGGCAAACCgg                                                                                                                                     >  1:946132/1‑151 (MQ=255)              cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGGAGGGGTAAGAATTGGCGCGCCCACCGGGGGGGCGGGGCCGGTGTGT‑CCGGGCGGGGGGA‑GTTCGGGGAAACCgg                                                                                                                                     >  1:1055647/1‑151 (MQ=255)              cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGGAGGGGTAAGAAATGGCGCGCCCACCGGCGTGGCGGGCCCGGGGGGT‑CCCGGCCGGGTGA‑CTTTGGGGAAacccg                                                                                                                                     >  1:1021368/1‑149 (MQ=255)              cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTTGGGCGCCCACCGGCGTGGCGGGCGCGGGGTGT‑CCGGGGGGGGTGG‑CTTCGGGCAAccccg                                                                                                                                     >  1:217395/1‑149 (MQ=255)              cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGGGCGCCCCCCGGCGTGGCGGGCTCGGGGTTT‑GCGGGCGGGGTGA‑GTTCGGGCaaacagg                                                                                                                                     >  2:974815/1‑146 (MQ=255)              cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGGGCTCGGGGTGT‑CCCGGGGGGGTGA‑CTTCGGGGAACCCgg                                                                                                                                     >  1:195643/1‑151 (MQ=255)              cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TCGGCAGGGGTAAGAATTGGCGCGCCCCCCGGCGGGGCGGGCCCGGGGGGT‑CCGGGGGGGGTGA‑CTTCGGGCAACCCggg                                                                                                                                    >  1:898273/1‑150 (MQ=255)              cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAATGTCCCAGGGTCGGCAGGGGTAAAAATTGGCGCGCCCACCGGCGTGGGCTGGCCGGGTTGT‑CCTGGCGGGGAGACCTTCGGGAAAccgg                                                                                                                                      >  2:854398/1‑149 (MQ=255)              cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAATGTCCAGGGGTCGGGAGGGGTAAGAATTTGCGCGCCCACCGGCGTGGCGGGCTCGGTGTGG‑CCGGGGGGGGGGG‑CGTCCGGGAACCCgg                                                                                                                                     >  2:216124/1‑151 (MQ=255)              cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAATGTCCAGGGGTCGGCAGGTGTAAGGATTGGCGCGCCCACCCGCGTGGCCTGCTCGGTGTGT‑GCCGGGGGGGTTC‑CATTGGGCAAccccg                                                                                                                                     >  2:526544/1‑149 (MQ=255)                ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGGACCTTCTTTCAGGGGTGTGTAGGGGTGAGAATGTGGGCCCCCACCGGTGCTGCGTGCCCGGGGCGT‑CGCGGGGGGGTGA‑GGTTGGGGAACCCggga                                                                                                                                   >  2:455016/1‑149 (MQ=255)                 tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑aa                                                                                                                                 <  2:800009/151‑1 (MQ=255)                 tataTGAGCGGAAAACCAGCGGCGCGTCAGAGAGATAAGACTATGTATGGCATTCCCATGGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTTCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑aa                                                                                                                                 <  2:517269/151‑1 (MQ=255)                    atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATc                                                                                                                              <  2:312813/151‑1 (MQ=255)                           gAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGATCGGCAGGGGTAAGAAATTGCGCGCCCCCCGGTGTGGCGTGGTCGGGGAGG‑CCGGGCGGGGTGA‑CCTTGGGGAACACGGTGAAACCCgcgc                                                                                                                         >  2:1144441/1‑148 (MQ=255)                           gAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCATGTCCAGGGGTCCGCAGGGGTAAGAATTGGGGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGGTGA‑CTTCGGGCAACCCGGT‑GAATCCgcttc                                                                                                                        >  2:105951/1‑148 (MQ=255)                                                 agaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGGGGGT‑CCGGGCGGGGTGA‑CTTCGGGGAAACCGGG‑AAAACCGCTTCTGGGGGGGGAGGGGGCGCCCgg                                                                                                  >  1:627824/1‑150 (MQ=255)                                                  gaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGGGGGT‑CCGGGCGGGGTGA‑CTTCGGGGAACCCGGG‑AAATCCGCCGCTGGGGGGGGAGGGGCTGCCCGGCg                                                                                                >  1:940126/1‑151 (MQ=255)                                                   atatGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGCCGGTGTAAGAATTTGCGCGCCCACCCGCCGGGCCTGCGCCGGGGGG‑CCCGGCGGGGTGA‑CTTCGGGGAACCCCGT‑AAAACCCCTGCTGGGGGGGGAGGGGTGGGCCGGCgg                                                                                               >  2:791824/1‑150 (MQ=255)                                                   atatGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGTGCTCGGGGTGT‑CCGGGCGGGGTGA‑CTTCGGGCAACCCGGT‑AAATCCCCTGCTGGGGGGGGAGGGGCTGCCCGGCgg                                                                                               >  1:488359/1‑150 (MQ=255)                                                        aCTCAGTATGGCGGTCCCATTTTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCCCCGGCGTGGGGGGCTCGGTGTGTGCGGGGGGGGGGGA‑TTCCGGGCACCCCGG‑‑AAATTCGCTGCTGGGGGGGGAGGGGGTGCCCGGGGGGACgg                                                                                          >  2:140720/1‑151 (MQ=255)                                                        aCTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGGTGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGGGCTGCCCGGCGGGACgg                                                                                          >  1:946786/1‑151 (MQ=255)                                                          tgaGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACGGGCGTGGCGTGCGCGGTGTGT‑CCGGGGGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAGGGTGCTGCCCGGCGAGACGGAc                                                                                        <  2:852622/149‑1 (MQ=255)                                                            aGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGT‑CCGGGGGGGGGGA‑ATTCGGGGAACCCGGG‑AAATCCCCTGCTGGGGGGGGAGGGGGTGCCCgggggggcgggc                                                                                        >  1:689092/1‑139 (MQ=255)                                                                tGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgcg                                                                                  <  1:33459/151‑1 (MQ=255)                                                                     gTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCCCCGGCGTGGCGTGCTCGGTGTGT‑TCGGGCGGGATGA‑CCTCGGGCAACCCGGC‑AAATCCGCTGTTGGGGGCGAAGGCTCCGCTCGGCGAGCCCGACCTTGCGCTgc                                                                              >  2:1142248/1‑150 (MQ=255)                                                                       cccATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccg                                                                           >  1:174827/1‑151 (MQ=255)                                                                       cccATTGTCCAGGGTTCGGCAGGCGTAAGAATTGGCGCCCCCACCGGCGCGGCGCGCTCGGTGTGT‑CCGGGCGGGAGGA‑CCTCGGGCAACCCGAT‑AAATCCGCGGCTGGGGGCGAAGGTGATTCCCGGCGAGCCGGACCTTGCGCTccccg                                                                           >  2:490293/1‑151 (MQ=255)                                                                             gTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggcccg                                                                      >  1:952370/1‑150 (MQ=255)                                                                                     ttCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTc                                                             <  1:757783/151‑1 (MQ=255)                                                                                             gtgtAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctct                                                      >  1:34279/1‑150 (MQ=255)                                                                                                   gAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCa                                                >  1:195871/1‑150 (MQ=255)                                                                                                     aTTGACGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGA‑CTGCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCt                                             <  2:727146/151‑1 (MQ=255)                                                                                                          cgcgcCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGGCCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcagc                                        <  2:773862/151‑1 (MQ=255)                                                                                                          cgcgcCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcagc                                        <  1:39695/151‑1 (MQ=255)                                                                                                            cgcCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcagcag                                      <  2:404610/151‑1 (MQ=255)                                                                                                                     gcgtggcgtgCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGa                              <  1:1122962/150‑1 (MQ=255)                                                                                                                         ggcgtgCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAg                         >  2:573541/1‑151 (MQ=255)                                                                                                                           cgtgCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGa                        <  2:491921/150‑1 (MQ=255)                                                                                                                                    tgtgt‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCAc               <  1:443554/150‑1 (MQ=255)                                                                                                                                    ggtgt‑GCGGGCGGGATGA‑CTTCGGGGAACACGGT‑AAATCCGCTGCTAGGGGCGACGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCTGTTCATACTCTCCCGAACCTACAGCAGCTACCGGACGAAGACGCCTGCAc               <  2:927823/149‑1 (MQ=255)                                                                                                                                      tgt‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGCCCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCgg            <  2:611449/151‑1 (MQ=255)                                                                                                                                      tgt‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCgg            <  1:975226/151‑1 (MQ=255)                                                                                                                                                gggAAGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCTAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGCCGCCGGCACCGGTGGGCGgtt   >  2:1110111/1‑151 (MQ=255)                                                                                                                                                  gATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGtttt  <  2:585103/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                       AACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTT  >  NC_000913/523242‑523535

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