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%	130.8 / NA	36	S151P (TCC→CCC)	speA	biosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/16);  total (20/16)

ACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGAACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCG  >  NC_000913/3085334‑3085583                                                                                                                                |                                                                                                                           aCCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCa                                                                                                                 >  2:216803/1‑139 (MQ=255)      gtgCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGAccc                                                                                                            >  1:802608/1‑139 (MQ=255)       tgCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCa                                                                                                           >  1:237690/1‑139 (MQ=255)                          ttAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCgaatgaat                                                                                        >  1:118642/1‑139 (MQ=255)                                    gCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGATCCAGCg                                                                                                        >  2:245270/1‑113 (MQ=255)                                    gCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCg                                                                                                        <  1:245270/113‑1 (MQ=255)                                                  ggTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTg                                                                <  1:730800/139‑1 (MQ=255)                                                   gTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCga                                                                                              <  2:737154/108‑1 (MQ=255)                                                   gTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCga                                                                                              >  1:737154/1‑108 (MQ=255)                                                     cTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCTATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTg                                                             >  1:218326/1‑139 (MQ=255)                                                                ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGa                                                  <  1:290325/139‑1 (MQ=255)                                                                  gCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATc                                                >  1:206192/1‑139 (MQ=255)                                                                                      tCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATc                            <  2:206192/139‑1 (MQ=255)                                                                                       cATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCg                           >  2:56803/1‑139 (MQ=255)                                                                                           ccAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGt                       >  2:401772/1‑139 (MQ=255)                                                                                                  gtgCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAAGCAGGGACTAAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGcc                <  2:237690/139‑1 (MQ=255)                                                                                                   tgCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCTCCCgaat                                                                                            <  2:550113/62‑1 (MQ=255)                                                                                                   tgCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCgaat                                                                                            >  1:550113/1‑62 (MQ=255)                                                                                                   tgCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATc                            >  2:717568/1‑126 (MQ=255)                                                                                                   tgCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATc                            <  1:717568/126‑1 (MQ=255)                                                                                                    gCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCgg                                              <  2:564063/107‑1 (MQ=255)                                                                                                    gCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCgg                                              >  1:564063/1‑107 (MQ=255)                                                                                                        gTACTGCTATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGa          <  2:218326/139‑1 (MQ=255)                                                                                                         tACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGt                                                           >  1:837103/1‑89 (MQ=255)                                                                                                         tACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGt                                                           <  2:837103/89‑1 (MQ=255)                                                                                                         tACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGAt         <  1:56803/139‑1 (MQ=255)                                                                                                           cTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTc       >  2:104330/1‑139 (MQ=255)                                                                                                             gCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAAc                                                       <  1:647197/89‑1 (MQ=255)                                                                                                             gCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAAc                                                       >  2:647197/1‑89 (MQ=255)                                                                                                              ccATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGa                                                  >  2:155169/1‑93 (MQ=255)                                                                                                              ccATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGa                                                  <  1:155169/93‑1 (MQ=255)                                                                                                              ccATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATc                                                <  2:313185/95‑1 (MQ=255)                                                                                                              ccATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATc                                                >  1:313185/1‑95 (MQ=255)                                                                                                              ccATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCt    >  1:754164/1‑139 (MQ=255)                                                                                                                     cTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCg  >  2:310413/1‑134 (MQ=255)                                                                                                                     cTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCg  <  1:310413/134‑1 (MQ=255)                                                                                                                                |                                                                                                                           ACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGAACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCG  >  NC_000913/3085334‑3085583

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