Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,581,825	T→C	Q318R (CAG→CGG)	ypfG ←	DUF1176 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,581,825	0	T	C	100.0%	190.5 / NA	55	Q318R (CAG→CGG)	ypfG	DUF1176 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/24);  total (31/24)

CTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATCAGTTCCAGTTCATTCGTCTCCTGACCGTTGTTGAACGGCAAAC  >  NC_000913/2581684‑2581972                                                                                                                                              |                                                                                                                                                    cTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCa                                                                                                                                             >  2:747524/1‑150 (MQ=255)   tATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCaaa                                                                                                                                           <  1:180136/150‑1 (MQ=255)   tATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCaaa                                                                                                                                           <  1:376/150‑1 (MQ=255)    aTCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAgc                                                                                                                                         >  1:724585/1‑151 (MQ=255)      ccGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAgcgc                                                                                                                                       >  2:850376/1‑151 (MQ=255)                   tCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGa                                                                                                                           <  2:199219/150‑1 (MQ=255)                     ttATCCGGCTTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAAt                                                                                                                         <  1:488303/150‑1 (MQ=255)                     ttATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAAt                                                                                                                         >  2:612043/1‑150 (MQ=255)                       aTCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGcc                                                                                                                      >  1:765188/1‑151 (MQ=255)                        tCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGaa                                                                                                                          <  1:574216/146‑1 (MQ=255)                           ggCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACaa                                                                                                                  <  1:296449/151‑1 (MQ=255)                            gCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAAt                                                                                                                 >  1:771004/1‑151 (MQ=255)                             ccTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATc                                                                                                                >  2:763660/1‑151 (MQ=255)                             ccTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATc                                                                                                                <  1:165593/151‑1 (MQ=255)                              cTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATc                                                                                                                <  1:236908/150‑1 (MQ=255)                                 cATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCtt                                                                                                            <  1:123567/151‑1 (MQ=255)                                  aTGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCtt                                                                                                            >  1:758897/1‑150 (MQ=255)                                    gTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAAt                                                                                                         >  2:599257/1‑151 (MQ=255)                                         gCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCcgcg                                                                                                    >  1:224655/1‑151 (MQ=255)                                          ccATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCgcgc                                                                                                   >  1:209804/1‑151 (MQ=255)                                          ccATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCgcgc                                                                                                   >  2:843561/1‑151 (MQ=255)                                           cATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGcc                                                                                                  >  2:894592/1‑151 (MQ=255)                                                tGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCttt                                                                                              <  2:209945/150‑1 (MQ=255)                                                        gggTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTACGGt                                                                                     >  2:15448/1‑151 (MQ=255)                                                          gTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTc                                                                                    >  1:199870/1‑150 (MQ=255)                                                          gTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCa                                                                                   <  1:844098/151‑1 (MQ=255)                                                                  caACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTc                                                                            >  2:216185/1‑150 (MQ=255)                                                                   aaCGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAc                                                                          >  1:648481/1‑151 (MQ=255)                                                                     cGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcgc                                                                        >  1:173000/1‑151 (MQ=255)                                                                      gTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcgcg                                                                       >  2:909087/1‑151 (MQ=255)                                                                        gggCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGAt                                                                     <  1:773092/151‑1 (MQ=255)                                                                            cAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTtct                                                                 >  2:450328/1‑151 (MQ=255)                                                                             aGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTtctc                                                                <  1:562858/151‑1 (MQ=255)                                                                                 gTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCtcatca                                                            <  2:601423/151‑1 (MQ=255)                                                                                      gCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAAt                                                         <  1:765930/149‑1 (MQ=255)                                                                                        ccATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTg                                                     <  1:263326/151‑1 (MQ=255)                                                                                         cATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCg                                                                                                               >  2:243586/1‑92 (MQ=255)                                                                                         cATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCg                                                                                                               <  1:243612/92‑1 (MQ=255)                                                                                         cATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGc                                                    >  2:360684/1‑151 (MQ=255)                                                                                              cAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTtc                                                >  2:220303/1‑150 (MQ=255)                                                                                                   aTCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATcagt                                          <  1:887105/151‑1 (MQ=255)                                                                                                   aTCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATcagt                                          >  1:158/1‑151 (MQ=255)                                                                                                       cAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATcagttcca                                      <  2:882422/151‑1 (MQ=255)                                                                                                        aGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATcagttcc                                       >  2:243228/1‑149 (MQ=255)                                                                                                            gggTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATcagttccagttc                                  >  2:727764/1‑150 (MQ=255)                                                                                                             ggTTCTGCGGCGTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATCAGTTCCAGTTCAt                                <  1:879739/151‑1 (MQ=255)                                                                                                                   gcggcgTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATCAGTTCCAGTTCATTCGt                            <  1:407496/149‑1 (MQ=255)                                                                                                                    cggcgTAACGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATCAGTTCCAGTTCATTCGtctc                         >  2:568559/1‑151 (MQ=255)                                                                                                                           aCGCACCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATCAGTTCCAGTTCATTCGTCTCCTGACCg                  >  2:668184/1‑151 (MQ=255)                                                                                                                               aCCAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATCAGTTCCAGTTCATTCGTCTCCTGACCgttgt              <  1:44635/151‑1 (MQ=255)                                                                                                                                ccAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATCAGTTCCAGTTCATTCGTCTCCTGACCgttgtt             >  2:666063/1‑151 (MQ=255)                                                                                                                                 cAGGCGAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATCAGTTCCAGTTCATTCGTCTCCTGACCGttgttg            >  2:615353/1‑151 (MQ=255)                                                                                                                                      gAAAACGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATCAGTTCCAGTTCATTCGTCTCCTGACCGTTGTTGAACgg       >  1:125665/1‑151 (MQ=255)                                                                                                                                        aaaCGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATCAGTTCCAGTTCATTCGTCTCCTGACCGTTGTTGAACGGCa     <  2:567446/151‑1 (MQ=255)                                                                                                                                           cGCCGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATCAGTTCCAGTTCATTCGTCTCCTGACCGTTGTTGAACGGCAAAc  <  2:770954/151‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                    CTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATCAAATGTTGCGTTCATCAGTTCCAGTTCATTCGTCTCCTGACCGTTGTTGAACGGCAAAC  >  NC_000913/2581684‑2581972

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