Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,346,685	0	C	T	42.0%	51.3 / 22.5	50	E214E (GAG→GAA)	fumB	anaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base C (15/14);  new base T (21/0);  total (36/14)
Fisher's exact test for biased strand distribution p-value = 8.91e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAG  >  NC_000913/4346537‑4346833                                                                                                                                                     |                                                                                                                                                         cATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTtctc                                                                                                                                                        <  2:639162/151‑1 (MQ=255)  aTGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTcg                                                                                                                                                       <  1:446022/151‑1 (MQ=255)     ttcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTcgac                                                                                                                                                     <  2:341851/150‑1 (MQ=255)      tcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTcgacg                                                                                                                                                    <  1:625742/150‑1 (MQ=255)             cTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGTGCAGTACCGGGGGGACGCATTTTTTCGACGA‑GGAAggt                                                                                                                                           >  1:983811/1‑149 (MQ=255)              ttCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCCGCCGTAACGCGGGTACGCCTTTTCTCGACGA‑GGAAGgtt                                                                                                                                          >  1:1149962/1‑151 (MQ=255)              ttCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGGACCGCGGGGAAGCATTTTTTCGACGA‑GGAAggt                                                                                                                                           >  2:828213/1‑148 (MQ=255)                ccGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGttttt                                                                                                                                        <  1:313909/151‑1 (MQ=255)                  gTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGGTGCAGGACCGAGGGGACGCCATTTCTCGGCGA‑GGGAGgtttt                                                                                                                                        >  2:204161/1‑149 (MQ=255)                     ggCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGtt                                                                                                                                   <  2:355722/151‑1 (MQ=255)                      gCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGttt                                                                                                                                  <  2:414560/151‑1 (MQ=255)                       cAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTAACGAGGGGGCGCATTTTTTTCACGG‑GGAAGGTTTTTAGttt                                                                                                                                  >  2:635030/1‑150 (MQ=255)                        aGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑c                                                                                                                               <  2:720389/151‑1 (MQ=255)                        aGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAAGGCGCAGTACCGCGGGTACGCCTTTTTTCGCAGG‑GGAAGATTTTCAGTTTgt                                                                                                                                >  2:249638/1‑150 (MQ=255)                            cATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCggg                                                                                                                           <  2:535348/151‑1 (MQ=255)                               cGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCGGACCGAGGGGACGGATTTTCTTGCCGAGGGAGGTTTTTCAGTTGG‑Cggggg                                                                                                                          >  2:817056/1‑149 (MQ=255)                                  aaTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGGCAGTAACGAGGGGACGCATTTTTTCGACGG‑GGGAGGTTTTTCGTTTT‑GCggg                                                                                                                           >  2:136925/1‑145 (MQ=255)                                     aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTcagc‑a‑g                                                                                                                <  2:1121330/151‑1 (MQ=255)                                       tGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGTACCGAGGGGACGCATTTTTTCGACGA‑GGAAGTTTTTTAGTTTG‑CCGGGAGGCAGC‑A‑ggg                                                                                                              >  1:886786/1‑151 (MQ=255)                                         aGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCGTACCGAGGGGACGCATTTTTTCGGCGA‑GGAAGGTTTTTAGTTTT‑TCGGGGGGCAGC‑C‑ggg                                                                                                              >  1:200648/1‑149 (MQ=255)                                          gcgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGA‑GGAAGTTTTTTAGTTTG‑CCGGGGGGCAGC‑A‑GGGCtt                                                                                                           >  1:1019808/1‑151 (MQ=255)                                            gcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGGCCGGAACCAGGGGACGCATTTTTTCGACGG‑GGGAGGTTTTTAGTTTT‑GCGGGGGGCAGC‑C‑GGGCttt                                                                                                          >  2:644868/1‑150 (MQ=255)                                               tGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGt                                                                                                       <  2:387654/150‑1 (MQ=255)                                                 cTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCGGGGGGGCGGACCGGGGGGACGCGATTTTTTTACGA‑GGGAGGTTTTTAGGTTT‑TCGGGGGGCAGC‑AGGGGGTTTGGtt                                                                                                      >  2:707696/1‑150 (MQ=255)                                                      ttGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGGGCAGTACCGGGGGGACGCATTTTTTCGACGA‑GGAAGGTTTTCAGATTG‑GCGGGAGGCAGC‑A‑GGGCTTTGGttt                                                                                                     >  1:866074/1‑145 (MQ=255)                                                      ttGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGA‑GGGAGGTTTTCAGTTTG‑CCGGGGGTCAGC‑A‑GGGCGTTGGTTTCCTgg                                                                                                >  1:681209/1‑150 (MQ=255)                                                        gACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGGCGCAGTACCGGGGGGAAGGATTTTTTTGACGA‑GGGAGGTTTTTAGGTTT‑CCGGGGGTCAGC‑A‑GGGGTTTGGgtttttgg                                                                                                >  2:370282/1‑143 (MQ=255)                                                         aCGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTaga                                                                                            <  1:930216/151‑1 (MQ=255)                                                                        gttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGGTTAGGACCGCGGGGACGCCCTTTTTCCACGA‑GGGAGGTTTTTAGTTTT‑TCGGGAGGCAGC‑A‑GGGCTTTTGGTTTCTGGTTGGGGTAAGTTttgtt                                                                               >  1:357546/1‑149 (MQ=255)                                                                             cGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCGGACCCAGGGGACGCCTTTTTTCGACGG‑GGGAGGTTTTTTGTTTG‑CCGGGGGGGAGC‑A‑GGGCTTTGGTTTTCTGGGAGAGATACGTTTTTTTTGCa                                                                           >  2:540246/1‑148 (MQ=255)                                                                                aGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCGGGGGGCAGACTCCTGGGGAAGCATTTTCTCGTCGG‑GGGAGTTTTTGCGTTTT‑CCGGGAGTCAGA‑C‑GGGGGTTGGGTTTCTCGTATAGATAAGTTTTTTTGGcggcagc                                                                      >  1:1047885/1‑144 (MQ=255)                                                                                         gCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGA‑GGGAGTTTTTTAGTTTG‑CCGGGGGGCAGC‑A‑GGGGTTTGGGTTCCTGGGAGAGATACGGTTTGTTGGGAGAGCCGCCGCCttt                                                             >  1:719149/1‑150 (MQ=255)                                                                                          ccAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCAAACCCAGGGAAAACATTTTTTCGTCCA‑AGAAGTTTTTTAGGTTT‑CCGGGGGGCCAG‑A‑GGGGGTTGGGTTCCGTGTAGAGAGAAATTTTTTTGGCggagggcccgccct                                                              >  1:1091091/1‑138 (MQ=255)                                                                                           cAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCCAAACCCCGGGAAAGCATTTTTTCGACGA‑GGGAGTTTTTCAGCTTG‑CCGGGGGTCAGC‑A‑GGGGTTTTGTTTCCTTGTAGGGATACGTTTTTTTGGCCGAGGCGCCGCCctt                                                             >  1:52740/1‑148 (MQ=255)                                                                                           cAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgc                                                          <  1:986539/151‑1 (MQ=255)                                                                                           cAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGTACCGAGGGGACGCATTTTTTCGACGA‑GGAAGGTTTTTAGTTTG‑CCGGGGGGCAGC‑A‑GGGCTTTTGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgc                                                          >  1:942601/1‑151 (MQ=255)                                                                                                acaAACGCGATATGGTACGGCGGGGCGGGCGCAGGAACCAGGGGACGCCTTTTTTCGTCGA‑GGGAGTTTTTTTGATTGCCCGGGGGTAAG‑‑A‑GGGCGTTTGTTTCCTGGGAGGGATAAGTTTTGTTGGGAGAGCCGCCGCCCTTtgcgac                                                       >  2:839544/1‑149 (MQ=255)                                                                                                acaAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTTCACGA‑GGGAGGTTTTTAGTTTG‑CCGGGGGTCAGC‑A‑GGGCTTTTGGTTTCTGGTAGAGATACGTTTTGTTGGGAGGGGCGCCGCCctttcccacgc                                                     >  2:440683/1‑143 (MQ=255)                                                                                                 caAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgca                                                    <  1:154144/151‑1 (MQ=255)                                                                                                   aaCGCGATATGGTACGGCGGGCAGGCTCCAGTACCGAGGGGACGCAATTTCTTGACGC‑GGCAGGTTTTTAGGTTG‑CCGGGGGTCAGCAA‑GGGCTTTGGTTT‑CTGGCAGAGATACGGTTTTTTTGCAGGGGCGCCGCCCTTcgcaac                                                       >  2:261684/1‑146 (MQ=255)                                                                                                     cgcgATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGA‑GGGAGTTTTTCAGTTTG‑CCGGGGGGCAGC‑A‑GGGGTTTTGTTTTCTGGGAGAGGTACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAgg                                                >  2:372485/1‑151 (MQ=255)                                                                                                           atGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAAttt                                           <  1:1021836/150‑1 (MQ=255)                                                                                                              gTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGGTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTAc                                       >  1:675787/1‑151 (MQ=255)                                                                                                                 cggcggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTc                                     >  2:948269/1‑150 (MQ=255)                                                                                                                                      gAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg               >  2:850337/1‑151 (MQ=255)                                                                                                                                            aCGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGAtt         >  2:874505/1‑151 (MQ=255)                                                                                                                                            aCGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGAtt         >  2:983185/1‑151 (MQ=255)                                                                                                                                              gCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTg       >  1:942617/1‑151 (MQ=255)                                                                                                                                                 ttttCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTgcgc    >  1:1007398/1‑151 (MQ=255)                                                                                                                                                 ttttCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTgcgc    >  1:1024602/1‑151 (MQ=255)                                                                                                                                                   ttCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTgcgc    >  1:42318/1‑149 (MQ=255)                                                                                                                                                    tctcGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCgcag  <  2:172693/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                         CATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAG  >  NC_000913/4346537‑4346833

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