Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,827,351	G→A	*320* (TAG→TAA)	srlE →	glucitol/sorbitol‑specific enzyme IIB component of PTS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,827,351	0	G	A	100.0%	179.0 / NA	55	*320* (TAG→TAA)	srlE	glucitol/sorbitol‑specific enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/32);  total (23/32)

GTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTA  >  NC_000913/2827217‑2827477                                                                                                                                        |                                                                                                                               gtcggtTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGc                                                                                                                            >  1:102089/1‑139 (MQ=255) gtcggtTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGc                                                                                                                            <  1:566812/139‑1 (MQ=255)            gCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                                 <  2:36358/139‑1 (MQ=255)            gCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                                 <  1:654566/139‑1 (MQ=255)                     aGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGttattta                                                                                                        <  2:1545117/139‑1 (MQ=255)                          gTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGa                                                                                                   <  2:442981/139‑1 (MQ=255)                          gTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGa                                                                                                   <  2:1093603/139‑1 (MQ=255)                               gACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGaccacc                                                                                              <  2:102089/139‑1 (MQ=255)                               gACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGaccacc                                                                                              <  1:1454721/139‑1 (MQ=255)                                acacGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGAccacca                                                                                             <  2:1720613/139‑1 (MQ=255)                                    ggTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCt                                                                                                                           <  2:195370/105‑1 (MQ=255)                                    ggTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCt                                                                                                                           >  1:195370/1‑105 (MQ=255)                                              ggTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGgc                                                                               >  2:1209437/1‑139 (MQ=255)                                                       tCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGc                                                                                                                            <  1:1065611/86‑1 (MQ=255)                                                       tCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGc                                                                                                                            >  2:1065611/1‑86 (MQ=255)                                                            aCTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGa                                                                 >  2:501619/1‑139 (MQ=255)                                                              tGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGc                                                                                                                            <  2:1263900/78‑1 (MQ=255)                                                              tGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGc                                                                                                                            >  1:1263900/1‑78 (MQ=255)                                                              tGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATATATCAGACCACCATCACCCGTATCGGCGCGAg                                                                          >  2:1124401/1‑128 (MQ=255)                                                              aggTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAg                                                                          <  1:1124401/127‑1 (MQ=255)                                                                gTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGc                                                                                                                            >  1:189894/1‑76 (MQ=255)                                                                gTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGc                                                                                                                            <  2:189894/76‑1 (MQ=255)                                                                gTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGcc                                                             >  2:1723756/1‑139 (MQ=255)                                                                 tGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGc                                                              >  1:685015/1‑137 (MQ=255)                                                                 tGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGc                                                              <  2:685015/137‑1 (MQ=255)                                                                  ggtttgttttttATCCGGCTCGCCGACTGTAC‑TGATCGCCTGGTGTGTCTCCGGTGTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCt                                                           <  2:16363/132‑1 (MQ=255)                                                                  gAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCt                                                           <  2:1415188/139‑1 (MQ=255)                                                                     ccGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCt                                                           <  1:475572/136‑1 (MQ=255)                                                                     ccGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCt                                                           >  2:475572/1‑136 (MQ=255)                                                                     ccGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAg                                                        >  2:284504/1‑139 (MQ=255)                                                                             tAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGa                                                <  1:284504/139‑1 (MQ=255)                                                                             tAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGa                                                <  1:1209437/139‑1 (MQ=255)                                                                                    cgcgcCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCtcat                                         <  1:501619/139‑1 (MQ=255)                                                                                                 c‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAg                                                 <  2:556161/118‑1 (MQ=255)                                                                                                 c‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAg                                                 >  1:556161/1‑118 (MQ=255)                                                                                                     aTCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACa                                                                                                                     >  2:1573714/1‑47 (MQ=255)                                                                                                     aTCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACa                                                                                                                     <  1:1573714/47‑1 (MQ=255)                                                                                                      tCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGgcg                        <  1:1382956/139‑1 (MQ=255)                                                                                                       cGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGaccac                                                                                               <  2:1253992/67‑1 (MQ=255)                                                                                                       cGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGaccac                                                                                               >  1:1253992/1‑67 (MQ=255)                                                                                                          cTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGa                                                                           >  2:931102/1‑84 (MQ=255)                                                                                                          cTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGa                                                                           <  1:931102/84‑1 (MQ=255)                                                                                                          cTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGcc                    >  1:1301549/1‑139 (MQ=255)                                                                                                                  tctcCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGAcc            <  2:1383185/139‑1 (MQ=255)                                                                                                                    tcCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTc          >  2:1370219/1‑139 (MQ=255)                                                                                                                         ttttATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGac                                                                                                  <  2:1619880/46‑1 (MQ=255)                                                                                                                         ttttATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGac                                                                                                  >  1:1619880/1‑46 (MQ=255)                                                                                                                          tttATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAagag    <  2:1301549/139‑1 (MQ=255)                                                                                                                            tatctatcAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTa  >  1:843540/1‑139 (MQ=255)                                                                                                                                   cAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGc                                                              <  2:1035748/72‑1 (MQ=255)                                                                                                                                   cAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGc                                                              >  2:1529962/1‑72 (MQ=255)                                                                                                                                   cAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGc                                                              <  2:1022827/72‑1 (MQ=255)                                                                                                                                   cAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGc                                                              <  1:1529962/72‑1 (MQ=255)                                                                                                                                   cAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGc                                                              >  1:1035748/1‑72 (MQ=255)                                                                                                                                   cAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGc                                                              >  1:1022827/1‑72 (MQ=255)                                                                                                                                        |                                                                                                                               GTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAC‑TGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTA  >  NC_000913/2827217‑2827477

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