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%	200.5 / NA	62	*176* (TAG→TAA)	gadE	gad regulon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/30);  total (32/30)

AAATGAAAGCAATATCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAGTCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCACCAG  >  NC_000913/3658758‑3659021                                                                                                                                        |                                                                                                                                 aaatGAAAGCAATATCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCa                                                                                                                               <  1:399217/139‑1 (MQ=255)             tatCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCaa                                                                                                                   >  1:299291/1‑139 (MQ=255)                               tatTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtggctggc                                                                                                 <  2:299291/139‑1 (MQ=255)                                               ttAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAg                                                                                 <  2:218062/139‑1 (MQ=255)                                                                cAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTAcc                                                                >  2:116312/1‑139 (MQ=255)                                                                      tATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGt                                                          <  2:777597/139‑1 (MQ=255)                                                                              aGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATg                                                                              >  2:93354/1‑111 (MQ=255)                                                                              aGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATg                                                                              <  1:93354/111‑1 (MQ=255)                                                                              aGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTAcccc                                                              <  2:729220/127‑1 (MQ=255)                                                                              aGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTAcccc                                                              >  1:729220/1‑127 (MQ=255)                                                                               gCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTaa                                                 >  2:390758/1‑139 (MQ=255)                                                                               gCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTaa                                                 >  1:719282/1‑139 (MQ=255)                                                                                cTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGc                                                                                   >  1:219940/1‑104 (MQ=255)                                                                                cTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGc                                                                                   <  2:219940/104‑1 (MQ=255)                                                                                  gAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAgg                                              >  1:666347/1‑139 (MQ=255)                                                                                        aaGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTg                                                         <  2:154982/122‑1 (MQ=255)                                                                                        aaGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGc                                        <  2:719282/139‑1 (MQ=255)                                                                                        aaGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCACCGGTg                                                         >  1:154982/1‑122 (MQ=255)                                                                                         aGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAAtt                                                                                                           <  1:333539/71‑1 (MQ=255)                                                                                         aGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAAtt                                                                                                           >  2:333539/1‑71 (MQ=255)                                                                                                tGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGa                                >  1:173751/1‑139 (MQ=255)                                                                                                tGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGa                                >  2:385530/1‑139 (MQ=255)                                                                                                tGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGa                                >  1:622837/1‑139 (MQ=255)                                                                                                      aCATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTAt                          >  2:341410/1‑139 (MQ=255)                                                                                                       cATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTAtt                         <  1:344847/139‑1 (MQ=255)                                                                                                        aTCGTCACTCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTc                        >  2:45519/1‑139 (MQ=255)                                                                                                        aTCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTc                        >  2:396711/1‑139 (MQ=255)                                                                                                        aTCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTc                        >  2:810110/1‑139 (MQ=255)                                                                                                        aTCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTc                        >  2:318973/1‑139 (MQ=255)                                                                                                         tCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTcc                       <  1:311505/139‑1 (MQ=255)                                                                                                            tCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCt                    <  1:810110/139‑1 (MQ=255)                                                                                                            tCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCt                    <  2:294309/139‑1 (MQ=255)                                                                                                            tCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCt                    <  2:666347/139‑1 (MQ=255)                                                                                                               cccTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGt                 <  1:385530/139‑1 (MQ=255)                                                                                                                ccTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTa                <  1:346944/139‑1 (MQ=255)                                                                                                                 cTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTAt               <  2:173751/139‑1 (MQ=255)                                                                                                                    ggTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCAt                                                                    >  1:347270/1‑83 (MQ=255)                                                                                                                    ggTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCAt                                                                    <  2:347270/83‑1 (MQ=255)                                                                                                                         cacaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTa                >  1:502254/1‑130 (MQ=255)                                                                                                                         cacaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTa                <  2:502254/130‑1 (MQ=255)                                                                                                                         cacaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATcc       <  1:493684/139‑1 (MQ=255)                                                                                                                         cacaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATcc       >  1:428512/1‑139 (MQ=255)                                                                                                                         cacaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATcc       >  2:493684/1‑139 (MQ=255)                                                                                                                          acaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTg                   >  1:507660/1‑126 (MQ=255)                                                                                                                          acaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTg                   <  2:507660/126‑1 (MQ=255)                                                                                                                              cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGt                                                          <  1:402537/83‑1 (MQ=255)                                                                                                                              cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGt                                                          >  2:402537/1‑83 (MQ=255)                                                                                                                              cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACt                                                   >  1:421995/1‑90 (MQ=255)                                                                                                                              cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACt                                                   <  2:421995/90‑1 (MQ=255)                                                                                                                              cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTccc                      <  1:472890/119‑1 (MQ=255)                                                                                                                              cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTccc                      >  2:472890/1‑119 (MQ=255)                                                                                                                              cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCACCAg  >  1:604582/1‑139 (MQ=255)                                                                                                                                tATTTTTAATCAGGACATAAGCAACTGAAATTGAtgg                                                                                                      <  1:434128/37‑1 (MQ=255)                                                                                                                                tATTTTTAATCAGGACATAAGCAACTGAAATTGAtgg                                                                                                      >  2:434128/1‑37 (MQ=255)                                                                                                                                tATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                           >  2:55342/1‑48 (MQ=255)                                                                                                                                tATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                           <  1:55342/48‑1 (MQ=255)                                                                                                                                tATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTAtt                         <  2:737525/114‑1 (MQ=255)                                                                                                                                tATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTAtt                         >  1:737525/1‑114 (MQ=255)                                                                                                                                tATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTa          >  2:458595/1‑129 (MQ=255)                                                                                                                                tATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTa          <  1:458595/129‑1 (MQ=255)                                                                                                                                 aTTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                           <  2:249177/111‑1 (MQ=255)                                                                                                                                 aTTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                           >  1:249177/1‑111 (MQ=255)                                                                                                                                        |                                                                                                                                 AAATGAAAGCAATATCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAGTCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCACCAG  >  NC_000913/3658758‑3659021

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