Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,227,170	T→C	H34R (CAC→CGC)	yohF ←	putative oxidoreductase with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,227,170	0	T	C	100.0%	195.5 / NA	55	H34R (CAC→CGC)	yohF	putative oxidoreductase with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/30);  total (25/30)

CCCAGCCGTTGAATGAGTTTCTCCAGCGCCAGTGCCCCTTCTGGTAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGTGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTA  >  NC_000913/2227035‑2227302                                                                                                                                        |                                                                                                                                     cccAGCCGTTGAATGAGTTTCTCCAGCGCCAGTGCCCCTTCTGGTAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCa                                                                                                                                  >  1:947167/1‑140 (MQ=255)    agtcGTTGAATGAGTTTCTCCAGCGCCAGTGCCCCTTCTGGTAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGt                                                                                                                               <  1:765096/137‑1 (MQ=255)       cGTTGAATGAGTTTCTCCAGCGCCAGTGCCCCTTCTGGTAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAAt                                                                                                                            <  2:648618/140‑1 (MQ=255)          tGAATGAGTTTCTCCAGCGCCAGTGCCCCTTCTGGTAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATAcc                                                                                                                         <  2:940165/140‑1 (MQ=255)                      tcCAGCGCCAGTGCCCCTTCTGGTAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAAccc                                                                                                             >  2:734228/1‑140 (MQ=255)                          gcgcCAGTGCCCCTTCTGGTAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCctgc                                                                                                         >  2:149295/1‑140 (MQ=255)                          gcgcCAGTGCCCCTTCTGGTAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCctgc                                                                                                         >  1:779519/1‑140 (MQ=255)                           cgcCAGTGCCCCTTCTGGTAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCctgct                                                                                                        >  1:795284/1‑140 (MQ=255)                             ccAGTGCCCCTTCTGGTAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCtgctgc                                                                                                      <  1:587899/140‑1 (MQ=255)                                gTGCCCCTTCTGGTAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGcc                                                                                                   >  2:324522/1‑140 (MQ=255)                                      cTTCTGGTAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCGATATCAAACCCCTGCTGCGCCAGtaat                                                                                             <  2:887746/140‑1 (MQ=255)                                           ggTAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAAcgc                                                                                        >  2:787831/1‑140 (MQ=255)                                           ggTAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAAcgc                                                                                        >  2:661558/1‑140 (MQ=255)                                            gTAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAAcgcg                                                                                       <  1:734228/140‑1 (MQ=255)                                             tAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCctgct                                                                                                        <  1:914526/122‑1 (MQ=255)                                             tAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCctgct                                                                                                        >  2:914526/1‑122 (MQ=255)                                              aGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCa                                                                                     <  1:671002/140‑1 (MQ=255)                                                  tGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCActct                                                                                 >  1:443425/1‑140 (MQ=255)                                                  tGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCActct                                                                                 >  1:935228/1‑140 (MQ=255)                                                             aGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATccc                                                                      <  1:787831/140‑1 (MQ=255)                                                             aGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATccc                                                                      <  1:661558/140‑1 (MQ=255)                                                                tGCACGATCTCCGCACGTACGCCGTGGCTAACTACCCCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGa                                                                   <  1:66732/140‑1 (MQ=255)                                                                   aCGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATc                                                                >  1:915028/1‑140 (MQ=255)                                                                    cGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCg                                                               <  1:196451/140‑1 (MQ=255)                                                                    cGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCg                                                               <  2:3911/140‑1 (MQ=255)                                                                       tctcCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAg                                                            >  1:14855/1‑140 (MQ=255)                                                                              aCGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGtaa                                                     <  2:935228/140‑1 (MQ=255)                                                                               cGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGtaa                                                     >  1:932820/1‑139 (MQ=255)                                                                               cGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGtaa                                                     <  2:932820/139‑1 (MQ=255)                                                                                 tACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAgg                                                           >  1:783114/1‑131 (MQ=255)                                                                                 tACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAgg                                                           <  2:783114/131‑1 (MQ=255)                                                                                 tACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGtaataa                                                  >  2:126405/1‑140 (MQ=255)                                                                                   cGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATc                                                >  1:930394/1‑140 (MQ=255)                                                                                       gTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCaa                                            >  1:795055/1‑140 (MQ=255)                                                                                         ggCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAAcc                                          >  1:391928/1‑140 (MQ=255)                                                                                          gCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCt                                         >  2:421379/1‑140 (MQ=255)                                                                                              aCTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCg                                               <  2:687675/130‑1 (MQ=255)                                                                                              aCTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCg                                               >  1:687675/1‑130 (MQ=255)                                                                                                  ccTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATc                                 <  2:559422/140‑1 (MQ=255)                                                                                                    tCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGa                               <  2:14855/140‑1 (MQ=255)                                                                                                       cgcgcgGTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGtt                            <  1:421379/140‑1 (MQ=255)                                                                                                             gTATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCAc                      >  2:439411/1‑140 (MQ=255)                                                                                                              tATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACt                     <  2:391928/140‑1 (MQ=255)                                                                                                              tATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACt                     <  1:149295/140‑1 (MQ=255)                                                                                                              tATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACt                     <  2:915028/140‑1 (MQ=255)                                                                                                              tATCTTTTGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACt                     <  1:439411/140‑1 (MQ=255)                                                                                                                  ttttGCCCCTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAAc                 <  2:795055/140‑1 (MQ=255)                                                                                                                        cccTTCTTCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGaa           >  1:854139/1‑140 (MQ=255)                                                                                                                          cttcttCATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATa         <  2:854139/140‑1 (MQ=255)                                                                                                                           ttcttcATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATaa        >  1:507496/1‑140 (MQ=255)                                                                                                                            tcttcATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATaaa       >  1:777502/1‑140 (MQ=255)                                                                                                                              ttcATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACg     <  2:507496/140‑1 (MQ=255)                                                                                                                              ttcATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACg     <  2:482805/140‑1 (MQ=255)                                                                                                                               tcATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGt    <  2:795284/140‑1 (MQ=255)                                                                                                                               tcATCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGt    <  1:126405/140‑1 (MQ=255)                                                                                                                                 aTCTGAGCGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTa  <  2:443425/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     CCCAGCCGTTGAATGAGTTTCTCCAGCGCCAGTGCCCCTTCTGGTAGATTGCCGAGATCCAGCTGCACGATCTCCGCACGTACGCCGTGGCTAACTACCTCACGCGCGGTATCTTTTGCCCCTTCTTCATCTGAGTGCCAGGTAATACCAATATCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTA  >  NC_000913/2227035‑2227302

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