Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,299,250	0	A	G	41.4%	18.3 / 25.3	31	G39G (GGT→GGC)	fdhF	formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base A (6/11);  new base G (6/6);  total (12/19)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTACCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAA  >  NC_000913/4299002‑4299483                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          ttttGCATTACATAGTGGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTAc                                                                                                                                                                                                                                          <  1:143296/250‑1 (MQ=255)           cATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTACCCTggttggtt                                                                                                                                                                                                                               <  1:326246/251‑1 (MQ=255)             tAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTGCCCTGGTTGGtttt                                                                                                                                                                                                                             >  1:330422/1‑251 (MQ=255)                       cGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTGCCCTGGTTGGTTTTCCCCTGCGCt                                                                                                                                                                                                                   >  1:163197/1‑251 (MQ=255)                                               cGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTACCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCg                                                                                                                                                                                           >  2:122622/1‑251 (MQ=255)                                                    gtcTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTGCCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATc                                                                                                                                                                                      <  1:240694/251‑1 (MQ=255)                                                        ggATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTACCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATcg                                                                                                                                                                                     <  1:320951/248‑1 (MQ=255)                                                          aTGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTACCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCgacgac                                                                                                                                                                                >  2:258171/1‑251 (MQ=255)                                                             gCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTGCCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAg                                                                                                                                                                             <  1:70621/251‑1 (MQ=255)                                                               aTCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCACCGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTACCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGt                                                                                                                                                                           <  1:192944/251‑1 (MQ=255)                                                                          actgcttttgggtGGCGCTAAGGTTCTCGGCAAGGTAATGGAGTGCCTCAGCCCAGGAAACATGGTCGAGTTTGCGGCCACGCTGGCGACGGAGCATGGGGGTTTTCAGTCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCGTCAGACACTGGGTTCCCTGTTTGGGTTTCCCCTGCGCTGCCTCCGCCCGTACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCa                                                                                                                                                                 <  2:376088/238‑1 (MQ=255)                                                                             tctcTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTACCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAAcc                                                                                                                                                              <  1:281599/250‑1 (MQ=255)                                                                                                  gcTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTACCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAAcc                                                                                                                                        <  1:309332/251‑1 (MQ=255)                                                                                                                 tCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTACCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTc                                                                                                                         >  1:164104/1‑251 (MQ=255)                                                                                                                                           tCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTGCCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTaatcaaat                                                                                               >  1:250356/1‑251 (MQ=255)                                                                                                                                               gTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTGCCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACg                                                                                           >  2:236407/1‑251 (MQ=255)                                                                                                                                                       ccACGCTGGCGACGGATCATGTGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTGCCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATAcgcg                                                                                   <  2:38583/251‑1 (MQ=255)                                                                                                                                                                   cGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTGCCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACa                                                                       <  1:361209/251‑1 (MQ=255)                                                                                                                                                                        cATGGGGGTTTTCAGGAGGGGTGTCAGGATCTGGGTATCGTTAAGGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTTCCCTGGTTGGTTTTCCCCGGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAAGAGGGGCAAACCGTGACGACTTTTTTCAGCGGTCTCGCTCCAGTTAATCGAATCACGCATACGCGCTCTCGACTACAGTaa                                                                   <  2:130154/250‑2 (MQ=255)                                                                                                                                                                                           ggggTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTACCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACAttt                                               <  1:42100/251‑1 (MQ=255)                                                                                                                                                                                                       tGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTACCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTc                                    >  1:208281/1‑250 (MQ=255)                                                                                                                                                                                                        gggTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTGCCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTgggg                                      <  1:236407/247‑1 (MQ=255)                                                                                                                                                                                                            aTCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTACCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTccc                              >  2:95321/1‑251 (MQ=255)                                                                                                                                                                                                              cGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTACCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGGATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCt                             <  1:20514/250‑1 (MQ=255)                                                                                                                                                                                                                        gTCCCAGCCATAATAACCCTTCAGACACCGGGTACCATGGTTGGTTTTCCCCTGCGCTGCCTCAGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACggg                  >  2:124903/1‑251 (MQ=255)                                                                                                                                                                                                                            cAGCCATAATAACCCTTCAGACACAGGGTGCCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACt              <  2:250356/251‑1 (MQ=255)                                                                                                                                                                                                                            cAGCCATAATAACCCTTCAGACACAGGGTACCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACt              <  1:362562/251‑1 (MQ=255)                                                                                                                                                                                                                                cATAATAACCCTTCAGACACAGGGTACCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACtattt          <  1:95321/251‑1 (MQ=255)                                                                                                                                                                                                                                  taataaCCCTTCAGACACAGGGTGCCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACtatttat        >  1:3963/1‑251 (MQ=255)                                                                                                                                                                                                                                       aCCCTTCAGACACAGGGTGCCCTGGATGGTTTTCACCTGCGCTGCCTCCGCCCGGACGCTTTTGCCGTTATCGACGACAAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCTTCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCaa    >  2:358627/1‑250 (MQ=255)                                                                                                                                                                                                                                        cccTTCAGACACAGGGTACCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCaaaa  <  1:258171/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          TTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATAACCCTTCAGACACAGGGTACCCTGGTTGGTTTTCCCCTGCGCTGCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAA  >  NC_000913/4299002‑4299483

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