Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,085,935	C→A	intergenic (‑113/‑81)	fdoG ← / → fdhD	formate dehydrogenase O subunit alpha/sulfurtransferase for molybdenum cofactor sulfuration

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,085,935	0	C	A	100.0%	187.8 / NA	55	intergenic (‑113/‑81)	fdoG/fdhD	formate dehydrogenase O subunit alpha/sulfurtransferase for molybdenum cofactor sulfuration
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (29/26);  total (29/26)

CTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAA  >  NC_000913/4085802‑4086056                                                                                                                                      |                                                                                                                          cTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAAtt                                                                                                                     >  1:506946/1‑140 (MQ=255)  tGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAAttt                                                                                                                    >  1:239314/1‑140 (MQ=255)      ttCTACTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCt                                                                                                                <  1:282653/140‑1 (MQ=255)        ctgctgACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCtcc                                                                                                              >  1:275972/1‑140 (MQ=255)              aCCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAg                                                                                                        >  2:293922/1‑140 (MQ=255)              aCCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAg                                                                                                        >  2:207833/1‑140 (MQ=255)              aCCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAg                                                                                                        >  1:527498/1‑140 (MQ=255)                 tGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGaa                                                                                                     <  1:207833/140‑1 (MQ=255)                 tGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGaa                                                                                                     <  2:506946/140‑1 (MQ=255)                     tGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGt                                                                                                 <  1:661179/140‑1 (MQ=255)                     tGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGt                                                                                                 <  1:350010/140‑1 (MQ=255)                           gCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTaa                                                                                           >  1:234458/1‑140 (MQ=255)                           gCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTaa                                                                                           >  2:478057/1‑140 (MQ=255)                             tCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACa                                                                                         <  1:293922/140‑1 (MQ=255)                             tCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAACGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGa                                                                                                      >  2:582660/1‑127 (MQ=255)                             tCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAACGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGa                                                                                                      <  1:582660/127‑1 (MQ=255)                               cTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGa                                                                                       <  2:239314/140‑1 (MQ=255)                                  gTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATa                                                                                    >  1:347487/1‑140 (MQ=255)                                      cGACATTGTCACGTCCGATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCAt                                                                                >  2:668192/1‑140 (MQ=255)                                      cGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACa                                                                                         >  2:34454/1‑131 (MQ=255)                                      cGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACa                                                                                         <  1:34454/131‑1 (MQ=255)                                      cGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCAt                                                                                >  1:270871/1‑140 (MQ=255)                                         cATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTTACAGAATACCATAAt                                                                             >  2:754922/1‑140 (MQ=255)                                            tGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGa                                                                                       <  1:362850/127‑1 (MQ=255)                                            tGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGa                                                                                       >  2:362850/1‑127 (MQ=255)                                             gTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTg                                                                         >  1:351493/1‑140 (MQ=255)                                             gTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTg                                                                         >  2:45185/1‑140 (MQ=255)                                              tCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTgg                                                                        <  1:45185/140‑1 (MQ=255)                                                   tCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGgtgtgt                                                                   <  2:275972/140‑1 (MQ=255)                                                   tCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGgtgtgt                                                                   >  2:29953/1‑140 (MQ=255)                                                   acccATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGgtgtgt                                                                   <  1:448386/139‑1 (MQ=255)                                                        tttACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGgtgtgt                                                                   <  1:659086/135‑1 (MQ=255)                                                        tttACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGgtgtgt                                                                   >  2:659086/1‑135 (MQ=255)                                                        tttACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGtt                                                              >  1:710263/1‑140 (MQ=255)                                                             aTTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTAt                                                         <  1:483315/140‑1 (MQ=255)                                                             aTTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTAt                                                         <  2:347487/140‑1 (MQ=255)                                                             aTTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTAt                                                         <  2:270871/140‑1 (MQ=255)                                                             aTTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTAt                                                         <  1:270061/140‑1 (MQ=255)                                                              ttCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATc                                                        <  2:527498/140‑1 (MQ=255)                                                                         cGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTAt                                                         >  2:79305/1‑128 (MQ=255)                                                                         cGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTAt                                                         <  1:79305/128‑1 (MQ=255)                                                                               cAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGa                                       <  2:351493/140‑1 (MQ=255)                                                                                  ggagtgggagtTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGaaaa                                    >  2:37946/1‑140 (MQ=255)                                                                                     gcgggagtTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGaaaaaaa                                 >  2:443450/3‑140 (MQ=255)                                                                                      tgggagtTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAac                                >  1:57789/1‑140 (MQ=255)                                                                                      tgggagtTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAac                                >  1:216034/1‑140 (MQ=255)                                                                                       gggagtTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAaca                               <  2:216034/140‑1 (MQ=255)                                                                                              aTTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGaa                        >  1:426932/1‑140 (MQ=255)                                                                                              aTTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGaa                        >  1:324128/1‑140 (MQ=255)                                                                                                    cTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGaaa                  >  1:494082/1‑140 (MQ=255)                                                                                                                 gCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGAc     <  2:234458/140‑1 (MQ=255)                                                                                                                 gCCGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGAc     <  1:668192/140‑1 (MQ=255)                                                                                                                  ccGCAAAGGTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACa    <  1:478057/140‑1 (MQ=255)                                                                                                                          gTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACaaa  <  1:655492/134‑1 (MQ=255)                                                                                                                          gTCCAGAATTGGATCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACaaa  >  2:655492/1‑134 (MQ=255)                                                                                                                                      |                                                                                                                          CTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAA  >  NC_000913/4085802‑4086056

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