Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,085,461	A→G	S151P (TCC→CCC)	speA ←	biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,085,461	0	A	G	100.0%	222.7 / NA	63	S151P (TCC→CCC)	speA	biosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/30);  total (33/30)

CAGATAGACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGAACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGCACGT  >  NC_000913/3085327‑3085588                                                                                                                                       |                                                                                                                                cAGATAGACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCg                                                                                                                             <  1:454971/139‑1 (MQ=255)  aGATAGACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCgg                                                                                                                            <  2:328685/139‑1 (MQ=255)           ttGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGAc                                                                                                                   >  2:105803/1‑139 (MQ=255)           ttGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGAc                                                                                                                   >  2:847100/1‑139 (MQ=255)             gtgCCCCATCTTCTCGCCAATTAATGTCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGAccc                                                                                                                 >  2:21633/1‑139 (MQ=255)             gtgCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGAccc                                                                                                                 >  2:1265112/1‑139 (MQ=255)              tgCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCa                                                                                                                >  2:1042008/1‑139 (MQ=255)              tgCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCa                                                                                                                >  1:1154706/1‑139 (MQ=255)              tgCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCa                                                                                                                >  1:704650/1‑139 (MQ=255)              tgCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCa                                                                                                                >  2:701099/1‑139 (MQ=255)               gCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTcc                                                                                                                       <  2:183591/131‑1 (MQ=255)               gCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTcc                                                                                                                       >  1:183591/1‑131 (MQ=255)                             ccAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCCGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCgaat                                                                                                 >  1:1074770/1‑139 (MQ=255)                                aTTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCgaatgaa                                                                                              >  1:888025/1‑139 (MQ=255)                                 ttAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCgg                                                                                                            <  2:352431/124‑1 (MQ=255)                                 ttAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCgg                                                                                                            >  1:352431/1‑124 (MQ=255)                                  tAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATc                                                                                            <  1:1042008/139‑1 (MQ=255)                                  tAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATc                                                                                            <  2:1074770/139‑1 (MQ=255)                                   aaTGTCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCa                                                                                           <  1:21633/139‑1 (MQ=255)                                          ggCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATc                                                                                            <  1:414956/131‑1 (MQ=255)                                          ggCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATc                                                                                            >  2:414956/1‑131 (MQ=255)                                          ggCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTc                                                                                    >  2:1248549/1‑139 (MQ=255)                                                    ttCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                          <  1:701099/139‑1 (MQ=255)                                                        cgGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGt                                                                      <  1:105803/139‑1 (MQ=255)                                                              ttATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGt                                                                <  1:1265112/139‑1 (MQ=255)                                                              ttATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGt                                                                <  1:111661/139‑1 (MQ=255)                                                                     cGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGCCCGGCTTCCAGACCCAGCgg                                                                                                            >  2:836493/1‑88 (MQ=255)                                                                     cGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCgg                                                                                                            <  1:836493/88‑1 (MQ=255)                                                                     cGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCa                                                                                           >  1:190326/1‑105 (MQ=255)                                                                     cGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCa                                                                                           <  2:190326/105‑1 (MQ=255)                                                                       ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCaa                                                                                  >  2:1029258/1‑112 (MQ=255)                                                                       ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCaa                                                                                  <  1:1029258/112‑1 (MQ=255)                                                                       ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGt                                                                >  1:1404350/1‑130 (MQ=255)                                                                       ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGt                                                                <  2:1404350/130‑1 (MQ=255)                                                                       ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGa                                                       >  1:756213/1‑139 (MQ=255)                                                                       ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGa                                                       <  2:756213/139‑1 (MQ=255)                                                                       ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGa                                                       >  1:1376811/1‑139 (MQ=255)                                                                        tGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGAt                                                      <  1:847100/139‑1 (MQ=255)                                                                        tGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGAt                                                      <  2:1154706/139‑1 (MQ=255)                                                                        tGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGAt                                                      <  1:1248549/139‑1 (MQ=255)                                                                        tGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGAt                                                      <  2:704650/139‑1 (MQ=255)                                                                         gCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                          >  1:1031387/1‑118 (MQ=255)                                                                         gCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                          <  2:1031387/118‑1 (MQ=255)                                                                            gACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCACTCAcg                                                                             >  2:238340/1‑112 (MQ=255)                                                                            gACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCAcg                                                                             <  1:238340/112‑1 (MQ=255)                                                                                  gACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACa                                            >  1:1268103/1‑139 (MQ=255)                                                                                         cGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGt                                     >  2:523165/1‑139 (MQ=255)                                                                                         cGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGt                                     >  2:1869/1‑139 (MQ=255)                                                                                            gTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGt                                                                >  2:937695/1‑109 (MQ=255)                                                                                            gTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGt                                                                <  1:937695/109‑1 (MQ=255)                                                                                               aTGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAg                                                                                          >  1:722997/1‑80 (MQ=255)                                                                                               aTGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAg                                                                                          <  2:722997/80‑1 (MQ=255)                                                                                                     gCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTtat                         <  2:888025/139‑1 (MQ=255)                                                                                                      cATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTtata                        >  2:810422/1‑139 (MQ=255)                                                                                                          tgCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGATATAGCCg                    >  2:202538/1‑139 (MQ=255)                                                                                                                tACTGGCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGt                                                                                                           <  2:526440/46‑1 (MQ=255)                                                                                                                tACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGt                                                                                                           >  1:526440/1‑46 (MQ=255)                                                                                                                tACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGAt              <  2:1268103/139‑1 (MQ=255)                                                                                                                    gCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTccc          >  2:376991/1‑139 (MQ=255)                                                                                                                            cTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGCACGt  >  1:414223/1‑139 (MQ=255)                                                                                                                               ggCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGa                                                       >  1:1330694/1‑83 (MQ=255)                                                                                                                               ggCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGa                                                       <  2:1330694/83‑1 (MQ=255)                                                                                                                                    tGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                          <  2:563539/59‑1 (MQ=255)                                                                                                                                       |                                                                                                                                CAGATAGACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGAACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGCACGT  >  NC_000913/3085327‑3085588

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