Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,381,728	0	A	C	46.2%	70.8 / 10.1	52	G197G (GGT→GGG)	frdA	anaerobic fumarate reductase catalytic and NAD/flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (22/6);  new base C (0/24);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 1.60e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTG  >  NC_000913/4381585‑4381874                                                                                                                                                 |                                                                                                                                                   gTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCAc                                                                                                                                             <  1:570811/151‑1 (MQ=255) gTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCAc                                                                                                                                             >  1:298344/1‑151 (MQ=255)   catcaGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCg                                                                                                                                           >  2:999843/1‑151 (MQ=255)    atcaGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGc                                                                                                                                          <  2:368250/151‑1 (MQ=255)     tcaGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGcc                                                                                                                                         >  2:199930/1‑151 (MQ=255)     tcaGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGcc                                                                                                                                         >  1:447148/1‑151 (MQ=255)       aGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAg                                                                                                                                       <  2:192861/151‑1 (MQ=255)                gAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGAAGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATa                                                                                                                               >  1:109956/1‑150 (MQ=255)                     tGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGAcc                                                                                                                         >  2:1051309/1‑151 (MQ=255)                       gcgggcgggttggttGATCCTGCACGAATCCCATGTCACGGGGCGGACCGCGGTGCTTGAGCGCCCTCCCCACGCCGCCCGCGGCGACGGGGCCGCCTTTGG‑TTTTTTAACAAAAAACGGCCCCCGCCCCGCCAGTAGCCATAACGACcgc                                                                                                                       <  1:144327/144‑1 (MQ=255)                           aCCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGtt                                                                                                                    >  1:161660/1‑150 (MQ=255)                            cgggTTGGGTGATATTGAACGAATTCCATGTCCCGCGGGGGAACGCCGTGGTTTGGCCCCATCCCCATCCCGTCACCGGTAACGATCCCGCCGTTGG‑TGTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGtt                                                                                                                    <  2:630933/147‑1 (MQ=255)                              ggTGGGGTGATACTGAACGAATTCCATGTCCCGCAGCGGACCGCCGTGGTTTAGCGCCATCCCCACCCCGTCCCCGATAACGATGCCGCCGTTGG‑TGTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAg                                                                                                                  <  2:33015/149‑1 (MQ=255)                                  gggTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCAACGCCAGTAGCCATAACGACCGCGTTAGCACgg                                                                                                             <  2:961243/150‑1 (MQ=255)                                   ggTGATCCTAAACGAATTCCATGTCCCGCAGCGGAACGCGGGGGCTTAGCCCCATCCCCACCCCGTCACCGGTAACGAGGCCGCCGTTGG‑TGTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGAt                                                                                                           <  2:620283/151‑1 (MQ=255)                                        tACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGc                                                                                                       >  1:909646/1‑150 (MQ=255)                                          ctaaccaaaTCCCATGCCCCGCAGGGGACCGCGGTGGTTTGCCCCCATCCCCACCCCTCCCCCGGAAACGATGCCGCCGTGGG‑TTTTTAAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCa                                                                                                      <  1:614847/142‑1 (MQ=255)                                             aaCGAATTCCATGTCCCGCAGGGGACCGCCGTGGTTTAGCGCCATCCCCATACCGTCCCCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCa                                                                                                   <  2:1019124/149‑1 (MQ=255)                                             aaCGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGc                                                                                                 >  2:485267/1‑151 (MQ=255)                                              aCGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCg                                                                                                >  1:925530/1‑151 (MQ=255)                                                         gTCACGCGGCGGAGCGCCGTGGTTTAGCGCAACCCCCATACCTCCACGGGAACCGAGGCCGCCTTTGG‑GGTTTTAACGAAAAACGCCCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTg                                                                                              <  2:778976/142‑1 (MQ=255)                                                                   ggACCGCGGGGGTTTAGCGCCACCCCCATCCCGCCCCGGGTAACGAGCCCGCCGTTGG‑TGTTGAAACAAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCAtt                                                                            <  2:141749/150‑1 (MQ=255)                                                                   ggAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTg                                                                           >  2:1040030/1‑151 (MQ=255)                                                                       cGCCGTGGCTTGGCGCCACCCCCATCCCGTCCCCGGTAACGATGCCGCCGTGGG‑TGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTg                                                                           <  2:167921/147‑1 (MQ=255)                                                                        gCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGAGCCTTCCATCATGTTCATTGCTAcc                                                                      >  1:502564/1‑151 (MQ=255)                                                                         cgggggCTTAGCGCCATCCCCATCCCGTCCCGGGTACCGATGCCCCCTTGGG‑TTTTGTAAGGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTa                                                                        <  2:634470/144‑1 (MQ=255)                                                                         ccGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTAcc                                                                      >  2:1168508/1‑150 (MQ=255)                                                                           ggggttgggcgcCATCCCCATCCCGCCCCGGGAAGCGATCCCGCCGTTGG‑TGTTGAAACAAAAACCGCCCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTa                                                                        <  2:570792/138‑1 (MQ=255)                                                                           gTGGTTTAGCCCCATCCCCATCCCGTCCCCGTTAACGATGCCGCCGTTGG‑TGTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAg                                                                    <  2:466853/150‑1 (MQ=255)                                                                           gTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAgg                                                                   >  2:414279/1‑151 (MQ=255)                                                                             gggttAGCGCCATCCCCAACCCCTCCCGGGAAACGCTGCCGGCGGTGG‑TGTTGGATGGATAAACGCGCCCGCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAgg                                                                   <  1:650459/146‑1 (MQ=255)                                                                                        atacccataccGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGAc                                                      >  1:789706/1‑151 (MQ=255)                                                                                          ccccttcccGTCCCCGGTAACGA‑GCCCGCGTTGGTTTTTGAAACAAAAACCCCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACca                                                    <  2:792623/144‑1 (MQ=255)                                                                                            caatcccGTCCCCGGTAACGAGGCCGCCGTTGG‑TTTTGAAACAAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACcatc                                                  <  2:882385/146‑1 (MQ=255)                                                                                                 cctcccccGGGAACGAGGCCGCGGTGGG‑TGTTTTAACAATAACCGCACCCCCCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCatcatc                                               <  2:875159/143‑1 (MQ=255)                                                                                                       ccGGTACCGCTCCCGCCGTTGG‑TTTTGTAACAAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGaa                                       <  2:70523/151‑1 (MQ=255)                                                                                                         ggTAGCGATGCCGCCGTTGG‑TGTTTAAACAATAACCGCCCCCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAat                                      <  1:810921/150‑1 (MQ=255)                                                                                                               gagcccgccgTTGG‑TGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAg                                <  2:588683/146‑1 (MQ=255)                                                                                                                   ccgccgTTGG‑TTTTGTAACGATAACCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCAc                             <  2:940277/149‑1 (MQ=255)                                                                                                                   ccgccgTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACg                            >  1:947435/1‑150 (MQ=255)                                                                                                                     gccgTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGaa                          >  2:246711/1‑150 (MQ=255)                                                                                                                      cctttGG‑TTTTGAAACAAAAAGCGGGACCATCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAAt                        <  2:685491/148‑1 (MQ=255)                                                                                                                      ccgTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAAt                        >  1:489005/1‑151 (MQ=255)                                                                                                                        gTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGt                      >  1:1181278/1‑151 (MQ=255)                                                                                                                         ttGG‑TTTTGTAACAAAAACCGCGCCCCCCTCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAAt                        <  1:883593/148‑1 (MQ=255)                                                                                                                         ttGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGtt                     >  1:1243098/1‑151 (MQ=255)                                                                                                                            ggttttaaaacaaaaaaCGCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGt                      <  1:1181483/134‑1 (MQ=255)                                                                                                                               gttatAACAATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGt                  <  2:656155/149‑1 (MQ=255)                                                                                                                                     aCGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCaaaa             <  2:634415/148‑1 (MQ=255)                                                                                                                                         tAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTg        >  1:943405/1‑149 (MQ=255)                                                                                                                                            aCGCGCCCCGGCCCGCCAGTAGCCATAACGACCGCGGTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATc    <  2:934922/150‑1 (MQ=255)                                                                                                                                             cgcgACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTg  >  1:1149512/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                   GTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGG‑TGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTG  >  NC_000913/4381585‑4381874

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