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_0009132,979,9430CT98.3% 188.8 / NA 60A308V (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 (1/0);  new base T (32/27);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTTCCGTTCACCGTCAGCCTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGCTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCAAAAATGACGCCCTGTGCGCCCTGT  >  NC_000913/2979810‑2980075
                                                                                                                                     |                                                                                                                                    
gCGGTTCCGTTCACCGTCAGCCTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGa                                                                                                                                 <  2:1533788/139‑1 (MQ=255)
 cGGTTCCGTTCACCGTCAGCCTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGa                                                                                                                                 <  1:1025328/138‑1 (MQ=255)
 cGGTTCCGTTCACCGTCAGCCTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGa                                                                                                                                 >  2:1025328/1‑138 (MQ=255)
              cGTCAGCCTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGa                                                                                                                   >  1:733707/1‑139 (MQ=255)
                  aGCCTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGc                                                                                                               <  2:955827/139‑1 (MQ=255)
                  aGCCTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGc                                                                                                               <  2:469983/139‑1 (MQ=255)
                     cTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATc                                                                                                            <  2:733707/139‑1 (MQ=255)
                        aTCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACa                                                                                                                 <  1:1388702/131‑1 (MQ=255)
                        aTCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACa                                                                                                                 >  2:1388702/1‑131 (MQ=255)
                            gCCCCCTGCACCGCCCGACATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCa                                                                                                     >  2:1284926/1‑139 (MQ=255)
                                  tGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATGAAAGCGACAGCATCCTCGGCATGGAt                                                                                                >  2:1962727/1‑138 (MQ=255)
                                  tGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATGAAAGCGACAGCATCCTCGGCATGGAt                                                                                                <  1:1962727/138‑1 (MQ=255)
                                       cGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATcgccgc                                                                                          >  2:754221/1‑139 (MQ=255)
                                        gCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATcgccg                                                                                           <  1:1216115/137‑1 (MQ=255)
                                        gCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATcgccg                                                                                           >  2:1216115/1‑137 (MQ=255)
                                        gCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCcgcg                                                                                         >  1:1153632/1‑139 (MQ=255)
                                        gCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCcgcg                                                                                         >  1:324796/1‑139 (MQ=255)
                                         cccGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATcgccg                                                                                           >  1:540359/1‑136 (MQ=255)
                                         cccGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATcgccg                                                                                           >  1:674718/1‑136 (MQ=255)
                                         cccGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATcgccg                                                                                           <  2:540359/136‑1 (MQ=255)
                                         cccGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATcgccg                                                                                           <  2:674718/136‑1 (MQ=255)
                                              catcaGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTAc                                                                                                                                   >  2:1746677/1‑91 (MQ=255)
                                              catcaGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTAc                                                                                                                                   <  1:1746677/91‑1 (MQ=255)
                                                 caGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAac                                                                                <  2:1361612/139‑1 (MQ=255)
                                                   gcgcTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAacac                                                                              <  1:1284926/139‑1 (MQ=255)
                                                       tGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGc                                                                          >  1:1039867/1‑139 (MQ=255)
                                                       tGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGc                                                                          >  1:1981181/1‑139 (MQ=255)
                                                        ggTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCa                                                                         <  2:594990/139‑1 (MQ=255)
                                                        ggTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCa                                                                         <  2:1214772/139‑1 (MQ=255)
                                                         gTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCaa                                                                        >  1:427513/1‑139 (MQ=255)
                                                                  tttAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTgcg                                                               <  2:1981181/139‑1 (MQ=255)
                                                                   ttAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTgcgc                                                              <  2:1039867/139‑1 (MQ=255)
                                                                                   aGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTa                                              >  1:1339875/1‑139 (MQ=255)
                                                                                    gCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTaa                                             <  2:1153632/139‑1 (MQ=255)
                                                                                       gggTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAAcc                                          >  1:67807/1‑139 (MQ=255)
                                                                                        ggTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCaaa                                                                       >  1:1164152/1‑109 (MQ=255)
                                                                                        ggTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCaaa                                                                       <  2:1164152/109‑1 (MQ=255)
                                                                                                      gTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTgcg                                                               <  1:4680/103‑1 (MQ=255)
                                                                                                      gTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTgcg                                                               >  2:4680/1‑103 (MQ=255)
                                                                                                      gTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGcc                           >  1:127550/1‑139 (MQ=255)
                                                                                                      gTCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGcc                           >  1:1767430/1‑139 (MQ=255)
                                                                                                       tCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGAt                                                                                                >  2:152948/1‑69 (MQ=255)
                                                                                                       tCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGAt                                                                                                <  1:152948/69‑1 (MQ=255)
                                                                                                       tCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTgcg                                                               >  1:614366/1‑102 (MQ=255)
                                                                                                       tCACTTCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTgcg                                                               <  2:614366/102‑1 (MQ=255)
                                                                                                           ttCTCTTGACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCaaaaa                      >  2:332396/1‑139 (MQ=255)
                                                                                                                  gACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATc                                                                                                            <  2:124994/46‑1 (MQ=255)
                                                                                                                  gACGCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATc                                                                                                            >  1:124994/1‑46 (MQ=255)
                                                                                                                     gCTATTTTGTCGTCAGTTCCGACAGCATAAAAGCGTCAGCATCCTCGGCATAGCTCGCCGAGGTAACAAACACAGGCAAAACACTACGCTCTTCTGGCACCAGAAAACGAAATTCTGTGCAGCCAAAAATTATGCCAtg            >  1:518250/1‑139 (MQ=255)
                                                                                                                     gCTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCAAAAATGACGCCCtg            >  2:1767037/1‑139 (MQ=255)
                                                                                                                      cTATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCAAAAATGACGCCCtgt           <  2:1767430/139‑1 (MQ=255)
                                                                                                                       tATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAacac                                                                              >  2:1807253/1‑71 (MQ=255)
                                                                                                                       tATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAacac                                                                              <  1:1807253/71‑1 (MQ=255)
                                                                                                                       tATTTTGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCAAAAATGACGCCCtgtg          >  1:1250418/1‑139 (MQ=255)
                                                                                                                         ttttGTCGTCAGCTACGAAAGCATCAAAGTGACCGCATCCTCGGCATGGATCGCCGCGGTATAAAACACAGGCCAAACACTGCGCTCTCCTGGCACCAGCAAACCACTTTCTGTGCAGCGAAAAATCACACCCTGTgcg        >  1:1317096/1‑139 (MQ=255)
                                                                                                                           ttGTCGTCAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCAAAAATGACGCCCTGTGCGcc      <  2:324796/139‑1 (MQ=255)
                                                                                                                            tGTCGTCCGTTACGACAGCATAAAAGCGACAGCATCCTCCGCATGGATCGCCGCGGTATCAAACACACGCAAAACACTCCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGAAGCCAAAAATGACGCCCTGTGCGccc     >  1:1932594/1‑139 (MQ=255)
                                                                                                                               cgtcAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCAAAAATGACGCCCTGTGCGCCCTGt  >  2:848165/1‑139 (MQ=255)
                                                                                                                                 tcAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTgcg                                                               >  2:305058/1‑76 (MQ=255)
                                                                                                                                 tcAGTTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTgcg                                                               <  1:305058/76‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                    
GCGGTTCCGTTCACCGTCAGCCTGATCCGCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAAGCGGGGTTACCGAAACTGGTCACTTCTCTTGACGCTATTTTGTCGTCAGCTACGACAGCATAAAAGCGACAGCATCCTCGGCATGGATCGCCGCGGTATCAAACACAGGCAAAACACTGCGCTCTTCTGGCACCAGTAAACCAATTTCTGTGCAGCCAAAAATGACGCCCTGTGCGCCCTGT  >  NC_000913/2979810‑2980075

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