Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,348,633	0	T	G	44.9%	72.9 / 11.9	49	G40G (GGT→GGG)	yfaE	ferredoxin involved with ribonucleotide reductase diferric‑tyrosyl radical (Y*) cofactor maintenance
Reads supporting (aligned to +/- strand):  ref base T (5/22);  new base G (22/0);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 1.64e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCGTTAG‑CCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAG  >  NC_000913/2348485‑2348779                                                                                                                                                     |                                                                                                                                                       acCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGTTa                                                                                                                                                      <  2:803419/151‑1 (MQ=255)   cgacgaTTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGGAGGGTACt                                                                                                                                                    >  1:359245/1‑151 (MQ=255)     acgaTTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGTTACTGc                                                                                                                                                  <  1:260519/151‑1 (MQ=255)         tttGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTTGGGACCCGTGTCGCGAAGGGTaattcg                                                                                                                                                 >  2:149187/1‑143 (MQ=255)           tGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGGGTACCAGGGGCGCGAAGGTTACTGCGGCGcc                                                                                                                                            >  1:161338/1‑151 (MQ=255)           tGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGGTACTGCGGcccc                                                                                                                                            >  2:366759/1‑151 (MQ=255)                aaCTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCg                                                                                                                                       <  1:801959/151‑1 (MQ=255)                   ttCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTTGCGAAGGGTACTGCGGGCCCTGTCGcac                                                                                                                                    >  1:800146/1‑151 (MQ=255)                     ccAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGca                                                                                                                                     <  1:710821/148‑1 (MQ=255)                     ccAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCacac                                                                                                                                  <  1:342102/151‑1 (MQ=255)                       aGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGGAGGGTACTGCGGCTCCTGTCGcac                                                                                                                                    >  1:15681/1‑147 (MQ=255)                        gCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGGTACTGCGGCTCCTGTCGCACACGcc                                                                                                                               >  1:164825/1‑151 (MQ=255)                          tctGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAATGTCGCGGAGGGTACTGCGGCTCCTGTCGCACACGCCt                                                                                                                              >  1:860315/1‑150 (MQ=255)                           ctGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGCTGGGTACCAGTGGCGCGGAGGGTACTGCGGCTCCCGTCGCCCACAcc                                                                                                                               >  2:377601/1‑148 (MQ=255)                             gATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGGTGAGTACCAGGGACGGGAAGGGTACTGCGGCACCCTGAGCACACGc                                                                                                                                >  1:708678/1‑145 (MQ=255)                               tGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGATTGTTAACAGCGTTGCGAAGGGTATTGCGGGTccctgcccccccac                                                                                                                                >  1:712243/1‑131 (MQ=255)                                   ccGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTTCCAGTGTCGCGAAGGGTTATGCGGCTCCTGTCGCCCACGCCTGGTTGCAgg                                                                                                                     >  1:470947/1‑150 (MQ=255)                                    cgcgTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTc                                                                                                                   <  1:130259/151‑1 (MQ=255)                                          aCCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTTAA‑Gtt                                                                                                             <  1:417864/150‑1 (MQ=255)                                                   aTCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTTAGTAACAGTGTCGCGGAGGGTACTGCGGCTCCCGTCGCACACGCCCGGGTGCAGGGCAA‑GGTGAATGGGtt                                                                                                    >  1:509390/1‑150 (MQ=255)                                                    tCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGGGTACCAGTGTCGCGAAGGGTACTGCGGCTCCTTTCGCCCACGCCCGGTTGCAGGGCAA‑GGTGACTGGGtt                                                                                                    >  2:294112/1‑149 (MQ=255)                                                                      tgtgCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGGGTACCACTGTCCCGAAGGTTACTGCGGGGCCTTTTGCCCACGCCCCGTTGCGGGGCAA‑ACTTACTGGGTTACCCG‑AACCGTT‑G‑CCttttt                                                                             >  2:474492/1‑149 (MQ=255)                                                                            aGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGGTAATGCGGCCCCTGTCTCCCACGCCTGGGTGCAGGTCAAGGTTGACT‑GATT‑GCCG‑GACCCGTAG‑CCCTTTTTTAGcc                                                                       >  1:47855/1‑151 (MQ=255)                                                                            aGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGATTAGTAACAAGGGGGGGGAGGGTACCGGGGGTCCCTGCGCCCACAGCCGGGTGCTGGGCCA‑GTTGACTGGATT‑TACG‑GAGCGTTAG‑CCCTTATTAAggc                                                                       >  2:457576/1‑149 (MQ=255)                                                                                        cttcccttcTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGGACCAGTGGCGCGAAGGTTTCTTCGGCTCCCGTCGCCCACGCCCGGTTGCAGGGCAA‑GATGACTGGATT‑GCCGAAACCGTTAG‑CC‑TTATTTAGCCGGGGGGAAtttt                                                           >  1:479227/1‑151 (MQ=255)                                                                                         ttcccttcTGGCGGCGCTGGAATCCCACAATGTGGCGGTTTAGTACCAGTGTCGCGAAGGGTACTGCGGCCCCCGTCGCACACGCCCGGGTGCAGGTCAA‑GTTGACTGGGTT‑TCCG‑GACCGGTAG‑CCTTTATTTAGCCGGGGGGAAttttt                                                          >  2:445398/1‑150 (MQ=255)                                                                                          tcccttcTGGCGGCGCTGGAATCCCACAATGTGGCGGTTTGGTACCAGTGTCGCGGAGGGTACTTCGGCTCCCTTCGCACACGCCTGGGTGCAGGGCAA‑GTTTACTGGATT‑GCCG‑GACCGTTTGCCCTTTTTTCGCCGGGGGGAAAtttt                                                           >  2:425283/1‑150 (MQ=255)                                                                                                 tGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGGCGCGGAGGGTACTTCGGCTCCTGTTGCACACGCCCGGTTGCAGGTCAA‑ATTTACTGGATT‑TCCG‑AACCGTTTG‑CCTTTATTTAGCCGGGGGGAATTTTTTCCtgtt                                                   >  1:482884/1‑150 (MQ=255)                                                                                                  ggggccgtgggAACCCCCCAAGGGGCCGTTTGAGTCCCATGGCCGCGAGGTTTCCTGGGGTCCCTGTCGCCCCCGCCGGGTTGAAGGTAAA‑GTTGCTTGGATT‑GCCG‑AACGGTT‑G‑CCTTTATTCAGCCGGGGGAAATTTTGCCCtgttgt                                                 <  2:111846/141‑1 (MQ=255)                                                                                                  ggcggcGCTGGAATCCCACAATGTGGCGGTTGAGTATCAGTTTCGCGGAGGGTACTGCGGCCCCTGTCGGACACGGCCGGTTGCAGGGCAA‑GTTGGCTGGCTT‑GCCG‑AACCCCTTG‑CCTTTTTTTTTCCGGGGGGAAATTTGCCCcgttgt                                                 >  2:229122/1‑151 (MQ=255)                                                                                                  ggcggcGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCGTTAG‑CCTTTATTCAGCCGGGGGAAATTTTGCCCtgtt                                                   <  1:337628/149‑1 (MQ=255)                                                                                                  ggcggcGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCGTTAG‑CCTTTATTCAGCCGGGGGAAATTTTGCCCtgtt                                                   <  2:345424/149‑1 (MQ=255)                                                                                                    cggcgCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTTGCGGAGGGTAATGCGGCTCCTGTCGCCCACCCCTTGTTGCAGGGCAA‑GTTGTCTGGGTT‑GCCG‑GACCGGTTG‑ACTTTTTTCAGCCGGGGGAAATTTTTCCCTgttgtt                                                >  1:464298/1‑150 (MQ=255)                                                                                                       cgcTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCGTTAG‑CCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCg                                            <  1:294138/151‑1 (MQ=255)                                                                                                        gcTGGAATCCCACAATGTGGCGGTTTAGTACCAGTGTCGCGGAGGGTACTGCGGCTCCTGTCGCACACGCCTGGGTGCAGGGCAA‑ATTGACTGGGTT‑GCCG‑AAACGTTAG‑CCCTTTTTTAGCCGGGGGGAATTTTTCCCCGTTGTTGcccg                                           >  2:183640/1‑149 (MQ=255)                                                                                                         cTGGAATCCCACAATGTGGCGGTTGAGGACCAGTGTCGCGCAGGTTAAATCGGCCCCTTGCGCAAGCGCCAGGTTGCAGGGCAA‑AATGACTGGATT‑GCCG‑AAACCGTAA‑CCTTTATTCAGGCGGGGGGAATTTTTCCCCGTTGTTTCCggg                                          >  1:676065/1‑151 (MQ=255)                                                                                                          tGGAATCCCACAATGTGGCGGTTGAGGACCAGTGTCGCGGAGGGTACTGCGGCCCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGGTT‑GCCG‑CACCGTTAG‑CCTTTTTTTAACCGGGGGGAATTTTGCCCTGTTGTTGCCggg                                          >  2:704946/1‑150 (MQ=255)                                                                                                           ggAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCGTTAG‑CCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGc                                         <  1:456764/150‑1 (MQ=255)                                                                                                             aaTCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCGTTAG‑CCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCaaa                                      <  2:74001/151‑1 (MQ=255)                                                                                                                 ccACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCGTTAG‑CCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGc                                  <  1:445445/151‑1 (MQ=255)                                                                                                                    caATGTGGCGGTTGAGTACCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCGTTAG‑CCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGat                               <  2:320144/151‑1 (MQ=255)                                                                                                                            cGGTTGAGTACCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCGTTAG‑CCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAAt                       <  2:668416/151‑1 (MQ=255)                                                                                                                                tGAGTACCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCTTTAG‑CCTTTATTCAGCCGGGTGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCga                    <  2:279212/150‑1 (MQ=255)                                                                                                                                         gtgtCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCGTTAG‑CCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAAtt          <  1:135465/151‑1 (MQ=255)                                                                                                                                         gtgtCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCGTTAG‑CCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAAt           <  1:310244/150‑1 (MQ=255)                                                                                                                                          tgtCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCGTTAG‑CCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAAtt          <  2:621634/150‑1 (MQ=255)                                                                                                                                             cgcgAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCGTTAG‑CCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGt       <  1:158321/150‑1 (MQ=255)                                                                                                                                              gcgAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCGTTAG‑CCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTg     <  2:692625/151‑1 (MQ=255)                                                                                                                                                  aGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCGTTAG‑CCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAg  >  2:889816/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                       ACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGAGTACCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAA‑GTTGACTGGATT‑GCCG‑AACCGTTAG‑CCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAG  >  NC_000913/2348485‑2348779

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: ‑