Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,453 G→T 10.3% W317L (TGG→TTG) ‡ ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,4530GT10.3% 273.5 / 10.1 87W317L (TGG→TTG) ‡ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base G (40/38);  new base T (5/4);  total (45/42)
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.36e-01

CAAGCACAATGGCCGGTTCATGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTT  >  minE/189383‑189520
                                                                      |                                                                   
cAAGCACAATGGCCGGTTCATGCGGTGATCACCAACTCCACCTATGATTGCTTGCTCTACAACACCGACTg                                                                     >  1:614252/1‑71 (MQ=255)
  aGCACAATGGCCGGTTCATGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGa                                                                   <  1:1781864/71‑1 (MQ=255)
   gCACAATGGCCGGTTCATGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGAt                                                                  >  1:96285/1‑71 (MQ=255)
    cacaATGGCCGGTTCATGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGAt                                                                  >  1:3254704/1‑70 (MQ=255)
     acaATGGCCGGTTCATGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATc                                                                 >  1:3245451/1‑70 (MQ=255)
      caATGGCCGGTTCATGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCaa                                                               <  1:3846992/71‑1 (MQ=255)
      caATGGCCGGTTCATGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCa                                                                >  1:2117228/1‑70 (MQ=255)
      caATGGCCGGTTCATGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCa                                                                >  1:1124031/1‑70 (MQ=255)
       aaTGGCCGGTTCATGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCaaa                                                              >  1:3158935/1‑71 (MQ=255)
         tGGCCGGTTCATGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACa                                                            <  1:368761/71‑1 (MQ=255)
               gTTCATGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGc                                                       <  1:1536412/70‑1 (MQ=255)
                ttCATGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTg                                                     <  1:743585/71‑1 (MQ=255)
                 tCATGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTg                                                     >  1:2940996/1‑70 (MQ=255)
                  cATGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGa                                                   >  1:3651934/1‑71 (MQ=255)
                    tGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCaaa                                                              <  1:3611367/58‑1 (MQ=255)
                    tGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCaaa                                                              <  1:3684712/58‑1 (MQ=255)
                    tGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCaaa                                                              <  1:1923082/58‑1 (MQ=255)
                       ggTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACa                                                            <  1:3825692/57‑1 (MQ=255)
                       ggTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTcc                                              >  1:426246/1‑71 (MQ=255)
                        gTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTccc                                             >  1:1472104/1‑71 (MQ=255)
                        gTGATCACCAACTCCACATATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTccc                                             >  1:1073877/1‑71 (MQ=255)
                         tGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCg                                            <  1:1479137/71‑1 (MQ=255)
                         tGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCg                                            >  1:2854688/1‑71 (MQ=255)
                           aTCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGt                                           <  1:18964/70‑1 (MQ=255)
                            tCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTc                                          >  1:1779957/1‑70 (MQ=255)
                            tCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCg                                         <  1:2126380/71‑1 (MQ=255)
                            tCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCg                                         <  1:2492462/71‑1 (MQ=255)
                             cacaAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAg                                                           >  1:3764650/1‑52 (MQ=255)
                             cACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGa                                        >  1:518877/1‑71 (MQ=255)
                              aCCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGAt                                       >  1:1186609/1‑71 (MQ=255)
                              aCCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGAt                                       <  1:2671379/71‑1 (MQ=255)
                               ccAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGa                                        <  1:1395885/69‑1 (MQ=255)
                               ccAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGAtt                                      <  1:3814791/71‑1 (MQ=255)
                                cAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGAtt                                      >  1:108453/1‑70 (MQ=255)
                                cAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTc                                     <  1:1216009/71‑1 (MQ=255)
                                cAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTc                                     <  1:3592842/71‑1 (MQ=255)
                                 aaCTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCaaa                                                              >  1:3764887/1‑45 (MQ=255)
                                 aaCTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGt                                                >  1:2198093/1‑59 (MQ=255)
                                 aaCTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTc                                          >  1:568710/1‑65 (MQ=255)
                                 aaCTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCa                                    >  1:2619085/1‑71 (MQ=255)
                                 aaCTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCa                                    >  1:110506/1‑71 (MQ=255)
                                  aCTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCAc                                   <  1:3885216/71‑1 (MQ=255)
                                     ccACCTATGATGGCTTGCTCTACAACACCGACTGGATCa                                                                <  1:3482720/39‑1 (MQ=255)
                                     ccACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACtt                                 >  1:1136506/1‑70 (MQ=255)
                                      cACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTccc                                             <  1:1761842/57‑1 (MQ=255)
                                      cACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTc                                >  1:1987341/1‑70 (MQ=255)
                                      cACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTc                                <  1:124671/70‑1 (MQ=255)
                                      cACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTc                                <  1:1565380/70‑1 (MQ=255)
                                        ccTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGt                                           >  1:2022580/1‑57 (MQ=255)
                                         cTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGAtt                            >  1:789431/1‑71 (MQ=255)
                                         cTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGAtt                            >  1:3353835/1‑71 (MQ=255)
                                         cTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGAt                             >  1:647505/1‑70 (MQ=255)
                                           atgatgGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGa                              >  1:2435795/1‑67 (MQ=255)
                                              atgGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTcc                                              >  1:1345790/1‑48 (MQ=255)
                                               tgGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCt                      <  1:2713246/71‑1 (MQ=255)
                                                 gCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTgg                    <  1:168469/71‑1 (MQ=255)
                                                 gCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTg                     <  1:2679645/70‑1 (MQ=255)
                                                  cTTGCTCTACAACACCGACTGGATCAAACAGACGCTgg                                                    <  1:2525929/38‑1 (MQ=255)
                                                  cTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTggg                   <  1:3873571/71‑1 (MQ=255)
                                                  cTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTcgg                   <  1:2089225/71‑1 (MQ=255)
                                                    tGCTGTACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTg                 >  1:3499844/1‑71 (MQ=255)
                                                     gCTGTACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGc                >  1:87838/1‑71 (MQ=255)
                                                       tctACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGa                              >  1:1899268/1‑55 (MQ=255)
                                                         tACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGcctt             <  1:3242697/70‑3 (MQ=255)
                                                         tACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTa            <  1:3537515/71‑1 (MQ=255)
                                                         tACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTa            >  1:2404752/1‑71 (MQ=255)
                                                         tACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTa            >  1:3923733/1‑71 (MQ=255)
                                                         tACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTa            <  1:3480662/71‑1 (MQ=255)
                                                         tACAACACCGACTGGATCAAACAGACTCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCg              <  1:1280254/69‑1 (MQ=255)
                                                         tACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCg              <  1:1557950/69‑1 (MQ=255)
                                                         tACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCg              <  1:1474286/69‑1 (MQ=255)
                                                         tACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTa            >  1:476107/1‑71 (MQ=255)
                                                         tACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTa            >  1:540775/1‑71 (MQ=255)
                                                          acaacaCCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGcctt             <  1:673216/69‑3 (MQ=255)
                                                          acaacaCCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTa            >  1:2144761/1‑70 (MQ=255)
                                                          acaacaCCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTa            >  1:1123100/1‑70 (MQ=255)
                                                           caacaCCGACTTCATCAACAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTaca          >  1:3742465/1‑71 (MQ=255)
                                                              cacCGACTGGATCAAACAGACGCTGGATGTCCCGTCGAtt                                      <  1:424376/40‑1 (MQ=255)
                                                               acagACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCa      >  1:669736/4‑71 (MQ=255)
                                                               acCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCa      >  1:3811930/1‑71 (MQ=255)
                                                               acCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACATAGATTCt                          >  1:3890559/1‑51 (MQ=255)
                                                                ccGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTa            <  1:59962/64‑1 (MQ=255)
                                                                  gACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCAttt   <  1:3002159/71‑1 (MQ=255)
                                                                  gACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCAttt   <  1:757113/71‑1 (MQ=255)
                                                                   aCTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCAtttt  >  1:2482060/1‑71 (MQ=255)
                                                                   aCTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCAtttt  >  1:3379048/1‑71 (MQ=255)
                                                                    cTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCAtttt  <  1:3861866/70‑1 (MQ=255)
                                                                     tGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTaca          <  1:1512994/61‑1 (MQ=255)
                                                                      ggATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                  >  1:781594/1‑52 (MQ=255)
                                                                      |                                                                   
CAAGCACAATGGCCGGTTCATGCGGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTT  >  minE/189383‑189520

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: