Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,086,618	A→G	intergenic (+216/‑88)	yqgC → / → metK	uncharacterized protein/S‑adenosylmethionine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,086,618	0	A	G	100.0%	203.5 / NA	55	intergenic (+216/‑88)	yqgC/metK	uncharacterized protein/S‑adenosylmethionine synthetase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/22);  total (33/22)

AGGAAAAATTTCCAGCCACGTTTTTAACAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCT  >  NC_000913/3086477‑3086765                                                                                                                                              |                                                                                                                                                    aGGAAAAATTTCCAGCCACGTTTTTAACAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGc                                                                                                                                               >  1:13098/1‑148 (MQ=255)    aaaaaTTTCCAGCCACGTTTTTAACAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTc                                                                                                                                         >  1:322194/1‑151 (MQ=255)     aaaaTTTCCAGCCACGTTTTTAACAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCg                                                                                                                                        <  1:81856/151‑1 (MQ=255)       aaTTTCCAGCCACGTTTTTAACAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTg                                                                                                                                      >  1:46701/1‑151 (MQ=255)        aTTTCCAGCCACGTTTTTAACAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTG‑TTTTGACCTGCAGAGTCGTGGt                                                                                                                                    <  2:326712/151‑1 (MQ=255)             cAGCCACGTTTTTAACAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAgg                                                                                                                                 >  2:343255/1‑150 (MQ=255)                        ttAACAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACcac                                                                                                                     >  1:236471/1‑151 (MQ=255)                         tAACAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACcac                                                                                                                     <  2:159792/150‑1 (MQ=255)                          aaCAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACcaca                                                                                                                    <  1:373937/150‑1 (MQ=255)                            cAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGa                                                                                                                  >  1:395555/1‑150 (MQ=255)                               aaTGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAAt                                                                                                              >  1:380781/1‑151 (MQ=255)                               aaTGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAAt                                                                                                              <  1:441014/151‑1 (MQ=255)                                   agaCACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCcac                                                                                                          <  2:578418/151‑1 (MQ=255)                                       acGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCcacaca                                                                                                       >  1:631701/1‑150 (MQ=255)                                             tcaaaaaaaGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGttt                                                                                                >  1:331021/3‑151 (MQ=255)                                             cAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGttt                                                                                                >  2:12230/1‑151 (MQ=255)                                               aaaaaaaGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTg                                                                                               >  2:658123/1‑150 (MQ=255)                                                 aaaaaGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGc                                                                                            <  2:172395/151‑1 (MQ=255)                                                   aaaGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTa                                                                                          <  1:437621/151‑1 (MQ=255)                                                       tGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAAcc                                                                                       >  1:658585/1‑150 (MQ=255)                                                         gAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCaa                                                                                     >  1:576285/1‑150 (MQ=255)                                                            aTAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATtc                                                                                 >  1:277662/1‑151 (MQ=255)                                                              aGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATtctc                                                                               <  1:537544/151‑1 (MQ=255)                                                              aGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATtctc                                                                               <  1:394129/151‑1 (MQ=255)                                                               gggTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTctct                                                                              >  2:650304/1‑151 (MQ=255)                                                                ggTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCtt                                                                             <  1:126940/151‑1 (MQ=255)                                                                    aagaaTTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAgg                                                                         <  1:67822/151‑1 (MQ=255)                                                                     agaaTTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGt                                                                        <  2:536837/151‑1 (MQ=255)                                                                      gaaTTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTg                                                                       >  2:460674/1‑151 (MQ=255)                                                                         ttGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGat                                                                     >  2:518252/1‑150 (MQ=255)                                                                         ttGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGat                                                                     >  1:581374/1‑150 (MQ=255)                                                                                aaaaTAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTaaa                                                               <  2:485706/149‑1 (MQ=255)                                                                                   aTAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATgg                                                          <  2:194780/151‑1 (MQ=255)                                                                                        cATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCaaa                                                      >  2:383340/1‑150 (MQ=255)                                                                                         aTCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAac                                                    >  2:255300/1‑151 (MQ=255)                                                                                          accAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGATTTTGACGTGCAGAGTCGTGGTAGGATCCGCTAAGACAGAACATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCaaaa                                                     >  1:662329/2‑149 (MQ=255)                                                                                           ccAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAacac                                                  >  1:456686/1‑151 (MQ=255)                                                                                             aGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCt                                                <  1:615348/151‑1 (MQ=255)                                                                                                tGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCtttt                                             >  2:411501/1‑151 (MQ=255)                                                                                                 gTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCtttt                                             <  1:336804/150‑1 (MQ=255)                                                                                                       ccatccatACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTcc                                      >  1:265913/1‑151 (MQ=255)                                                                                                        catccatACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCg                                     >  2:427326/1‑151 (MQ=255)                                                                                                        catccatACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCg                                     <  2:145225/151‑1 (MQ=255)                                                                                                             atACCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTcc                                <  2:442748/151‑1 (MQ=255)                                                                                                               aCCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCg                               <  2:236472/150‑1 (MQ=255)                                                                                                               aCCGATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGt                              <  2:631711/151‑1 (MQ=255)                                                                                                                  gATTAACACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGtctc                           >  2:500114/1‑151 (MQ=255)                                                                                                                      aaCACTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGa                        >  1:514892/1‑150 (MQ=255)                                                                                                                        cacTCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAgg                     <  1:518249/151‑1 (MQ=255)                                                                                                                           tCAGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGc                   >  2:524285/1‑150 (MQ=255)                                                                                                                             aGACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATc                >  1:119073/1‑151 (MQ=255)                                                                                                                              gACTGCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATcc               >  1:659756/1‑151 (MQ=255)                                                                                                                                  gCCAGTGTTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGAc           >  2:543250/1‑151 (MQ=255)                                                                                                                                      gtgtTTTTGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACaaa        >  2:39775/1‑150 (MQ=255)                                                                                                                                           tttGACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCt  >  1:641650/1‑151 (MQ=255)                                                                                                                                              |                                                                                                                                                    AGGAAAAATTTCCAGCCACGTTTTTAACAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCT  >  NC_000913/3086477‑3086765

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