Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	325,394	0	G	A	48.6%	2.2 / 45.3	35	intergenic (‑30/+183)	ykgH/betA	uncharacterized protein YkgH/choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (11/7);  new base A (7/10);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGC  >  NC_000913/325160‑325637                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    tGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAAt                                                                                                                                                                                                                                     >  1:30961/1‑251 (MQ=255)     aTCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTc                                                                                                                                                                                                                                 <  1:123637/251‑1 (MQ=255)       cacaGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCt                                                                                                                                                                                                                                <  2:133349/250‑1 (MQ=255)          aGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAGTCTGTCGCAGATAAGGTTGTATTAATAGTCTGTAt                                                                                                                                                                                                                            >  2:257260/1‑251 (MQ=255)           gCTAACAAGTCTACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATc                                                                                                                                                                                                                           >  1:77420/1‑251 (MQ=255)                           ggCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGtt                                                                                                                                                                                                            <  1:257260/250‑1 (MQ=255)                              gTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAAGCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAat                                                                                                                                                                                                        <  1:284770/251‑1 (MQ=255)                                      gTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCaa                                                                                                                                                                                                <  2:64871/251‑1 (MQ=255)                                          aTTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCaaaaa                                                                                                                                                                                             <  2:238061/250‑1 (MQ=255)                                                       aaaTCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGaa                                                                                                                                                                               <  1:244518/251‑1 (MQ=255)                                                       aaaTCACTTTTAGATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTTATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAATGTTGTATTAATAGTCTTGATCATGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGaa                                                                                                                                                                               <  1:32155/251‑1 (MQ=255)                                                                    atTGTCTTTATTGATATACAGAAACAACAGCACCCTAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGa                                                                                                                                                                   >  1:281891/1‑250 (MQ=255)                                                                        tCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCAt                                                                                                                                                              >  1:82048/1‑251 (MQ=255)                                                                                       aGAAACCACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAgg                                                                                                                                               >  1:250623/1‑250 (MQ=255)                                                                                        gAAACAACAGCACCATAGCCATACATATAGCCATAATAAACAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAgtg                                                                                                                                              >  1:151605/1‑251 (MQ=255)                                                                                          aacaacAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGtt                                                                                                                                            >  1:90928/1‑251 (MQ=255)                                                                                                                 tataGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCctc                                                                                                                     >  2:52939/1‑251 (MQ=255)                                                                                                                  ataGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCctct                                                                                                                    <  1:214711/251‑1 (MQ=255)                                                                                                                              aaaaGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTa                                                                                                         >  1:141433/1‑250 (MQ=255)                                                                                                                                      aaCACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCaa                                                                                                 <  1:126445/250‑1 (MQ=255)                                                                                                                                          cGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGa                                                                                            <  1:52939/251‑1 (MQ=255)                                                                                                                                                  tttCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTAGTAATAGTCTGTATCAGGAATGTTCGGGGTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTc                                                                                     <  1:12139/250‑1 (MQ=255)                                                                                                                                                       tcAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGt                                                                               >  1:190858/1‑251 (MQ=255)                                                                                                                                                                           aTCAGATCGCTATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAt                                                           >  2:155618/1‑251 (MQ=255)                                                                                                                                                                           aTCAGATCGATATCCTGTTTGCTTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATTAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAATGATTCGGTCCCCCTCTATGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAt                                                           >  2:223592/1‑251 (MQ=255)                                                                                                                                                                                      atCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGctct                                                >  2:136809/1‑251 (MQ=255)                                                                                                                                                                                        ccTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCtctcc                                              >  2:176894/1‑250 (MQ=255)                                                                                                                                                                                                 tttGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATg                                                                                                                                                                                 <  1:236315/111‑1 (MQ=255)                                                                                                                                                                                                 tttGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATg                                                                                                                                                                                 >  2:236315/1‑111 (MQ=255)                                                                                                                                                                                                  ttGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGGATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCGCTGGAACGGTAACTCTCAATCTGATCTGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCCCCGGCAttc                                     <  2:100442/250‑3 (MQ=255)                                                                                                                                                                                                   tGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATccc                                    >  1:31481/1‑250 (MQ=255)                                                                                                                                                                                                       cACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCAttt                                <  1:54052/250‑1 (MQ=255)                                                                                                                                                                                                                     ccAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATcc                  <  1:41126/250‑1 (MQ=255)                                                                                                                                                                                                                           aaTAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTg           >  1:255412/1‑251 (MQ=255)                                                                                                                                                                                                                                     gTCGCAAATAAGGTTGGATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAAAACTGGATATGAAGCATGGGAGTTACCTCAGTGGTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCCCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGc  <  1:285782/250‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    TGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGC  >  NC_000913/325160‑325637

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