Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	787,153	A→G	L148P (CTG→CCG)	gpmA ←	2,3‑bisphosphoglycerate‑dependent phosphoglycerate mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	787,153	0	A	G	100.0%	178.5 / NA	50	L148P (CTG→CCG)	gpmA	2,3‑bisphosphoglycerate‑dependent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/27);  total (23/27)

TCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCAGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCTCATCATCTTTAGTCAGTTCCGGCGGAGTCACTGCAAAACCACGACGCCACTGTTTCACCTGCTCGTCGCCATACTT  >  NC_000913/787019‑787280                                                                                                                                       |                                                                                                                                 tCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGTGCCAGGCTTTCCGTCAGCGGCGGttct                                                                                                                             >  2:961639/1‑140 (MQ=255) tCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGttct                                                                                                                             <  1:246216/140‑1 (MQ=255)        cgcACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTTTCGc                                                                                                                      >  1:842523/1‑140 (MQ=255)        cgcACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGc                                                                                                                      >  2:619557/1‑140 (MQ=255)         gcACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCt                                                                                                                     <  2:93222/140‑1 (MQ=255)          cACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑c                                                                                                                   >  1:38084/1‑140 (MQ=255)              tAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGtt                                                                                                               <  1:388000/140‑1 (MQ=255)                   agTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTcgcg                                                                                                          >  1:83089/1‑140 (MQ=255)                      tACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTaa                                                                                                       <  1:330310/140‑1 (MQ=255)                       aCCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTcg                                                                                                            >  2:1274263/1‑134 (MQ=255)                       aCCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTcg                                                                                                            <  1:1274263/134‑1 (MQ=255)                       aCCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTaaa                                                                                                      >  1:218275/1‑139 (MQ=255)                        ccGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGTGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAAcg                                                                                                     <  1:961639/140‑1 (MQ=255)                         cGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAAcgc                                                                                                    >  1:1110317/1‑140 (MQ=255)                             tgCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCCCT‑CAGTTTCGCGGAACGCCGaa                                                                                                >  1:122926/1‑139 (MQ=255)                             tgCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGGTCTTTTTCGCC‑CCGTTTTGCGTAACGCgggg                                                                                                >  2:756769/1‑138 (MQ=255)                              gCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATc                                                                                               <  1:1035507/140‑1 (MQ=255)                               cAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTTTTTTTCGCT‑CCGTTTTTCGTAACGCggggcg                                                                                              >  2:330310/1‑136 (MQ=255)                               cAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCCCT‑CAGTTTTCCCGAACGCGGATCg                                                                                              >  1:832770/1‑140 (MQ=255)                                   gatgatCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGAc                                                                                          <  2:346824/140‑1 (MQ=255)                                   gatgatCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGAc                                                                                          <  2:38084/140‑1 (MQ=255)                                   gatgatCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGAc                                                                                          <  1:619557/140‑1 (MQ=255)                                   gatgatCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGAc                                                                                          <  2:148302/140‑1 (MQ=255)                                      gatCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTTTTTTTCGCTCCAGTTT‑GCGCTACGCGGATTGTGAcccc                                                                                       >  2:902521/1‑139 (MQ=255)                                       atCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTTTTTTTCCCT‑CCCTTTTCCGGAAAGCGGATCGGGAccccg                                                                                      >  1:1277512/1‑138 (MQ=255)                                         cACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGAt                                                                                    <  1:902521/140‑1 (MQ=255)                                                tCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCtc                                                                             <  2:939066/140‑1 (MQ=255)                                                 cACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTTTTTTTCTCT‑CAGTTTCGCGTAACGCGGATCGGGACCCCGATAAAGCtcc                                                                            >  2:381166/1‑139 (MQ=255)                                                  aCCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTTTTTTTCGCT‑CCGTTTCTCGTAACGCGGATCGTGACCCCGGTAAAGcccct                                                                           >  1:840024/1‑136 (MQ=255)                                                  aCCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTTTCGCT‑CCGTTTTTCGTAACGCGGAGCGGGGCCCGGAGAACGCtcct                                                                           >  2:1309022/1‑138 (MQ=255)                                                       tctTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCTCATCATCt                                                                      <  1:640842/140‑1 (MQ=255)                                                               cgcgGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCTCATCATCTTTAGTCAg                                                              <  1:381166/140‑1 (MQ=255)                                                                 cgGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCTCATCATCTTTAGTCAGtt                                                            <  1:530771/140‑1 (MQ=255)                                                                           gTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCTCATCATCTTTAGTCAGTTCCGGCGGAGt                                                  <  2:1110317/140‑1 (MQ=255)                                                                                 ttCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTTTCGCT‑CCGTTTTGCGTAACGCGGGTCGTGACCCGGGTAACGCTCATCATCTTTTGTCAGTTCCGGCGGGGGCACTGc                                            >  2:1139599/1‑140 (MQ=255)                                                                                    cAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCTCATCATCTTTAGTCAGTTCCGGCGGAGTCACTGCaaa                                         <  2:95535/140‑1 (MQ=255)                                                                                      gTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTTTCGCT‑CAGTTTTGCGTAAAGCCGATCGTGACCCCGGTAAAGCTCCTCATCTTTTTTTAGTTCCGGGGGGGGCACTGCaaaaa                                       >  2:233433/1‑139 (MQ=255)                                                                                      gTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCTCATCATCTTTAGTCAGTTCCGGCGGAGTCACTGCaaaaa                                       >  2:722702/1‑139 (MQ=255)                                                                                              tCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCTCATCATCTTTAGTCAGTTCCGGCGGAGTCACTGCAAAACCACGACGc                               <  2:1056579/140‑1 (MQ=255)                                                                                               cACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATTGTGACCCGGATAACGCTCATCATCTTTAGTCAGTTCCGGCGGAGTCACTGCAAAACCACGACGcc                              <  1:188030/140‑1 (MQ=255)                                                                                               cACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCTCATCATCTTTAGTCAGTTCCGGCGGAGTCACTGCAAAACCACGACGcc                              <  2:1102043/140‑1 (MQ=255)                                                                                               cACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCTCATCATCTTTAGTCAGTTCCGGCGGAGTCACTGCAAAACCACGACGcc                              <  2:998165/140‑1 (MQ=255)                                                                                               cACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTCCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCTCATCATCTTTAGTCAGTTCCGGCGGAGTCACTGCAAAACCACGACGcc                              <  2:1297921/140‑1 (MQ=255)                                                                                                 cgcgGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CCGTTTTGCCTAACGCGGGTCGGGACCCGGGTAAAGCTCATCATCTTTTGTCAGTTCCCGCCGGGTCACTGCAAAAACACGACGcccc                            >  2:862788/1‑138 (MQ=255)                                                                                                  gcgGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTTTTTTTCGCT‑CAGTTTTGCGTAAAACGGATCGTGACCCCGGTAAAGCTCATCATCTTTTGTTAGTTCCCGCGGAGGCACTGCAAAAACACCACGCCAcc                           >  2:367856/1‑139 (MQ=255)                                                                                                  gcgGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCTCATCATCTTTAGTCAGTTCCGGCGGAGTCACTGCAAAACCACGACGCCACt                           <  1:298365/140‑1 (MQ=255)                                                                                                   cgGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTTTCGCT‑CCGTTTCGCGGAAAGCGGGGCGTGACCCCGGGAACGCTCCTCATCTTTTGTCAGTTTCCGCGGGGGCACTGCAAAACCCCCACGCCCCtc                          >  2:188677/1‑139 (MQ=255)                                                                                                         aTGGTCAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTTTCTCT‑CAGTTTCGCGGAACGCGGATCGTGACCCCGGTAACGCTCATCATCTTTTGTCAGTTTCGGCGGAGTCACTGCAAAAACACGACGCCACTGTTTCAc                    >  2:1154089/1‑140 (MQ=255)                                                                                                              cAGCGCCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCTCATCATCTTTAGTCAGTTCCGGCGGAGTCACTGCAAAACCACGACGCCACTGTTTCACCTGCt               <  2:194705/140‑1 (MQ=255)                                                                                                                 cgcCAGGCTTTCCGTCAGCGGCGGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCTCATCATCTTTAGTCAGTTCCGGCGGAGTCACTGCAAAACCACGACGCCACTGTTTCACCTGCtcgt            <  2:83089/140‑1 (MQ=255)                                                                                                                      ggCTTTCCGTCAGCGGCGGTTCTTTCCCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCTCATCATCTTTAGTCAGTTCCGGCGGAGTCACTGCAAAACCACGACGCCACTGTTTCACCTGCTCGTCGCCa       <  1:1154089/140‑1 (MQ=255)                                                                                                                           tCCGTCAGCGGCGGGTCTTTTTCGCC‑CAGTTTCGCGTAAAGCGGATCGTGACCCCGGTAACGCTCATCATCTTTTGTCAGTTCCGGCGGGGGCACTGCAAAACCACCACGCCACTGTTTCACCTGCTCGTCGCCCTACtt  >  2:595353/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                 TCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGTAAGGGATCACGCGGTCAATGGTCAGCGCCAGGCTTTCCGTCAGCGGCAGTTCTTTCTCGCT‑CAGTTTCGCGTAACGCGGATCGTGACCCGGATAACGCTCATCATCTTTAGTCAGTTCCGGCGGAGTCACTGCAAAACCACGACGCCACTGTTTCACCTGCTCGTCGCCATACTT  >  NC_000913/787019‑787280

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