Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,411,689	0	G	T	79.4%	58.0 / 11.6	35	E13* (GAA→TAA)	yhdJ	DNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base T (12/15);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCG  >  NC_000913/3411461‑3411935                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       aCATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACGCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACgg                                                                                                                                                                                                                                  <  1:85534/251‑1 (MQ=255)                 cacCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGaa                                                                                                                                                                                                                  <  1:338026/251‑1 (MQ=255)                     cGTGGTAACCAGACCAGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACtt                                                                                                                                                                                                               >  1:377671/1‑250 (MQ=255)                         gTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTaaaa                                                                                                                                                                                                           >  1:139522/1‑250 (MQ=255)                                             gtgataaTGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATGCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACGCGGCATGGGGAGGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGAGGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTGAGGGTGGTGCGTTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCg                                                                                                                                                                                      <  2:463501/250‑1 (MQ=255)                                                              ccGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGAc                                                                                                                                                                     >  2:59998/1‑251 (MQ=255)                                                               cGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAGGCGCCTTTTTTATAGGTGTCACAAAGTGAGTGACCATGATAACAGGATGTGACCCGACCCGGTTTGGTAATGAAGCTAAGCCCATTATTCACGGTGAGGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGAcc                                                                                                                                                                    <  1:327252/251‑1 (MQ=255)                                                                   gTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGAGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACccacc                                                                                                                                                                >  2:310372/1‑251 (MQ=255)                                                                     aCGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACccacca                                                                                                                                                               <  1:245627/250‑1 (MQ=255)                                                                      cGCTGCGTAAAAAATTGAAAAAATACGGCAGTAACTAATTCAGGTTAGCTAAAGGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCata                                                                                                                                                             <  1:425053/251‑1 (MQ=255)                                                                              aaaaaaTTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCg                                                                                                                                                     >  2:11267/1‑251 (MQ=255)                                                                                              aCGGCATGAACTAATGCAGCTTAGCTAAATGCTTAATTAAAAAGGCGCTACTCGGCATGGGGAAGCGACTTTTTTATAGGTGTCACAAAGGGAGTGACCATGGGAACAGGAGGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATACACGGTGATGCTCTTGCCAAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATg                                                                                                                                      <  1:252304/250‑1 (MQ=255)                                                                                                ggCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTc                                                                                                                                   >  2:121323/1‑251 (MQ=255)                                                                                                 gCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCTTTTTTTATAGGTGTCACAAAGGGAGTGACAATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCt                                                                                                                                  <  1:418618/251‑1 (MQ=255)                                                                                                   aTGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTg                                                                                                                                 >  1:121453/1‑250 (MQ=255)                                                                                                   aTGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTg                                                                                                                                 >  1:108965/1‑250 (MQ=255)                                                                                                     gAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACACGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATc                                                                                                                              >  2:366325/1‑251 (MQ=255)                                                                                                                ggTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCCCGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGaa                                                                                                                   <  1:271635/251‑1 (MQ=255)                                                                                                                 gTTAGCTAAATGCTTGATTAAAAAGGCGATACTCGGCAGGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACGGACCCGGTTTGGTAAAAAAGCTAAAACCATTATTAACGGGGATGCGCTTGCCGACCTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGCTCGAAGCCTGGaaa                                                                                                                  <  2:317248/251‑1 (MQ=255)                                                                                                                    aGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGCGTGACAATGAGAACAGGATGTGAACCGACACGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCACGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAaaga                                                                                                                >  2:381877/1‑250 (MQ=255)                                                                                                                     gCTAAATGCTTAAATAAAAAGGCGCTACTCGGCTTGGGGAGGCGCCTTTTTTATAGGTTTCACAAAGGGAGTGACCAGGAGAACAGTATGTGATCCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGACCTTAAAAAGACCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAaagaag                                                                                                              <  1:366325/251‑1 (MQ=255)                                                                                                                        aaaggtttGATTAAAAAGTCGCTCCTCGGAAGGGGGAAGCGCTTTTTTTATAGGTGTCAAAAAGGAAGTGACCAGGAGAACAGGCTGTGAACCGACCCGGTTTGGTAATTAGGCAAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATc                                                                                                           <  1:14240/245‑1 (MQ=255)                                                                                                                                               cgagcgcagggggAAGCTCTTTTTTTATAGGGTTCACAAAGGGAGTGACGATGATAACAGGATGTGAACCGACCCGTTTTGGTAATTAAGCTAAGACCATTATTCTCGGTGATGCGCTGGCTGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAAGTTTGATGGTCTGGTCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGGTTGaa                                                                                    <  2:468028/242‑1 (MQ=255)                                                                                                                                                     cATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGCACCGACCCGGTTTGGCAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGaa                                                                                    >  1:468028/1‑245 (MQ=255)                                                                                                                                                            aaGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGGTCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGt                                                                       <  1:366975/251‑1 (MQ=255)                                                                                                                                                            aaGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGACGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGt                                                                       >  2:123890/1‑251 (MQ=255)                                                                                                                                                                                tCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGaaaa                                                     >  1:401406/1‑249 (MQ=255)                                                                                                                                                                                   aaaaGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAg                                                <  1:470904/250‑1 (MQ=255)                                                                                                                                                                                                              gtgAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTAc                      <  1:68143/250‑1 (MQ=255)                                                                                                                                                                                                                gcaCCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGCACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGAAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACAGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACgg                    >  2:279643/3‑250 (MQ=255)                                                                                                                                                                                                                 aaCCGACCCGGTTTGGTAAGTAAGCAAAGACAATTATTCACGGTGATGCGCTTGCCAAACTTAAAAAGACCACCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGGGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGaa                  <  1:48134/251‑1 (MQ=255)                                                                                                                                                                                                                   ccGACCCGGTTTGTTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGaaaa                <  1:152944/251‑1 (MQ=255)                                                                                                                                                                                                                    cGACCCGGTGTGGTAATTAAGTGAAGACAATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTTTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTTTTTGAAGTGATTGCAGAGTGCGACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAAc               <  1:35644/251‑1 (MQ=255)                                                                                                                                                                                                                         cGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGcc          >  1:442375/1‑251 (MQ=255)                                                                                                                                                                                                                                  aaTGAAGCTAAGCCCATTATTCACGGTTATGCTTTTGCCGAACTTAAAACGACCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATACCATCGGTCAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCg  <  1:109458/250‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       ACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCG  >  NC_000913/3411461‑3411935

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