Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,499 C→G 11.1% A332A (GCC→GCG)  ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,4990CG11.1% 167.5 / 10.2 63A332A (GCC→GCG) ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base C (27/29);  new base G (5/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

ATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTG  >  minE/189429‑189566
                                                                      |                                                                    
atgGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGcc                                                                      >  1:1373105/1‑71 (MQ=255)
 tgGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCt                                                                     <  1:4418760/71‑1 (MQ=255)
   gCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTg                                                                    <  1:614753/70‑1 (MQ=255)
     ttGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                                 >  1:2674363/1‑71 (MQ=255)
     ttGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                                 >  1:2428289/1‑71 (MQ=255)
      tGCTGTACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTg                                                                >  1:130997/1‑71 (MQ=255)
      tGCTGTACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTg                                                                >  1:69485/1‑71 (MQ=255)
      tGCTGTACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTg                                                                >  1:636499/1‑71 (MQ=255)
       gCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGc                                                               >  1:2076418/1‑71 (MQ=255)
        ctctACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGcc                                                              >  1:2328671/1‑71 (MQ=255)
        ctctACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGcc                                                              <  1:5310417/71‑1 (MQ=255)
          ctACAACACCGACTGGATCAAACAGCCGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCg                                                             >  1:2405386/1‑70 (MQ=255)
             caacaCCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTaca                                                         >  1:19933/1‑71 (MQ=255)
             caacaCCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTaca                                                         >  1:2258878/1‑71 (MQ=255)
             caacaCCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTaca                                                         >  1:485429/1‑71 (MQ=255)
              aacaCCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTacac                                                        >  1:2438796/1‑71 (MQ=255)
               acacCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                                 <  1:4547776/61‑1 (MQ=255)
               acacCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                                 <  1:3466722/61‑1 (MQ=255)
               acacCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                                 <  1:5215632/61‑1 (MQ=255)
               acacCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                                 <  1:302618/61‑1 (MQ=255)
                cacCGACTGGACCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACAccc                                                      <  1:5215005/71‑1 (MQ=255)
                 acCGTCTGGATCAAAGAGACGCTGGATGTCCCGTCGATTGACTTCGATTGTGCCTGGGTGCCGTACACCCa                                                     <  1:5465383/71‑1 (MQ=255)
                   cGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTAGATTCTGCCTGGGTGCCGTACACCCA‑tt                                                  >  1:7307/1‑71 (MQ=255)
                            cAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGa                                         >  1:3387964/1‑71 (MQ=255)
                             aaaCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTgg                                                                   >  1:801927/1‑45 (MQ=255)
                             aaaCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTgg                                                                   >  1:4588438/1‑45 (MQ=255)
                             aaaCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTgg                                                                   >  1:2494008/1‑45 (MQ=255)
                             aaaCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTgg                                                                   >  1:2431789/1‑45 (MQ=255)
                               aCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTacac                                                        >  1:2049515/1‑54 (MQ=255)
                               aCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTaca                                                         >  1:4164848/1‑53 (MQ=255)
                                cAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTa                                     >  1:1274597/1‑71 (MQ=255)
                                cAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTa                                     <  1:3199739/71‑1 (MQ=255)
                                   aCGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGAt                                        >  1:4213022/1‑65 (MQ=255)
                                    cGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCa                                  <  1:1872445/70‑1 (MQ=255)
                                    cGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAg                                 >  1:317984/1‑71 (MQ=255)
                                         gATGTCCCGTCGATTCACTTTGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAGGGTaa                            <  1:5687877/71‑1 (MQ=255)
                                         gATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCTATCTACCAGGGTaa                            <  1:5057101/71‑1 (MQ=255)
                                         gATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCTATCTACCAGGGTaa                            <  1:3289155/71‑1 (MQ=255)
                                         gATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAGGGTaa                            <  1:4871286/71‑1 (MQ=255)
                                           tGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAGGGTaaaa                          <  1:862260/71‑1 (MQ=255)
                                             tCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAGGGTAAAAg                         <  1:1131259/70‑1 (MQ=255)
                                                cGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCa                                  <  1:4289751/58‑1 (MQ=255)
                                                cGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCa                                  <  1:4283153/58‑1 (MQ=255)
                                                cGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCa                                  <  1:684431/58‑1 (MQ=255)
                                                  tCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTAcc                                   <  1:1856686/55‑1 (MQ=255)
                                                  tCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTAcc                                   <  1:4846197/55‑1 (MQ=255)
                                                  tCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTAcc                                   <  1:359900/55‑1 (MQ=255)
                                                    gATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAGGGTAAAAGTGGTa                    >  1:4615643/1‑68 (MQ=255)
                                                    gATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAGGGTAAAAGTGGTAt                   >  1:792084/1‑69 (MQ=255)
                                                    gATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAGGGTAAAAGTGGTAt                   >  1:2351111/1‑69 (MQ=255)
                                                    gATTCACTTCGATTCTGCCTAGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAGGGTAAAAGTGGTATGa                 <  1:4226539/71‑1 (MQ=255)
                                                       tCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAgcg              >  1:2957342/1‑71 (MQ=255)
                                                        cACTTTGACTCCGCGTTGGTGCCTTACACCAACTTCTCA‑CCGATTTACGAAGGTAAATGCGGTATGAgcg              <  1:2040140/70‑1 (MQ=255)
                                                        cACTTTGACTCCGCGTGGGTGCCTTACACCAACTTCTCA‑CCGATTTACGAAGGTAAATGCGGTATGAgcgg             <  1:609115/71‑1 (MQ=255)
                                                         aCTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAGGGTaa                            <  1:1966890/55‑1 (MQ=255)
                                                         aCTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAGGGTaa                            <  1:5018399/55‑1 (MQ=255)
                                                          cTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTAc                                    <  1:3561015/46‑1 (MQ=255)
                                                             cGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGAt                                        >  1:2358286/1‑39 (MQ=255)
                                                             cGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGAt                                        >  1:3369274/1‑39 (MQ=255)
                                                             cGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCa                                  <  1:2805639/45‑1 (MQ=255)
                                                             cGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAgg                                >  1:5009936/1‑47 (MQ=255)
                                                             cGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAgg                                >  1:4457879/1‑47 (MQ=255)
                                                             cGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAg                                 >  1:4993203/1‑46 (MQ=255)
                                                              gATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCg       <  1:5686929/71‑1 (MQ=255)
                                                                     ccTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTg  <  1:4371691/69‑1 (MQ=255)
                                                                      |                                                                    
ATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCA‑TTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTG  >  minE/189429‑189566

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