Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,082,219	T→C	T181A (ACG→GCG)	fdoH ←	formate dehydrogenase‑O, Fe‑S subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,082,219	0	T	C	100.0%	189.7 / NA	54	T181A (ACG→GCG)	fdoH	formate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/24);  total (30/24)

CAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGTTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACT  >  NC_000913/4082076‑4082365                                                                                                                                                |                                                                                                                                                   cagcaCGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAg                                                                                                                                             <  1:513984/151‑1 (MQ=255)  agcaCGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAgg                                                                                                                                            <  1:564733/151‑1 (MQ=255)         aCATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCa                                                                                                                                      <  2:189837/150‑1 (MQ=255)                  gtgtACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTg                                                                                                                             <  1:906499/150‑1 (MQ=255)                   tgtACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGt                                                                                                                           >  2:262421/1‑151 (MQ=255)                   tgtACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGt                                                                                                                           >  2:182902/1‑151 (MQ=255)                    gtACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGtt                                                                                                                          >  2:118199/1‑151 (MQ=255)                          ccGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGc                                                                                                                     <  2:380156/150‑1 (MQ=255)                           cGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGccc                                                                                                                   >  1:873711/1‑151 (MQ=255)                                  tGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAaccacc                                                                                                            >  2:682303/1‑151 (MQ=255)                                   gCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAaccacc                                                                                                            >  1:42231/1‑150 (MQ=255)                                   gCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAaccacc                                                                                                            <  2:545469/150‑1 (MQ=255)                                       gATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCAcgcg                                                                                                       <  2:828172/151‑1 (MQ=255)                                         tCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGATGGTTCTTGCCCAACCACCACGCGGt                                                                                                     >  2:349435/1‑151 (MQ=255)                                          cttACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTc                                                                                                    >  1:254225/3‑151 (MQ=255)                                           gTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTc                                                                                                    >  2:557576/1‑150 (MQ=255)                                            tACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCTCTCGCCCGCCAGCGTTTTCACCGACTCTTTCGTACCGACGTGAACCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAacc                                                                                                               <  2:609760/138‑1 (MQ=255)                                             aCAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCaa                                                                                                  <  2:512524/150‑1 (MQ=255)                                                    cGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGc                                                                                          >  2:367878/1‑151 (MQ=255)                                                     gCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCg                                                                                         >  1:814888/1‑151 (MQ=255)                                                       aTTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTa                                                                                       <  1:358009/151‑1 (MQ=255)                                                                  cgcgGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGAc                                                                            >  2:167971/1‑151 (MQ=255)                                                                   gcgGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGAcg                                                                           <  2:435776/151‑1 (MQ=255)                                                                   gcgGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGAcg                                                                           >  1:179482/1‑151 (MQ=255)                                                                       gTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGt                                                                       <  2:177239/151‑1 (MQ=255)                                                                           tcagctcagcAACGCGCTCGCTTGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGt                                                                    <  2:664002/150‑1 (MQ=255)                                                                                 cagcAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCg                                                              <  1:465059/150‑1 (MQ=255)                                                                                     aaCGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGtt                                                          <  2:807385/150‑1 (MQ=255)                                                                                       cgcgCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTg                                                         >  1:862954/1‑149 (MQ=255)                                                                                          gctcgctcGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGgcg                                                    >  2:833988/1‑151 (MQ=255)                                                                                            tcgctcgcCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGcgcg                                                  >  2:279537/1‑151 (MQ=255)                                                                                              gctcgcCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCgg                                                 >  2:87280/1‑150 (MQ=255)                                                                                                 cgcCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCCGaa                                               >  2:867118/1‑149 (MQ=255)                                                                                                 cgcCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGGTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTTAGGCCCGGAATg                                             >  2:197263/1‑151 (MQ=255)                                                                                                   ccAGCGTTTTCATCGACTCTTTCGTACCGAAGTGTATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTTGTTCTTGCCCATCCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTTTCTTCCGGGTTGAGGCGCGGAATGTc                                           <  2:370480/151‑1 (MQ=255)                                                                                                   ccAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGt                                            <  2:358458/150‑1 (MQ=255)                                                                                                    cAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCg                                          >  2:304203/1‑151 (MQ=255)                                                                                                    cAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCg                                          >  2:294651/1‑151 (MQ=255)                                                                                                    cAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCg                                          >  1:611173/1‑151 (MQ=255)                                                                                                     aGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGa                                         <  2:5077/151‑1 (MQ=255)                                                                                                          tttCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACgg                                     <  2:553570/150‑1 (MQ=255)                                                                                                                cGACTCTTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGccc                              <  2:79688/151‑1 (MQ=255)                                                                                                                    tctTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCgca                           >  1:138292/1‑150 (MQ=255)                                                                                                                     ctTTCGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCgcaat                         >  2:230545/1‑151 (MQ=255)                                                                                                                         cGTACCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCgcaatgca                      >  1:443406/1‑150 (MQ=255)                                                                                                                            aCCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATaa                  <  1:629335/151‑1 (MQ=255)                                                                                                                            aCCGAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATa                   <  1:137188/150‑1 (MQ=255)                                                                                                                               gAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCg               >  1:857005/1‑151 (MQ=255)                                                                                                                               gAAGTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCg               <  1:864970/151‑1 (MQ=255)                                                                                                                                  gTGAATCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAg            <  1:846309/151‑1 (MQ=255)                                                                                                                                     aaTCGCGCCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCCATGCAATAACCGCAGc           >  1:218542/1‑149 (MQ=255)                                                                                                                                         gcgcCCGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGc     >  2:761009/1‑151 (MQ=255)                                                                                                                                            cccGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCAc   >  1:506936/1‑150 (MQ=255)                                                                                                                                            cccGCTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACt  >  1:547741/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                   CAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGCGCCCGTTGGGCAGGTCTTCACGCAGGCTGGTTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACT  >  NC_000913/4082076‑4082365

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