Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,658,893	G→A	*176* (TAG→TAA)	gadE →	gad regulon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,658,893	0	G	A	100.0%	169.0 / NA	52	*176* (TAG→TAA)	gadE	gad regulon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/28);  total (24/28)

ACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAGTCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCAC  >  NC_000913/3658780‑3659018                                                                                                                  |                                                                                                                              aCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGa                                                                                                      >  1:713423/1‑139 (MQ=255)     tCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtggc                                                                                                  >  1:920445/1‑139 (MQ=255)     tCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtggc                                                                                                  <  2:624192/139‑1 (MQ=255)              cGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGc                                                                                                                  <  2:938048/114‑1 (MQ=255)              cGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGc                                                                                                                  >  1:938048/1‑114 (MQ=255)              cGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGAc                                                                                         <  2:493265/139‑1 (MQ=255)               gCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                        <  1:272737/139‑1 (MQ=255)                caacaaACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGa                                                                                       <  1:844225/139‑1 (MQ=255)                        aTTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCa                                                                               <  2:713423/139‑1 (MQ=255)                                 tCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGa                                                                             >  2:898140/1‑132 (MQ=255)                                 tCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGa                                                                             <  1:898140/132‑1 (MQ=255)                                                 aTGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTg                                                      <  1:7857/139‑1 (MQ=255)                                                   gTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                    <  2:755899/139‑1 (MQ=255)                                                   gTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTCAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                    <  2:817221/139‑1 (MQ=255)                                                        aGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                        <  2:815197/98‑1 (MQ=255)                                                        aGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                        >  1:815197/1‑98 (MQ=255)                                                         gCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATg                                                                           >  2:860389/1‑110 (MQ=255)                                                         gCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATg                                                                           <  1:860389/110‑1 (MQ=255)                                                         gCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                    <  2:161365/133‑1 (MQ=255)                                                         gCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                    >  1:161365/1‑133 (MQ=255)                                                         gCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTaa                                              >  1:818082/1‑139 (MQ=255)                                                              aaCTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGgag                                         <  2:642642/139‑1 (MQ=255)                                                                       gTATGAGCGACATCGTCACACTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTaaaa                                >  1:785059/1‑139 (MQ=255)                                                                         aTGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAg                                                                                      <  1:345751/83‑1 (MQ=255)                                                                         aTGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAg                                                                                      >  2:345751/1‑83 (MQ=255)                                                                         aTGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACg                              >  2:483065/1‑139 (MQ=255)                                                                          tGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGAACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                        >  1:253241/1‑80 (MQ=255)                                                                          tGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGAACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                        <  2:253241/80‑1 (MQ=255)                                                                            aGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAAGTGAAATTGATGGCTGGCATGACg                                                                                        >  2:708548/1‑78 (MQ=255)                                                                            aGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAATTTGATGGCTGGCATGACg                                                                                        <  1:708548/78‑1 (MQ=255)                                                                             gCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGga                                          >  1:407690/1‑123 (MQ=255)                                                                             gCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGga                                          <  2:407690/123‑1 (MQ=255)                                                                               gACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                        <  2:785059/139‑1 (MQ=255)                                                                                          ccTGGGTATCACATCTTATTTTTAATCAGGACATAAGc                                                                                                                  >  1:591776/1‑38 (MQ=38)                                                                                          ccTGGGTATCACATCTTATTTTTAATCAGGACATAAGc                                                                                                                  <  2:591776/38‑1 (MQ=38)                                                                                              ggTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCa                                                                               <  1:182782/69‑1 (MQ=255)                                                                                              ggTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCa                                                                               >  2:182782/1‑69 (MQ=255)                                                                                               gTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                    <  1:737985/95‑1 (MQ=255)                                                                                               gTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                    >  2:737985/1‑95 (MQ=255)                                                                                                tATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                    <  2:444616/94‑1 (MQ=255)                                                                                                tATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                    >  1:444616/1‑94 (MQ=255)                                                                                                   cacaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATg                                                                           >  2:346880/1‑68 (MQ=255)                                                                                                   cacaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATg                                                                           <  1:346880/68‑1 (MQ=255)                                                                                                   cacaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                        >  2:14793/1‑119 (MQ=255)                                                                                                   cacaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                        <  1:14793/119‑1 (MQ=255)                                                                                                   cacaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTcccc                  <  2:654304/125‑1 (MQ=255)                                                                                                   cacaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTcccc                  >  1:654304/1‑125 (MQ=255)                                                                                                   cacaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATcc    >  2:453413/1‑139 (MQ=255)                                                                                                     caTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATccac  <  2:920445/139‑1 (MQ=255)                                                                                                        cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTAcccc                                                           <  1:14225/79‑1 (MQ=255)                                                                                                        cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTAcccc                                                           >  2:14225/1‑79 (MQ=255)                                                                                                              tttAATCAGGACATAGGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTg                                                      >  2:717641/1‑78 (MQ=255)                                                                                                                  |                                                                                                                              ACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAGTCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCAC  >  NC_000913/3658780‑3659018

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