Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,086,523	A→G	T170A (ACT→GCT)	fdhD →	sulfurtransferase for molybdenum cofactor sulfuration

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,086,523	0	A	G	98.2%	196.5 / ‑4.3	57	T170A (ACT→GCT)	fdhD	sulfurtransferase for molybdenum cofactor sulfuration
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (28/28);  total (29/28)
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.27e-01

CGTGGAGCAACTTAATGACATCGGAAAACCGGTGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGAAAGCTGGCAGCAAGGTGCGGTAC  >  NC_000913/4086390‑4086658                                                                                                                                      |                                                                                                                                        cgTGGAGCAACTTAATGACATCGGAAAACCGGTGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAAGTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCAcg                                                                                                                                   >  2:675588/1‑140 (MQ=255)          aCTTAATGACATCGGAAAACCGGTGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCt                                                                                                                          <  2:698870/140‑1 (MQ=255)                     tCGGAAAACCGGTGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCa                                                                                                               >  2:791911/1‑140 (MQ=255)                      cGGAAAACCGGTGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCAt                                                                                                              >  1:487289/1‑140 (MQ=255)                              cGGTGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGaa                                                                                                      >  2:165285/1‑140 (MQ=255)                                gTGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACt                                                                                                    >  1:940497/1‑140 (MQ=255)                                gTGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACt                                                                                                    >  1:668696/1‑140 (MQ=255)                                gTGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACt                                                                                                    >  1:358645/1‑140 (MQ=255)                                    aGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTc                                                                                                <  2:668696/140‑1 (MQ=255)                                      ccGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTcgg                                                                                              >  2:866318/1‑140 (MQ=255)                                       cGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTcggc                                                                                             <  2:667771/140‑1 (MQ=255)                                           aCCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGc                                                                                         >  1:170992/1‑140 (MQ=255)                                           aCCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGc                                                                                         <  1:628274/140‑1 (MQ=255)                                             cGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCAt                                                                                       <  1:675588/140‑1 (MQ=255)                                              gTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATg                                                                                      <  1:869864/140‑1 (MQ=255)                                                  aCCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGa                                                                                  <  1:620062/140‑1 (MQ=255)                                                    ccAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACg                                                                                >  1:500876/1‑140 (MQ=255)                                                    ccAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACg                                                                                >  1:208751/1‑140 (MQ=255)                                                      aGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTg                                                                              <  1:309227/140‑1 (MQ=255)                                                       gACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTgg                                                                             >  2:1091937/1‑140 (MQ=255)                                                        aCGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTggg                                                                            >  1:514783/1‑140 (MQ=255)                                                        aCGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTggg                                                                            >  1:174565/1‑140 (MQ=255)                                                           tttGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCg                                                                         >  2:535258/1‑140 (MQ=255)                                                           tttGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCg                                                                         <  1:1091937/140‑1 (MQ=255)                                                              gATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCa                                                                      >  2:334623/1‑140 (MQ=255)                                                              gATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCa                                                                      >  1:946286/1‑140 (MQ=255)                                                                tctcAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCAtg                                                                    <  1:165285/140‑1 (MQ=255)                                                                  tcAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCAtgtg                                                                  >  2:545451/1‑140 (MQ=255)                                                                         aCTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTgg                                                           >  2:899132/1‑140 (MQ=255)                                                                          cTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGa                                                          <  1:791911/140‑1 (MQ=255)                                                                            ggATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACa                                                        <  1:535258/140‑1 (MQ=255)                                                                                         cGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGgccg                                           <  1:866318/140‑1 (MQ=255)                                                                                            cATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGccgccg                                        >  1:404203/1‑140 (MQ=255)                                                                                              tctcAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGt                                      <  1:899132/140‑1 (MQ=255)                                                                                              tctcAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGt                                      >  2:559204/1‑140 (MQ=255)                                                                                              tctcAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGt                                      <  1:559204/140‑1 (MQ=255)                                                                                               ctcAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTc                                     >  1:715607/1‑140 (MQ=255)                                                                                                 cAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTcac                                   >  2:738486/1‑140 (MQ=255)                                                                                                   aCGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGt                                      <  1:800397/135‑1 (MQ=255)                                                                                                   aCGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGt                                      >  2:800397/1‑135 (MQ=255)                                                                                                   aCGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACaa                                 >  2:912756/1‑140 (MQ=255)                                                                                                   aCGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACaa                                 >  2:267663/1‑140 (MQ=255)                                                                                                           cAGCCAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGa                         <  1:334623/140‑1 (MQ=255)                                                                                                               cAGTGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGAAAGc                     <  2:715607/140‑1 (MQ=255)                                                                                                                  tGGGGCAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGAAAGCTgg                  >  2:950777/1‑140 (MQ=255)                                                                                                                       cAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGAAAGCTGgcagca             <  1:545451/140‑1 (MQ=255)                                                                                                                       cAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGAAAGCTGgcagca             <  1:738486/140‑1 (MQ=255)                                                                                                                       cAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGAAAGCTGgcagca             <  1:912756/140‑1 (MQ=255)                                                                                                                       cAACTGACTGGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGAAAGCTGgcagca             <  2:946286/140‑1 (MQ=255)                                                                                                                            gactgGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGAAAGCTGGCAGCAAGGTg        <  1:651704/140‑1 (MQ=255)                                                                                                                            gactgGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGAAAGCTGGCAGCAAGGTg        <  1:754259/140‑1 (MQ=255)                                                                                                                            gactgGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGAAAGCTGGCAGCAAGGTg        <  2:824284/140‑1 (MQ=255)                                                                                                                            gactgGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGAAAGCTGGCAGCAAGGTg        >  1:824284/1‑140 (MQ=255)                                                                                                                              ctgGTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGAAAGCTGGCAGCAAGGTGCg      >  2:319451/1‑140 (MQ=255)                                                                                                                                 gTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGAAAGCTGGCAGCAAGGTGCGGTa   <  2:174565/140‑1 (MQ=255)                                                                                                                                 gTTGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGAAAGCTGGCAGCAAGGTGCGGTa   <  1:267663/140‑1 (MQ=255)                                                                                                                                  ttGTGCTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGAAAGCTGGCAGCAAGGTGCGGTAc  <  1:950777/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                        CGTGGAGCAACTTAATGACATCGGAAAACCGGTGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCTGGACAAACTGTTAGGCCGCCGGTCACAAGAAGGGGAAAGCTGGCAGCAAGGTGCGGTAC  >  NC_000913/4086390‑4086658

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