Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,902,543	G→A	G216G (GGG→GGA)	ygcE →	putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,902,543	0	G	A	100.0%	203.5 / NA	64	G216G (GGG→GGA)	ygcE	putative kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/31);  total (33/31)

TTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCAAACATTG  >  NC_000913/2902418‑2902675                                                                                                                              |                                                                                                                                     tttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGGGACACATTTCCCGGCAGGACTGCCGGCTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTgg                                                                                                                         >  2:1661581/1‑139 (MQ=255) tttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTgg                                                                                                                         >  2:122610/1‑139 (MQ=255)             aTGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAAcgg                                                                                                             >  1:695136/1‑139 (MQ=255)                    gCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTa                                                                                                                             <  2:126788/116‑1 (MQ=255)                    gCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTa                                                                                                                             >  1:126788/1‑116 (MQ=255)                         gTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTtctc                                                                                                 <  2:406491/139‑1 (MQ=255)                           ccTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                                                               <  1:1154426/139‑1 (MQ=255)                                    tatCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTgg                                                                                             <  2:1093434/132‑1 (MQ=255)                                    tatCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTgg                                                                                             >  1:1093434/1‑132 (MQ=255)                                           ccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCg                                                                                  >  2:1075099/1‑136 (MQ=255)                                           ccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTAGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGTATTACTGGATCATGAAACGGCGGTAATTTCTCTTGGCACTTACATCg                                                                                  <  1:1075099/136‑1 (MQ=255)                                              caAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTtc                                                                                                   <  1:940044/116‑1 (MQ=255)                                              caAGCCGCACTGGCGACAAATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTtc                                                                                                   >  2:940044/1‑116 (MQ=255)                                                      aCTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAAc                                                                    >  1:106137/1‑139 (MQ=255)                                                             acacATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAg                                                             >  1:662075/1‑139 (MQ=255)                                                             acacATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAg                                                             <  1:122610/139‑1 (MQ=255)                                                               acaTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCa                                                           <  2:308416/139‑1 (MQ=255)                                                                    tCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGAt                                                                                                                               <  1:629471/66‑1 (MQ=255)                                                                    tCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGAt                                                                                                                               >  2:629471/1‑66 (MQ=255)                                                                    tCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATGGCATTGATGATGAACGGCAAAGCACTGcc                                                      >  2:1091183/1‑139 (MQ=255)                                                                     cccGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCg                                                     <  1:1066049/139‑1 (MQ=255)                                                                         gCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACAt                                                                                    <  1:1505312/104‑1 (MQ=255)                                                                         gCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACAt                                                                                    >  2:1505312/1‑104 (MQ=255)                                                                            ggACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATc                                              <  1:1091183/139‑1 (MQ=255)                                                                             gACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATcc                                             >  1:1627955/1‑139 (MQ=255)                                                                                 gCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTa                                         <  2:662075/139‑1 (MQ=255)                                                                                        gtttgtACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACt                                  >  1:113541/1‑139 (MQ=255)                                                                                           tgtACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTa                                                                                                                             <  1:1598144/45‑1 (MQ=255)                                                                                           tgtACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTa                                                                                                                             >  2:1598144/1‑45 (MQ=255)                                                                                           tgtACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGc                               <  2:113541/139‑1 (MQ=255)                                                                                              accaccAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTa                                         <  2:793691/126‑1 (MQ=255)                                                                                              accaccAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTa                                         >  1:793691/1‑126 (MQ=255)                                                                                                caccaGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGga                                                                                                                        >  1:657491/1‑45 (MQ=255)                                                                                                caccaGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGga                                                                                                                        <  2:657491/45‑1 (MQ=255)                                                                                                  ccaGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaa                                                   >  1:390572/1‑112 (MQ=255)                                                                                                  ccaGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaa                                                   <  2:390572/112‑1 (MQ=255)                                                                                                   caGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTg                                 <  2:1577259/129‑1 (MQ=255)                                                                                                   caGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTg                                 >  1:1577259/1‑129 (MQ=255)                                                                                                     gTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACAt                                                                                    >  2:1232246/1‑76 (MQ=255)                                                                                                     gTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACAt                                                                                    <  1:1232246/76‑1 (MQ=255)                                                                                                     gTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCa                                                                >  1:1184232/1‑96 (MQ=255)                                                                                                     gTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCa                                                                <  2:1184232/96‑1 (MQ=255)                                                                                                     gTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGtc                     >  1:1564075/1‑139 (MQ=255)                                                                                                        aTAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATc                                              <  2:1251125/111‑1 (MQ=255)                                                                                                        aTAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATc                                              >  1:1251125/1‑111 (MQ=255)                                                                                                           aaCCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGaa                                                                     <  1:1281271/85‑1 (MQ=255)                                                                                                           aaCCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGaa                                                                     >  2:1281271/1‑85 (MQ=255)                                                                                                            aCCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGaa                                                                     >  2:1624840/1‑84 (MQ=255)                                                                                                            aCCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGaa                                                                     <  1:1624840/84‑1 (MQ=255)                                                                                                             ccGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCa                                                                                <  1:1075182/72‑1 (MQ=255)                                                                                                             ccGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCa                                                                                >  2:1075182/1‑72 (MQ=255)                                                                                                                gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGaaa                                                                                                                >  2:940591/1‑37 (MQ=255)                                                                                                                gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGaaa                                                                                                                <  1:940591/37‑1 (MQ=255)                                                                                                                gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGc          <  2:695136/139‑1 (MQ=255)                                                                                                                 tAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCa         >  1:47560/1‑139 (MQ=255)                                                                                                                       ctctGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                                                               <  1:674922/47‑1 (MQ=255)                                                                                                                       ctctGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                                                               >  2:674922/1‑47 (MQ=255)                                                                                                                       ctctGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaa                                                   >  2:1468777/1‑91 (MQ=255)                                                                                                                       ctctGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaa                                                   <  1:1468777/91‑1 (MQ=255)                                                                                                                        tctGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCAAACATTg  >  2:522733/1‑139 (MQ=255)                                                                                                                          tGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCAt                                                                               <  2:1244767/60‑1 (MQ=255)                                                                                                                          tGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCAt                                                                               >  1:1244767/1‑60 (MQ=255)                                                                                                                          tGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCAAACATTg  >  2:1416230/1‑137 (MQ=255)                                                                                                                          tGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCAAACATTg  <  1:1416230/137‑1 (MQ=255)                                                                                                                           gggagCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCAAACATTg  >  2:1165626/1‑136 (MQ=255)                                                                                                                              |                                                                                                                                     TTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCAAACATTG  >  NC_000913/2902418‑2902675

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: ‑