Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,508 G→T 11.3% P335P (CCG→CCT ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,5080GT11.3% 162.1 / 5.4 53P335P (CCG→CCTldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base G (18/29);  new base T (2/4);  total (20/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

CTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTG  >  minE/189439‑189566
                                                                     |                                                           
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)
      cacCGACTGGACCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACAccc                                                      <  1:5215005/71‑1 (MQ=255)
       acCGTCTGGATCAAAGAGACGCTGGATGTCCCGTCGATTGACTTCGATTGTGCCTGGGTGCCGTACACCCa                                                     <  1:5465383/71‑1 (MQ=255)
         cGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTAGATTCTGCCTGGGTGCCGTACACCCAtt                                                  >  1:7307/1‑71 (MQ=255)
                  cAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGa                                         >  1:3387964/1‑71 (MQ=255)
                     aCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTacac                                                        >  1:2049515/1‑54 (MQ=255)
                     aCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTaca                                                         >  1:4164848/1‑53 (MQ=255)
                      cAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTa                                     <  1:3199739/71‑1 (MQ=255)
                      cAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTa                                     >  1:1274597/1‑71 (MQ=255)
                         aCGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGAt                                        >  1:4213022/1‑65 (MQ=255)
                          cGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCa                                  <  1:1872445/70‑1 (MQ=255)
                          cGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAg                                 >  1:317984/1‑71 (MQ=255)
                               gATGTCCCGTCGATTCACTTTGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaa                            <  1:5687877/71‑1 (MQ=255)
                               gATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCTATCTACCAGGGTaa                            <  1:5057101/71‑1 (MQ=255)
                               gATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCTATCTACCAGGGTaa                            <  1:3289155/71‑1 (MQ=255)
                               gATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaa                            <  1:4871286/71‑1 (MQ=255)
                                 tGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaaaa                          <  1:862260/71‑1 (MQ=255)
                                   tCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAg                         <  1:1131259/70‑1 (MQ=255)
                                      cGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCa                                  <  1:684431/58‑1 (MQ=255)
                                      cGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCa                                  <  1:4289751/58‑1 (MQ=255)
                                      cGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCa                                  <  1:4283153/58‑1 (MQ=255)
                                        tCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTAcc                                   <  1:4846197/55‑1 (MQ=255)
                                        tCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTAcc                                   <  1:1856686/55‑1 (MQ=255)
                                        tCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTAcc                                   <  1:359900/55‑1 (MQ=255)
                                          gATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTa                    >  1:4615643/1‑68 (MQ=255)
                                          gATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTAt                   >  1:2351111/1‑69 (MQ=255)
                                          gATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTAt                   >  1:792084/1‑69 (MQ=255)
                                          gATTCACTTCGATTCTGCCTAGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGa                 <  1:4226539/71‑1 (MQ=255)
                                             tCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAgcg              >  1:2957342/1‑71 (MQ=255)
                                              cACTTTGACTCCGCGTTGGTGCCTTACACCAACTTCTCACCGATTTACGAAGGTAAATGCGGTATGAgcg              <  1:2040140/70‑1 (MQ=255)
                                              cACTTTGACTCCGCGTGGGTGCCTTACACCAACTTCTCACCGATTTACGAAGGTAAATGCGGTATGAgcgg             <  1:609115/71‑1 (MQ=255)
                                               aCTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaa                            <  1:1966890/55‑1 (MQ=255)
                                               aCTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaa                            <  1:5018399/55‑1 (MQ=255)
                                                cTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTAc                                    <  1:3561015/46‑1 (MQ=255)
                                                   cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGAt                                        >  1:3369274/1‑39 (MQ=255)
                                                   cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGAt                                        >  1:2358286/1‑39 (MQ=255)
                                                   cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCa                                  <  1:2805639/45‑1 (MQ=255)
                                                   cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAgg                                >  1:4457879/1‑47 (MQ=255)
                                                   cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAgg                                >  1:5009936/1‑47 (MQ=255)
                                                   cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAg                                 >  1:4993203/1‑46 (MQ=255)
                                                    gATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCg       <  1:5686929/71‑1 (MQ=255)
                                                           ccTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTg  <  1:4371691/69‑1 (MQ=255)
                                                               ggTGCGGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTgg                      <  1:3907938/45‑1 (MQ=255)
                                                               ggTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTgg                      <  1:4903872/45‑1 (MQ=255)
                                                                gTGCCTTACACCAACTTCTCACCGATTTACGAAGGTAAATGCGGTATGAGCGGTGGCCgtgt    <  1:1284978/62‑1 (MQ=255)
                                                                gTGCCTTACACCAACTTCTCACCGATTTACGAAGGTAAATGCGGTATGAGCGGTGGCCgtgt    <  1:3328815/62‑1 (MQ=255)
                                                                gTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGt                     <  1:2329540/45‑1 (MQ=255)
                                                                gTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGt                     <  1:4528008/45‑1 (MQ=255)
                                                                 tGCCGTACACCCATTTTCATCCGATCTACCAGGGTaa                            <  1:29654/37‑1 (MQ=255)
                                                                   ccGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTgg                      <  1:2079447/41‑1 (MQ=255)
                                                                     gTACACCCATTTTCATCCGATCTACCAGGGTaaaa                          >  1:4486317/1‑35 (MQ=255)
                                                                     |                                                           
CTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTG  >  minE/189439‑189566

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