Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,085,818	C→T	A32T (GCC→ACC)	speA ←	arginine decarboxylase, biosynthetic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,085,818	0	C	T	100.0%	175.9 / NA	53	A32T (GCC→ACC)	speA	arginine decarboxylase, biosynthetic
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/25);  total (28/25)

CTGCGCGAGATCGACGCGAGCTTCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTA  >  NC_000913/3085690‑3085950                                                                                                                                 |                                                                                                                                     cTGCGCGAGATCGACGCGAGCTGCCTGGACGTCCGGGGCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGc                                                                                                                           <  1:221818/140‑1 (MQ=255)    cgcgAGATCGACGCGAGCTTCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCa                                                                                                                        >  1:80435/1‑140 (MQ=255)    cgcgAGATCGACGCGAGCTTCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCa                                                                                                                        >  2:591127/1‑140 (MQ=255)    cgcgAGATCGACGCGAGCTTCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCa                                                                                                                        >  2:334144/1‑140 (MQ=255)    cgcgAGATCGACGCGAGCTTCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCa                                                                                                                        >  1:492225/1‑140 (MQ=255)    cgcgAGATCGACGCGAGCTTCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTATAAGTACCCAGTATGTTGCTGTTTTCCTGGGAGCTCa                                                                                                                        >  1:234457/1‑140 (MQ=255)         gaTCGACGCGAGCTCCCGGGACGTCCGGGTCCGGGCCCCCGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCa                                                                                                                   <  2:667637/140‑1 (MQ=255)           tCGACGCGAGCTTCCGGGACGTCCGGGTCCGGGGACACGCTAATGTGGCCCAGCTCGTTAGCGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCACCTTGCTGGTTTCCTGGGAGCTCATTGCAAc                                                                                                                 <  1:642545/140‑1 (MQ=255)            cGACGCGAGCTTCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAAcc                                                                                                                >  1:67500/1‑140 (MQ=255)              aCGCGAGCTTCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAAcctc                                                                                                              >  1:482545/1‑140 (MQ=255)                     cTTCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATg                                                                                                       >  1:618925/1‑140 (MQ=255)                     cTTCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATg                                                                                                       <  2:67500/140‑1 (MQ=255)                       tCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGa                                                                                                     <  2:618925/140‑1 (MQ=255)                       tCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGa                                                                                                     >  2:543946/1‑140 (MQ=255)                       tCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGa                                                                                                     <  2:511767/140‑1 (MQ=255)                       tCCGGGACGTCCGGGTCAGGGCGCACGCTAATGTGGCCCTGCTCGTTATCGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGa                                                                                                     <  2:243610/140‑1 (MQ=255)                          gggACGTCCGGGTCCAGGCACACGCTAATGTGGCCCAACTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCg                                                                                                  <  1:228292/140‑1 (MQ=255)                           ggACGTCCGGGTCCGCGCACACCCTTATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCgt                                                                                                 <  1:170693/140‑1 (MQ=255)                               gtccgggtccggGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCATGGGAGCGCATGGCGACCTCCAGCATGGAGCgtagta                                                                                             >  1:322933/1‑140 (MQ=255)                                                    aaTGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACg                                                                        <  1:334144/140‑1 (MQ=255)                                                          gCCCAGCTGGTTAACGGCATAGTAGTTATTTCCGCTCCATGCAATATTATAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCa                                                                  <  2:322933/140‑1 (MQ=255)                                                          gCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCa                                                                  <  2:571964/140‑1 (MQ=255)                                                           cccAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCaa                                                                 <  2:80435/140‑1 (MQ=255)                                                           cccAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCa                                                                  <  1:41403/139‑1 (MQ=255)                                                           cccAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCa                                                                  >  2:41403/1‑139 (MQ=255)                                                               gCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAAccc                                                             >  2:527804/1‑140 (MQ=255)                                                                cTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCa                                                            >  2:156797/1‑140 (MQ=255)                                                                     tAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCCCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGAc                                                       <  2:431304/140‑1 (MQ=255)                                                                      aaCGTCATAGTAGTTATGGCCCCACCAGGCAATATTGTATGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACa                                                      <  2:234457/140‑1 (MQ=255)                                                                            aTAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATgtcgt                                                <  2:482545/140‑1 (MQ=255)                                                                             tagtagTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATgtcgtc                                               >  2:46841/1‑140 (MQ=255)                                                                                     aTTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAg                                       >  1:383889/1‑140 (MQ=255)                                                                                     aTTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAg                                       >  2:449551/1‑140 (MQ=255)                                                                                     aTTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCGGACGAAGGCAAGCCAATAGACATGTCGTCAGACATAg                                       >  1:234954/1‑140 (MQ=255)                                                                                          cccACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAAc                                  >  2:125415/1‑140 (MQ=255)                                                                                            caccaGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCt                                >  2:400296/1‑140 (MQ=255)                                                                                             accaGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTc                               >  2:384077/1‑140 (MQ=255)                                                                                                aGGCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCaaa                            <  2:72589/140‑1 (MQ=255)                                                                                                  gCAATATTGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAatt                          <  1:400296/140‑1 (MQ=255)                                                                                                         tGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAattat                        >  2:35580/1‑135 (MQ=255)                                                                                                         tGTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAattat                        <  1:35580/135‑1 (MQ=255)                                                                                                          gTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGCCATGTCGTCAGCCATAGCGAACCTCAAAttattttatt                  <  2:492225/140‑1 (MQ=255)                                                                                                          gTAAGTACGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAAttattttatt                  <  1:384077/140‑1 (MQ=255)                                                                                                                aCGCAGCATCTTGCTGGTTTCCTGGTAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAgtgt            >  1:560301/1‑140 (MQ=255)                                                                                                                aCGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGTAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAgtgt            >  2:616924/1‑140 (MQ=255)                                                                                                                aCGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAgtgt            >  2:448102/1‑140 (MQ=255)                                                                                                                aCGCAGCATCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAgtgt            >  1:335451/1‑140 (MQ=255)                                                                                                                  gcagcaTCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTaa          >  2:398039/1‑140 (MQ=255)                                                                                                                  gcagcaTCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTaa          >  1:555228/1‑140 (MQ=255)                                                                                                                    cagaTCTTGCTGGTTTCCTGGGTGCTCATTGAAACTTCCCGCATGGAGAGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTaaaa        <  2:214174/137‑1 (MQ=255)                                                                                                                     gcaTCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTaa                <  1:354522/131‑1 (MQ=255)                                                                                                                     gcaTCTTGCTGGTTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTaa                >  2:354522/1‑131 (MQ=255)                                                                                                                          ttGCTGGTTTCCTGGGAACTCACTGCAACCTCCTGTATGCAGCGTAGTCCACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTa  <  2:176385/140‑1 (MQ=255)                                                                                                                                 |                                                                                                                                     CTGCGCGAGATCGACGCGAGCTTCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTA  >  NC_000913/3085690‑3085950

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