Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,865,867	0	T	G	42.9%	74.7 / 13.5	49	V256G (GTT→GGT)	ygbN	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (6/22);  new base G (21/0);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 7.74e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGTAGAAGTACTGGAACAGATGCAACTGGCTCCGGCCAGTGAGGAAGGCGCGACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGCTCTAATGATTGCGCT‑GGTGCTGGCATTCTGGTTATTGGCTTTACGTCGCGGCTGGAGCTTACAACACACCAGCGACA  >  NC_000913/2865721‑2866016                                                                                                                                                   |                                                                                                                                                         gTCAGTAGAAGTACTGGAACAGATGCAACTGGCTCCGGCCAGTGAGGAAGGCGCGACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTcc                                                                                                                                                      <  1:747207/151‑1 (MQ=255)   cAGTAGAAGTACTGGAACAGATGCAACTGGCTCCGGCCAGTGAGGAAGGCGCGACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCg                                                                                                                                                     <  2:560565/150‑1 (MQ=255)    aGTAGAAGTACTGGAACAGATGCAACTGGCTCCGGCCAGTGAGGAAGGCGCGACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGGTTcgtc                                                                                                                                                    >  2:824348/1‑147 (MQ=255)       aGAAGTACTGGAACAGATGCAACTGGCTCCGGCCAGTGAGGAAGGCGCGACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGgtt                                                                                                                                                        >  1:597375/1‑143 (MQ=255)                  aaCAGATGCAACTGGCTCCGGCCAGTGAGGAAGGCGCGACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATgccgcc                                                                                                                                     <  1:502147/151‑1 (MQ=255)                  aaCAGATGCAACTGGCTCCGGCCAGTGAGGAAGGCGCGACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGGACGGGTTTCCCAACACCGGTaccgcc                                                                                                                                     >  1:334832/1‑151 (MQ=255)                              tGGCTCCGGCCAGTGAGGAAGGCGCGACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGGACGGGTTTCGCAACACAAGTGCCGCCCTCGCATccccc                                                                                                                         >  2:214797/1‑150 (MQ=255)                               ggCTCCGGCCAGTGAGGAAGGCGCGACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATcccc                                                                                                                          <  1:836454/149‑1 (MQ=255)                               ggCTCCGGCCAGTGAGGAAGGCGCGACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCt                                                                                                                         <  1:976849/150‑1 (MQ=255)                               ggCTCCGGCCAGTGAGGAAGGCGCGACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGGACGGGTTTCGCAAAACCAATTCCGCCCTCGCcccccc                                                                                                                          >  2:923891/1‑149 (MQ=255)                                   ccGGCCAGTGAGGAAGGCGCGACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGGGCGGGTTTCGCCACACCGAAGCCGCCCtcccctcccc                                                                                                                          >  2:731873/1‑145 (MQ=255)                                       ccAGTGAGGAAGGCGCGACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTAc                                                                                                               <  1:524414/151‑1 (MQ=255)                                              ggAAGGCGCGACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAg                                                                                                        <  2:606355/151‑1 (MQ=255)                                                aaGGCGCGACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCt                                                                                                      <  1:437671/151‑1 (MQ=255)                                                    cgcgACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATc                                                                                                 <  1:887925/151‑1 (MQ=255)                                                        aCAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGGACGGGTTCCGCAACACTCATGCCGCCTTCGCATACCCTC‑CTTGGGACGCTACAGCTG‑a                                                                                                   >  2:564039/1‑145 (MQ=255)                                                          aaaaTTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGGTTTCGCCACACCGATGCCGCCCTCCCATCCCCCG‑CTTTGTGCGCTACAGCTG‑ATCGGc                                                                                              >  1:605707/1‑148 (MQ=255)                                                          aaaaTTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGGACGGGTTTCGCCACACCCATGCCGCCCTCGCCTCCCCCC‑CTTTGTACGCTACAGCCG‑ATCgg                                                                                               >  1:724910/1‑147 (MQ=255)                                                            aaTTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGGACGGGGTTCGCCACACTGATGCCGCCCTCCCACCCCCCGCCTTGGT‑CGCCCCAGCTG‑ATCGGCTcccc                                                                                         >  2:474865/1‑151 (MQ=255)                                                             aTTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGCGAGGTTTTCGCCACACCGAAGCCCCCCTTCCCTCCCCCC‑CTTGGTACGCTCCAGCAG‑CTCGGCTccc                                                                                          >  2:595736/1‑147 (MQ=255)                                                               tAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCaa                                                                                       <  2:484612/150‑1 (MQ=255)                                                               tAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGGACGGGTTCCCCCACACTCATGCCGCCCTCTCCTCCCCCC‑CTTTGGACGCCTCAGCTG‑ATCGGCTcccc                                                                                         >  2:596260/1‑148 (MQ=255)                                                                 aGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATgg                                                                                    <  2:306650/151‑1 (MQ=255)                                                                 aGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGGAAGGGGTTCCCCACACCGATGCCCCCCTCGCCTCCCCCC‑CTTTGTACGCCCCACCCG‑CTCGGCTccccca                                                                                       >  2:968942/1‑146 (MQ=255)                                                                       aaaataaaTCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGctc                                                                              <  1:890391/151‑1 (MQ=255)                                                                          ataaatCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGGAAGGGTTTCGCCACACTGATGCCGCCCTCGCATCCCCCC‑CTTTGTACGCTACAGCTG‑ATCGGCTCCCCAATGGGAGctct                                                                             >  1:536045/1‑149 (MQ=255)                                                                              atCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGCTCTAATGAt                                                                       <  2:227403/151‑1 (MQ=255)                                                                                   ccGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGGAAGGGTTTCCCCACCCACTTGCCGCCCTCCCATCCCCCT‑CTTGGTGCGCCACCACCG‑AGCGGCCCCCCCATGGGAGCTCTTATGAtt                                                                      >  1:806697/1‑147 (MQ=255)                                                                                    cGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGCTCTAATGATTGCGCt                                                                 <  1:760152/151‑1 (MQ=255)                                                                                        cGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGGGAGGTTTTCCCCACACTTCTGCCGCCCTCGCATCCCCCG‑CTTGGGACGCTAAACCTG‑GTCGGGCCCCCCAAGGGAGCTCTTATTATTTGGCT‑gggg                                                            >  1:131263/1‑149 (MQ=255)                                                                                                 ggTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGGACGGGTTTCGCCACACCGATGCCGCCCTCCCCTCCCCCG‑CCTGGTACGCCACCGCCG‑ATCGGCCCACCCAACGGAGGTCTTATGATTTTGCT‑GGGGGTgggctt                                                    >  1:327321/1‑146 (MQ=255)                                                                                                 ggTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGGACGGGTTCCGCAACACCGATGCCGCCCTCGCCTCCCCCG‑CCTTGGGCGCCACAGCCG‑GTCCGGTCCCCCATGGTAGCTCTTATGATTTCGCT‑GGGGGTGGcgt                                                     >  2:836143/1‑147 (MQ=255)                                                                                                   tCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGGACGGGTTTCGCCACACTTATGCCCCCCTCGCATCCCCCG‑CTTGGTACGCTACAGCTG‑ATCGGCCCCCCCATGGGAGCCCTAATTATTGTGCT‑GGGGGTGGCGTTCt                                                  >  2:79034/1‑150 (MQ=255)                                                                                                    cACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGCTCTAATGATTGCGCT‑GGTGCTGGCATTCTg                                                 <  2:93158/150‑1 (MQ=255)                                                                                                         cGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGCTCTAATGATTGCGCT‑GGTGCTGGCATTCTGGTTa                                             <  2:94726/149‑1 (MQ=255)                                                                                                         cGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGCTCTAATGATTGCGCT‑GGTGCTGGCATTCTGGTTAtt                                           <  2:779400/151‑1 (MQ=255)                                                                                                          gCTAATTGTTATTCCTATCGCGATTATCATGGCGGGGATGGGTTTCGCAAAACAGATGCCGCCCTCGCCTCCCCCT‑CTTTGTAAGCTACAACTG‑ATCGGCTCACCCATAGTAGATCTAATGATTTTGGT‑GGTGTTGGGCTTCTTGTtttt                                           >  1:138917/1‑150 (MQ=255)                                                                                                           cTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGGTTTCGCAACAATTATGCCCCCCTCGCCTCCCCCG‑CTTGGGACGCTACAGCTG‑ATCCGCTCACCCACGGGAGCTCTAATTATTTCGCT‑GGGGCTGGGGTTCTgggtttt                                           >  2:511478/1‑144 (MQ=255)                                                                                                                   tATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGCTCTAATGATTGCGCT‑GGTGCTGGCATTCTGGTTATTGGCTTTACg                                  <  1:981608/150‑1 (MQ=255)                                                                                                                          aTCGCGATTATCATGGCGGGGTCGGGTTTCGCCAAACAGAGGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTGCAGCTG‑AAGGGCTCACCAAAAATAGCTCTAATAATTATGCT‑GGTGCTGGGATTCTGGTTATTGGCTTTTTGTCTCgg                            >  2:14442/1‑149 (MQ=255)                                                                                                                           tCGCGATTATCATGGCGGGTACGGGTTTCGCCACACCGATGCCCCCTTCCCATCCCCCG‑CTTTGTTCGCTACAGCCG‑GTTGGCTCCCCCATCGGAGCTCTTATTATTTTGCT‑GGTGGTGGGCTTCttgttttttgtttttctt                                 >  2:138695/1‑139 (MQ=255)                                                                                                                            cgcgATTATCATGGCGGGTAACGGTTCCGCAACCCCACCGCCCCCCTCTCATCCCCCT‑CTTGGTACGGTTAAGCTGAACCGGCTACCCAATGG‑AGCTCTCATTATTGCGCTGGGTGTCGGGTTCCTGGTT‑TTTGCTTTTCGTcgc                              >  2:841900/1‑145 (MQ=255)                                                                                                                             gcgATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGCTCTAATGATTGCGCT‑GGTGCTGGCATTCTGGTTATTGGCTTTACGTCGCGGCTGGa                       <  2:845550/151‑1 (MQ=255)                                                                                                                                  tATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGCTCTAATGATTGCGCT‑GGTGCTGGCATTCTGGTTATTGGCTTTACGTCGCGGCTGGAGCtt                   <  2:815043/150‑1 (MQ=255)                                                                                                                                        ggCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGCTCTAATGATTGCGCT‑GGTGCTGGCATTCTGGTTATTGGCTTTACGTCGCGGCTGGAGCTTacaactc            <  1:464908/151‑3 (MQ=255)                                                                                                                                          cGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGCTCTAATGATTGCGCT‑GGTGCTGGCATTCTGGTTATTGGCTTTACGTCGCGGCTGGAGCTTACAacacac          >  1:866917/1‑151 (MQ=255)                                                                                                                                             gTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGCTCTAATGATTGCGCT‑GGTGCTGGCATTCTGGTTATTGGCTTTACGTCGCGGCTGGAGCTTACAACACAcc         >  1:470328/1‑149 (MQ=255)                                                                                                                                               aCGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGCTCTAATGATTGCGCT‑GGTGCTGGCATTCTGGTTATTGGCTTTACGTCGCGGCTGGAGCTTACAACACACCAGCg     >  1:435072/1‑151 (MQ=255)                                                                                                                                                cGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGCTCTAATGATTGCGCT‑GGTGCTGGCATTCTGGTTATTGGCTTTACGTCGCGGCTGGAGCTTACAACACACCAGCg     >  1:385433/1‑150 (MQ=255)                                                                                                                                                  gTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGCTCTAATGATTGCGCT‑GGTGCTGGCATTCTGGTTATTGGCTTTACGTCGCGGCTGGAGCTTACAACACACCAGCGACa  <  1:731954/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                         GTCAGTAGAAGTACTGGAACAGATGCAACTGGCTCCGGCCAGTGAGGAAGGCGCGACAAAATTAAGCGATAAAATAAATCCACCGGGCGTCGCGCTGGTCACCTCGCTAATTGTTATTCCTATCGCGATTATCATGGCGGGTACGGTTTCCGCAACACTGATGCCGCCTTCGCATCCCCTG‑CTTGGTACGCTACAGCTG‑ATCGGCTCACCAATGGTAGCTCTAATGATTGCGCT‑GGTGCTGGCATTCTGGTTATTGGCTTTACGTCGCGGCTGGAGCTTACAACACACCAGCGACA  >  NC_000913/2865721‑2866016

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