Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	761,133	G→T	D810Y (GAT→TAT)	sucA →	2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	761,133	0	G	T	100.0%	156.9 / NA	49	D810Y (GAT→TAT)	sucA	2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (22/27);  total (22/27)

TCAGGCGCTGCGCGGGATGCGTCGTCCGCTGGTCGTGATGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTGATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAACACGATGTCGCCATTGTGCGTATCGAGCAACTCTACCCGTTCCCGCATAAA  >  NC_000913_3_bme_pgi/760994‑761267                                                                                                                                            |                                                                                                                                       tcAGGCGCTGCGCGGGATGCGTCGTCCGCTGGTCGTGATGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAgg                                                                                                                               >  1:186831/1‑149 (MQ=255)   aGGCGCTGCGCGGGATGCGTCGTCCGCTGGTCGTGATGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGc                                                                                                                             <  2:646010/149‑1 (MQ=255)      cgcTGCGCGGGATGCGTCGTCCGCTGGTCGTGATGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTg                                                                                                                          <  2:644879/149‑1 (MQ=255)      cgcTGCGCGGGATGCGTCGTCCGCTGGTCGTGATGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTg                                                                                                                          <  1:376334/149‑1 (MQ=255)        cTGTGCGGGATGCGTCGTCCGCTGGGCGTGATGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGaa                                                                                                                        <  1:68361/149‑1 (MQ=255)        cTGCGCGGGATGCGTCGTCCGCTGGTCGTGATGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGaa                                                                                                                        <  1:213021/149‑1 (MQ=255)            gcgGGATGCGTCGTCCACTGGTCGTGATGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTGTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAgcgc                                                                                                                    >  1:601198/1‑149 (MQ=255)                 aTGCGTCGTCCGCCGGTCGTGATGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGt                                                                                                               <  1:321403/149‑1 (MQ=255)                  tGCGTCGTCCGCTGGTCGTGATGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTg                                                                                                              <  2:279049/149‑1 (MQ=255)                  tGCGTCGTCCGCTGGTCGTGACGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTg                                                                                                              <  2:278628/149‑1 (MQ=255)                      tcgtcCGCTGGTCGTGATGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGAtgt                                                                                                          >  2:576045/1‑149 (MQ=255)                           cGCTGGTCGTGATGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCt                                                                                                     >  1:235992/1‑149 (MQ=255)                            gCTGGTCGTGATGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCATCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTg                                                                                                    <  1:145947/149‑1 (MQ=255)                             cTGGTCGGGATGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTgg                                                                                                   <  1:295768/149‑1 (MQ=255)                                     gATGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGttt                                                                                           >  2:312667/1‑149 (MQ=255)                                           gCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACg                                                                                     >  2:402414/1‑149 (MQ=255)                                                aaTCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATTCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACctg                                                                                <  1:457142/149‑1 (MQ=255)                                                    ccTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCGAGCCATCGGTGAAATCGACGAGATTTATCCGAAGGGCATGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTgg                                                                            >  2:624782/1‑149 (MQ=255)                                                    ccTGCTGCGTCATCCGCTGGCGGTATCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTgg                                                                            >  1:598694/1‑149 (MQ=255)                                                      tgctgcGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTTACCGAAGGGCGTGAAGCGCGTGGTGTTGTGTTCTGGTAAGGTTTATTGGGACCTGCTggca                                                                          >  1:98117/1‑147 (MQ=255)                                                       gctgcGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAAc                                                                         >  2:460614/1‑149 (MQ=255)                                                         tgcGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAg                                                                       <  2:443637/149‑1 (MQ=255)                                                          gcGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCAGAAGGGAATGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGc                                                                      <  1:481227/149‑1 (MQ=255)                                                               aTCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCAGCCAGCCATCGGTGAAAACGACGAGCATTATCCGAAGGGCGCGAAACGCGTAGTGATGTGATCTGGTCACGTTTCTTCCGACATGCCGCAACAGcgtcgt                                                                 >  2:421733/1‑149 (MQ=255)                                                                tCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGcgtcgaa                                                                >  1:424574/1‑147 (MQ=255)                                                                  cGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGa                                                                                    >  2:60283/1‑127 (MQ=255)                                                                  cGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGa                                                                                    <  1:60283/127‑1 (MQ=255)                                                                    cTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGaa                                                            <  1:402414/149‑1 (MQ=255)                                                                    cTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGaa                                                            <  1:498118/149‑1 (MQ=255)                                                                    cTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACCGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGCACACCATCGTAACaa                                                            <  2:181362/149‑1 (MQ=255)                                                                               cAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAACACg                                                 >  1:521283/1‑149 (MQ=255)                                                                                   cTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGACATAGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAACACGATGt                                             >  2:323869/1‑149 (MQ=255)                                                                                      gaagaaCTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAATACgttatagc                                          >  2:518996/1‑142 (MQ=255)                                                                                          aaCTGGCGAACGGCACCTTCCTTCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAACACGATGTCGCCAtt                                      <  1:424708/149‑1 (MQ=255)                                                                                            cTGCCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCTACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAACACGATGTCGCCATtgt                                    <  2:364833/149‑1 (MQ=255)                                                                                             tGGCGAACGGCACCTTCCTGCCAGCAATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAAAAGAGTCGCAAGAACAATCAACACGATGTCGCCAGtgtg                                   >  2:605978/1‑149 (MQ=255)                                                                                              tgctaaCGGCATCTTCCTTCCAGCCATGGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGTAACAGCGTCGTAAGAACAATCAACACGATGTCGCCATTGTGc                                  <  2:539291/145‑1 (MQ=255)                                                                                              ggCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAACACGATGTCGCCATTGTGc                                  >  2:41824/1‑149 (MQ=255)                                                                                                cGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAACACGATGTCGCCATTGTGCGt                                <  1:460614/149‑1 (MQ=255)                                                                                                  aaCGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAACACGATGTCGCCATTGTGCGTAt                              <  1:312667/149‑1 (MQ=255)                                                                                                     ggCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAACACGATGTCGCCATTGTGCGTATCGa                           <  2:9099/149‑1 (MQ=255)                                                                                                          cTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAACACGATGTCGCCATTGTGCGTATCGAGCAAc                      >  1:55850/1‑149 (MQ=255)                                                                                                          cTTCCTGCCAGCCAGCGGGGAAATCGACGAGCTTTATGCGAAGGTCTGGAAGCGCTTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAATAATCAACACGATGTCGCCATTGTGCGTATCGAGCAAc                      <  2:521283/149‑1 (MQ=255)                                                                                                               tGCCAGCCATCGGTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAACACGATGTCGCCATTGTGCGTATCGAGCAACTCTAc                 >  1:427848/1‑149 (MQ=255)                                                                                                                        tCGGGGAGATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTAGTAAGAACAATCAACACGATGTCGCCATTGTGCGTATCGAGCAACTCTACCCGTTCCCg        <  2:202870/149‑1 (MQ=255)                                                                                                                          ggTGAAATCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAACACGATGTCGCCATTGTGCGTATAGAGCAACTCTACCCGTTCCCGCa      >  1:657066/1‑149 (MQ=255)                                                                                                                          ggTGAAAGCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAACACGATGTCGCCATTGTGCGTACCGAGCAACTCTACCCGTTCCCGCa      <  2:448455/149‑1 (MQ=255)                                                                                                                           gTGAAATCGACGAGCTTTATACGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAACACGATGTCGCCATTGTGCGTATCGAACAACTCTACCCGTTCCCGCAt     >  1:593488/1‑149 (MQ=255)                                                                                                                              aaaTCGACGAGCTTTATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAACACGATGTCGCCATTGTGCGTATCGAGCAACTCTACCCGTTCCCGCATaaa  <  1:41824/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                       TCAGGCGCTGCGCGGGATGCGTCGTCCGCTGGTCGTGATGTCGCCGAAATCCCTGCTGCGTCATCCGCTGGCGGTTTCCAGCCTCGAAGAACTGGCGAACGGCACCTTCCTGCCAGCCATCGGTGAAATCGACGAGCTTGATCCGAAGGGCGTGAAGCGCGTAGTGATGTGTTCTGGTAAGGTTTATTACGACCTGCTGGAACAGCGTCGTAAGAACAATCAACACGATGTCGCCATTGTGCGTATCGAGCAACTCTACCCGTTCCCGCATAAA  >  NC_000913_3_bme_pgi/760994‑761267

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑