Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	951,454	C→T	E701K (GAA→AAA)	pflB ←	pyruvate formate‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	951,454	0	C	T	100.0%	174.7 / NA	50	E701K (GAA→AAA)	pflB	pyruvate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/27);  total (23/27)

TCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGT  >  NC_000913/951325‑951579                                                                                                                                  |                                                                                                                              tCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCt                                                                                                                     >  1:279363/1‑140 (MQ=255)  cAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCtt                                                                                                                    <  1:196790/140‑1 (MQ=255)   aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTc                                                                                                                   >  2:634658/1‑140 (MQ=255)     cGAGTTGAAACGTACTGTGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCtt                                                                                                                    >  2:425352/1‑137 (MQ=255)     cGAGTTGAAACGTACTGTGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCtt                                                                                                                    <  1:425352/137‑1 (MQ=255)     cGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCgt                                                                                                                 >  2:423284/1‑140 (MQ=255)     cGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCgt                                                                                                                 >  1:419958/1‑140 (MQ=255)          tGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGtggtgg                                                                                                            <  1:634658/140‑1 (MQ=255)                gTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTaa                                                                                                      <  1:159008/140‑1 (MQ=255)                  aCTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAAcc                                                                                                    >  2:9073/1‑140 (MQ=255)                  aCTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAAcc                                                                                                    <  2:632961/140‑1 (MQ=255)                   cTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAAcca                                                                                                   <  1:217966/140‑1 (MQ=255)                      cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAAccatcc                                                                                                >  1:173404/1‑140 (MQ=255)                       gTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAAccatcca                                                                                               <  2:721967/140‑1 (MQ=255)                        tAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAAccatcca                                                                                               >  1:486634/1‑139 (MQ=255)                        tAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAAccatcca                                                                                               <  2:486634/139‑1 (MQ=255)                            cAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGa                                                                                          <  1:423284/140‑1 (MQ=255)                                 acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagc                                                                                     <  1:9073/140‑1 (MQ=255)                                     ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAgg                                                                                 <  2:737090/140‑1 (MQ=255)                                         ggTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTgg                                                                             >  2:346694/1‑140 (MQ=255)                                         ggTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTgg                                                                             >  1:596679/1‑140 (MQ=255)                                          gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag                                                                                      <  2:827029/130‑1 (MQ=255)                                          gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag                                                                                      >  1:827029/1‑130 (MQ=255)                                          gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGt                                                                            >  2:711307/1‑140 (MQ=255)                                                  cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGaa                                                                    >  2:367137/1‑140 (MQ=255)                                                    gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACt                                                                  >  1:863704/1‑140 (MQ=255)                                                         ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcgtc                                                             >  2:189332/1‑140 (MQ=255)                                                            tCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcg                                                               >  1:374750/1‑135 (MQ=255)                                                            tCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcg                                                               <  2:374750/135‑1 (MQ=255)                                                                  ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAcc                                                    >  1:370610/1‑140 (MQ=255)                                                                           gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGtt                                           >  2:291188/1‑140 (MQ=255)                                                                            cgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTc                                          >  2:237945/1‑140 (MQ=255)                                                                              tCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                                        <  2:596679/140‑1 (MQ=255)                                                                                 agcagcATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAAc                                     <  2:419958/140‑1 (MQ=255)                                                                                         ttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAgaga                             <  1:291188/140‑1 (MQ=255)                                                                                         ttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAgaga                             <  2:5078/140‑1 (MQ=255)                                                                                         ttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAgaga                             <  1:237945/140‑1 (MQ=255)                                                                                              ggTTCATCACGTTAACGTTCAGGTGCTGTCCAGCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtg                        <  1:346694/140‑1 (MQ=255)                                                                                               gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtgt                       <  2:863704/140‑1 (MQ=255)                                                                                                ttCATCACGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTa                      >  1:758466/1‑140 (MQ=255)                                                                                                   atcaCGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGga                   >  1:624371/1‑140 (MQ=255)                                                                                                    tcaCGTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGgag                  <  2:173404/140‑1 (MQ=255)                                                                                                        gTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGCCCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTTCCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATAc              <  2:305387/140‑1 (MQ=255)                                                                                                        gTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATAc              <  1:632082/140‑1 (MQ=255)                                                                                                        gTTAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATAc              <  2:279363/140‑1 (MQ=255)                                                                                                         ttAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATAcc             <  2:758466/140‑1 (MQ=255)                                                                                                          tAACGTTCAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCa            >  2:496034/1‑140 (MQ=255)                                                                                                                 cAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAACCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGAt                <  2:155576/129‑1 (MQ=255)                                                                                                                 cAGGTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAACCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGAt                >  1:155576/1‑129 (MQ=255)                                                                                                                    gTGCTGACCACCTTTGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGt  <  2:370610/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                              TCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGT  >  NC_000913/951325‑951579

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