Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,897,050	+G	coding (277/666 nt)	yieH →	phosphoenolpyruvate and 6‑phosphogluconate phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,897,046	1	.	G	100.0%	239.6 / NA	66	E91E (GAA→GAG)	yieH	phosphoenolpyruvate and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (31/35);  total (31/35)

GTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAA‑GGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGG  >  NC_000913/3896926‑3897177                                                                                                                          |                                                                                                                                    gTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCt                                                                                                                       <  1:827331/136‑1 (MQ=255) gTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCt                                                                                                                       >  2:827331/1‑136 (MQ=255)  tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTg                                                                                                                   >  2:194802/1‑139 (MQ=255)  tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTg                                                                                                                   <  1:831166/139‑1 (MQ=255)  tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTg                                                                                                                   >  1:115009/1‑139 (MQ=255)            tCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTg                                                                                                                   <  1:1024193/129‑1 (MQ=255)            tCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                         >  2:1428787/1‑139 (MQ=255)            tCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCACGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTg                                                                                                                   >  2:1024193/1‑129 (MQ=255)                 atTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTg                                                                                                                   >  1:221666/1‑124 (MQ=255)                 atTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTg                                                                                                                   <  2:221666/124‑1 (MQ=255)                 atTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTgcg                                                                                                    >  2:1332202/1‑139 (MQ=255)                 atTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTgcg                                                                                                    >  1:1346224/1‑139 (MQ=255)                       tCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGTGATCACTGCGCCAAtg                                                                                              <  1:250541/139‑1 (MQ=255)                        cccTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgt                                                                                             <  2:258835/139‑1 (MQ=255)                          cTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgt                                                                                           >  2:1016225/1‑139 (MQ=255)                          cTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgt                                                                                           <  2:1346224/139‑1 (MQ=255)                           ttGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          <  1:194802/139‑1 (MQ=255)                           ttGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          <  1:1016225/139‑1 (MQ=255)                               aCATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCg                                                                                                             >  1:1246244/1‑116 (MQ=255)                               aCATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCg                                                                                                             <  2:1246244/116‑1 (MQ=255)                               aCATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGt                                                                                      >  2:641639/1‑139 (MQ=255)                                   ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                      <  2:132588/87‑1 (MQ=255)                                   ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                      >  1:132588/1‑87 (MQ=255)                                   ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                                                 <  2:837701/124‑1 (MQ=255)                                   ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                                                 >  1:837701/1‑124 (MQ=255)                                   ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          <  2:1047745/131‑1 (MQ=255)                                   ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          >  1:1047745/1‑131 (MQ=255)                                    gtgtTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                         <  2:1203118/115‑1 (MQ=255)                                    gtgtTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                         >  1:1203118/1‑115 (MQ=255)                                       ttACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTg                                                                                                                   >  1:1100265/1‑102 (MQ=255)                                       ttACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTg                                                                                                                   <  2:1100265/102‑1 (MQ=255)                                        tACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTaa                                                                                >  1:99202/1‑136 (MQ=255)                                        tACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTaa                                                                                <  2:99202/136‑1 (MQ=255)                                            ttAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCa                                                                         <  2:115009/139‑1 (MQ=255)                                               gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGCAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                         >  2:504012/1‑104 (MQ=255)                                               gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                         <  1:504012/104‑1 (MQ=255)                                               gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                                                 >  2:674832/1‑112 (MQ=255)                                               gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                                                 <  1:674832/112‑1 (MQ=255)                                               gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAAt                                                                      >  1:1281962/1‑139 (MQ=255)                                                       aGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGAt                                                                                                           <  1:700906/94‑1 (MQ=255)                                                       aGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGAt                                                                                                           >  2:700906/1‑94 (MQ=255)                                                         aaGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGTTGGAGCGCTcc                                                                                                                     <  1:279044/82‑1 (MQ=255)                                                         aaGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGTTGGAGCGCTcc                                                                                                                     >  2:279044/1‑82 (MQ=255)                                                         aaGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgc                                                            >  1:692592/1‑139 (MQ=255)                                                                       tACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCt                                                                                                                       >  1:213266/1‑66 (MQ=255)                                                                       tACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCt                                                                                                                       <  2:213266/66‑1 (MQ=255)                                                                       tACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                              <  1:1332202/139‑1 (MQ=255)                                                                         ccGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCaa                                            <  1:1428787/139‑1 (MQ=255)                                                                         ccGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAACCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCaa                                            <  2:119876/139‑1 (MQ=255)                                                                              cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                         <  1:1327994/73‑1 (MQ=255)                                                                              cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                         >  2:1327994/1‑73 (MQ=255)                                                                              cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          >  1:319650/1‑88 (MQ=255)                                                                              cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          <  2:319650/88‑1 (MQ=255)                                                                              cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGa                                       >  1:1029921/1‑139 (MQ=255)                                                                                   gTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATg                                  <  1:489375/139‑1 (MQ=255)                                                                                         cGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCac                            <  1:641639/139‑1 (MQ=255)                                                                                          ggCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCact                           <  2:692592/139‑1 (MQ=255)                                                                                               ttCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          >  1:164538/1‑71 (MQ=255)                                                                                               ttCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          <  2:164538/71‑1 (MQ=255)                                                                                               ttCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACa                                                                  <  2:765426/95‑1 (MQ=255)                                                                                               ttCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACa                                                                  >  1:765426/1‑95 (MQ=255)                                                                                               ttCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTc                      >  1:305629/1‑139 (MQ=255)                                                                                                   aTTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCgg                  <  2:1281962/139‑1 (MQ=255)                                                                                                     tCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGAt                <  2:305629/139‑1 (MQ=255)                                                                                                         aaCTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAAc            >  1:36615/1‑139 (MQ=255)                                                                                                              gAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAt                                                                                               >  2:196742/1‑51 (MQ=255)                                                                                                              gAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAt                                                                                               <  1:196742/51‑1 (MQ=255)                                                                                                                   cATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCgg  <  2:1029921/139‑1 (MQ=255)                                                                                                                          |                                                                                                                                    GTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAA‑GGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGG  >  NC_000913/3896926‑3897177

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