Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,897,439	G→A	*222* (TAG→TAA)	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,439	0	G	A	100.0%	206.5 / NA	65	*222* (TAG→TAA)	yieH	phosphoenolpyruvate and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/28);  total (37/28)

ACGCAGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGT  >  NC_000913/3897303‑3897573                                                                                                                                         |                                                                                                                                       aCGCAGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAAtt                                                                                                                                      >  2:1041031/1‑139 (MQ=255)                  ttCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCAc                                                                                                                     >  2:296660/1‑139 (MQ=255)                    cTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACtt                                                                                                                   >  2:849416/1‑139 (MQ=255)                      aCTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTAc                                                                                                                 <  2:415251/139‑1 (MQ=255)                          cTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCg                                                                                                             >  1:459062/1‑139 (MQ=255)                                       cacaATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGgcgc                                                                                                <  1:411048/139‑1 (MQ=255)                                          aaTAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCa                                                                                             <  1:934993/139‑1 (MQ=255)                                            tAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTccc                                                                                                                          <  2:770922/108‑1 (MQ=255)                                            tAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTccc                                                                                                                          >  1:770922/1‑108 (MQ=255)                                              aGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAgg                                                                                         >  2:358804/1‑139 (MQ=255)                                                ccGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGgta                                                                                       <  1:296660/139‑1 (MQ=255)                                                   aTCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTg                                                                                                    >  1:57569/1‑123 (MQ=255)                                                   aTCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTg                                                                                                    <  2:57569/123‑1 (MQ=255)                                                    tCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCa                                                                                             >  1:223214/1‑129 (MQ=255)                                                    tCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCa                                                                                             <  2:223214/129‑1 (MQ=255)                                                    tCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGgtaagta                                                                                   >  2:1082144/1‑139 (MQ=255)                                                      gTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTaaa                                                                                 >  1:131055/1‑139 (MQ=255)                                                       ttCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAg                                                                                <  1:137481/139‑1 (MQ=255)                                                       ttCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAg                                                                                >  2:884260/1‑139 (MQ=255)                                                           cccGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAg                                                                            <  1:884260/139‑1 (MQ=255)                                                                aaGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGa                                                                       >  1:85018/1‑139 (MQ=255)                                                                aaGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGa                                                                       <  1:358804/139‑1 (MQ=255)                                                                 aGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCt                                                                                               <  1:471747/114‑1 (MQ=255)                                                                 aGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCt                                                                                               >  2:471747/1‑114 (MQ=255)                                                                     accaccTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGt                                                                  <  2:1093018/139‑1 (MQ=255)                                                                      ccaccTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTc                                                                 >  1:771381/1‑139 (MQ=255)                                                                               cccATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACg                                                        <  2:131055/139‑1 (MQ=255)                                                                               cccATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACg                                                        >  1:736048/1‑139 (MQ=255)                                                                                    cTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTAc                                                                                                                 <  1:652683/77‑1 (MQ=255)                                                                                    cTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTAc                                                                                                                 >  2:652683/1‑77 (MQ=255)                                                                                    cTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCggg                                                   >  2:951414/1‑139 (MQ=255)                                                                                    cTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCggg                                                   >  1:30503/1‑139 (MQ=255)                                                                                    cTTTCGCAGTTACCTGAACTGTGGAAAGCGCGCGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCggg                                                   >  2:494080/1‑139 (MQ=255)                                                                                       tCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTaa                                                >  2:106927/1‑139 (MQ=255)                                                                                        cGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAAt                                               >  1:1049207/1‑139 (MQ=255)                                                                                           aggtaactgcaCTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACa                                            >  2:191797/11‑139 (MQ=255)                                                                                             gtaacTGCACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATc                                          <  1:191797/135‑1 (MQ=255)                                                                                                ccTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCt                                                                                               >  1:19015/1‑83 (MQ=255)                                                                                                ccTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCt                                                                                               <  2:19015/83‑1 (MQ=255)                                                                                                       tgtgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGgtaagta                                                                                   <  2:320930/88‑1 (MQ=255)                                                                                                       tgtgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGgtaagta                                                                                   >  1:320930/1‑88 (MQ=255)                                                                                                       tgtgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTaaa                                                                                 <  1:181121/90‑1 (MQ=255)                                                                                                       tgtgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTaaa                                                                                 >  2:181121/1‑90 (MQ=255)                                                                                                         tgGAAAGCGCGTGTTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGt                                                                  <  2:1026371/103‑1 (MQ=255)                                                                                                         tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGt                                                                  >  1:1026371/1‑103 (MQ=255)                                                                                                         tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTAtt                              >  2:630120/1‑139 (MQ=255)                                                                                                         tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTAtt                              >  2:108221/1‑139 (MQ=255)                                                                                                         tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTAtt                              >  2:631229/1‑139 (MQ=255)                                                                                                         tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTAtt                              >  1:244791/1‑139 (MQ=255)                                                                                                          ggAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGat                                                                      <  2:1043761/98‑1 (MQ=255)                                                                                                          ggAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGat                                                                      >  1:1043761/1‑98 (MQ=255)                                                                                                            aaaGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTg                           <  1:631229/139‑1 (MQ=255)                                                                                                             aaGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGc                          <  1:630120/139‑1 (MQ=255)                                                                                                              aGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGgcgc                                                                                                <  2:694828/68‑1 (MQ=255)                                                                                                              aGCGCGTGGTTGGGATAATACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGgcgt                                                                                                >  1:694828/1‑67 (MQ=255)                                                                                                                cgcgTGGTTGGGATATTACGGCATAATTCTTCACACTCGCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAGGGGAGTTTGATATGTCTGTTTCACg                                                        >  2:549372/1‑106 (MQ=255)                                                                                                                cgcgTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACg                                                        <  1:549372/106‑1 (MQ=255)                                                                                                                 gcgTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTaaa                                                                                 >  2:533029/1‑80 (MQ=255)                                                                                                                 gcgTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTaaa                                                                                 <  1:533029/80‑1 (MQ=255)                                                                                                                   gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGtttt                    >  2:5313/1‑139 (MQ=255)                                                                                                                                tACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAAtt       >  2:367520/1‑139 (MQ=255)                                                                                                                                 aCGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGat                                                                      >  2:1011822/1‑75 (MQ=255)                                                                                                                                 aCGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGat                                                                      <  1:1011822/75‑1 (MQ=255)                                                                                                                                  cGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTgg     <  1:108221/139‑1 (MQ=255)                                                                                                                                     cATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGgtgt  <  1:367520/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       ACGCAGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGT  >  NC_000913/3897303‑3897573

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