Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,979,943	C→T	A308V (GCT→GTT)  *231* (TAG→TAA)	lysR → ygeA ←	transcriptional activator of lysA, autorepressor Asp/Glu_racemase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,979,943	0	C	T	100.0%	162.3 / NA	51	A308V (GCT→GTT)  *231* (TAG→TAA)	lysR ygeA	transcriptional activator of lysA, autorepressor Asp/Glu_racemase family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/30);  total (21/30)

GCGGTTCCGTTCACCGTCAGCCTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGCTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCAAAAATGACGCCCTGTGCGCC  >  NC_000913/2979810‑2980071                                                                                                                                      |                                                                                                                                 gCGGTTCCGTTCACCGTCAGCCTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGa                                                                                                                             <  1:1029708/139‑1 (MQ=255) gCGGTTCCGTTCACCGTCAGCCTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGa                                                                                                                             <  1:135661/139‑1 (MQ=255)       ccGTTCACCGTCAGCCTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCAt                                                                                                                       >  2:313437/1‑139 (MQ=255)                 tCAGCCTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACa                                                                                                             >  1:856231/1‑139 (MQ=255)                 tCAGCCTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACa                                                                                                             >  1:1231316/1‑139 (MQ=255)                   aGCCTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGc                                                                                                           <  2:1058379/139‑1 (MQ=255)                     ccTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCAt                                                                                                         <  1:737899/139‑1 (MQ=255)                     ccTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCAt                                                                                                         <  1:1042307/139‑1 (MQ=255)                       tGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATcc                                                                                                       <  1:313437/139‑1 (MQ=255)                        gATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCt                                                                                                      <  1:335463/139‑1 (MQ=255)                               ccccTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCAt                                                                                                                       <  2:288047/115‑1 (MQ=255)                               ccccTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCAt                                                                                                                       >  1:288047/1‑115 (MQ=255)                               ccccTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATg                                                                                               <  2:137214/139‑1 (MQ=255)                                  cTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGAt                                                                                            <  2:664773/139‑1 (MQ=255)                                  cTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTGACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGAt                                                                                            <  2:103491/139‑1 (MQ=255)                                   tGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCg                                                                                                                >  2:159332/1‑118 (MQ=255)                                   tGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCg                                                                                                                <  1:159332/118‑1 (MQ=255)                                       ccGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAACATAAAAGCGACAGCATCCTCGGCATGGATcgccg                                                                                       <  1:828114/139‑1 (MQ=255)                                          cccGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCgg                                                                                    >  2:509502/1‑139 (MQ=255)                                           ccGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGt                                                                                   >  2:589890/1‑139 (MQ=255)                                                  caGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAac                                                                            <  1:389164/139‑1 (MQ=255)                                                  caGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTATTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAac                                                                            <  2:940805/139‑1 (MQ=255)                                                                   tttAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTgcg                                                           <  1:589890/139‑1 (MQ=255)                                                                   tttAGTCGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTgcg                                                           <  2:396855/139‑1 (MQ=255)                                                                    ttAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTgcgc                                                          <  2:1231316/139‑1 (MQ=255)                                                                          ggCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATa                                                                                                                      >  1:91185/1‑73 (MQ=255)                                                                          ggCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATa                                                                                                                      <  2:91185/73‑1 (MQ=255)                                                                          ggCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAg                                                                                                            >  2:558915/1‑83 (MQ=255)                                                                          ggCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAg                                                                                                            <  1:558915/83‑1 (MQ=255)                                                                          ggCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCtcttct                                                    >  1:1019358/1‑139 (MQ=255)                                                                          ggCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCtcttct                                                    >  1:379247/1‑139 (MQ=255)                                                                                      cGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGc                                                                                                           <  1:1031917/72‑1 (MQ=255)                                                                                      cGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGc                                                                                                           >  2:1031917/1‑72 (MQ=255)                                                                                         ggTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCt                                                                                                      <  2:1229315/74‑1 (MQ=255)                                                                                         ggTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCt                                                                                                      >  1:1229315/1‑74 (MQ=255)                                                                                          gTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGa                                                                                                                             >  1:568428/1‑50 (MQ=255)                                                                                          gTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGa                                                                                                                             <  2:568428/50‑1 (MQ=255)                                                                                           ttACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAAt                                   >  2:557143/1‑139 (MQ=255)                                                                                                   aCTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATcgccg                                                                                       >  1:795715/1‑79 (MQ=255)                                                                                                   aCTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATcgccg                                                                                       <  2:795715/79‑1 (MQ=255)                                                                                                         cACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCaa                     <  1:272187/139‑1 (MQ=255)                                                                                                         cACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGATATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCaa                     <  2:856231/139‑1 (MQ=255)                                                                                                                 ttGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGAt                                                                                            <  1:706/60‑1 (MQ=255)                                                                                                                 ttGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGAt                                                                                            >  2:706/1‑60 (MQ=255)                                                                                                                 ttGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCAAAAATGACg             >  1:481096/1‑139 (MQ=255)                                                                                                                  tGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCAAAAATGACGc            <  1:378619/139‑1 (MQ=255)                                                                                                                   gACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCAAAAATGACGcc           >  2:1297071/1‑139 (MQ=255)                                                                                                                      gCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATcgccg                                                                                       >  1:636254/1‑60 (MQ=255)                                                                                                                      gCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATcgccg                                                                                       <  2:636254/60‑1 (MQ=255)                                                                                                                        tATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCAAAAATGACGCCCtgtg      >  2:945028/1‑139 (MQ=255)                                                                                                                            ttGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCAAAAATGACGCCCTGTGCGcc  <  1:1297071/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                 GCGGTTCCGTTCACCGTCAGCCTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGCTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCAAAAATGACGCCCTGTGCGCC  >  NC_000913/2979810‑2980071

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: ‑