Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,475 C→G 11.5% V324V (GTC→GTG ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,4750CG11.5% 213.9 / 15.2 80V324V (GTC→GTGldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base C (36/34);  new base G (5/4);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

GGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGT  >  minE/189406‑189544
                                                                     |                                                                     
ggTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTcc                                                                      >  1:426246/1‑71 (MQ=255)
 gTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTccc                                                                     >  1:1472104/1‑71 (MQ=255)
 gTGATCACCAACTCCACATATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTccc                                                                     >  1:1073877/1‑71 (MQ=255)
  tGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCg                                                                    >  1:2854688/1‑71 (MQ=255)
  tGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCg                                                                    <  1:1479137/71‑1 (MQ=255)
    aTCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGt                                                                   <  1:18964/70‑1 (MQ=255)
     tCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTc                                                                  >  1:1779957/1‑70 (MQ=255)
     tCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCg                                                                 <  1:2492462/71‑1 (MQ=255)
     tCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCg                                                                 <  1:2126380/71‑1 (MQ=255)
      cACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGa                                                                >  1:518877/1‑71 (MQ=255)
       aCCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGAt                                                               <  1:2671379/71‑1 (MQ=255)
       aCCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGAt                                                               >  1:1186609/1‑71 (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)
          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)
              ccACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACtt                                                         >  1:1136506/1‑70 (MQ=255)
               cACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTccc                                                                     <  1:1761842/57‑1 (MQ=255)
               cACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTc                                                        <  1:1565380/70‑1 (MQ=255)
               cACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTc                                                        <  1:124671/70‑1 (MQ=255)
               cACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTc                                                        >  1:1987341/1‑70 (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)
                           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:2404752/1‑71 (MQ=255)
                                  tACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTa                                    >  1:3923733/1‑71 (MQ=255)
                                  tACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTa                                    <  1:3537515/71‑1 (MQ=255)
                                  tACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTa                                    <  1:3480662/71‑1 (MQ=255)
                                  tACAACACCGACTGGATCAAACAGACTCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCg                                      <  1:1280254/69‑1 (MQ=255)
                                  tACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCg                                      <  1:1474286/69‑1 (MQ=255)
                                  tACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCg                                      <  1:1557950/69‑1 (MQ=255)
                                  tACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTa                                    >  1:540775/1‑71 (MQ=255)
                                  tACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTa                                    >  1:476107/1‑71 (MQ=255)
                                   acaacaCCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGcctt                                     <  1:673216/69‑3 (MQ=255)
                                   acaacaCCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTa                                    >  1:1123100/1‑70 (MQ=255)
                                   acaacaCCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTa                                    >  1:2144761/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)
                                                gATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACAcc                                >  1:755731/1‑61 (MQ=255)
                                                gATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCa                        <  1:1265227/69‑1 (MQ=255)
                                                gATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCa                        <  1:2709764/69‑1 (MQ=255)
                                                gATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCAt                       <  1:1294762/70‑1 (MQ=255)
                                                   cAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTaca                                  >  1:1569762/1‑56 (MQ=255)
                                                   cAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTaca                                  >  1:416130/1‑56 (MQ=255)
                                                     aaCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATc                 <  1:2252234/71‑1 (MQ=255)
                                                     aaCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATc                 >  1:382467/1‑71 (MQ=255)
                                                      aCAGACGCTGGATGTCCCGTCGATTGACTTCGATTCTGCCTGGGTGCCGTacac                                 >  1:1430428/1‑54 (MQ=255)
                                                      aCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCt                <  1:2792172/71‑1 (MQ=255)
                                                      aCAGACGCTGGATGTCCCGTCGATTCACTTCGATTATGCCTGGGTGCCGTACACCCATTTTCAta                      >  1:1193992/1‑64 (MQ=255)
                                                       cAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCg                    >  1:2710033/1‑66 (MQ=255)
                                                          aCGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTAc              >  1:483846/1‑69 (MQ=255)
                                                            gCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTgg                                            >  1:2573694/1‑37 (MQ=255)
                                                            gCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCa                        <  1:872736/57‑1 (MQ=255)
                                                                gATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTa       >  1:651933/1‑70 (MQ=255)
                                                                gATGTCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaaa     <  1:667737/71‑1 (MQ=255)
                                                                 aTGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaaa     <  1:3153832/71‑1 (MQ=255)
                                                                    tCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGt  <  1:2497198/71‑1 (MQ=255)
                                                                     |                                                                     
GGTGATCACCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGT  >  minE/189406‑189544

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: