Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,353,121	C→T	*310* (TAG→TAA)	yhcC ←	putative Fe‑S oxidoreductase, Radical SAM superfamily protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,353,121	0	C	T	100.0%	212.4 / NA	65	*310* (TAG→TAA)	yhcC	putative Fe‑S oxidoreductase, Radical SAM superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/30);  total (35/30)

GTCAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACCTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCGGG  >  NC_000913/3352987‑3353250                                                                                                                                         |                                                                                                                                  gTCAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACt                                                                                                                               <  1:1429028/139‑1 (MQ=255)     aCATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGt                                                                                                                           >  2:28254/1‑139 (MQ=255)                cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGcc                                                                                                                >  1:1061770/1‑139 (MQ=255)                                       ttAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATc                                                                                           >  2:725243/1‑139 (MQ=255)                                       ttAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATccc                                                                                         >  2:934349/1‑139 (MQ=255)                                          gggAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGt                                                                                      >  2:1282211/1‑139 (MQ=255)                                               ttCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGcc                                                                                               <  2:898007/125‑1 (MQ=255)                                               ttCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGcc                                                                                               >  1:898007/1‑125 (MQ=255)                                               ttCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGcc                                                                                 >  2:711443/1‑139 (MQ=255)                                               ttCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGcc                                                                                 >  1:387649/1‑139 (MQ=255)                                                tCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCa                                                                                >  2:925669/1‑139 (MQ=255)                                                  cTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATg                                                                              <  2:618440/139‑1 (MQ=255)                                                     cACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttc                                                                           <  2:387649/139‑1 (MQ=255)                                                      aCGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttca                                                                          >  1:1481281/1‑139 (MQ=255)                                                          cAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcattca                                                                      <  1:934349/139‑1 (MQ=255)                                                             ccTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGAt                                                                   >  2:1112082/1‑139 (MQ=255)                                                                 ccccAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAgg                                                                                                                     <  1:1550642/85‑1 (MQ=255)                                                                 ccccAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAgg                                                                                                                     >  2:1550642/1‑85 (MQ=255)                                                                 ccccAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCg                                                               >  2:1684232/1‑139 (MQ=255)                                                                  cccAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGa                                                              >  2:260559/1‑139 (MQ=255)                                                                      ataataCCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGt                                                                                      >  2:105775/1‑111 (MQ=255)                                                                      ataataCCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGt                                                                                      <  1:105775/111‑1 (MQ=255)                                                                          taCCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTc                                                      <  2:1481281/139‑1 (MQ=255)                                                                            cccAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGa                                                    >  1:961541/1‑139 (MQ=255)                                                                            cccAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGACCCAGCGCCGATCCCTGTACGCCATGttcat                                                                         <  1:761374/118‑1 (MQ=255)                                                                            cccAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGACCCAGCGCCGATCCCTGTACGCCATGttcat                                                                         >  2:761374/1‑118 (MQ=255)                                                                                 ttCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                      >  1:1678503/1‑68 (MQ=255)                                                                                 ttCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                      <  2:1678503/68‑1 (MQ=255)                                                                                 ttCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcattca                                                                      >  1:287614/1‑116 (MQ=255)                                                                                 ttCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcattca                                                                      <  2:287614/116‑1 (MQ=255)                                                                                      aCAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcattca                                                                      >  1:435862/1‑111 (MQ=255)                                                                                      aCAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcattca                                                                      <  2:435862/111‑1 (MQ=255)                                                                                       cAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCa                                                                                >  2:334745/1‑100 (MQ=255)                                                                                       cAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCa                                                                                <  1:334745/100‑1 (MQ=255)                                                                                         gTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCa                                       <  1:711443/139‑1 (MQ=255)                                                                                          ttACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATg                                                                              <  1:1599375/99‑1 (MQ=255)                                                                                          ttACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATg                                                                              >  2:1599375/1‑99 (MQ=255)                                                                                           tACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAg                                                         <  2:1388050/119‑1 (MQ=255)                                                                                           tACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAg                                                         >  1:1388050/1‑119 (MQ=255)                                                                                            aCGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATc                                                                                           >  2:259967/1‑84 (MQ=255)                                                                                            aCGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATc                                                                                           <  1:259967/84‑1 (MQ=255)                                                                                              gTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGAt                                  <  1:1112082/139‑1 (MQ=255)                                                                                                      t‑‑tGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcac                          <  2:1061770/139‑1 (MQ=255)                                                                                                          gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcacca                        >  1:561275/1‑139 (MQ=255)                                                                                                          gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGca                           <  1:407961/136‑1 (MQ=255)                                                                                                          gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGca                           >  2:407961/1‑136 (MQ=255)                                                                                                           tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAgg                                                                                                                     >  2:1338257/1‑45 (MQ=255)                                                                                                           tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAgg                                                                                                                     <  1:1338257/45‑1 (MQ=255)                                                                                                           tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAgcg                                                                                                 >  2:860398/1‑65 (MQ=255)                                                                                                           tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAgcg                                                                                                 <  1:860398/65‑1 (MQ=255)                                                                                                            gtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                            >  2:1489344/1‑69 (MQ=255)                                                                                                            gtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                            <  1:1489344/69‑1 (MQ=255)                                                                                                            gtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACcaca                      <  2:961541/139‑1 (MQ=255)                                                                                                             tgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAg                     >  1:1288893/1‑139 (MQ=255)                                                                                                              gagaAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGc                    <  1:725243/139‑1 (MQ=255)                                                                                                                     gCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGAt                                  <  2:1143273/118‑1 (MQ=255)                                                                                                                     gCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGAt                                  >  1:1143273/1‑118 (MQ=255)                                                                                                                      cGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGagcaag            >  1:766370/1‑139 (MQ=255)                                                                                                                       ggAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGa                                                    <  1:1528048/98‑1 (MQ=255)                                                                                                                       ggAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGa                                                    >  2:1528048/1‑98 (MQ=255)                                                                                                                        gAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGagcaagca          <  2:561275/139‑1 (MQ=255)                                                                                                                        gAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGagcaagca          <  1:1282211/139‑1 (MQ=255)                                                                                                                           ttCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGa                                                    <  1:269615/94‑1 (MQ=255)                                                                                                                           ttCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGa                                                    >  2:269615/1‑94 (MQ=255)                                                                                                                                cTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCggg  >  2:462112/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                  GTCAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAAT‑‑TGTGTGAGAAATGCGGAATTCCCCTCAACGACCTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCGGG  >  NC_000913/3352987‑3353250

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