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

CCAGATAGACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGAACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGCACGTTTGAAC  >  NC_000913/3085326‑3085594                                                                                                                                        |                                                                                                                                      ccAGATAGACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGAcc                                                                                                                                    >  2:896012/1‑139 (MQ=255) ccAGATAGACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGAcc                                                                                                                                    <  1:484885/139‑1 (MQ=255)          ccTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATGAACTCGGCTTTGGGACCGGCTTCCAg                                                                                                                           >  2:591824/1‑139 (MQ=255)              gtgCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGAccc                                                                                                                       <  2:710776/139‑1 (MQ=255)                        ttctCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGc                                                                                                             >  1:855042/1‑139 (MQ=255)                               cAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCgaatg                                                                                                      >  2:301809/1‑139 (MQ=255)                                  ttAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCgaatgaat                                                                                                   >  1:626972/1‑139 (MQ=255)                                                 tatatTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCg                                                                                                                                   >  2:405061/1‑92 (MQ=255)                                                 tatatTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCg                                                                                                                                   <  1:405061/92‑1 (MQ=255)                                                     ttCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCCGCTTTGGGACCGGCTACCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                                <  1:591824/139‑1 (MQ=255)                                                        gcgGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGcggcgg                                                                             >  1:266378/1‑139 (MQ=255)                                                         cgGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGAcc                                                                                                                                    >  1:458899/1‑83 (MQ=255)                                                         cgGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGAcc                                                                                                                                    <  2:458899/83‑1 (MQ=255)                                                          ggTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCaa                                                                                        <  2:845750/126‑1 (MQ=255)                                                          ggTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCaa                                                                                        >  1:845750/1‑126 (MQ=255)                                                                tataACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAg                                                                                                <  1:395390/112‑1 (MQ=255)                                                                tataACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAg                                                                                                >  2:395390/1‑112 (MQ=255)                                                                        ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGa                                                             >  2:492202/1‑139 (MQ=255)                                                                         tGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGAt                                                            <  2:626972/139‑1 (MQ=255)                                                                              aCGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGAt                                                       >  2:674700/1‑139 (MQ=255)                                                                               cGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATa                                                      <  1:896012/139‑1 (MQ=255)                                                                                   gACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCgg                                                                                                                  >  1:334085/1‑75 (MQ=255)                                                                                   gACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCgg                                                                                                                  <  2:334085/75‑1 (MQ=255)                                                                                        tACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACa                                                  >  1:282013/1‑134 (MQ=255)                                                                                        tACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACa                                                  <  2:282013/134‑1 (MQ=255)                                                                                          cGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGt                                                                      >  2:487080/1‑112 (MQ=255)                                                                                          cGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGt                                                                      <  1:487080/112‑1 (MQ=255)                                                                                                aTGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGc                                     >  1:760454/1‑139 (MQ=255)                                                                                                     aGCATGTGCCAGTACTGCCATTAACTCGGCTTTGGGACCGGCTTCCAGACACATCGGTTCGCCCGAATGAATCAGGGACGCAAGCACGCGGCGGTGCTGGTTAACTTTGTTCGGATAAACAAGGAAGTAAGCGCCGTta                                <  2:237161/139‑1 (MQ=255)                                                                                                     aGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTta                                <  2:760454/139‑1 (MQ=255)                                                                                                      gCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTtat                               <  2:866179/139‑1 (MQ=255)                                                                                                       cATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTtata                              >  1:818804/1‑139 (MQ=255)                                                                                                         tgtgCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGc                            <  1:674700/139‑1 (MQ=255)                                                                                                           tgCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCg                          >  2:130180/1‑139 (MQ=255)                                                                                                                  aCTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGAtt                   >  2:100247/1‑139 (MQ=255)                                                                                                                      ccATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCt               >  1:750507/1‑139 (MQ=255)                                                                                                                                  cTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCgaatg                                                                                                      >  2:379450/1‑40 (MQ=255)                                                                                                                                  cTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCgaatg                                                                                                      <  1:379450/40‑1 (MQ=255)                                                                                                                                  cTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGCACGTTTGaa   >  1:452939/1‑139 (MQ=255)                                                                                                                                   tttGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGCACGTTTGAAc  <  2:750507/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      CCAGATAGACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGAACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGCACGTTTGAAC  >  NC_000913/3085326‑3085594

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