Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,873,099	T→C	D80G (GAC→GGC)	ygbE ←	DUF3561 family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,873,099	0	T	C	100.0%	188.4 / NA	52	D80G (GAC→GGC)	ygbE	DUF3561 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/21);  total (31/21)

AGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGTCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAA  >  NC_000913/2872951‑2873243                                                                                                                                                     |                                                                                                                                                 aGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGccc                                                                                                                                                >  2:119807/1‑151 (MQ=255)    aCATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATc                                                                                                                                             <  2:169316/151‑1 (MQ=255)     cATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATc                                                                                                                                             >  1:541672/1‑150 (MQ=255)      aTAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCa                                                                                                                                            >  2:606148/1‑150 (MQ=255)        aaTCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCaatg                                                                                                                                         >  2:60544/1‑151 (MQ=255)        aaTCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCaatg                                                                                                                                         <  2:600282/151‑1 (MQ=255)           ctcATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCaatgaat                                                                                                                                      >  2:33431/1‑151 (MQ=255)                        aCGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATccc                                                                                                                         >  1:606612/1‑151 (MQ=255)                         cGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATcccc                                                                                                                        >  2:390443/1‑151 (MQ=255)                          ggTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCcca                                                                                                                       <  2:97636/151‑1 (MQ=255)                          ggTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCcca                                                                                                                       <  2:194285/151‑1 (MQ=255)                           gTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCccac                                                                                                                      >  2:315707/1‑151 (MQ=255)                               gAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCcaccac                                                                                                                   <  2:538750/150‑1 (MQ=255)                                aaTTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTa                                                                                                                 <  1:133498/151‑1 (MQ=255)                                  ttATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTAcc                                                                                                               >  2:234622/1‑151 (MQ=255)                                   tATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACcg                                                                                                              >  2:437038/1‑151 (MQ=255)                                   tATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACcg                                                                                                              >  2:810041/1‑151 (MQ=255)                                       aaCAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGAc                                                                                                          >  2:15669/1‑151 (MQ=255)                                       aaCAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGAc                                                                                                          >  2:150362/1‑151 (MQ=255)                                                gaCGGGAAAGTACAAATTAGCCCAGTAGCCCCATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAg                                                                                                 <  1:809985/151‑1 (MQ=255)                                                  cGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAgcg                                                                                               >  2:152245/1‑151 (MQ=255)                                                    ggAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGcc                                                                                             <  1:539530/151‑1 (MQ=255)                                                    ggAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGcc                                                                                             >  2:45017/1‑151 (MQ=255)                                                      aaaGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGcc                                                                                             >  2:350374/1‑149 (MQ=255)                                                          tACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCagaa                                                                                         >  1:802770/1‑149 (MQ=255)                                                             aaaTTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCagaaagaa                                                                                     >  2:665375/1‑150 (MQ=255)                                                                     ccAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGgt                                                                            >  1:174769/1‑151 (MQ=255)                                                                        gTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTa                                                                         <  2:429869/151‑1 (MQ=255)                                                                                aTAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAaaga                                                                 <  1:516113/151‑1 (MQ=255)                                                                                  aaaCAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAaagaag                                                               <  1:437004/151‑1 (MQ=255)                                                                                   aaCAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAaaga                                                                 >  1:196501/1‑148 (MQ=255)                                                                                     cAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAAc                                                            >  2:653229/1‑151 (MQ=255)                                                                                        tGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGc                                                         >  1:239361/1‑151 (MQ=255)                                                                                         gcgcCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCg                                                        >  1:203471/1‑151 (MQ=255)                                                                                         gcgcCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCg                                                        <  2:203494/151‑1 (MQ=255)                                                                                           gcCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCgtg                                                      >  1:46882/1‑151 (MQ=255)                                                                                             cAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGtt                                                    >  1:807308/1‑151 (MQ=255)                                                                                             cAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGtt                                                    >  2:289175/1‑151 (MQ=255)                                                                                             cAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGtt                                                    <  2:134892/151‑1 (MQ=255)                                                                                               aaCATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGtt                                                    >  1:6786/1‑149 (MQ=255)                                                                                               aaCATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTag                                                  <  1:178901/151‑1 (MQ=255)                                                                                                   tAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGcc                                              >  2:601038/1‑151 (MQ=255)                                                                                                   tAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGcc                                              >  1:752241/1‑151 (MQ=255)                                                                                                    aaTGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGcc                                              >  2:719303/1‑150 (MQ=255)                                                                                                       gCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATaa                                          <  1:360138/151‑1 (MQ=255)                                                                                                       gCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATa                                           <  2:518450/150‑1 (MQ=255)                                                                                                        ccTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATaaa                                         <  2:46904/151‑1 (MQ=255)                                                                                                           aCTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAAt                                        <  1:555864/149‑1 (MQ=255)                                                                                                                 aCCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCg                                >  1:325410/1‑151 (MQ=255)                                                                                                                                        gcgtagcTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAg         <  1:90311/151‑1 (MQ=255)                                                                                                                                            agcTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGgaga     <  1:33409/151‑1 (MQ=255)                                                                                                                                                ttCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATaa  <  1:289150/150‑1 (MQ=255)                                                                                                                                                ttCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATaa  >  2:157069/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                 AGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGTCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAA  >  NC_000913/2872951‑2873243

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