Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,263,580	A→G	C27C (TGT→TGC)	tdcD ←	propionate kinase/acetate kinase C, anaerobic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,263,580	0	A	G	100.0%	256.9 / NA	69	C27C (TGT→TGC)	tdcD	propionate kinase/acetate kinase C, anaerobic
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/33);  total (36/33)

GCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGA  >  NC_000913/3263464‑3263705                                                                                                                     |                                                                                                                              gCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACgg                                                                                                         >  1:383747/1‑139 (MQ=255)  ccTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGa                                                                                                        >  2:396144/1‑139 (MQ=255)    ttCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGaaa                                                                                                      >  2:814245/1‑139 (MQ=255)                   tAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACg                                                                                       <  2:315066/139‑1 (MQ=255)                      cTGTGGTGAGCCAGCGGTGCTGGCTCTCCTCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAAc                                                                                    >  1:1319697/1‑139 (MQ=255)                         tggtgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACcac                                                                                 <  1:1247773/139‑1 (MQ=255)                         tggtgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACcac                                                                                 >  1:938898/1‑139 (MQ=255)                          ggtgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACcaca                                                                                >  1:1421713/1‑139 (MQ=255)                             gAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCGGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGTACcac                                                                                 >  1:966125/1‑135 (MQ=255)                             gAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACcac                                                                                 <  2:966125/135‑1 (MQ=255)                              aGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTcgca                                                                                                                             <  1:1040183/90‑1 (MQ=255)                              aGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATGAATACTTcgca                                                                                                                             >  2:1040183/1‑90 (MQ=255)                               gCCAGCGGTGCTGGCTCTCCTCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTtaa                                                                           <  2:1319697/139‑1 (MQ=255)                                ccAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCt                                                                                                                       >  2:649980/1‑94 (MQ=255)                                ccAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCt                                                                                                                       <  1:649980/94‑1 (MQ=255)                                        gCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaaa                                                                  <  2:1421713/139‑1 (MQ=255)                                                 cccccATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTta                                                                            >  2:684686/1‑120 (MQ=255)                                                 cccccATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTta                                                                            <  1:684686/120‑1 (MQ=255)                                                 cccccATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGaa                                                         <  2:383747/139‑1 (MQ=255)                                                  ccccATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGACAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACTACAGTTAATAACCAAAAAAACCGGaaa                                                        >  2:1066728/1‑139 (MQ=255)                                                        ttACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAattcatt                                                  <  2:938898/139‑1 (MQ=255)                                                          aCGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAAttcattca                                                >  2:113772/1‑139 (MQ=255)                                                          aCGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAAttcattca                                                >  2:599594/1‑139 (MQ=255)                                                          aCGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAAttcattca                                                <  2:1393394/139‑1 (MQ=255)                                                          aCGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAAttcattca                                                >  2:3161/1‑139 (MQ=255)                                                            ggATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACTGGAAACTCATTCAtc                                              >  2:365188/1‑139 (MQ=255)                                                            ggATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCg                                                            <  1:696016/125‑1 (MQ=255)                                                            ggATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCg                                                            >  2:696016/1‑125 (MQ=255)                                                              aTAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACca                                                                                  <  1:104170/101‑1 (MQ=255)                                                              aTAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACca                                                                                  >  2:104170/1‑101 (MQ=255)                                                               tAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCAt                                                                                                                  >  1:986149/1‑68 (MQ=255)                                                               tAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCAt                                                                                                                  <  2:986149/68‑1 (MQ=255)                                                                 aGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCAtctc                                            <  2:999619/136‑1 (MQ=255)                                                                 aGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCAtctc                                            >  1:999619/1‑136 (MQ=255)                                                                          tttCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCt                                >  1:1430209/1‑139 (MQ=255)                                                                            tCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGa                                                                                                        <  1:129951/65‑1 (MQ=255)                                                                            tCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGa                                                                                                        >  2:129951/1‑65 (MQ=255)                                                                                      tACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGaa                                                                                     >  2:652993/1‑74 (MQ=255)                                                                                      tACCGTCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGaa                                                                                     <  1:652993/74‑1 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAg                                                                                                              <  1:588099/44‑1 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAg                                                                                                              >  2:588099/1‑44 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACg                                                                                                          >  2:81002/1‑48 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACg                                                                                                          <  1:81002/48‑1 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATc                                                                                              <  2:1031252/60‑1 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATc                                                                                              >  1:1031252/1‑60 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAAc                                                                                    <  1:336373/70‑1 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAAc                                                                                    >  2:336373/1‑70 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACca                                                                                  <  1:170388/72‑1 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACca                                                                                  >  2:170388/1‑72 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCTGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAAc               <  1:365188/139‑1 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAatt                                                      >  1:275699/1‑100 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAatt                                                      <  2:275699/100‑1 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCAtctc                                            >  1:540164/1‑110 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCAtctc                                            <  2:540164/110‑1 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATGCTGAGTTACGGATTAAAAc               >  2:986113/1‑139 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAAc               >  1:97168/1‑139 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAAc               >  1:607972/1‑139 (MQ=255)                                                                                           tCGGCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAAc               >  1:160767/1‑139 (MQ=255)                                                                                             ggCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAg             <  1:814245/139‑1 (MQ=255)                                                                                              gCAATACCTGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGt            <  2:1109431/139‑1 (MQ=255)                                                                                                      tGACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATAc    <  2:1430209/139‑1 (MQ=255)                                                                                                       gACATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACg   <  1:113772/139‑1 (MQ=255)                                                                                                        aCATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCtttt                                        >  1:435991/1‑101 (MQ=255)                                                                                                        aCATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCtttt                                        <  2:435991/101‑1 (MQ=255)                                                                                                        aCATTAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGa  <  1:3161/139‑1 (MQ=255)                                                                                                            tAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACg                                                                                       <  1:490910/50‑1 (MQ=255)                                                                                                            tAATACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACg                                                                                       >  2:490910/1‑50 (MQ=255)                                                                                                             aaTACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACg                                                                                       >  2:1440590/1‑49 (MQ=255)                                                                                                             aaTACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACg                                                                                       <  1:1440590/49‑1 (MQ=255)                                                                                                              aTACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTtc                                       >  1:454395/1‑96 (MQ=255)                                                                                                              aTACTTCGCAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTtc                                       <  2:454395/96‑1 (MQ=255)                                                                                                                     |                                                                                                                              GCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGA  >  NC_000913/3263464‑3263705

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