Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	523,335	0	T	G	58.6%	13.9 / 27.5	30	G25G (GGT→GGG)	rhsD	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (8/4);  new base G (17/0);  total (26/4)
Fisher's exact test for biased strand distribution p-value = 2.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.75e-05

GTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAT  >  NC_000913/523189‑523478                                                                                                                                                   |                                                                                                                                                  gTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTaa                                                                                                                                               <  1:719957/151‑1 (MQ=255)     tgtgtgTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGaa                                                                                                                                            >  2:257247/1‑150 (MQ=255)                ttattttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGGCGGGGTAAGAATTGGGGCGccc                                                                                                                                 >  2:477589/1‑150 (MQ=255)                     tttttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGcccccc                                                                                                                              >  1:292285/1‑148 (MQ=255)                      ttttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGAGCGTCAGGGAGATATGACTCAGTATGGCGGACCCATTGTCCAGGGTTCGGGGGGCGTAAGCATTTGCGCGCCCACCgggg                                                                                                                          >  1:220615/1‑149 (MQ=255)                         tttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCCGGGGTAAGAATTGGGGGGCCCCCCggcgg                                                                                                                         >  2:635379/1‑148 (MQ=255)                            gttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGCCGGGGTAAGAATTGGCGCGCCCACCggcgtgg                                                                                                                       >  2:506138/1‑148 (MQ=255)                               gCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCACTTTCCCGGGTTCGGCCTGGGTACGAATTGGCGCGCCCCCCgg                                                                                                                            >  2:46361/1‑140 (MQ=255)                                cTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCAGGTTCGGCCGGGGGAAGAATTTGCGCGCCCACCcgcgtgg                                                                                                                       >  1:435133/1‑144 (MQ=255)                                 tAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCCGGGTTCGGGAGGGGGAAGAATTGGCGCGCCCCCCggc                                                                                                                           >  2:487249/1‑139 (MQ=255)                                  aaCACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTTGGGGGCCCACCGGCGGGGCGTGGTCg                                                                                                               >  2:464813/1‑150 (MQ=255)                                            tatGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTTGGGCGCCCACCGGCGTGGCGGGCTCGGTGGGTCcggg                                                                                                    >  2:124553/1‑151 (MQ=255)                                              tGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTACCACTGTCCCGGGTTCGGCAGGGGGAAGAATTTGCGCGCCCCCCCGCCGGGCGCGCTCGGGGGGGCcgggcg                                                                                                  >  2:536918/1‑151 (MQ=255)                                                   aaTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCCCCGGCGTGGCGTGCTCCGGGTGTCcgggcggg                                                                                                >  1:162000/1‑148 (MQ=255)                                                   aaTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCCGCAGGGGTAAGAATTTGCGCGCCCCCCGGCGTGGCGGGCTCGGGGTGTTccgggggg                                                                                                >  2:595717/1‑148 (MQ=255)                                                    aTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcgggcggggt                                                                                              >  1:601059/1‑147 (MQ=255)                                                                 cACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAAGAATGGGGGGCCCACCCGGGTGGCGCGCTCCGGGGGTTCGGGGGGGGAGA‑ACT‑CGGGCAAacc                                                                             >  2:381736/1‑151 (MQ=255)                                                                    tttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTTGGCAGGGGTAAGAATTGGCGCGCCCCCCCGGGGGGCGTGGTCGGGGGGGCCGGGCGGGGTGACCGT‑CGGGAAAccgg                                                                            >  2:259990/1‑148 (MQ=255)                                                                          tGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGCGGGGTGA‑CTT‑CCGGCAACCCGGGAAAAcc                                                                    >  1:350057/1‑151 (MQ=255)                                                                                    aCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCCGCGTGGCGTGCTCGGGGTGTCCGGGGGGGGTGA‑CTT‑CGGGCAAACCGGGAAAACCGCTGCTgggg                                                          >  1:544657/1‑151 (MQ=255)                                                                                       aGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGCAGGGGGAAGAATTTGGGCGCCCACCGGCGGGGGGGGCTCGGGGGGTCCCGGGGGGGGGA‑ATTCCGGGAAACCCGGAAAACCCCCGCCTgggggg                                                        >  2:725367/1‑150 (MQ=255)                                                                                                                cAGTATGGCGGTCCCCTTGTCCCGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCCGCGTGGCGGGCCCGGCGTGTCCGGGCGGGATGA‑CTT‑CGGGCAACCCCGTAAATCCCCTGCTGGGGGCGAAGGGGCTGCCCGGCGAGACGGAc                               >  2:184942/1‑150 (MQ=255)                                                                                                                cAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGAcc                              <  2:612872/151‑1 (MQ=255)                                                                                                                               aTTGCCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTT‑CGGGCAACCCGGGAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGGCCTTGCGCTtcccggcc               >  2:269672/1‑151 (MQ=255)                                                                                                                                 tGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggccc              <  1:464834/150‑1 (MQ=255)                                                                                                                                    ccAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTg          <  1:725411/151‑1 (MQ=255)                                                                                                                                      aGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTT‑CGGGCAACCCGGTCAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGcc        >  2:374722/1‑151 (MQ=255)                                                                                                                                       gggTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTT‑CGGGCAACCCGGGAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCg       >  1:682303/1‑151 (MQ=255)                                                                                                                                           tCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGCGGGATGA‑CTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTc    >  1:339011/1‑150 (MQ=255)                                                                                                                                            cGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAt  >  1:386339/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAT  >  NC_000913/523189‑523478

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