Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,594,520	A→G	S531P (TCC→CCC)	ydeK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,594,520	0	A	G	100.0%	270.8 / NA	86	S531P (TCC→CCC)	ydeK	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (43/43);  total (43/43)

ATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAT  >  NC_000913/1594430‑1594609                                                                                           |                                                                                          aTACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGgggg                                                                                         >  1:3518099/1‑93 (MQ=255)   acacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCa                                                                                       <  1:2332557/93‑1 (MQ=255)    cacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCAt                                                                                      <  2:2170749/93‑1 (MQ=255)     acCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCAtt                                                                                     <  1:412390/93‑1 (MQ=255)      ccGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCAttt                                                                                    >  1:768373/1‑93 (MQ=255)        gACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTtcc                                                                                  >  1:1000454/1‑93 (MQ=255)          ccTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTtcctc                                                                                <  2:2894397/93‑1 (MQ=255)          ccTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTtcctc                                                                                >  2:2000027/1‑93 (MQ=255)            tGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCg                                                                              <  2:1759946/93‑1 (MQ=255)             gTAGCAATTGTTCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGa                                                                             <  2:2351069/93‑1 (MQ=255)             gTAGCAATTGTGCGTTAGGAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGa                                                                             <  1:1894441/93‑1 (MQ=255)               aGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATa                                                                           >  2:1579017/1‑93 (MQ=255)                gCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATaa                                                                          >  2:3842007/1‑93 (MQ=255)                 cAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATaaa                                                                         <  1:4543454/93‑1 (MQ=255)                  aaTTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAAt                                                                        >  2:2339795/1‑93 (MQ=255)                  aaTTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAAt                                                                        <  2:4183694/93‑1 (MQ=255)                    ttGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATcc                                                                      >  2:3852263/1‑93 (MQ=255)                      gtgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCAt                                                                    >  2:3661773/1‑93 (MQ=255)                           ttAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGccc                                                               <  1:2248427/93‑1 (MQ=255)                              gtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGtt                                                            >  1:1664361/1‑93 (MQ=255)                                 gatgatGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGaa                                                         >  2:3936550/1‑93 (MQ=255)                                 gatgatGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGaa                                                         >  2:3288407/1‑93 (MQ=255)                                 gatgatGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGaa                                                         >  1:1660071/1‑93 (MQ=255)                                 gatgatGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGaa                                                         >  2:4198333/1‑93 (MQ=255)                                    gatgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTa                                                      <  1:2650474/93‑1 (MQ=255)                                     atgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTag                                                     >  1:4533000/1‑92 (MQ=255)                                     atgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTaa                                                     <  2:795438/93‑1 (MQ=255)                                     atgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTaa                                                     <  2:3853916/93‑1 (MQ=255)                                      tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                    >  1:4236227/1‑93 (MQ=255)                                      agTCTTTAAATACCAGAGACTGTCCGTTGAAATATAAACTATCCCCTTCCCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                    <  2:4193399/92‑1 (MQ=255)                                        tCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGTTAAATCCATAGCCCGTTGAATTAAggg                                                  <  1:83713/93‑1 (MQ=255)                                        tCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAggg                                                  <  1:2689584/93‑1 (MQ=255)                                         cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGt                                                 <  2:4165925/93‑1 (MQ=255)                                          tttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTg                                                >  2:2110859/1‑93 (MQ=255)                                          tttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTg                                                >  2:3401122/1‑93 (MQ=255)                                          tttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTg                                                <  1:4641343/93‑1 (MQ=255)                                           ttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAgggggg                                               >  2:4385328/1‑93 (MQ=255)                                           ttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTgg                                               <  2:3982961/93‑1 (MQ=255)                                            tAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGc                                              <  1:4538296/93‑1 (MQ=255)                                            tAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGc                                              <  2:4462970/93‑1 (MQ=255)                                            tAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGc                                              <  1:4480705/93‑1 (MQ=255)                                              aaTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGcgc                                            >  2:1936737/1‑91 (MQ=255)                                              aaTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGcgc                                            <  2:2739655/93‑3 (MQ=255)                                               aTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGcc                                           >  2:43319/1‑93 (MQ=255)                                               aTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGcc                                           >  1:191015/1‑93 (MQ=255)                                                ttCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCa                                          <  2:1746635/93‑1 (MQ=255)                                                  ccAGAGCCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCAtt                                        <  2:392494/93‑1 (MQ=255)                                                   cAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTa                                       <  1:2260723/93‑1 (MQ=255)                                                   cAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTa                                       <  1:2680223/93‑1 (MQ=255)                                                    agagACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTGGAATTAAGGGTGGCGCCATTat                                      >  1:45645/1‑92 (MQ=255)                                                       gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAGATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                   <  1:2044224/93‑1 (MQ=255)                                                       gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                   >  2:4676393/1‑93 (MQ=255)                                                       gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                   >  1:2974936/1‑93 (MQ=255)                                                        aCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGAc                                  <  1:3962506/93‑1 (MQ=255)                                                         cTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                 <  2:3253978/93‑1 (MQ=255)                                                         cTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                 >  2:489082/1‑93 (MQ=255)                                                         cTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                 >  2:3566373/1‑93 (MQ=255)                                                         cTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                 >  1:1313277/1‑93 (MQ=255)                                                         cTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                 <  1:4067425/93‑1 (MQ=255)                                                          tGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAt                                >  2:3319087/1‑93 (MQ=255)                                                           gTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                               <  1:4373634/93‑1 (MQ=255)                                                            tCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAttt                              >  1:4084283/1‑93 (MQ=255)                                                               gTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACg                           >  1:1913978/1‑93 (MQ=255)                                                               gTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACg                           >  1:4600901/1‑93 (MQ=255)                                                                ttGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                          >  1:3736504/1‑93 (MQ=255)                                                                  gAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATa                        >  1:1138966/1‑93 (MQ=255)                                                                   aaaTATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGAtcc                       <  2:417440/93‑3 (MQ=255)                                                                   aaaTATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATAc                       <  1:1664774/93‑1 (MQ=255)                                                                    aaTATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATAcc                      <  1:2497056/93‑1 (MQ=255)                                                                     atatTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCg                     >  2:3221804/1‑93 (MQ=255)                                                                      tatTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGc                    <  2:3833627/93‑1 (MQ=255)                                                                      tatTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGc                    <  1:206226/93‑1 (MQ=255)                                                                            cTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaa              >  2:688162/1‑93 (MQ=255)                                                                              aTCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaaa            <  2:2603015/93‑1 (MQ=255)                                                                               tCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAt           >  1:3877447/1‑93 (MQ=255)                                                                                ccccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAtt          >  1:4328961/1‑93 (MQ=255)                                                                                ccccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAtt          <  1:4631823/93‑1 (MQ=255)                                                                                 cccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTa         >  2:1067318/1‑93 (MQ=255)                                                                                 cccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTa         >  2:3107321/1‑93 (MQ=255)                                                                                  ccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAc        >  1:3933511/1‑93 (MQ=255)                                                                                   ctgaCCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTTACATTTACGATACCGCTACCAAAATTAcc       <  2:3785434/90‑1 (MQ=255)                                                                                   cTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAcc       <  1:958856/93‑1 (MQ=255)                                                                                   cTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAcc       >  2:3957813/1‑93 (MQ=255)                                                                                   cTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCAATTACTGACATTTACGATACCGCTACCAAAATTAcc       >  1:1988760/1‑93 (MQ=255)                                                                                      aCCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGAc    >  1:639633/1‑93 (MQ=255)                                                                                       ccGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa   <  2:2321563/93‑1 (MQ=255)                                                                                       ccGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa   <  1:1095757/93‑1 (MQ=255)                                                                                        cGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAt  <  1:2565018/93‑1 (MQ=255)                                                                                           |                                                                                          ATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAT  >  NC_000913/1594430‑1594609

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